winbugs做mcmc

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

初学者提问，我想在sv-t的基础上加入代表假日前和假日后的哑元变量，请问这段代码哪部分出现错误？
model{
for(i in 1:n)
{y[i]~dt(mu[i],p[i],omega)
p[i]<-exp(-theta[i])
mu[i]<-alpha0+alpha[1]*HOL[i]^pre+alpha[2]*HOL[i]^post
}
HOL~dbin(p,n)
alpha0~dnorm(0,0.001)
alpha1~dnorm(0,0.001)
alpha2~dnorm(0,0.001)
theta[1]~dnorm(lambda,delta,tau2)
for(j in 2:n)
{theta[j]~dnorm(theta2[j],itau2)
theta2[j]<-lambda+phi*(theta[j-1]-lambda)}
delta[i]<-delta[1]*HOL[i]^pre+delta[2]*HOL[i]^post
delta1~dnorm(0,0.001)
delta2~dnorm(0,0.001)
phi<-2*phi1-1
tau<-sqrt(1/itau2)
lambda~dnorm(0,0.01)
itau2~dgamma(2.5,0.025)
phi1~dbeta(20,1.5)
omega~dchisqr(8)
}

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：winbugs WINBUG BUGS mcmc Win

taening 发表于 2019-12-4 09:13
初学者提问，我想在sv-t的基础上加入代表假日前和假日后的哑元变量，请问这段代码哪部分出现错误？
model{ ...

。。