[疑难杂症]哪位大神帮忙把这段代码注释一下?

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

毕业论文里用到马尔可夫链蒙特卡罗模拟方法模拟跳扩散模型的参数，打算用openbugs来进行模拟，但是看了一些别人的代码后，还是没理清楚openbugs的代码里一些符号和句子的意思。
比如下面这段代码，里面dpois，sumNumNeigh是什么意思？？能不能有哪位大神帮忙把这段代码注释一下，感激不尽！！！

#MODEL
model
{
for (i in 1:N) {
O ~ dpois(mu)
log(mu) <- log(E) + alpha0 + alpha1 * X/10 +
b[i
# Area-specific relative risk (for maps)
RR <- exp(alpha0 + alpha1 * X/10 + b)
}
# CAR prior distribution for random effects:
b[1:N] ~ car.normal(adj[], weights[], num[], tau)
for (k in 1:sumNumNeigh) {
weights[k] <- 1
}


# Other priors:
alpha0 ~ dflat()
alpha1 ~ dnorm(0, 1e-05)
tau ~ dgamma(0.5, 5e-04)
# prior on precision
sigma <- sqrt(1/tau)
# standard deviation
}
#DATA
list(N = 56, O = c(9, 39, 11, 9, 15, 8, 26, 7, 6,
20, 13, 5, 3, 8, 17, 9, 2, 7, 9, 7, 16, 31, 11, 7, 19, 15,
7, 10, 16, 11, 5, 3, 7, 8, 11, 9, 11, 8, 6, 4, 10, 8, 2,
6, 19, 3, 2, 3, 28, 6, 1, 1, 1, 1, 0, 0), E = c(1.4, 8.7,
3, 2.5, 4.3, 2.4, 8.1, 2.3, 2, 6.6, 4.4, 1.8, 1.1, 3.3, 7.8,
4.6, 1.1, 4.2, 5.5, 4.4, 10.5, 22.7, 8.8, 5.6, 15.5, 12.5,
6, 9, 14.4, 10.2, 4.8, 2.9, 7, 8.5, 12.3, 10.1, 12.7, 9.4,
7.2, 5.3, 18.8, 15.8, 4.3, 14.6, 50.7, 8.2, 5.6, 9.3, 88.7,
19.6, 3.4, 3.6, 5.7, 7, 4.2, 1.8), X = c(16, 16, 10, 24,
10, 24, 10, 7, 7, 16, 7, 16, 10, 24, 7, 16, 10, 7, 7, 10,
7, 16, 10, 7, 1, 1, 7, 7, 10, 10, 7, 24, 10, 7, 7, 0, 10,
1, 16, 0, 1, 16, 16, 0, 1, 7, 1, 1, 0, 1, 1, 0, 1, 1, 16,
10), num = c(3, 2, 1, 3, 3, 0, 5, 0, 5, 4, 0, 2, 3, 3, 2,
6, 6, 6, 5, 3, 3, 2, 4, 8, 3, 3, 4, 4, 11, 6, 7, 3, 4, 9,
4, 2, 4, 6, 3, 4, 5, 5, 4, 5, 4, 6, 6, 4, 9, 2, 4, 4, 4,
5, 6, 5), adj = c(19, 9, 5, 10, 7, 12, 28, 20, 18, 19, 12,
1, 17, 16, 13, 10, 2, 29, 23, 19, 17, 1, 22, 16, 7, 2, 5,
3, 19, 17, 7, 35, 32, 31, 29, 25, 29, 22, 21, 17, 10, 7,
29, 19, 16, 13, 9, 7, 56, 55, 33, 28, 20, 4, 17, 13, 9, 5,
1, 56, 18, 4, 50, 29, 16, 16, 10, 39, 34, 29, 9, 56, 55,
48, 47, 44, 31, 30, 27, 29, 26, 15, 43, 29, 25, 56, 32, 31,
24, 45, 33, 18, 4, 50, 43, 34, 26, 25, 23, 21, 17, 16, 15,
9, 55, 45, 44, 42, 38, 24, 47, 46, 35, 32, 27, 24, 14, 31,
27, 14, 55, 45, 28, 18, 54, 52, 51, 43, 42, 40, 39, 29, 23,
46, 37, 31, 14, 41, 37, 46, 41, 36, 35, 54, 51, 49, 44, 42,
30, 40, 34, 23, 52, 49, 39, 34, 53, 49, 46, 37, 36, 51, 43,
38, 34, 30, 42, 34, 29, 26, 49, 48, 38, 30, 24, 55, 33, 30,
28, 53, 47, 41, 37, 35, 31, 53, 49, 48, 46, 31, 24, 49, 47,
44, 24, 54, 53, 52, 48, 47, 44, 41, 40, 38, 29, 21, 54, 42,
38, 34, 54, 49, 40, 34, 49, 47, 46, 41, 52, 51, 49, 38, 34,
56, 45, 33, 30, 24, 18, 55, 27, 24, 20, 18), sumNumNeigh = 234)
#INITIAL VALUES
list(tau = 1, alpha0 = 0, alpha1 = 0, b = c(0, 0,
0, 0, 0, NA, 0, NA, 0, 0, NA, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0))

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：疑难杂症 openbugs 蒙特卡罗模拟 model 马尔可夫链 蒙特卡罗 毕业论文 模型

我跟你的问题一样啊，winBUGS指导全是英文的，根本看不懂。你是学什么专业的，可否探讨一下？

snowwhite@HXH 发表于 2014-6-12 17:27
我跟你的问题一样啊，winBUGS指导全是英文的，根本看不懂。你是学什么专业的，可否探讨一下？

学金融的，这两天看这个东西看得懊恼死了

吃不饱的鱼 发表于 2014-6-12 21:03
学金融的，这两天看这个东西看得懊恼死了

dpois是服从泊松分布的意思，sum是求和的意思。我是学环境工程的，你现在这个软件学的怎么样了啊，我懂一点点了，可是好多还是闹不明白

Nicolle 发表于 2014-6-15 11:03
Page 209,Disease Mapping with WinBUGS and MLwiN

MODEL ckls{
#Hyperprior
dt<-1/365;
alphastar~dnorm(0,0,1);
betastar~dnorm(0,0,1);
rho~dunif(0.5,1.5);
isigma2~dgamma(2.5,0.025);
sigma2<-1/isigam2;
Y0itheta2<-1/(dt*sigma2);
Y0~dnorm(alphastar,Y0itheta2);
Ymean[1]<-(alphastar+betastar*Y0);
Yitheta[1]<-1/(dt*sigma2*pow(Y0,2*rho));
#Prior Distribution of rs
for (i in 1:k){
Ymean[i+1]<-(alphastar+beatstar*Y);
Yitheta[i+1]<-1/(dt*sigma2*pow(Y,2*rho));
Y[i+1]~dnorm(Ymean[i+1],Yitheta[i+1]);
}
}
烦大神帮我翻译下这段代码吧，不胜感激