Simulation and Calibration of a Fully Bayesian Marked Multidimensional Hawkes Process with Dissimilar Decays
|
|
- Horatio Jenkins
- 5 years ago
- Views:
Transcription
1 Simulation and Calibration of a Fully Bayesian Marked Multidimensional Hawkes Process with Dissimilar Decays Kar Wai Lim, Young Lee, Leif Hanlen, Hongbiao Zhao Australian National University Data61/CSIRO Xiamen University 1 Proceedings of the 8th Asian Conference in Machine Learning (ACML) / 19
2 Talk Outline Introduction on Hawkes Processes Simulation of Hawkes Processes Bayesian Inference for Hawkes 2 / 19
3 Background Poisson distributions Commonly used to model the number of times an event occurs in an interval of time or space. Textbook example: the number of cars passing an intersection in half an hour. Figure: (left) Probability mass functions (right) Observed histogram Additive Property: If X Poi(λ1), Y Poi(λ 2), then X + Y Poi(λ 1 + λ 2) 3 / 19
4 Background (Homogeneous) Poisson process It is a stochastic process that keep track of the running counts of an event over time (and space). For example, the number of cars passing an intersection is an evolution of counts with time: Call the evolution of counts as the counting process N(t) and the times of an event happening the event times t i. 4 / 19
5 Background Properties of Poisson process The counting process starts at zero: N(t = 0) = 0. Parameterised by the expected number of events per unit time, e.g. λ = 3 vehicles per minute. The counting process N(t) at time t follows Poi(λt). (number of events observed until time t) The difference (also called increment) in counting processes N(t) N(s) Poi(λ(t s)) t > s Superposition property: If N(t) PP(λ1), M(t) PP(λ 2), then N(t) + M(t) PP(λ 1 + λ 2) 5 / 19
6 Background What if some events are more frequent at certain times? More cars during peak hours! Instead of constant intensity, allow the intensity to vary with time: λ(t) becomes a function of time. 6 / 19
7 Background Extension: Inhomogeneous Poisson process Example λ(t): Piecewise linear; Piecewise polynomial; Cyclical functions such as sine curve. Figure: Generated data 7 / 19
8 Background Properties of inhomogeneous Poisson process (IPP) The counting process starts at zero: N(t = 0) = 0. Parameterised by intensity function λ(t). t The counting process at time t follows Poi( λ(u) du). 0 The difference (also called increment) in counting processes ( t ) N(t) N(s) Poi λ(u) du t > s s Superposition property still holds: If N(t) IPP(λ1(t)), M(t) IPP(λ 2(t)), then N(t) + M(t) IPP(λ 1(t) + λ 2(t)) 8 / 19
9 Hawkes Processes Hawkes process is a point process in which an occurrence of an event triggers future events (self-excitation). Our formulation of Hawkes (univariate): λ(t) = µ(t) + N(T ) i=1: t>t i Y i e δ(t ti ) Decaying background intensity: µ(t) = k + Y (0) e δ t Random self excitations: Y i i.i.d. Gamma Terminology: ti, i = 1,..., N(T ) is a sequence of non-negative random variables such that t i < t i+1, known as event times. i = t i t i 1 is called the inter-arrival time. 9 / 19
10 Multivariate Hawkes Captures multiple event types for which the events mutually excite one another. Our formulation (Bivariate Hawkes): λ 1 (t) = µ 1 (t) + N 1 (t) Y 1 1,j e δ1 1 t + N 2 (t) Y 2 1,j e δ2 1 t j=1: t t 1 j j=1: t t 2 j λ 2 (t) = µ 2 (t) + N 1 (t) Y 1 2,j e δ1 2 t + N 2 (t) Y 2 2,j e δ2 2 t j=1: t t 1 j j=1: t t 2 j where λ 1 (t) and λ 2 (t) are the intensity functions for events 1 and 2, respectively. Note that the decay parameters δ are different for each process. 10 / 19
11 Illustration of Multivariate 0 Hawkes Time Intensity Time Dimension 1 Dimension 2 Dimension 3 Counting Process Time 11 / 19
12 Detour: Stationarity of Hawkes process Due to self-excitation property, a Hawkes process is only stable (stationary) when certain condition is satisfied. The intensity process λ(t) explodes if this condition is not satisfied: Causing chain reactions: intensity increases more future events further increases in intensity... We present a theoretical result on the expected stationary intensities for our Hawkes formulation. [Extension of Hawkes (1971) and Bacry et al. (2015)] 12 / 19
13 Outline Introduction on Hawkes Processes Simulation of Hawkes Processes Bayesian Inference for Hawkes 13 / 19
14 Simulation of Hawkes Processes There are three categories of simulation methods. 1. Inverse Sampling (Ozaki, 1979) Derives cdf (cumulative distribution function) of inter-arrival times, then performs inverse sampling. Cdf cannot be inverted directly so approximation is needed. 2. Thinning (Lewis and Shedler, 1979; Ogata, 1981) Simulate samples from a Poisson process and then thin the samples. Akin to a rejection sampler. 3. Cluster method (Brix & Kendall, 2002; Møller & Rasmussen, 2005) Recast Hawkes using a Poisson cluster representation. Each observed event generates an IPP. Superposition of all of them forms a Hawkes process. Notable mention: exact sampler of Dassios & Zhao (2013) Performs inverse sampling without approximation by decomposing a variable into two need to satisfy a Markovian constraint. Our method: exploits superposition theory and first order statistics for efficient sampling. 14 / 19
15 Our Simulation Method in One Slide Illustration with bivariate Hawkes λ 1(t) = µ 1(t) + N 1 (t) Y 1 1,j e δ1 1 t + N 2 (t) Y 2 1,j e δ2 1 t j=1: t t 1 j j=1: t t 2 j λ 2(t) = µ 2(t) + N 1 (t) Y 1 2,j e δ1 2 t + N 2 (t) Y 2 2,j e δ2 2 t j=1: t t 1 j j=1: t t 2 j A Hawkes with intensity λ 1 (t) is a superposition of IPP (with intensities µ 1 etc). Inter-arrival times (a i, b i, c i...) for these IPP can be sampled easily. We show that the inter-arrival time i for a Hawkes process is a first order statistics of these inter-arrival times: i = min{a i, b i, c i...} No need to resort to approximation or satisfy Markovian constraint. 2 Note: efficient caching can be performed if the Hawkes is Markov. 15 / 19
16 Simulation Statistics We compare the simulated statistics against theoretical expectations (over 1 million simulation paths): Intensity Dimension 1 Dimension 2 Dimension Time Table 1: Comparison of simulated mean intensities against their theore ten fixed time t. The simulated mean intensities are comput simulated intensities over one million sample paths. The perce the simulated means against the theoretical expectations are m the correctness of the proposed simulation algorithm. [Maybe 2 for space] Process m = 1 Process m = 2 Pro Time Sim. Expt. %Diff. Sim. Expt. %Diff. Sim. E In addition, through the simulation of univariate Hawkes processe parison against the inverse sampling method of Ozaki (1979) and the cl of Brix and Kendall (2002). Since implementation on these methods we able, we implement them based on the algorithms outlined in the resp might not be fully optimised, however, the same can be said to the im own sampler. The simulation settings (subscripts m and i suppressed) Figure: Plot of simulated mean intensities vs the theoretical stationary average intensities of the three-dimensional Hawkes processes. Verifies that our algorithm and implementation is correct. See paper for other results / 19
17 Outline Introduction on Hawkes Processes Simulation of Hawkes Processes Bayesian Inference for Hawkes 17 / 19
18 Bayesian Inference in One Slide Fully Gibbs sampling achieved by Auxiliary variables augmentation introduce additional parameters called branching structures allow decoupling of existing parameters. Adaptive rejection sampling (ARS) for variables that do not have known posterior distributions, we show conditions for which the posteriors are log-concave, and sample via ARS. On simulated data, we demonstrate that the parameters learned using Bayesian inference is accurate and superior to MLE: See paper for application on modelling Dark Networks. 18 / 19
19 Summary Theoretical result on expected stationary intensities Simulation of multivariate Hawkes with superposition theory and first order statistics Bayesian inference on Hawkes with auxiliary variable augmentation and adaptive rejection sampling 19 / 19
Simulation and Calibration of a Fully Bayesian Marked Multidimensional Hawkes Process with Dissimilar Decays
JMLR: Workshop and Conference Proceedings 63:238 262, 2016 ACML 2016 Simulation and Calibration of a Fully Bayesian Marked Multidimensional Hawkes Process with Dissimilar Decays Kar Wai Lim Data61 / CSIRO
More informationThe exponential distribution and the Poisson process
The exponential distribution and the Poisson process 1-1 Exponential Distribution: Basic Facts PDF f(t) = { λe λt, t 0 0, t < 0 CDF Pr{T t) = 0 t λe λu du = 1 e λt (t 0) Mean E[T] = 1 λ Variance Var[T]
More informationRENEWAL PROCESSES AND POISSON PROCESSES
1 RENEWAL PROCESSES AND POISSON PROCESSES Andrea Bobbio Anno Accademico 1997-1998 Renewal and Poisson Processes 2 Renewal Processes A renewal process is a point process characterized by the fact that the
More informationManual for SOA Exam MLC.
Chapter 10. Poisson processes. Section 10.5. Nonhomogenous Poisson processes Extract from: Arcones Fall 2009 Edition, available at http://www.actexmadriver.com/ 1/14 Nonhomogenous Poisson processes Definition
More informationTemporal point processes: the conditional intensity function
Temporal point processes: the conditional intensity function Jakob Gulddahl Rasmussen December 21, 2009 Contents 1 Introduction 2 2 Evolutionary point processes 2 2.1 Evolutionarity..............................
More informationLocally stationary Hawkes processes
Locally stationary Hawkes processes François Roue LTCI, CNRS, Télécom ParisTech, Université Paris-Saclay Talk based on Roue, von Sachs, and Sansonnet [2016] 1 / 40 Outline Introduction Non-stationary Hawkes
More informationExponential Distribution and Poisson Process
Exponential Distribution and Poisson Process Stochastic Processes - Lecture Notes Fatih Cavdur to accompany Introduction to Probability Models by Sheldon M. Ross Fall 215 Outline Introduction Exponential
More informationA Dynamic Contagion Process with Applications to Finance & Insurance
A Dynamic Contagion Process with Applications to Finance & Insurance Angelos Dassios Department of Statistics London School of Economics Angelos Dassios, Hongbiao Zhao (LSE) A Dynamic Contagion Process
More informationLet (Ω, F) be a measureable space. A filtration in discrete time is a sequence of. F s F t
2.2 Filtrations Let (Ω, F) be a measureable space. A filtration in discrete time is a sequence of σ algebras {F t } such that F t F and F t F t+1 for all t = 0, 1,.... In continuous time, the second condition
More informationPROBABILITY DISTRIBUTIONS
Review of PROBABILITY DISTRIBUTIONS Hideaki Shimazaki, Ph.D. http://goo.gl/visng Poisson process 1 Probability distribution Probability that a (continuous) random variable X is in (x,x+dx). ( ) P x < X
More informationFast MCMC sampling for Markov jump processes and extensions
Fast MCMC sampling for Markov jump processes and extensions Vinayak Rao and Yee Whye Teh Rao: Department of Statistical Science, Duke University Teh: Department of Statistics, Oxford University Work done
More informationECE 313 Probability with Engineering Applications Fall 2000
Exponential random variables Exponential random variables arise in studies of waiting times, service times, etc X is called an exponential random variable with parameter λ if its pdf is given by f(u) =
More informationContinuous-time Markov Chains
Continuous-time Markov Chains Gonzalo Mateos Dept. of ECE and Goergen Institute for Data Science University of Rochester gmateosb@ece.rochester.edu http://www.ece.rochester.edu/~gmateosb/ October 23, 2017
More informationLikelihood Function for Multivariate Hawkes Processes
Lielihood Function for Multivariate Hawes Processes Yuanda Chen January, 6 Abstract In this article we discuss the lielihood function for an M-variate Hawes process and derive the explicit formula for
More informationHigh-frequency data modelling using Hawkes processes
High-frequency data modelling using Hawkes processes Valérie Chavez-Demoulin 1 joint work J.A McGill 1 Faculty of Business and Economics, University of Lausanne, Switzerland Boulder, April 2016 Boulder,
More informationHigh-frequency data modelling using Hawkes processes
Valérie Chavez-Demoulin joint work with High-frequency A.C. Davison data modelling and using A.J. Hawkes McNeil processes(2005), J.A EVT2013 McGill 1 /(201 High-frequency data modelling using Hawkes processes
More informationPoissonprocessandderivationofBellmanequations
APPENDIX B PoissonprocessandderivationofBellmanequations 1 Poisson process Let us first define the exponential distribution Definition B1 A continuous random variable X is said to have an exponential distribution
More informationarxiv: v1 [q-fin.rm] 27 Jun 2017
Risk Model Based on General Compound Hawkes Process Anatoliy Swishchuk 1 2 arxiv:1706.09038v1 [q-fin.rm] 27 Jun 2017 Abstract: In this paper, we introduce a new model for the risk process based on general
More informationPoisson Processes. Particles arriving over time at a particle detector. Several ways to describe most common model.
Poisson Processes Particles arriving over time at a particle detector. Several ways to describe most common model. Approach 1: a) numbers of particles arriving in an interval has Poisson distribution,
More informationSTA 4273H: Statistical Machine Learning
STA 4273H: Statistical Machine Learning Russ Salakhutdinov Department of Computer Science! Department of Statistical Sciences! rsalakhu@cs.toronto.edu! h0p://www.cs.utoronto.ca/~rsalakhu/ Lecture 7 Approximate
More informationABC methods for phase-type distributions with applications in insurance risk problems
ABC methods for phase-type with applications problems Concepcion Ausin, Department of Statistics, Universidad Carlos III de Madrid Joint work with: Pedro Galeano, Universidad Carlos III de Madrid Simon
More informationAN EM ALGORITHM FOR HAWKES PROCESS
AN EM ALGORITHM FOR HAWKES PROCESS Peter F. Halpin new york university December 17, 2012 Correspondence should be sent to Dr. Peter F. Halpin 246 Greene Street, Office 316E New York, NY 10003-6677 E-Mail:
More informationReview. DS GA 1002 Statistical and Mathematical Models. Carlos Fernandez-Granda
Review DS GA 1002 Statistical and Mathematical Models http://www.cims.nyu.edu/~cfgranda/pages/dsga1002_fall16 Carlos Fernandez-Granda Probability and statistics Probability: Framework for dealing with
More informationMetropolis-Hastings Algorithm
Strength of the Gibbs sampler Metropolis-Hastings Algorithm Easy algorithm to think about. Exploits the factorization properties of the joint probability distribution. No difficult choices to be made to
More informationSection 7.5: Nonstationary Poisson Processes
Section 7.5: Nonstationary Poisson Processes Discrete-Event Simulation: A First Course c 2006 Pearson Ed., Inc. 0-13-142917-5 Discrete-Event Simulation: A First Course Section 7.5: Nonstationary Poisson
More information1 A Tutorial on Hawkes Processes
1 A Tutorial on Hawkes Processes for Events in Social Media arxiv:1708.06401v2 [stat.ml] 9 Oct 2017 Marian-Andrei Rizoiu, The Australian National University; Data61, CSIRO Young Lee, Data61, CSIRO; The
More informationCOPYRIGHTED MATERIAL CONTENTS. Preface Preface to the First Edition
Preface Preface to the First Edition xi xiii 1 Basic Probability Theory 1 1.1 Introduction 1 1.2 Sample Spaces and Events 3 1.3 The Axioms of Probability 7 1.4 Finite Sample Spaces and Combinatorics 15
More information3. Poisson Processes (12/09/12, see Adult and Baby Ross)
3. Poisson Processes (12/09/12, see Adult and Baby Ross) Exponential Distribution Poisson Processes Poisson and Exponential Relationship Generalizations 1 Exponential Distribution Definition: The continuous
More informationSimulation, Estimation and Applications of Hawkes Processes
Simulation, Estimation and Applications of Hawkes Processes A PROJECT SUBMITTED TO THE FACULTY OF THE GRADUATE SCHOOL OF THE UNIVERSITY OF MINNESOTA BY Katarzyna Obral IN PARTIAL FULFILLMENT OF THE REQUIREMENTS
More informationLecture 20. Poisson Processes. Text: A Course in Probability by Weiss STAT 225 Introduction to Probability Models March 26, 2014
Lecture 20 Text: A Course in Probability by Weiss 12.1 STAT 225 Introduction to Probability Models March 26, 2014 Whitney Huang Purdue University 20.1 Agenda 1 2 20.2 For a specified event that occurs
More informationSolutions For Stochastic Process Final Exam
Solutions For Stochastic Process Final Exam (a) λ BMW = 20 0% = 2 X BMW Poisson(2) Let N t be the number of BMWs which have passes during [0, t] Then the probability in question is P (N ) = P (N = 0) =
More informationChapter 2. Poisson Processes. Prof. Shun-Ren Yang Department of Computer Science, National Tsing Hua University, Taiwan
Chapter 2. Poisson Processes Prof. Shun-Ren Yang Department of Computer Science, National Tsing Hua University, Taiwan Outline Introduction to Poisson Processes Definition of arrival process Definition
More information(b) What is the variance of the time until the second customer arrives, starting empty, assuming that we measure time in minutes?
IEOR 3106: Introduction to Operations Research: Stochastic Models Fall 2006, Professor Whitt SOLUTIONS to Final Exam Chapters 4-7 and 10 in Ross, Tuesday, December 19, 4:10pm-7:00pm Open Book: but only
More informationLecture 2: Repetition of probability theory and statistics
Algorithms for Uncertainty Quantification SS8, IN2345 Tobias Neckel Scientific Computing in Computer Science TUM Lecture 2: Repetition of probability theory and statistics Concept of Building Block: Prerequisites:
More informationNonparametric Bayesian Methods - Lecture I
Nonparametric Bayesian Methods - Lecture I Harry van Zanten Korteweg-de Vries Institute for Mathematics CRiSM Masterclass, April 4-6, 2016 Overview of the lectures I Intro to nonparametric Bayesian statistics
More informationarxiv: v1 [stat.co] 2 Nov 2017
Binary Bouncy Particle Sampler arxiv:1711.922v1 [stat.co] 2 Nov 217 Ari Pakman Department of Statistics Center for Theoretical Neuroscience Grossman Center for the Statistics of Mind Columbia University
More informationThe story of the film so far... Mathematics for Informatics 4a. Continuous-time Markov processes. Counting processes
The story of the film so far... Mathematics for Informatics 4a José Figueroa-O Farrill Lecture 19 28 March 2012 We have been studying stochastic processes; i.e., systems whose time evolution has an element
More informationRenewal theory and its applications
Renewal theory and its applications Stella Kapodistria and Jacques Resing September 11th, 212 ISP Definition of a Renewal process Renewal theory and its applications If we substitute the Exponentially
More informationTemporal Point Processes the Conditional Intensity Function
Temporal Point Processes the Conditional Intensity Function Jakob G. Rasmussen Department of Mathematics Aalborg University Denmark February 8, 2010 1/10 Temporal point processes A temporal point process
More informationCTDL-Positive Stable Frailty Model
CTDL-Positive Stable Frailty Model M. Blagojevic 1, G. MacKenzie 2 1 Department of Mathematics, Keele University, Staffordshire ST5 5BG,UK and 2 Centre of Biostatistics, University of Limerick, Ireland
More informationFactorization of Seperable and Patterned Covariance Matrices for Gibbs Sampling
Monte Carlo Methods Appl, Vol 6, No 3 (2000), pp 205 210 c VSP 2000 Factorization of Seperable and Patterned Covariance Matrices for Gibbs Sampling Daniel B Rowe H & SS, 228-77 California Institute of
More informationModelling Trades-Through in a Limit Order Book Using Hawkes Processes
Vol. 6, 212-22 June 14, 212 http://dx.doi.org/1.518/economics-ejournal.ja.212-22 Modelling Trades-Through in a Limit Order Book Using Hawkes Processes Ioane Muni Toke Ecole Centrale Paris and University
More informationStochastic Processes
Stochastic Processes Stochastic Process Non Formal Definition: Non formal: A stochastic process (random process) is the opposite of a deterministic process such as one defined by a differential equation.
More informationBayesian Methods for Machine Learning
Bayesian Methods for Machine Learning CS 584: Big Data Analytics Material adapted from Radford Neal s tutorial (http://ftp.cs.utoronto.ca/pub/radford/bayes-tut.pdf), Zoubin Ghahramni (http://hunch.net/~coms-4771/zoubin_ghahramani_bayesian_learning.pdf),
More informationIEOR 6711: Stochastic Models I, Fall 2003, Professor Whitt. Solutions to Final Exam: Thursday, December 18.
IEOR 6711: Stochastic Models I, Fall 23, Professor Whitt Solutions to Final Exam: Thursday, December 18. Below are six questions with several parts. Do as much as you can. Show your work. 1. Two-Pump Gas
More informationprocesses Sai Xiao, Athanasios Kottas and Bruno Sansó Abstract
Nonparametric Bayesian modeling and inference for renewal processes Sai Xiao, Athanasios Kottas and Bruno Sansó Abstract We propose a flexible approach to modeling and inference for renewal processes.
More informationSimulation of truncated normal variables. Christian P. Robert LSTA, Université Pierre et Marie Curie, Paris
Simulation of truncated normal variables Christian P. Robert LSTA, Université Pierre et Marie Curie, Paris Abstract arxiv:0907.4010v1 [stat.co] 23 Jul 2009 We provide in this paper simulation algorithms
More informationMultivariate Hawkes Processes and Their Simulations
Multivariate Hawkes Processes and Their Simulations Yuanda Chen September, 2016 Abstract In this article we will extend our discussion to the multivariate Hawkes processes, which are mutually exciting
More informationContinuous Time Processes
page 102 Chapter 7 Continuous Time Processes 7.1 Introduction In a continuous time stochastic process (with discrete state space), a change of state can occur at any time instant. The associated point
More informationInference for non-stationary, non-gaussian, Irregularly-Sampled Processes
Inference for non-stationary, non-gaussian, Irregularly-Sampled Processes Robert Wolpert Duke Department of Statistical Science wolpert@stat.duke.edu SAMSI ASTRO Opening WS Aug 22-26 2016 Aug 24 Hamner
More informationBayesian Methods with Monte Carlo Markov Chains II
Bayesian Methods with Monte Carlo Markov Chains II Henry Horng-Shing Lu Institute of Statistics National Chiao Tung University hslu@stat.nctu.edu.tw http://tigpbp.iis.sinica.edu.tw/courses.htm 1 Part 3
More informationStat 535 C - Statistical Computing & Monte Carlo Methods. Arnaud Doucet.
Stat 535 C - Statistical Computing & Monte Carlo Methods Arnaud Doucet Email: arnaud@cs.ubc.ca 1 1.1 Outline Introduction to Markov chain Monte Carlo The Gibbs Sampler Examples Overview of the Lecture
More informationHawkes Processes and their Applications in Finance and Insurance
Hawkes Processes and their Applications in Finance and Insurance Anatoliy Swishchuk University of Calgary Calgary, Alberta, Canada Hawks Seminar Talk Dept. of Math. & Stat. Calgary, Canada May 9th, 2018
More informationStat 451 Lecture Notes Markov Chain Monte Carlo. Ryan Martin UIC
Stat 451 Lecture Notes 07 12 Markov Chain Monte Carlo Ryan Martin UIC www.math.uic.edu/~rgmartin 1 Based on Chapters 8 9 in Givens & Hoeting, Chapters 25 27 in Lange 2 Updated: April 4, 2016 1 / 42 Outline
More informationerrors every 1 hour unless he falls asleep, in which case he just reports the total errors
I. First Definition of a Poisson Process A. Definition: Poisson process A Poisson Process {X(t), t 0} with intensity λ > 0 is a counting process with the following properties. INDEPENDENT INCREMENTS. For
More informationData analysis and stochastic modeling
Data analysis and stochastic modeling Lecture 7 An introduction to queueing theory Guillaume Gravier guillaume.gravier@irisa.fr with a lot of help from Paul Jensen s course http://www.me.utexas.edu/ jensen/ormm/instruction/powerpoint/or_models_09/14_queuing.ppt
More informationApproximate Bayesian Computation and Particle Filters
Approximate Bayesian Computation and Particle Filters Dennis Prangle Reading University 5th February 2014 Introduction Talk is mostly a literature review A few comments on my own ongoing research See Jasra
More informationA Bayesian Approach to Phylogenetics
A Bayesian Approach to Phylogenetics Niklas Wahlberg Based largely on slides by Paul Lewis (www.eeb.uconn.edu) An Introduction to Bayesian Phylogenetics Bayesian inference in general Markov chain Monte
More informationStatistics 222, Spatial Statistics. Outline for the day: 1. Problems and code from last lecture. 2. Likelihood. 3. MLE. 4. Simulation.
Statistics 222, Spatial Statistics. Outline for the day: 1. Problems and code from last lecture. 2. Likelihood. 3. MLE. 4. Simulation. 1. Questions and code from last time. The difference between ETAS
More informationPetr Volf. Model for Difference of Two Series of Poisson-like Count Data
Petr Volf Institute of Information Theory and Automation Academy of Sciences of the Czech Republic Pod vodárenskou věží 4, 182 8 Praha 8 e-mail: volf@utia.cas.cz Model for Difference of Two Series of Poisson-like
More informationSystem Simulation Part II: Mathematical and Statistical Models Chapter 5: Statistical Models
System Simulation Part II: Mathematical and Statistical Models Chapter 5: Statistical Models Fatih Cavdur fatihcavdur@uludag.edu.tr March 20, 2012 Introduction Introduction The world of the model-builder
More informationLecture 4a: Continuous-Time Markov Chain Models
Lecture 4a: Continuous-Time Markov Chain Models Continuous-time Markov chains are stochastic processes whose time is continuous, t [0, ), but the random variables are discrete. Prominent examples of continuous-time
More informationSystem Simulation Part II: Mathematical and Statistical Models Chapter 5: Statistical Models
System Simulation Part II: Mathematical and Statistical Models Chapter 5: Statistical Models Fatih Cavdur fatihcavdur@uludag.edu.tr March 29, 2014 Introduction Introduction The world of the model-builder
More informationReading: Karlin and Taylor Ch. 5 Resnick Ch. 3. A renewal process is a generalization of the Poisson point process.
Renewal Processes Wednesday, December 16, 2015 1:02 PM Reading: Karlin and Taylor Ch. 5 Resnick Ch. 3 A renewal process is a generalization of the Poisson point process. The Poisson point process is completely
More informationModeling Recurrent Events in Panel Data Using Mixed Poisson Models
Modeling Recurrent Events in Panel Data Using Mixed Poisson Models V. Savani and A. Zhigljavsky Abstract This paper reviews the applicability of the mixed Poisson process as a model for recurrent events
More informationBrief introduction to Markov Chain Monte Carlo
Brief introduction to Department of Probability and Mathematical Statistics seminar Stochastic modeling in economics and finance November 7, 2011 Brief introduction to Content 1 and motivation Classical
More informationNonparametric Bayesian Methods (Gaussian Processes)
[70240413 Statistical Machine Learning, Spring, 2015] Nonparametric Bayesian Methods (Gaussian Processes) Jun Zhu dcszj@mail.tsinghua.edu.cn http://bigml.cs.tsinghua.edu.cn/~jun State Key Lab of Intelligent
More informationLearning with Temporal Point Processes
Learning with Temporal Point Processes t Manuel Gomez Rodriguez MPI for Software Systems Isabel Valera MPI for Intelligent Systems Slides/references: http://learning.mpi-sws.org/tpp-icml18 ICML TUTORIAL,
More informationComputational statistics
Computational statistics Markov Chain Monte Carlo methods Thierry Denœux March 2017 Thierry Denœux Computational statistics March 2017 1 / 71 Contents of this chapter When a target density f can be evaluated
More informationOptional Stopping Theorem Let X be a martingale and T be a stopping time such
Plan Counting, Renewal, and Point Processes 0. Finish FDR Example 1. The Basic Renewal Process 2. The Poisson Process Revisited 3. Variants and Extensions 4. Point Processes Reading: G&S: 7.1 7.3, 7.10
More informationProbabilistic Graphical Models Lecture 17: Markov chain Monte Carlo
Probabilistic Graphical Models Lecture 17: Markov chain Monte Carlo Andrew Gordon Wilson www.cs.cmu.edu/~andrewgw Carnegie Mellon University March 18, 2015 1 / 45 Resources and Attribution Image credits,
More informationRecap. Probability, stochastic processes, Markov chains. ELEC-C7210 Modeling and analysis of communication networks
Recap Probability, stochastic processes, Markov chains ELEC-C7210 Modeling and analysis of communication networks 1 Recap: Probability theory important distributions Discrete distributions Geometric distribution
More informationSome Applications of Hawkes Processes for Order Book Modelling
Some Applications of Hawkes Processes for Order Book Modelling Ioane Muni Toke Ecole Centrale Paris Chair of Quantitative Finance First Workshop on Quantitative Finance and Economics International Christian
More informationMCMC algorithms for fitting Bayesian models
MCMC algorithms for fitting Bayesian models p. 1/1 MCMC algorithms for fitting Bayesian models Sudipto Banerjee sudiptob@biostat.umn.edu University of Minnesota MCMC algorithms for fitting Bayesian models
More information6 Markov Chain Monte Carlo (MCMC)
6 Markov Chain Monte Carlo (MCMC) The underlying idea in MCMC is to replace the iid samples of basic MC methods, with dependent samples from an ergodic Markov chain, whose limiting (stationary) distribution
More information1 Probability and Random Variables
1 Probability and Random Variables The models that you have seen thus far are deterministic models. For any time t, there is a unique solution X(t). On the other hand, stochastic models will result in
More informationGlossary. The ISI glossary of statistical terms provides definitions in a number of different languages:
Glossary The ISI glossary of statistical terms provides definitions in a number of different languages: http://isi.cbs.nl/glossary/index.htm Adjusted r 2 Adjusted R squared measures the proportion of the
More informationQ = (c) Assuming that Ricoh has been working continuously for 7 days, what is the probability that it will remain working at least 8 more days?
IEOR 4106: Introduction to Operations Research: Stochastic Models Spring 2005, Professor Whitt, Second Midterm Exam Chapters 5-6 in Ross, Thursday, March 31, 11:00am-1:00pm Open Book: but only the Ross
More informationModeling and Performance Analysis with Discrete-Event Simulation
Simulation Modeling and Performance Analysis with Discrete-Event Simulation Chapter 9 Input Modeling Contents Data Collection Identifying the Distribution with Data Parameter Estimation Goodness-of-Fit
More informationPreliminary Statistics Lecture 2: Probability Theory (Outline) prelimsoas.webs.com
1 School of Oriental and African Studies September 2015 Department of Economics Preliminary Statistics Lecture 2: Probability Theory (Outline) prelimsoas.webs.com Gujarati D. Basic Econometrics, Appendix
More informationStat 426 : Homework 1.
Stat 426 : Homework 1. Moulinath Banerjee University of Michigan Announcement: The homework carries 120 points and contributes 10 points to the total grade. (1) A geometric random variable W takes values
More informationBayesian Inference for Discretely Sampled Diffusion Processes: A New MCMC Based Approach to Inference
Bayesian Inference for Discretely Sampled Diffusion Processes: A New MCMC Based Approach to Inference Osnat Stramer 1 and Matthew Bognar 1 Department of Statistics and Actuarial Science, University of
More informationLecture 2 APPLICATION OF EXREME VALUE THEORY TO CLIMATE CHANGE. Rick Katz
1 Lecture 2 APPLICATION OF EXREME VALUE THEORY TO CLIMATE CHANGE Rick Katz Institute for Study of Society and Environment National Center for Atmospheric Research Boulder, CO USA email: rwk@ucar.edu Home
More informationMarkov Chain Monte Carlo methods
Markov Chain Monte Carlo methods Tomas McKelvey and Lennart Svensson Signal Processing Group Department of Signals and Systems Chalmers University of Technology, Sweden November 26, 2012 Today s learning
More informationBayesian Regression Linear and Logistic Regression
When we want more than point estimates Bayesian Regression Linear and Logistic Regression Nicole Beckage Ordinary Least Squares Regression and Lasso Regression return only point estimates But what if we
More informationModels for Multivariate Panel Count Data
Semiparametric Models for Multivariate Panel Count Data KyungMann Kim University of Wisconsin-Madison kmkim@biostat.wisc.edu 2 April 2015 Outline 1 Introduction 2 3 4 Panel Count Data Motivation Previous
More informationProbability Methods in Civil Engineering Prof. Rajib Maity Department of Civil Engineering Indian Institute of Technology, Kharagpur
Probability Methods in Civil Engineering Prof. Rajib Maity Department of Civil Engineering Indian Institute of Technology, Kharagpur Lecture No. # 12 Probability Distribution of Continuous RVs (Contd.)
More informationTractable Nonparametric Bayesian Inference in Poisson Processes with Gaussian Process Intensities
Tractable Nonparametric Bayesian Inference in Poisson Processes with Gaussian Process Intensities Ryan Prescott Adams rpa23@cam.ac.uk Cavendish Laboratory, University of Cambridge, Cambridge CB3 HE, UK
More informationBayesian Inference in GLMs. Frequentists typically base inferences on MLEs, asymptotic confidence
Bayesian Inference in GLMs Frequentists typically base inferences on MLEs, asymptotic confidence limits, and log-likelihood ratio tests Bayesians base inferences on the posterior distribution of the unknowns
More informationCDA5530: Performance Models of Computers and Networks. Chapter 3: Review of Practical
CDA5530: Performance Models of Computers and Networks Chapter 3: Review of Practical Stochastic Processes Definition Stochastic ti process X = {X(t), t T} is a collection of random variables (rvs); one
More informationFigure 10.1: Recording when the event E occurs
10 Poisson Processes Let T R be an interval. A family of random variables {X(t) ; t T} is called a continuous time stochastic process. We often consider T = [0, 1] and T = [0, ). As X(t) is a random variable
More informationExample: physical systems. If the state space. Example: speech recognition. Context can be. Example: epidemics. Suppose each infected
4. Markov Chains A discrete time process {X n,n = 0,1,2,...} with discrete state space X n {0,1,2,...} is a Markov chain if it has the Markov property: P[X n+1 =j X n =i,x n 1 =i n 1,...,X 0 =i 0 ] = P[X
More informationBivariate Degradation Modeling Based on Gamma Process
Bivariate Degradation Modeling Based on Gamma Process Jinglun Zhou Zhengqiang Pan Member IAENG and Quan Sun Abstract Many highly reliable products have two or more performance characteristics (PCs). The
More informationStatistics of Stochastic Processes
Prof. Dr. J. Franke All of Statistics 4.1 Statistics of Stochastic Processes discrete time: sequence of r.v...., X 1, X 0, X 1, X 2,... X t R d in general. Here: d = 1. continuous time: random function
More informationDynamic models. Dependent data The AR(p) model The MA(q) model Hidden Markov models. 6 Dynamic models
6 Dependent data The AR(p) model The MA(q) model Hidden Markov models Dependent data Dependent data Huge portion of real-life data involving dependent datapoints Example (Capture-recapture) capture histories
More information13: Variational inference II
10-708: Probabilistic Graphical Models, Spring 2015 13: Variational inference II Lecturer: Eric P. Xing Scribes: Ronghuo Zheng, Zhiting Hu, Yuntian Deng 1 Introduction We started to talk about variational
More informationSTAT 302 Introduction to Probability Learning Outcomes. Textbook: A First Course in Probability by Sheldon Ross, 8 th ed.
STAT 302 Introduction to Probability Learning Outcomes Textbook: A First Course in Probability by Sheldon Ross, 8 th ed. Chapter 1: Combinatorial Analysis Demonstrate the ability to solve combinatorial
More informationTime-varying failure rate for system reliability analysis in large-scale railway risk assessment simulation
Time-varying failure rate for system reliability analysis in large-scale railway risk assessment simulation H. Zhang, E. Cutright & T. Giras Center of Rail Safety-Critical Excellence, University of Virginia,
More informationSpectra of some self-exciting and mutually exciting point processes
Biometrika (1971), 58, 1, p. 83 83 Printed in Great Britain Spectra of some self-exciting and mutually exciting point processes BY ALAN G. HAWKES University of Durham SUMMAKY In recent years methods of
More informationAre financial markets becoming systemically more unstable?
Are financial markets becoming systemically more unstable? Lucio Maria Calcagnile joint work with Giacomo Bormetti, Michele Treccani, Stefano Marmi, Fabrizio Lillo www.list-group.com www.quantlab.it Parma,
More information