Winbugs编程错误，求高手帮忙！

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

用winbugs软件运行到compile一步时，系统提示出现错误：array index is greater than array upper bound for R，但是如果只运行少量数据的话，就没有问题，希望大神们能帮我看一下程序有什么问题？怎么修改？万分感谢！！！
程序如下：
#liner regression for storks
model
{   
    for (i in 1:N)
        {   R ~ dnorm(mu,tau);
               mu <- intercept+slope*F;
            }
#prior for regression parameters
     intercept ~ dnorm(0,0.01);
     slope ~ dnorm(0,1000000);
#prior for precision parameter
  tau<- 1/tau1
  tau1~ dgamma(0.001,1000);
  #monitor the standard deviation
     sigma <- 1/sqrt(tau);

}

#data
list(N=60,R=c(-0.847335574，0.948132808，-0.015846396，0.493692142，1.503485125，0.467815568，-0.695609989，0.149998574，0.644530852，1.15767255，-0.904794484，0.695833328，0.839755242，0.254549023，-0.136739408，-1.13130779，-1.824518521，-0.165091925，-1.24270106，-0.516097238，-0.528051249，-2.195099911，-1.777527034，-0.175859617，-0.099581458，-0.654167732，-0.475567247，0.740174819，1.274970851，1.020761815，-0.468460018，-1.606178888，-1.386916986，0.910187771，0.227306949，0.650029346，0.404339696，-0.051155664，0.709363233，0.446425078，-0.286782117，0.317214063，-0.76660388，-0.295401085，-3.151480534，0.298469754，0.614632738，0.158743346，0.939005109，0.772503736，0.217853371，0.690332948，-1.148056227，-0.048242996，-0.888220296，0.381752332，-0.502714645，-0.452516044，0.898220844，0.200989919
),F=c(4.761932083，-0.303139501，-5.680798502，2.488715605，5.386920786，-0.049376452，-0.567056975，-5.339513943，7.652004282，0.115893176，4.298669179，-0.017239189，-13.60983163，-4.675317071，-3.964222479，-4.476225451，3.294683558，-6.023917828，-9.562518986，2.033881966，5.029464485，3.815298024，-9.279730114，14.15285344，-4.848883999，0.740226674，-5.285373406，-4.875563895，-5.052250501，-5.525350274，10.28183451，0.61890236，-0.279151254，0.530995914，3.420865422，-5.914043854，-0.818986884，-3.904567992，-4.077687837，-1.604102033，3.40317264，6.939586313，6.053081825，5.113201609，3.805389006，-8.562896448，-2.688869978，-10.30931221，-1.166470048，-3.421944698，3.870552995，-5.662569127，1.324117924，-5.743600092，-6.104322312，9.384305751，-9.698998775，-1.965484573，-8.586614937，1.810348892))
#intal value
list(intercept=0,slope=0,tau=0.01)

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：winbugs WINBUG BUGS Win bug 程序 软件

model
{   
    for (i in 1:N)
        {     R[i] ~ dnorm(mu[i],tau);
               mu[i] <- intercept+slope*F[i];
            }
#prior for regression parameters
     intercept ~ dnorm(0,0.01);
     slope ~ dnorm(0,1000000);
#prior for precision parameter
  tau<- 1/tau1
  tau1~ dgamma(0.001,1000);
  #monitor the standard deviation
     sigma <- 1/sqrt(tau);
}

你好，你会用WinBUGS吗？我想弄结构方程模型的，但不太会用这个软件，能否帮忙下，重谢

jnmq 发表于 2014-2-24 00:07
你好，你会用WinBUGS吗？我想弄结构方程模型的，但不太会用这个软件，能否帮忙下，重谢

你好，你的winbugs编程问题都解决了吗？想请教一下。

应该是因为样本数据不够

你试试重新粘贴数据，我之前也是这样，重新粘贴后可以啦

你好，我是学医学的，有贝叶斯统计建模问题（主要是WINBUGS）软件请教你可以么？qq 1400236514，我可以交学费。