winbugs程序问题

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

#logistic regression for shorks
#breeding success as a function of both rainfall and temperature in May
model
{
  for(i in 1:N)
    {
     nbsuccess[i]~dbin(p[i],nbpairs[i])
     logit(p[i])<-intercept+slope.temperature*(temperature[i]-mean(temperature[]))/(sd(temperature[]))
       +slope.rainfall*(rainfall[i]-mean(rainfall[]))/(sd(rainfall[]))
    }
   #priors for regression parameters
   intercept~dnorm(0.0,0.001)
   slope.temperature~dnorm(0.0,0.001)
   slope.rainfall~dnorm(0.0,0.001)
}
#data
list(nbsuccess=c(151,105,73,107,113,87,77,108,118,122,112,120,122,89,69,71,53
41,53,31,35,14,18),nbpairs=c(173,164,103,113,122,112,98,121,132,136,133,137,145,117,
90,80,67,54,58,39,42,23,23),temperature=c(15.1,13.3,15.3,13.3,14.6,15.6,13.1,13.1,15,11.7,15.3,14.4,14.4,12.7,11.7,11.9,15.9,13.4,14,13.9,12.9,15.1,13),
rainfall=c(67,52,88,61,32,36,72,43,92,32,86,28,57,55,66,26,28,96,48,90,86,78,87),
N=23)
#initial values
list(intercept=0,slope.temperature=0,slope.rainfall=0,tau=1)
   

这个程序不知道为什么在加载数据时老是提示：expected right parenthesis?是个逻辑回归的例子，按理说是对的啊，，我是PPT上看到的，人家是可以运行的。望高手帮助下，谢谢。

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：winbugs WINBUG BUGS bug Win function success 程序