楼主: duweiqi46
1637 1

Winbugs中能不能这样实现循环,第一次用,希望有人帮我看看, [推广有奖]

  • 0关注
  • 0粉丝

VIP1

小学生

64%

还不是VIP/贵宾

-

威望
0
论坛币
1549 个
通用积分
0
学术水平
0 点
热心指数
0 点
信用等级
0 点
经验
160 点
帖子
3
精华
0
在线时间
7 小时
注册时间
2017-2-21
最后登录
2017-12-27

楼主
duweiqi46 发表于 2017-4-27 16:05:06 |AI写论文

+2 论坛币
k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群
赵安豆老师微信:zhaoandou666

经管之家联合CDA

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

感谢您参与论坛问题回答

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

+2 论坛币
model
{
    for(i in 1:n) {
                            for(j in 1:m) {
                                                           N[i, j]~dnorm(u[i], sigma[i])
                                                      }
                     }
             u[i]<-5.03+beta*log(theta[i])
            sigma[i]<-1/(0.052+beta*log(theta[i])-beta1*log(theta[i])
            theta[i]<-(300/s[i])
       }

       beta~dnorm(1.67, 0.00001)
       beta1~dnorm(1.87, 0.00001)
}

list(N=structure(.Data=c(5.603, 5.544, 5.528, 5.630, 5.594, 5.540, 5.581, 5.548, 5.426, 5.567, 5.554, 5.627, 5.630, 5.596, 5.626, 5.028, 5.074, 5.016, 4.894, 4.993, 5.071, 5.024, 5.035, 4.954, 5.039, 5.098, 5.057, 5.092, 5.082, 5.005, 4.784, 4.842, 4.776, 4.813, 4.813, 4.860, 4.798, 4.776, 4.758, 4.770, 4.755, 4.837, 4.736, 4.842, 4.796,  4.477, 4.400, 4.426, 4.462, 4.592, 4.411, 4.447, 4.402, 4.665, 4.475, 4.458, 4.551, 4.551, 4.525, 4.641), .Dim=c(4, 15)), s=c(200, 300, 350, 400),
     c=4, m=15)
list(beta=2, beta1=2)

二维码

扫码加我 拉你入群

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

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

关键词:winbugs WINBUG BUGS bug Win

沙发
lvpinpin 发表于 2017-10-27 09:10:24
你的问题解决了吗,你这个model运行时出现“expected variable name ”,我也遇到这问题,这是什么原因

您需要登录后才可以回帖 登录 | 我要注册

本版微信群
加好友,备注jltj
拉您入交流群
GMT+8, 2026-1-3 19:46