winbugs做模拟

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

请教大家个问题,我想做个模拟:v(t)=v*(2.25+z*3.33),z是服从正态分布的随机数,为求得v(t)的值,采用winbugs做MCMC模拟,但不知错在哪里?急需指教!十分感谢!

model
{
for (t in 1:n) {
z~dnorm(2.25,3.33);
p<-0.3*z+0.25;
v[t]<-478.55*exp(p);
}
}
list(z=c(1,2,1,1,3,6,4,2,5,7,1,1,6,0,,5,4,1,0,1,7),
n=20)

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：winbugs WINBUG BUGS bug Win 模拟 winbugs

程序如下

model
{
for (t in 1:n) {
z~dnorm(2.25,3.33);
p<-0.3*z+0.25;
v[t]<-478.55*exp(p);
}
}
list(z=c(1,2,1,1,3,6,4,2,5,7,1,1,6,0,,5,4,1,0,1,7),
n=20)

不知道楼主的意思是是什么

1. model
   
2. {
   
3. for (t in 1:n) {
   
4. z[t]~dnorm(2.25,3.33);
   
5. p[t]<-0.3*z[t]+0.25;
   
6. v[t]<-478.55*exp(p[t]);
   
7. }
   
8. }
   
9. list(z=c(1,2,1,1,3,6,4,2,5,7,1,1,6,0,,5,4,1,0,1,7),
   
10. n=20)

复制代码

另外要注意dnorm的第二个参数,这个参数不是一般我们用的variance，而是precision,也就是variance的倒数

也在学winbugs

感觉你模型里的表达式有点问题。