hus.gibbs(SMPracticals)
hus.gibbs()所属R语言包:SMPracticals
Gibbs Sampler for Poisson Changepoint Model, Practical 11.6
Gibbs抽样泊松Changepoint的模型,实际11.6
译者:生物统计家园网 机器人LoveR
描述----------Description----------
This function implements a Gibbs sampler for the Poisson changepoint model applied to the HUS data used in Example 4.40 and Practical 11.6 of Davison (2003), which should be consulted for details.
此功能实现了使用的HUS数据的4.40 11.6戴维森(2003年),应详细咨询和实用的的泊松结构变动模型的Gibbs抽样。
用法----------Usage----------
hus.gibbs(init, y, R = 10, a1 = 1, a2 = 1, c = 0.01, d = 0.01)
参数----------Arguments----------
参数:init
Initial values for parameters
参数的初始值
参数:y
A series of Poisson counts
一系列的泊松计数
参数:R
Number of iterations of sampler
取样器的迭代数
参数:a1
Value of a hyperparameter
一个超参数的值
参数:a2
Value of a hyperparameter
一个超参数的值
参数:c
Value of a hyperparameter
一个超参数的值
参数:d
Value of a hyperparameter
一个超参数的值
Details
详细信息----------Details----------
This is provided simply so that readers spend less time typing. It is not intended to be robust and general code.
这是简单地让读者花更少的时间打字。它的目的不是是坚固和一般的代码。
值----------Value----------
A matrix of size R x 7, whose first five columns contain the values of the parameters for the iterations. Columns 6 and 7 contain the log likelihood and log prior for that iteration.
尺寸R×7,其第5列中包含的参数的值的迭代的矩阵。 6和第7列包含对数似然迭代前登录。
(作者)----------Author(s)----------
Anthony Davison (<code>anthony.davison@epfl.ch</code>)
参考文献----------References----------
实例----------Examples----------
## From Example 11.6:[#从例11.6:]
hus <- c(1,5,3,2,2,1,0,0,2,1,1,7,11,4,7,10,16,16,9,15)
system.time( gibbs.out <- hus.gibbs(c(5, 5, 1, 1, 2), hus, R=1000))
plot.ts(gibbs.out[,1], main="lambda1") # time series plot for lam1[时间序列图lam1]
plot.ts(gibbs.out[,2], main="lambda1") # time series plot for lam2[时间序列图lam2]
plot.ts(gibbs.out[,6], main="log lik") # and of log likelihood[和对数似然]
table(gibbs.out[,5]) # tabulate observed values of tau[制表观测值的头]
rm(hus)
转载请注明:出自 生物统计家园网(http://www.biostatistic.net)。
注:
注1:为了方便大家学习,本文档为生物统计家园网机器人LoveR翻译而成,仅供个人R语言学习参考使用,生物统计家园保留版权。
注2:由于是机器人自动翻译,难免有不准确之处,使用时仔细对照中、英文内容进行反复理解,可以帮助R语言的学习。
注3:如遇到不准确之处,请在本贴的后面进行回帖,我们会逐渐进行修订。
|