请选择 进入手机版 | 继续访问电脑版
楼主: 3200302
813 2

新手。在WINBUGS里可以运行,改格式用R2winbugs运行不了,显示i in1:N的N没有被定义 [推广有奖]

  • 0关注
  • 0粉丝

高中生

90%

还不是VIP/贵宾

-

威望
0
论坛币
0 个
通用积分
0
学术水平
0 点
热心指数
0 点
信用等级
0 点
经验
155 点
帖子
7
精华
0
在线时间
68 小时
注册时间
2021-1-11
最后登录
2022-6-16

3200302 学生认证  发表于 2021-1-19 17:25:53 |显示全部楼层 |坛友微信交流群
3论坛币
代码如下

model
    {
        for(i in 1:N) {
       H[i]~dnorm(mu[i],tau)
       mu[i]<-1.3+a*pow((1-exp(-b*D[i])),c)    }
     a~ dnorm(0.0, 1.0E-6)
     b~ dnorm(0.0, 1.0E-6)
    c~ dnorm(0.0, 1.0E-6)
    tau ~ dgamma(0.01, 0.01)
}

D=c(6.5,9.86,7.5,20.13,23.41,7.18,23.1,5.24,19.11,26.06,4.6,4.95,31.69,7.39,27.29,9.82,16.08,5.85,2.182,24.42,21.65,4.78,5.05,12.44,19.8,12.68,10.78,23,5.05,5.01,11.82,24.63,19.32,13.96,13.74),H=c(5.2,6.5,5.7,20,19.4,3.65,20.2,6.4,18.9,19.6,5.7,6.3,20.5,5.6,19.8,8.1,19.1,5.2,1.7,19.3,21.2,6.9,3.5,9.3,21.2,8.1,10.3,19.3,6.8,6.7,8.31,19,19,13.3,9.5),N=35)

data=list(D,H,N)
inits=function(){list(a=40,b=0.06,c=1.15)}
output=bugs(data,inits,model.file="output.bug",parameters=c("a","b","c"),n.chains =1,n.iter=5100,n.burnin=100,n.thin=1,bugs.directory="D/WinBUGS14/",codaPkg=FALSE)


winbugs的log:
model is syntactically correct
data loaded
compile(1)
variable N is not defined
command #Bugs:inits cannot be executed (is greyed out)

关键词:R2WinBUGS winbugs WINBUG BUGS bug
已有 1 人评分经验 论坛币 收起 理由
nuomin + 100 + 20 鼓励积极发帖讨论

总评分: 经验 + 100  论坛币 + 20   查看全部评分

3200302 学生认证  发表于 2021-1-19 18:12:02 |显示全部楼层 |坛友微信交流群
帮我看看这个代码有问题么

使用道具

3200302 学生认证  发表于 2021-1-22 22:26:39 |显示全部楼层 |坛友微信交流群
大佬们帮帮忙,必有重谢

使用道具

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

本版微信群
加好友,备注jltj
拉您入交流群

京ICP备16021002-2号 京B2-20170662号 京公网安备 11010802022788号 论坛法律顾问:王进律师 知识产权保护声明   免责及隐私声明

GMT+8, 2024-3-29 20:36