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 程序

两个错

1.少了,

  nbsuccess=c(151,105,73,107,113,87,77,108,118,122,112,120,122,89,69,71,53,

2.initial values多了tau=1,remove it.

  ==>list(intercept=0,slope.temperature=0,slope.rainfall=0)

epoh 发表于 2011-8-4 18:52
两个错1.少了,  nbsuccess=c(151,105,73,107,113,87,77,108,118,122,112,120,122,89,69,71,53, 2.initial v ...

厉害啊，谢谢你！，，，顺便问下你，有没有SV-N SV-t SV-M的winbugs代码啊？有的话能否发我一份jiemin1217@126.com，谢谢

epoh 发表于 2011-8-4 18:52
两个错1.少了,  nbsuccess=c(151,105,73,107,113,87,77,108,118,122,112,120,122,89,69,71,53, 2.initial v ...

hi,昨天那个程序啊，数据能加载了，但编译时说educational version cannot do this model，难道是winbugs版本问题？

pls Obtaining the key for unrestricted use
   http://www.mrc-bsu.cam.ac.uk/bugs/winbugs/contents.shtml

epoh 发表于 2011-8-5 18:07
pls Obtaining the key for unrestricted use
   http://www.mrc-bsu.cam.ac.uk/bugs/winbugs/contents.sh ...

网页是打开了，以前也打开过这个网页，但不知道怎么弄。你能详细点说下吗？ 谢谢。

winbugs有哪些版本？
  https://bbs.pinggu.org/forum.php? ... p;page=1#pid5975906

epoh 发表于 2011-8-5 20:03
winbugs有哪些版本？
  https://bbs.pinggu.org/forum.php?mod=viewthread&tid=816264&page=1#pid5975906

问题解决了。原来我博客上有相关资料。呵呵，也非常感谢你！

epoh 发表于 2011-8-5 20:03
winbugs有哪些版本？
  https://bbs.pinggu.org/forum.php?mod=viewthread&tid=816264&page=1#pid5975906

一般SV模型：
y[t]=sigma[t]*u[t],u[t]~N(0,1);
log theta[t]=mu+phi*(log theta[t-1]-mu)+tau^2*v[t],v[t]~N(0,1)
这里E[y]当然明白；
但是E[y^2]=exp[mu+tau^2/2*(1-phi^2)]是怎么出来的，我推了好久都没有推到出来，感觉很不常规。

还有，为什么winbugs里变量的方差都表示成倒数呢？这一点很不理解啊。
以上两点，尤其是第一点，望老师指点下。

Stochastic volatility model:

  1.Observation Equations

  2.State Equations

数学式子不好表达

详细请参考

sv-models.pdf

sv-models.pdf (483.31 KB)

2011-8-9 20:38:19 上传