Simulation of right and left truncated gamma distributions by mixtures

被引:32
作者
Philippe, A [1 ]
机构
[1] UNIV ROUEN,CNRS,UPRES A6085,F-76821 MONT ST AIGNAN,FRANCE
关键词
minimum acceptance probability; mixture distribution; accept-reject algorithm;
D O I
10.1023/A:1018534102043
中图分类号
TP301 [理论、方法];
学科分类号
081202 ;
摘要
We study the properties of truncated gamma distributions and we derive simulation algorithms which dominate the standard algorithms for these distributions. For the right truncated gamma distribution, an optimal accept-reject algorithm is based on the fact that its density can be expressed as an infinite mixture of beta distribution. For integer values of the parameters, the density of the left truncated distributions can be rewritten as a mixture which can be easily generated. We give an optimal accept-reject algorithm for the other values of the parameter. We compare the efficiency of our algorithm with the previous method and show the improvement in terms of minimum acceptance probability. The algorithm proposed here has an acceptance probability which is superior to e/4.
引用
收藏
页码:173 / 181
页数:9
相关论文
共 10 条
[1]  
Abramowitz M., 1964, HDB MATH FUNCTIONS
[2]  
Dagpunar J. S., 1978, Journal of Statistical Computation and Simulation, V8, P59, DOI 10.1080/00949657808810248
[3]   MAXIMUM LIKELIHOOD FROM INCOMPLETE DATA VIA EM ALGORITHM [J].
DEMPSTER, AP ;
LAIRD, NM ;
RUBIN, DB .
JOURNAL OF THE ROYAL STATISTICAL SOCIETY SERIES B-METHODOLOGICAL, 1977, 39 (01) :1-38
[4]  
DEVROYE L, 1985, NONUNIFORM RANDOM VA
[5]   BAYESIAN-ANALYSIS OF CONSTRAINED PARAMETER AND TRUNCATED DATA PROBLEMS USING GIBBS SAMPLING [J].
GELFAND, AE ;
SMITH, AFM ;
LEE, TM .
JOURNAL OF THE AMERICAN STATISTICAL ASSOCIATION, 1992, 87 (418) :523-532
[6]  
GEYER CJ, 1992, J R STAT SOC B, V54, P657
[7]  
Robert C., 1996, METHODES MONTE CARLO
[8]   SIMULATION OF TRUNCATED NORMAL VARIABLES [J].
ROBERT, CP .
STATISTICS AND COMPUTING, 1995, 5 (02) :121-125
[9]  
ROBERT CP, 1996, 9624 CREST
[10]  
[No title captured]