A modified adaptive accept-reject algorithm for univariate densities with bounded support

被引:5
作者
Botts, Carsten [1 ]
机构
[1] Williams Coll, Dept Math & Stat, Williamstown, MA 01267 USA
关键词
log-concave; T; -concave; -convex; inflection points; bayesian analysis;
D O I
10.1080/00949651003645326
中图分类号
TP39 [计算机的应用];
学科分类号
081203 ; 0835 ;
摘要
The need to simulate from a univariate density arises in several settings, particularly in Bayesian analysis. An especially efficient algorithm which can be used to sample from a univariate density, f(X), is the adaptive accept-reject algorithm. To implement the adaptive accept-reject algorithm, the user has to envelope T circle f(X), where T is some transformation such that the density g(x) proportional to T(-1) (alpha + beta x) is easy to sample from. Successfully enveloping T circle f(X), however, requires that the user identify the number and location of T circle f(X)'s inflection points. This is not always a trivial task. In this paper, we propose an adaptive accept-reject algorithm which relieves the user of precisely identifying the location of T circle f(X)'s inflection points. This new algorithm is shown to be efficient and can be used to sample from any density such that its support is bounded and its log is three-times differentiable.
引用
收藏
页码:1039 / 1053
页数:15
相关论文
共 11 条
[1]  
[Anonymous], APPL STAT, DOI DOI 10.2307/2347565
[2]  
[Anonymous], J ROY STAT SOC C APP, DOI DOI 10.2307/2986138
[3]  
[Anonymous], QUALITY TECHNOLOGY Q
[4]   Random variable generation using concavity properties of transformed densities [J].
Evans, M ;
Swartz, T .
JOURNAL OF COMPUTATIONAL AND GRAPHICAL STATISTICS, 1998, 7 (04) :514-528
[5]   A REJECTION TECHNIQUE FOR SAMPLING FROM T-CONCAVE DISTRIBUTIONS [J].
HORMANN, W .
ACM TRANSACTIONS ON MATHEMATICAL SOFTWARE, 1995, 21 (02) :182-193
[6]   Black-box algorithms for sampling from continuous distributions [J].
Leydold, Josef ;
Hormann, Wolfgang .
PROCEEDINGS OF THE 2006 WINTER SIMULATION CONFERENCE, VOLS 1-5, 2006, :129-136
[7]   A note on the recurrence relations between moments of order statistics from right truncated log-logistic distribution [J].
Lin, CT .
STATISTICAL PAPERS, 2000, 41 (01) :99-107
[8]   Adaptive rejection Metropolis sampling using Lagrange interpolation polynomials of degree 2 [J].
Meyer, Renate ;
Cai, Bo ;
Perron, Francois .
COMPUTATIONAL STATISTICS & DATA ANALYSIS, 2008, 52 (07) :3408-3423
[9]   SIMULATION OF TRUNCATED NORMAL VARIABLES [J].
ROBERT, CP .
STATISTICS AND COMPUTING, 1995, 5 (02) :121-125
[10]  
Wakefield J. C., 1991, Statistics and Computing, V1, P129, DOI [10.1007/BF01889987, DOI 10.1007/BF01889987]