Non-Parametric Bayesian Inference for Controlled Branching Processes Through MCMC Methods
|
|
- Charleen Harris
- 5 years ago
- Views:
Transcription
1 Non-Parametric Bayesian Inference for Controlled Branching Processes Through MCMC Methods M. González, R. Martínez, I. del Puerto, A. Ramos Department of Mathematics. University of Extremadura Spanish Branching Processes Group 18th International Conference in Computational Statistics Porto, August 2008 M. González
2 Contents 1 Controlled Branching Processes 2 3 Introducing the Method Developing the Method Simulated Example 4 Concluding Remarks References M. González
3 Contents 1 Controlled Branching Processes 2 3 Introducing the Method Developing the Method Simulated Example 4 Concluding Remarks References M. González
4 Contents 1 Controlled Branching Processes 2 3 Introducing the Method Developing the Method Simulated Example 4 Concluding Remarks References M. González
5 Contents 1 Controlled Branching Processes 2 3 Introducing the Method Developing the Method Simulated Example 4 Concluding Remarks References M. González
6 Branching Processes Inside the general context concerning Stochastic Models, Branching Processes Theory provides appropriate mathematical models for description of the probabilistic evolution of systems whose components (cell, particles, individuals in general), after certain life period, reproduce and die. Therefore, it can be applied in several fields (biology, demography, ecology, epidemiology, genetics, medicine,...). M. González
7 Branching Processes Example Z 0 = 1. Z n Z n+1 = j=1 X nj M. González
8 Branching Processes Example Z 0 = 1. Z n Z n+1 = j=1 X nj M. González
9 Branching Processes Example Z 0 = 1 Z 1 = 2. Z n Z n+1 = j=1 X nj M. González
10 Branching Processes Example Z 0 = 1 Z 1 = 2. Z n Z n+1 = j=1 X nj M. González
11 Branching Processes Example Z 0 = 1 Z 1 = 2 Z 2 = 7. Z n Z n+1 = j=1 X nj M. González
12 Branching Processes Example Z 0 = 1 Z 1 = 2 Z 2 = 7. Z n Z n+1 = j=1 X nj M. González
13 Branching Processes Example Z 0 = 1 Z 1 = 2 Z 2 = 7 Z 3 = 10. Z n Z n+1 = j=1 X nj M. González
14 Branching Processes Example Z 0 = 1 Z 1 = 2 Z 2 = 7 Z 3 = 10. Z n Z n+1 = j=1 X nj M. González
15 Branching Processes Example Z 0 = 1 Z 1 = 2 Z 2 = 7 Z 3 = 10. Z n Z n+1 = j=1 X nj M. González
16 Branching Processes Main Results for Galton Watson Branching Processes Let m = E[X 01 ] and σ 2 = Var[X 01 ] Extinction Problem If m 1 the process dies out with probability 1 If m > 1 there exists a positive probability of non-extinction Asymptotic behaviour Statistical Inference M. González
17 Branching Processes Main Results for Galton Watson Branching Processes Let m = E[X 01 ] and σ 2 = Var[X 01 ] Extinction Problem If m 1 the process dies out with probability 1 If m > 1 there exists a positive probability of non-extinction Asymptotic behaviour Statistical Inference M. González
18 Branching Processes Main Results for Galton Watson Branching Processes Let m = E[X 01 ] and σ 2 = Var[X 01 ] Extinction Problem If m 1 the process dies out with probability 1 If m > 1 there exists a positive probability of non-extinction Asymptotic behaviour Statistical Inference M. González
19 Branching Processes Main Results for Galton Watson Branching Processes Let m = E[X 01 ] and σ 2 = Var[X 01 ] Extinction Problem If m 1 the process dies out with probability 1 If m > 1 there exists a positive probability of non-extinction Asymptotic behaviour Statistical Inference M. González
20 Branching Processes Many monographs about the theory and applications about the branching processes have been published: Harris, T. (1963). The Theory of branching processes. Springer-Verlag. Jagers, P. (1975) Branching processes with Biological Applications, John Wiley and Sons, Inc. Asmussen, S. and Hering, H. (1983). Branching processes. Birkhäuser. Boston. Athreya, K.B. and Jagers, P. (1997). Classical and modern branching processes. Springer-Verlag. Kimmel, M. and Axelrod, D.E. (2002) Branching processes in Biology, Springer-Verlag New York, Inc. Haccou, P., Jagers, P., and Vatutin, V. (2005) Branching Processes: Variation, Growth, and Extinction of Populations. Cambridge University Press. M. González
21 Branching Processes A Controlled Branching Process is a discrete-time stochastic growth population model in which the individuals with reproductive capacity in each generation are controlled by some function φ. This branching model is well-suited for describing the probabilistic evolution of populations in which, for various reasons of an environmental, social or other nature, there is a mechanism that establishes the number of progenitors who take part in each generation. M. González
22 Branching Processes Mathematically: Controlled Branching Process Z 0 = N, Z n+1 = {Z n } n 0 φ(z n) i=1 X ni, n = 0, 1,... {X ni : i = 1, 2,..., n = 0, 1,...} are i.i.d. random variables. {p k : k S} Offspring Distribution m = E[X 01 ], σ 2 = Var[X 01 ] φ : R + R + is assumed to be integer-valued for integer-valued arguments Control Function M. González
23 Controlled Branching Processes Properties {Z n } n 0 is a Homogeneous Markov Chain Duality Extinction-Explosion: P(Z n 0) + P(Z n ) = 1 Main Topics Investigated Extinction Problem Sevast yanov and Zubkov (1974) Zubkov (1974) Molina, González and Mota (1998) Asymptotic Behaviour: Growth rates Bagley (1986) Molina, González and Mota (1998) González, Molina, del Puerto (2002, 2003, 2004, 2005a,b) M. González
24 Controlled Branching Processes Properties {Z n } n 0 is a Homogeneous Markov Chain Duality Extinction-Explosion: P(Z n 0) + P(Z n ) = 1 Main Topics Investigated Extinction Problem Sevast yanov and Zubkov (1974) Zubkov (1974) Molina, González and Mota (1998) Asymptotic Behaviour: Growth rates Bagley (1986) Molina, González and Mota (1998) González, Molina, del Puerto (2002, 2003, 2004, 2005a,b) M. González
25 Controlled Branching Processes Main Topics Investigated Statistical Inference Dion, J. P. and Essebbar, B. (1995). On the statistics of controlled branching processes. Lecture Notes in Statistics, 99: M. González, R. Martínez, I. Del Puerto (2004). Nonparametric estimation of the offspring distribution and the mean for a controlled branching process. Test, 13(2), M. González, R. Martínez, I. Del Puerto (2005). Estimation of the variance for a controlled branching process. Test, 14(1), T.N. Sriram, A. Bhattacharya, M. González, R. Martínez, I. Del Puerto (2007). Estimation of the offspring mean in a controlled branching process with a random control function. Stochastic Processes and their Applications, 117, M. González
26 Non-Parametric Framework Offspring Distribution: p = {p k : k S} S finite. Sample: The entire family tree up to the current generation or at least {X ki : i = 1,..., φ(z k ), k = 0, 1,..., n} Z n = {Z j (k): k S, j = 0,..., n} where Z j (k) = φ(z j ) i=1 I {X ji =k} = number of parents in the jth-generation which generate exactly k offspring Objetive: Make inference on p M. González
27 Non-Parametric Framework Offspring Distribution: p = {p k : k S} S finite. Sample: The entire family tree up to the current generation or at least {X ki : i = 1,..., φ(z k ), k = 0, 1,..., n} Z n = {Z j (k): k S, j = 0,..., n} where Z j (k) = φ(z j ) i=1 I {X ji =k} = number of parents in the jth-generation which generate exactly k offspring Objetive: Make inference on p M. González
28 Likelihood Function L(p Z n ) n j=0 p Z j(k) k k S Conjugate Class of Distributions: Dirichlet Family Prior Distribution: p D(α k : k S) Posterior Distribution: n p Z n D(α k + Z j (k): k S) j=0 M. González
29 Likelihood Function L(p Z n ) n j=0 p Z j(k) k k S Conjugate Class of Distributions: Dirichlet Family Prior Distribution: p D(α k : k S) Posterior Distribution: n p Z n D(α k + Z j (k): k S) j=0 M. González
30 Setting out the Problem In real problems it is difficult to observe the entire family tree {X ki : i = 1, 2,..., k = 0, 1,..., n} or even the random variables Z n = {Z j (k): k S, j = 0,..., n} Usual Sample Information Z n = {Z j : j = 0,..., n} A Solution We introduce an algorithm to approximate the distribution p Z n using Markov Chain Monte Carlo Methods M. González
31 Setting out the Problem In real problems it is difficult to observe the entire family tree {X ki : i = 1, 2,..., k = 0, 1,..., n} or even the random variables Z n = {Z j (k): k S, j = 0,..., n} Usual Sample Information Z n = {Z j : j = 0,..., n} A Solution We introduce an algorithm to approximate the distribution p Z n using Markov Chain Monte Carlo Methods M. González
32 Setting out the Problem In real problems it is difficult to observe the entire family tree {X ki : i = 1, 2,..., k = 0, 1,..., n} or even the random variables Z n = {Z j (k): k S, j = 0,..., n} Usual Sample Information Z n = {Z j : j = 0,..., n} A Solution We introduce an algorithm to approximate the distribution p Z n using Markov Chain Monte Carlo Methods M. González
33 Introducing the Method Developing the Method Simulated Example Gibbs Sampler: Introducing the Method Sample: Z n = {Z j : j = 0,..., n} The Problem Latent Variables: p Z n Z n = {Z j (k): k S, j = 0,..., n} Gibbs Sampler: p Z n, Z n Z n Z n, p M. González
34 Introducing the Method Developing the Method Simulated Example Gibbs Sampler: Introducing the Method Sample: Z n = {Z j : j = 0,..., n} The Problem Latent Variables: p Z n Z n = {Z j (k): k S, j = 0,..., n} Gibbs Sampler: p Z n, Z n Z n Z n, p M. González
35 Introducing the Method Developing the Method Simulated Example Gibbs Sampler: Introducing the Method Sample: Z n = {Z j : j = 0,..., n} The Problem Latent Variables: p Z n Z n = {Z j (k): k S, j = 0,..., n} Gibbs Sampler: p Z n, Z n Z n Z n, p M. González
36 Introducing the Method Developing the Method Simulated Example Gibbs Sampler: Introducing the Method Sample: Z n = {Z j : j = 0,..., n} The Problem Latent Variables: p Z n Z n = {Z j (k): k S, j = 0,..., n} Gibbs Sampler: p Z n, Z n Z n Z n, p M. González
37 Introducing the Method Developing the Method Simulated Example Gibbs Sampler: Introducing the Method First Conditional Distribution: p Z n, Zn n p Z n, Zn p Z n D(α k + Z j (k): k S) j=0 For j = 0,..., n φ(z j ) = k S Z j (k) Z j+1 = k S kz j (k) M. González
38 Introducing the Method Developing the Method Simulated Example Gibbs Sampler: Introducing the Method First Conditional Distribution: p Z n, Zn n p Z n, Zn p Z n D(α k + Z j (k): k S) j=0 For j = 0,..., n φ(z j ) = k S Z j (k) Z j+1 = k S kz j (k) M. González
39 Introducing the Method Developing the Method Simulated Example Gibbs Sampler: Introducing the Method First Conditional Distribution: p Z n, Zn n p Z n, Zn p Z n D(α k + Z j (k): k S) j=0 For j = 0,..., n φ(z j ) = k S Z j (k) Z j+1 = k S kz j (k) M. González
40 Introducing the Method Developing the Method Simulated Example Gibbs Sampler: Introducing the Method Second Conditional Distribution: Z n Zn, p n P(Z n Zn, p) = P(Z j (k): k S Z j, Z j+1, p) j=0 (Z j (k): k S) Z j, Z j+1, p is obtained from a Multinomial(φ(Z j ), p) normalized by considering the constraint Z j+1 = kz j (k) k S M. González
41 Introducing the Method Developing the Method Simulated Example Gibbs Sampler: Introducing the Method Second Conditional Distribution: Z n Zn, p n P(Z n Zn, p) = P(Z j (k): k S Z j, Z j+1, p) j=0 (Z j (k): k S) Z j, Z j+1, p is obtained from a Multinomial(φ(Z j ), p) normalized by considering the constraint Z j+1 = kz j (k) k S M. González
42 Introducing the Method Developing the Method Simulated Example Gibbs Sampler: Introducing the Method Second Conditional Distribution: Z n Z n, p p φ(z 0 ) Z 0 (k), k S Z 1 φ(z 1 ) Z 1 (k), k S Z 2 φ(z 2 )... Z n φ(z n ) Z n (k), k S Z n+1 φ(z j ) = k S Z j (k), Z j+1 = k S kz j (k) M. González
43 Introducing the Method Developing the Method Simulated Example Gibbs Sampler: Introducing the Method Second Conditional Distribution: Z n Z n, p p φ(z 0 ) Z 0 (k), k S Z 1 φ(z 1 ) Z 1 (k), k S Z 2 φ(z 2 )... Z n φ(z n ) Z n (k), k S Z n+1 φ(z j ) = k S Z j (k), Z j+1 = k S kz j (k) M. González
44 Introducing the Method Developing the Method Simulated Example Gibbs Sampler: Introducing the Method Second Conditional Distribution: Z n Z n, p p φ(z 0 ) Z 0 (k), k S Z 1 φ(z 1 ) Z 1 (k), k S Z 2 φ(z 2 )... Z n φ(z n ) Z n (k), k S Z n+1 φ(z j ) = k S Z j (k), Z j+1 = k S kz j (k) M. González
45 Introducing the Method Developing the Method Simulated Example Gibbs Sampler: Introducing the Method Second Conditional Distribution: Z n Z n, p p φ(z 0 ) Z 0 (k), k S Z 1 φ(z 1 ) Z 1 (k), k S Z 2 φ(z 2 )... Z n φ(z n ) Z n (k), k S Z n+1 φ(z j ) = k S Z j (k), Z j+1 = k S kz j (k) M. González
46 Introducing the Method Developing the Method Simulated Example Gibbs Sampler: Developing the Method Algorithm Fixed p (0) Do l = 1 Generate Z n (l) Z n Zn, p (l 1) Generate p (l) p Z n (l) Do l = l + 1 For a run of the sequence {p (l) } l 0, we choose Q + 1 vectors in the way {p (N), p (N+G),..., p (N+QG) }, where G is a batch size. The vectors {p (N), p (N+G),..., p (N+QG) } are considered independent samples from p Z n if G and N are large enough (Tierney (1994)). Since these vectors could be affected by the initial state p (0), we apply the algorithm T times, obtaining a final sample of length T(Q + 1). M. González
47 Introducing the Method Developing the Method Simulated Example Gibbs Sampler: Developing the Method Algorithm Fixed p (0) Do l = 1 Generate Z n (l) Z n Zn, p (l 1) Generate p (l) p Z n (l) Do l = l + 1 For a run of the sequence {p (l) } l 0, we choose Q + 1 vectors in the way {p (N), p (N+G),..., p (N+QG) }, where G is a batch size. The vectors {p (N), p (N+G),..., p (N+QG) } are considered independent samples from p Z n if G and N are large enough (Tierney (1994)). Since these vectors could be affected by the initial state p (0), we apply the algorithm T times, obtaining a final sample of length T(Q + 1). M. González
48 Introducing the Method Developing the Method Simulated Example Gibbs Sampler: Developing the Method Algorithm Fixed p (0) Do l = 1 Generate Z n (l) Z n Zn, p (l 1) Generate p (l) p Z n (l) Do l = l + 1 For a run of the sequence {p (l) } l 0, we choose Q + 1 vectors in the way {p (N), p (N+G),..., p (N+QG) }, where G is a batch size. The vectors {p (N), p (N+G),..., p (N+QG) } are considered independent samples from p Z n if G and N are large enough (Tierney (1994)). Since these vectors could be affected by the initial state p (0), we apply the algorithm T times, obtaining a final sample of length T(Q + 1). M. González
49 Introducing the Method Developing the Method Simulated Example Gibbs Sampler: Developing the Method Algorithm Fixed p (0) Do l = 1 Generate Z n (l) Z n Zn, p (l 1) Generate p (l) p Z n (l) Do l = l + 1 For a run of the sequence {p (l) } l 0, we choose Q + 1 vectors in the way {p (N), p (N+G),..., p (N+QG) }, where G is a batch size. The vectors {p (N), p (N+G),..., p (N+QG) } are considered independent samples from p Z n if G and N are large enough (Tierney (1994)). Since these vectors could be affected by the initial state p (0), we apply the algorithm T times, obtaining a final sample of length T(Q + 1). M. González
50 Introducing the Method Developing the Method Simulated Example Gibbs Sampler: Simulated Example Offspring Distribution: k p k Parameters: m = 1.08, σ 2 = Control function: φ(x) = 7 if x 7; x if 7 < x 20; 20 if x > 20 Simulated Data Controlled Branching Process Galton Watson Branching Process M. González
51 Introducing the Method Developing the Method Simulated Example Gibbs Sampler: Simulated Example Observed Data: n = 15 n Z n n Z n p D(1/2,..., 1/2) Selection of N, G, Q and T N = 5000, G = 100, Q = 49 and T = 100 Gelman-Rubin-Brooks diagnostic plots. Estimated potential scale reduction factor. Autocorrelation values. M. González
52 Introducing the Method Developing the Method Simulated Example Gibbs Sampler: Simulated Example Gelman-Rubin-Brooks diagnostic plots (CODA package for R) p 0 p 1 p 2 shrink factor median 97.5% shrink factor median 97.5% shrink factor median 97.5% last iteration in chain last iteration in chain last iteration in chain p 3 p 4 shrink factor median 97.5% shrink factor median 97.5% last iteration in chain last iteration in chain M. González
53 Introducing the Method Developing the Method Simulated Example Gibbs Sampler: Simulated Examples p D(1/2,..., 1/2) Selection of N, G, Q and T N = 5000, G = 100, Q = 49 and T = 100 Gelman-Rubin-Brooks diagnostic plots. Estimated potential scale reduction factor. Autocorrelation values. Computation Time: 60.10s for each chain on an Intel(R) Core(TM)2 Duo CPU T7500 running at 2.20GHz with 2038 MB RAM. M. González
54 Introducing the Method Developing the Method Simulated Example Gibbs Sampler: Simulated Example Sample Information: Z n N = 5000, G = 100, Q = 49 and T = 100 (Sample Size: 5000) p 0 hpd 95% p^0 p 0 hpd 95% p 3 hpd 95% p^3 p 0 hpd 95% Offspring Mean hpd 95% mm^ hpd 95% Density Density Density Algorithm s Efficiency MEAN SD MCSE TSSE M. González
55 Concluding Remarks Concluding Remarks References In a non-parametric Bayesian framework we can make inference on the offspring distribution of Controlled Branching Processes, and consequently on the rest of offspring parameters, without observing the entire family tree, but only considering the total number of individuals in each generation. We use a MCMC method (Gibbs sampler) and the statistical software and programming environment R, in order to give a "likely" approach to family trees. M. González
56 Concluding Remarks Concluding Remarks References In a non-parametric Bayesian framework we can make inference on the offspring distribution of Controlled Branching Processes, and consequently on the rest of offspring parameters, without observing the entire family tree, but only considering the total number of individuals in each generation. We use a MCMC method (Gibbs sampler) and the statistical software and programming environment R, in order to give a "likely" approach to family trees. M. González
57 References Concluding Remarks References Bagley, J. H. (1986). On the almost sure convergence of controlled branching processes. Journal of Applied Probability, 23: M. González, M. Molina, I. del Puerto (2002). On the class of controlled branching process with random control functions. Journal of Applied Probability, 39 (4), M. González, M. Molina, I. del Puerto (2003). On the geometric growth in controlled branching processes with random control function. Journal of Applied Probability, 40(4), M. González, M. Molina, I. del Puerto (2004). Limiting distribution for subcritical controlled branching processes with random control function. Statistics and Probability Letters, 67(3), M. González, M. Molina, I. del Puerto (2005a). Asymptotic behaviour of critical controlled branching process with random control function. Journal of Applied Probability, 42(2), M. González, M. Molina, I. del Puerto (2005b). On the L2-convergence of controlled branching processes with random control function. Bernoulli, 11(1), M. Molina, M. González, M. Mota (1998). Some theoretical results about superadditive controlled Galton-Watson branching processes. Proceedings of the International Conference Prague Stochastics98, 2, Sevast yanov, B. A. and Zubkov, A. (1974). Controlled branching processes. Theory of Probability and its Applications, 19: Tierney, L. (1994). Markov chains for exploring posterior distributions. Annals of Statistics, 22, Zubkov, A. M. (1974). Analogies between Galton-Watson processes and φ-branching processes. Theory of Probability and its Applications,19: M. González
58 Concluding Remarks References Thank you very much! M. González
RECENT RESULTS FOR SUPERCRITICAL CONTROLLED BRANCHING PROCESSES WITH CONTROL RANDOM FUNCTIONS
Pliska Stud. Math. Bulgar. 16 (2004), 43-54 STUDIA MATHEMATICA BULGARICA RECENT RESULTS FOR SUPERCRITICAL CONTROLLED BRANCHING PROCESSES WITH CONTROL RANDOM FUNCTIONS Miguel González, Manuel Molina, Inés
More informationESTIMATION OF THE OFFSPRING AND IMMIGRATION MEAN VECTORS FOR BISEXUAL BRANCHING PROCESSES WITH IMMIGRATION OF FEMALES AND MALES
Pliska Stud. Math. Bulgar. 20 (2011), 63 80 STUDIA MATHEMATICA BULGARICA ESTIMATION OF THE OFFSPRING AND IMMIGRATION MEAN VECTORS FOR BISEXUAL BRANCHING PROCESSES WITH IMMIGRATION OF FEMALES AND MALES
More informationA Note on Bisexual Galton-Watson Branching Processes with Immigration
E extracta mathematicae Vol. 16, Núm. 3, 361 365 (2001 A Note on Bisexual Galton-Watson Branching Processes with Immigration M. González, M. Molina, M. Mota Departamento de Matemáticas, Universidad de
More informationLimiting distribution for subcritical controlled branching processes with random control function
Available online at www.sciencedirect.com Statistics & Probability Letters 67 (2004 277 284 Limiting distribution for subcritical controlled branching processes with random control function M. Gonzalez,
More informationA NOTE ON THE ASYMPTOTIC BEHAVIOUR OF A PERIODIC MULTITYPE GALTON-WATSON BRANCHING PROCESS. M. González, R. Martínez, M. Mota
Serdica Math. J. 30 (2004), 483 494 A NOTE ON THE ASYMPTOTIC BEHAVIOUR OF A PERIODIC MULTITYPE GALTON-WATSON BRANCHING PROCESS M. González, R. Martínez, M. Mota Communicated by N. M. Yanev Abstract. In
More informationST 740: Markov Chain Monte Carlo
ST 740: Markov Chain Monte Carlo Alyson Wilson Department of Statistics North Carolina State University October 14, 2012 A. Wilson (NCSU Stsatistics) MCMC October 14, 2012 1 / 20 Convergence Diagnostics:
More informationPrerequisite: STATS 7 or STATS 8 or AP90 or (STATS 120A and STATS 120B and STATS 120C). AP90 with a minimum score of 3
University of California, Irvine 2017-2018 1 Statistics (STATS) Courses STATS 5. Seminar in Data Science. 1 Unit. An introduction to the field of Data Science; intended for entering freshman and transfers.
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 informationBayesian Networks in Educational Assessment
Bayesian Networks in Educational Assessment Estimating Parameters with MCMC Bayesian Inference: Expanding Our Context Roy Levy Arizona State University Roy.Levy@asu.edu 2017 Roy Levy MCMC 1 MCMC 2 Posterior
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 informationSTAT 425: Introduction to Bayesian Analysis
STAT 425: Introduction to Bayesian Analysis Marina Vannucci Rice University, USA Fall 2017 Marina Vannucci (Rice University, USA) Bayesian Analysis (Part 2) Fall 2017 1 / 19 Part 2: Markov chain Monte
More informationDiffusion Approximation of Branching Processes
UDC 519.218.23 Diffusion Approximation of Branching Processes N. Limnios Sorbonne University, Université de technologie de Compiègne LMAC Laboratory of Applied Mathematics of Compiègne - CS 60319-60203
More informationNon-Parametric Bayes
Non-Parametric Bayes Mark Schmidt UBC Machine Learning Reading Group January 2016 Current Hot Topics in Machine Learning Bayesian learning includes: Gaussian processes. Approximate inference. Bayesian
More information17 : Markov Chain Monte Carlo
10-708: Probabilistic Graphical Models, Spring 2015 17 : Markov Chain Monte Carlo Lecturer: Eric P. Xing Scribes: Heran Lin, Bin Deng, Yun Huang 1 Review of Monte Carlo Methods 1.1 Overview Monte Carlo
More informationBranching within branching: a general model for host-parasite co-evolution
Branching within branching: a general model for host-parasite co-evolution Gerold Alsmeyer (joint work with Sören Gröttrup) May 15, 2017 Gerold Alsmeyer Host-parasite co-evolution 1 of 26 1 Model 2 The
More informationBayesian Phylogenetics:
Bayesian Phylogenetics: an introduction Marc A. Suchard msuchard@ucla.edu UCLA Who is this man? How sure are you? The one true tree? Methods we ve learned so far try to find a single tree that best describes
More informationPARAMETER ESTIMATION: BAYESIAN APPROACH. These notes summarize the lectures on Bayesian parameter estimation.
PARAMETER ESTIMATION: BAYESIAN APPROACH. These notes summarize the lectures on Bayesian parameter estimation.. Beta Distribution We ll start by learning about the Beta distribution, since we end up using
More informationSequential Monte Carlo Methods for Bayesian Computation
Sequential Monte Carlo Methods for Bayesian Computation A. Doucet Kyoto Sept. 2012 A. Doucet (MLSS Sept. 2012) Sept. 2012 1 / 136 Motivating Example 1: Generic Bayesian Model Let X be a vector parameter
More information16 : Approximate Inference: Markov Chain Monte Carlo
10-708: Probabilistic Graphical Models 10-708, Spring 2017 16 : Approximate Inference: Markov Chain Monte Carlo Lecturer: Eric P. Xing Scribes: Yuan Yang, Chao-Ming Yen 1 Introduction As the target distribution
More informationBayesian Nonparametrics
Bayesian Nonparametrics Peter Orbanz Columbia University PARAMETERS AND PATTERNS Parameters P(X θ) = Probability[data pattern] 3 2 1 0 1 2 3 5 0 5 Inference idea data = underlying pattern + independent
More informationSTA 4273H: Statistical Machine Learning
STA 4273H: Statistical Machine Learning Russ Salakhutdinov Department of Statistics! rsalakhu@utstat.toronto.edu! http://www.utstat.utoronto.ca/~rsalakhu/ Sidney Smith Hall, Room 6002 Lecture 7 Approximate
More informationStatistical Inference for Stochastic Epidemic Models
Statistical Inference for Stochastic Epidemic Models George Streftaris 1 and Gavin J. Gibson 1 1 Department of Actuarial Mathematics & Statistics, Heriot-Watt University, Riccarton, Edinburgh EH14 4AS,
More information2015 IAA Colloquium in OSLO
2015 IAA Colloquium in OSLO The 2014 West African Ebola outbreak as a Branching Process [Working Paper] Guillaume Ominetti 1 Benjamin Schannes 2 June 8, 2015 Working Paper URL upon request 1 SCOR, 5 avenue
More informationNikolay M. YANEV. List of Publications
Nikolay M. YANEV. List of Publications 1. Branching stochastic processes with immigration. Bulletin de l Institut de Mathematiques (Acad. Bulg. Sci.), XV (1972), 71-88. (In Bulgarian). 2. On a class of
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 informationContents. Part I: Fundamentals of Bayesian Inference 1
Contents Preface xiii Part I: Fundamentals of Bayesian Inference 1 1 Probability and inference 3 1.1 The three steps of Bayesian data analysis 3 1.2 General notation for statistical inference 4 1.3 Bayesian
More informationBayesian inference for multivariate extreme value distributions
Bayesian inference for multivariate extreme value distributions Sebastian Engelke Clément Dombry, Marco Oesting Toronto, Fields Institute, May 4th, 2016 Main motivation For a parametric model Z F θ of
More informationMCMC and Gibbs Sampling. Kayhan Batmanghelich
MCMC and Gibbs Sampling Kayhan Batmanghelich 1 Approaches to inference l Exact inference algorithms l l l The elimination algorithm Message-passing algorithm (sum-product, belief propagation) The junction
More informationGaussian Mixture Model
Case Study : Document Retrieval MAP EM, Latent Dirichlet Allocation, Gibbs Sampling Machine Learning/Statistics for Big Data CSE599C/STAT59, University of Washington Emily Fox 0 Emily Fox February 5 th,
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 informationThe Bayesian Choice. Christian P. Robert. From Decision-Theoretic Foundations to Computational Implementation. Second Edition.
Christian P. Robert The Bayesian Choice From Decision-Theoretic Foundations to Computational Implementation Second Edition With 23 Illustrations ^Springer" Contents Preface to the Second Edition Preface
More informationReconstruction of individual patient data for meta analysis via Bayesian approach
Reconstruction of individual patient data for meta analysis via Bayesian approach Yusuke Yamaguchi, Wataru Sakamoto and Shingo Shirahata Graduate School of Engineering Science, Osaka University Masashi
More informationApproximate Bayesian Computation
Approximate Bayesian Computation Michael Gutmann https://sites.google.com/site/michaelgutmann University of Helsinki and Aalto University 1st December 2015 Content Two parts: 1. The basics of approximate
More informationMultilevel Statistical Models: 3 rd edition, 2003 Contents
Multilevel Statistical Models: 3 rd edition, 2003 Contents Preface Acknowledgements Notation Two and three level models. A general classification notation and diagram Glossary Chapter 1 An introduction
More informationBAYESIAN METHODS FOR VARIABLE SELECTION WITH APPLICATIONS TO HIGH-DIMENSIONAL DATA
BAYESIAN METHODS FOR VARIABLE SELECTION WITH APPLICATIONS TO HIGH-DIMENSIONAL DATA Intro: Course Outline and Brief Intro to Marina Vannucci Rice University, USA PASI-CIMAT 04/28-30/2010 Marina Vannucci
More informationMarkov Chain Monte Carlo
Markov Chain Monte Carlo (and Bayesian Mixture Models) David M. Blei Columbia University October 14, 2014 We have discussed probabilistic modeling, and have seen how the posterior distribution is the critical
More informationBayesian Inference. Chapter 1. Introduction and basic concepts
Bayesian Inference Chapter 1. Introduction and basic concepts M. Concepción Ausín Department of Statistics Universidad Carlos III de Madrid Master in Business Administration and Quantitative Methods Master
More informationNONLINEAR APPLICATIONS OF MARKOV CHAIN MONTE CARLO
NONLINEAR APPLICATIONS OF MARKOV CHAIN MONTE CARLO by Gregois Lee, B.Sc.(ANU), B.Sc.Hons(UTas) Submitted in fulfilment of the requirements for the Degree of Doctor of Philosophy Department of Mathematics
More informationPart 8: GLMs and Hierarchical LMs and GLMs
Part 8: GLMs and Hierarchical LMs and GLMs 1 Example: Song sparrow reproductive success Arcese et al., (1992) provide data on a sample from a population of 52 female song sparrows studied over the course
More informationApplicability of subsampling bootstrap methods in Markov chain Monte Carlo
Applicability of subsampling bootstrap methods in Markov chain Monte Carlo James M. Flegal Abstract Markov chain Monte Carlo (MCMC) methods allow exploration of intractable probability distributions by
More information(5) Multi-parameter models - Gibbs sampling. ST440/540: Applied Bayesian Analysis
Summarizing a posterior Given the data and prior the posterior is determined Summarizing the posterior gives parameter estimates, intervals, and hypothesis tests Most of these computations are integrals
More informationMarkov Chain Monte Carlo
Markov Chain Monte Carlo Recall: To compute the expectation E ( h(y ) ) we use the approximation E(h(Y )) 1 n n h(y ) t=1 with Y (1),..., Y (n) h(y). Thus our aim is to sample Y (1),..., Y (n) from f(y).
More informationA note on stochastic context-free grammars, termination and the EM-algorithm
A note on stochastic context-free grammars, termination and the EM-algorithm Niels Richard Hansen Department of Applied Mathematics and Statistics, University of Copenhagen, Universitetsparken 5, 2100
More informationBayesian Inference and MCMC
Bayesian Inference and MCMC Aryan Arbabi Partly based on MCMC slides from CSC412 Fall 2018 1 / 18 Bayesian Inference - Motivation Consider we have a data set D = {x 1,..., x n }. E.g each x i can be the
More informationMarkov chain Monte Carlo
Markov chain Monte Carlo Feng Li feng.li@cufe.edu.cn School of Statistics and Mathematics Central University of Finance and Economics Revised on April 24, 2017 Today we are going to learn... 1 Markov Chains
More informationAsymptotic Behavior of a Controlled Branching Process with Continuous State Space
Asymptotic Behavior of a Controlled Branching Process with Continuous State Space I. Rahimov Department of Mathematical Sciences, KFUPM, Dhahran, Saudi Arabia ABSTRACT In the paper a modification of the
More informationLinear-fractional branching processes with countably many types
Branching processes and and their applications Badajoz, April 11-13, 2012 Serik Sagitov Chalmers University and University of Gothenburg Linear-fractional branching processes with countably many types
More informationeqr094: Hierarchical MCMC for Bayesian System Reliability
eqr094: Hierarchical MCMC for Bayesian System Reliability Alyson G. Wilson Statistical Sciences Group, Los Alamos National Laboratory P.O. Box 1663, MS F600 Los Alamos, NM 87545 USA Phone: 505-667-9167
More informationResistance Growth of Branching Random Networks
Peking University Oct.25, 2018, Chengdu Joint work with Yueyun Hu (U. Paris 13) and Shen Lin (U. Paris 6), supported by NSFC Grant No. 11528101 (2016-2017) for Research Cooperation with Oversea Investigators
More informationBayes: All uncertainty is described using probability.
Bayes: All uncertainty is described using probability. Let w be the data and θ be any unknown quantities. Likelihood. The probability model π(w θ) has θ fixed and w varying. The likelihood L(θ; w) is π(w
More informationA quick introduction to Markov chains and Markov chain Monte Carlo (revised version)
A quick introduction to Markov chains and Markov chain Monte Carlo (revised version) Rasmus Waagepetersen Institute of Mathematical Sciences Aalborg University 1 Introduction These notes are intended to
More informationSampling Methods (11/30/04)
CS281A/Stat241A: Statistical Learning Theory Sampling Methods (11/30/04) Lecturer: Michael I. Jordan Scribe: Jaspal S. Sandhu 1 Gibbs Sampling Figure 1: Undirected and directed graphs, respectively, with
More information27 : Distributed Monte Carlo Markov Chain. 1 Recap of MCMC and Naive Parallel Gibbs Sampling
10-708: Probabilistic Graphical Models 10-708, Spring 2014 27 : Distributed Monte Carlo Markov Chain Lecturer: Eric P. Xing Scribes: Pengtao Xie, Khoa Luu In this scribe, we are going to review the Parallel
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 informationSIMILAR MARKOV CHAINS
SIMILAR MARKOV CHAINS by Phil Pollett The University of Queensland MAIN REFERENCES Convergence of Markov transition probabilities and their spectral properties 1. Vere-Jones, D. Geometric ergodicity in
More informationAsymptotics and Simulation of Heavy-Tailed Processes
Asymptotics and Simulation of Heavy-Tailed Processes Department of Mathematics Stockholm, Sweden Workshop on Heavy-tailed Distributions and Extreme Value Theory ISI Kolkata January 14-17, 2013 Outline
More informationLecturer: David Blei Lecture #3 Scribes: Jordan Boyd-Graber and Francisco Pereira October 1, 2007
COS 597C: Bayesian Nonparametrics Lecturer: David Blei Lecture # Scribes: Jordan Boyd-Graber and Francisco Pereira October, 7 Gibbs Sampling with a DP First, let s recapitulate the model that we re using.
More informationMarkov Chain Monte Carlo in Practice
Markov Chain Monte Carlo in Practice Edited by W.R. Gilks Medical Research Council Biostatistics Unit Cambridge UK S. Richardson French National Institute for Health and Medical Research Vilejuif France
More information4 Branching Processes
4 Branching Processes Organise by generations: Discrete time. If P(no offspring) 0 there is a probability that the process will die out. Let X = number of offspring of an individual p(x) = P(X = x) = offspring
More informationFitting Narrow Emission Lines in X-ray Spectra
Outline Fitting Narrow Emission Lines in X-ray Spectra Taeyoung Park Department of Statistics, University of Pittsburgh October 11, 2007 Outline of Presentation Outline This talk has three components:
More informationModern Discrete Probability Branching processes
Modern Discrete Probability IV - Branching processes Review Sébastien Roch UW Madison Mathematics November 15, 2014 1 Basic definitions 2 3 4 Galton-Watson branching processes I Definition A Galton-Watson
More informationStat 516, Homework 1
Stat 516, Homework 1 Due date: October 7 1. Consider an urn with n distinct balls numbered 1,..., n. We sample balls from the urn with replacement. Let N be the number of draws until we encounter a ball
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 informationMolecular Epidemiology Workshop: Bayesian Data Analysis
Molecular Epidemiology Workshop: Bayesian Data Analysis Jay Taylor and Ananias Escalante School of Mathematical and Statistical Sciences Center for Evolutionary Medicine and Informatics Arizona State University
More informationCASE STUDY: EXTINCTION OF FAMILY NAMES
CASE STUDY: EXTINCTION OF FAMILY NAMES The idea that families die out originated in antiquity, particilarly since the establishment of patrilineality (a common kinship system in which an individual s family
More informationIntroduction to Bayesian Statistics and Markov Chain Monte Carlo Estimation. EPSY 905: Multivariate Analysis Spring 2016 Lecture #10: April 6, 2016
Introduction to Bayesian Statistics and Markov Chain Monte Carlo Estimation EPSY 905: Multivariate Analysis Spring 2016 Lecture #10: April 6, 2016 EPSY 905: Intro to Bayesian and MCMC Today s Class An
More informationMonte Carlo Methods. Leon Gu CSD, CMU
Monte Carlo Methods Leon Gu CSD, CMU Approximate Inference EM: y-observed variables; x-hidden variables; θ-parameters; E-step: q(x) = p(x y, θ t 1 ) M-step: θ t = arg max E q(x) [log p(y, x θ)] θ Monte
More informationLIMIT THEOREMS FOR NON-CRITICAL BRANCHING PROCESSES WITH CONTINUOUS STATE SPACE. S. Kurbanov
Serdica Math. J. 34 (2008), 483 488 LIMIT THEOREMS FOR NON-CRITICAL BRANCHING PROCESSES WITH CONTINUOUS STATE SPACE S. Kurbanov Communicated by N. Yanev Abstract. In the paper a modification of the branching
More informationLocal consistency of Markov chain Monte Carlo methods
Ann Inst Stat Math (2014) 66:63 74 DOI 10.1007/s10463-013-0403-3 Local consistency of Markov chain Monte Carlo methods Kengo Kamatani Received: 12 January 2012 / Revised: 8 March 2013 / Published online:
More informationFAV i R This paper is produced mechanically as part of FAViR. See for more information.
Bayesian Claim Severity Part 2 Mixed Exponentials with Trend, Censoring, and Truncation By Benedict Escoto FAV i R This paper is produced mechanically as part of FAViR. See http://www.favir.net for more
More informationThe Particle Filter. PD Dr. Rudolph Triebel Computer Vision Group. Machine Learning for Computer Vision
The Particle Filter Non-parametric implementation of Bayes filter Represents the belief (posterior) random state samples. by a set of This representation is approximate. Can represent distributions that
More informationMachine Learning for Data Science (CS4786) Lecture 24
Machine Learning for Data Science (CS4786) Lecture 24 Graphical Models: Approximate Inference Course Webpage : http://www.cs.cornell.edu/courses/cs4786/2016sp/ BELIEF PROPAGATION OR MESSAGE PASSING Each
More informationBayesian Quadrature: Model-based Approximate Integration. David Duvenaud University of Cambridge
Bayesian Quadrature: Model-based Approimate Integration David Duvenaud University of Cambridge The Quadrature Problem ˆ We want to estimate an integral Z = f ()p()d ˆ Most computational problems in inference
More information2 Inference for Multinomial Distribution
Markov Chain Monte Carlo Methods Part III: Statistical Concepts By K.B.Athreya, Mohan Delampady and T.Krishnan 1 Introduction In parts I and II of this series it was shown how Markov chain Monte Carlo
More informationRecent Advances in Bayesian Inference Techniques
Recent Advances in Bayesian Inference Techniques Christopher M. Bishop Microsoft Research, Cambridge, U.K. research.microsoft.com/~cmbishop SIAM Conference on Data Mining, April 2004 Abstract Bayesian
More informationMarkov Chain Monte Carlo Methods
Markov Chain Monte Carlo Methods Sargur Srihari srihari@cedar.buffalo.edu 1 Topics Limitations of Likelihood Weighting Gibbs Sampling Algorithm Markov Chains Gibbs Sampling Revisited A broader class of
More informationModelling and forecasting of offshore wind power fluctuations with Markov-Switching models
Modelling and forecasting of offshore wind power fluctuations with Markov-Switching models 02433 - Hidden Markov Models Pierre-Julien Trombe, Martin Wæver Pedersen, Henrik Madsen Course week 10 MWP, compiled
More informationMarkov Chain Monte Carlo (MCMC) and Model Evaluation. August 15, 2017
Markov Chain Monte Carlo (MCMC) and Model Evaluation August 15, 2017 Frequentist Linking Frequentist and Bayesian Statistics How can we estimate model parameters and what does it imply? Want to find the
More informationApproximate Bayesian Computation: a simulation based approach to inference
Approximate Bayesian Computation: a simulation based approach to inference Richard Wilkinson Simon Tavaré 2 Department of Probability and Statistics University of Sheffield 2 Department of Applied Mathematics
More informationBayesian Statistical Methods. Jeff Gill. Department of Political Science, University of Florida
Bayesian Statistical Methods Jeff Gill Department of Political Science, University of Florida 234 Anderson Hall, PO Box 117325, Gainesville, FL 32611-7325 Voice: 352-392-0262x272, Fax: 352-392-8127, Email:
More informationOutline. Clustering. Capturing Unobserved Heterogeneity in the Austrian Labor Market Using Finite Mixtures of Markov Chain Models
Capturing Unobserved Heterogeneity in the Austrian Labor Market Using Finite Mixtures of Markov Chain Models Collaboration with Rudolf Winter-Ebmer, Department of Economics, Johannes Kepler University
More informationSTAT 518 Intro Student Presentation
STAT 518 Intro Student Presentation Wen Wei Loh April 11, 2013 Title of paper Radford M. Neal [1999] Bayesian Statistics, 6: 475-501, 1999 What the paper is about Regression and Classification Flexible
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 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 informationTutorial on Probabilistic Programming with PyMC3
185.A83 Machine Learning for Health Informatics 2017S, VU, 2.0 h, 3.0 ECTS Tutorial 02-04.04.2017 Tutorial on Probabilistic Programming with PyMC3 florian.endel@tuwien.ac.at http://hci-kdd.org/machine-learning-for-health-informatics-course
More informationINTRODUCTION TO BAYESIAN STATISTICS
INTRODUCTION TO BAYESIAN STATISTICS Sarat C. Dass Department of Statistics & Probability Department of Computer Science & Engineering Michigan State University TOPICS The Bayesian Framework Different Types
More informationText Mining for Economics and Finance Latent Dirichlet Allocation
Text Mining for Economics and Finance Latent Dirichlet Allocation Stephen Hansen Text Mining Lecture 5 1 / 45 Introduction Recall we are interested in mixed-membership modeling, but that the plsi model
More informationPOSTERIOR ANALYSIS OF THE MULTIPLICATIVE HETEROSCEDASTICITY MODEL
COMMUN. STATIST. THEORY METH., 30(5), 855 874 (2001) POSTERIOR ANALYSIS OF THE MULTIPLICATIVE HETEROSCEDASTICITY MODEL Hisashi Tanizaki and Xingyuan Zhang Faculty of Economics, Kobe University, Kobe 657-8501,
More informationDecomposition of supercritical branching processes with countably many types
Pomorie 26 June 2012 Serik Sagitov and Altynay Shaimerdenova 20 slides Chalmers University and University of Gothenburg Al-Farabi Kazakh National University Decomposition of supercritical branching processes
More informationCondensed Table of Contents for Introduction to Stochastic Search and Optimization: Estimation, Simulation, and Control by J. C.
Condensed Table of Contents for Introduction to Stochastic Search and Optimization: Estimation, Simulation, and Control by J. C. Spall John Wiley and Sons, Inc., 2003 Preface... xiii 1. Stochastic Search
More information28 : Approximate Inference - Distributed MCMC
10-708: Probabilistic Graphical Models, Spring 2015 28 : Approximate Inference - Distributed MCMC Lecturer: Avinava Dubey Scribes: Hakim Sidahmed, Aman Gupta 1 Introduction For many interesting problems,
More informationAdvanced Statistical Modelling
Markov chain Monte Carlo (MCMC) Methods and Their Applications in Bayesian Statistics School of Technology and Business Studies/Statistics Dalarna University Borlänge, Sweden. Feb. 05, 2014. Outlines 1
More informationBayesian Methods in Multilevel Regression
Bayesian Methods in Multilevel Regression Joop Hox MuLOG, 15 september 2000 mcmc What is Statistics?! Statistics is about uncertainty To err is human, to forgive divine, but to include errors in your design
More informationDynamic Approaches: The Hidden Markov Model
Dynamic Approaches: The Hidden Markov Model Davide Bacciu Dipartimento di Informatica Università di Pisa bacciu@di.unipi.it Machine Learning: Neural Networks and Advanced Models (AA2) Inference as Message
More informationPREPRINT 2007:25. Multitype Galton-Watson processes escaping extinction SERIK SAGITOV MARIA CONCEIÇÃO SERRA
PREPRINT 2007:25 Multitype Galton-Watson processes escaping extinction SERIK SAGITOV MARIA CONCEIÇÃO SERRA Department of Mathematical Sciences Division of Mathematical Statistics CHALMERS UNIVERSITY OF
More information14 Branching processes
4 BRANCHING PROCESSES 6 4 Branching processes In this chapter we will consider a rom model for population growth in the absence of spatial or any other resource constraints. So, consider a population of
More informationBayesian GLMs and Metropolis-Hastings Algorithm
Bayesian GLMs and Metropolis-Hastings Algorithm We have seen that with conjugate or semi-conjugate prior distributions the Gibbs sampler can be used to sample from the posterior distribution. In situations,
More informationGeneral Construction of Irreversible Kernel in Markov Chain Monte Carlo
General Construction of Irreversible Kernel in Markov Chain Monte Carlo Metropolis heat bath Suwa Todo Department of Applied Physics, The University of Tokyo Department of Physics, Boston University (from
More informationTopic Modelling and Latent Dirichlet Allocation
Topic Modelling and Latent Dirichlet Allocation Stephen Clark (with thanks to Mark Gales for some of the slides) Lent 2013 Machine Learning for Language Processing: Lecture 7 MPhil in Advanced Computer
More informationMCMC: Markov Chain Monte Carlo
I529: Machine Learning in Bioinformatics (Spring 2013) MCMC: Markov Chain Monte Carlo Yuzhen Ye School of Informatics and Computing Indiana University, Bloomington Spring 2013 Contents Review of Markov
More information