Abstract
SUMMARY Stochastic processes with independent increments play a central role in Bayesian nonparametric inference. The distributions of the increments of these processes, aside from fixed points of discontinuity, are infinitely divisible and their Laplace and/or Fourier transforms in the Lévy representation are usually known. Conventional Bayesian inference in this context has been limited largely to providing point estimates of the random quantities of interest, although Markov chain Monte Carlo methods have been used to obtain a fuller analysis in the context of Dirichlet process priors. In this paper, we propose and implement a general method for simulating infinitely divisible random variates when their Fourier or Laplace transforms are available in the Lévy representation. Theoretical justification is established by proving a convergence theorem that is a ‘sampling form’ of a classical theorem in probability. The results provide a method for implementing Bayesian nonparametric inference by using a wide range of stochastic processes as priors.
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
More From: Journal of the Royal Statistical Society Series B: Statistical Methodology
Disclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.