困扰了数天的问题，R2Winbugs求助！附代码

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

在R中运行了：
bugs.sim<-bugs(data,inits=NULL,parameters,noyear,n.chains=1,n.iter=40000,n.burnin=1000,bugs.directory="D:/winbugs14/WinBUGS14/",working.directory="C:/BUGS",debug=TRUE)
之后，在winbugs的Log中得到如下信息：
display(log)
check(C:/BUGS/noyear.bug.txt)
model is syntactically correct
data(C:/BUGS/data.txt)
undefined variable
compile(1)
gen.inits()
command #Bugs:gen.inits cannot be executed (is greyed out)
thin.updater(39)
update(26)
command #Bugs:update cannot be executed (is greyed out)
set(beta)
command #Bugs:set cannot be executed (is greyed out)
set(p1)
command #Bugs:set cannot be executed (is greyed out)
set(p2)
command #Bugs:set cannot be executed (is greyed out)
set(p3)
command #Bugs:set cannot be executed (is greyed out)
set(deviance)
command #Bugs:set cannot be executed (is greyed out)
dic.set()
command #Bugs:dic.set cannot be executed (is greyed out)
update(1000)
command #Bugs:update cannot be executed (is greyed out)
coda(*,C:/BUGS/coda)
command #Bugs:coda cannot be executed (is greyed out)
stats(*)
command #Bugs:stats cannot be executed (is greyed out)
dic.stats()

DIC
history(*,C:/BUGS/history.odc)
command #Bugs:history cannot be executed (is greyed out)
save(C:/BUGS/log.odc)
save(C:/BUGS/log.txt)


请问是哪里出了问题，是给的原始数据的问题吗。。？但是检查了好几遍，实在是不知道问题在哪里了。。请各位仁兄帮忙！十分感谢！！

代码在2楼

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：R2WinBUGS winbugs WINBUG BUGS bug

为方便解答，给出代码：

1. model{
   
2.   
   
3.   for (y in 1:N.years.sites){  
   
4.     for (d in 1:response.day){  
   
5.       events[d,y] ~ dpois(p[d,y])  # raw data is {0/1} for each day
   
6.       p[d,y] <- max(p0[d,y],.001) # trick to avoid crashing from low values
   
7.       p0[d,y] <- a[d] * risk[d,y]
   
8.       risk[d,y] <- exp(beta[1]*temp.cum[d,y] + beta[2]*frost7[d,y] +
   
9.                          beta[3]*precip7[d,y] )   # risk function that can include any covariates of interest, time-dependent or otherwise
   
10.       
    
11.       risk2[d,y] <- min(risk[d,y], 100)
    
12.      # events.pred[d,y] <- 1 - pow(exp(-sum(a[1 : d])), risk2[d,y] )
    
13.     }}}
    
14. 
    
15. for (d in 1:N.days){  #
    
16.   a[d] <- p1/ (1+exp((2.2/p3) * (p2-d)))     # baseline function used in this study (could be modeled otherwise for different applications, or left non-parametric)
    
17.   
    
18.   #p1 <- 1 # max
    
19.   #p2 <- 200 # point of inflection
    
20.   #p3 <- 20 # time between .1 and .5 and .5 and .9
    
21.   
    
22. }  
    
23. p1~dlnorm(p1.sp, p1.tau)
    
24. p3~dlnorm(p3.sp, p3.tau)
    
25. 
    
26. p2~dlnorm(p2.sp, p2.tau)
    
27. 
    
28. 
    
29. ## Priors (BUGS normal distributions use precision instead of variance, where precision = 1/variance)
    
30. p1.sp ~ dnorm(0,.001)
    
31. p2.sp ~ dnorm(3,.001)
    
32. p3.sp ~ dnorm(50,.001)
    
33. 
    
34. p1.tau ~ dgamma(.01,.01)
    
35. p2.tau ~ dgamma(.01,.01)
    
36. p3.tau ~ dgamma(.01,.01)
    
37. 
    
38. # Overall species baseline curve
    
39. for (d in 1:N.days){    a.sp[d] <- p1.sp/ (1+exp((2.2/p3.sp) * (p2.sp-d)))   }   
    
40. 
    
41. 
    
42. for(pp in 1:n.beta){
    
43.   beta[pp] ~ dnorm(beta0[pp], beta.tau[pp])
    
44. }
    
45. 
    
46. for(pp in 1:n.beta){
    
47.   beta0[pp]~ dnorm(0,.001)  
    
48.   beta.tau[pp] <- pow(beta.sigma[pp], -2)
    
49.   beta.sigma[pp] ~ dunif(0,10)
    
50. }
    
51. 
    
52. } # end model

复制代码

数据则是用csv格式做的，有多个文档，已经转换成vector；
data<-list("event","light","prep","temp.cum","N.sites","response.day","n.beta","N.days")
parameters<-c("beta","p1","p2","p3")
inits<-function(){list()}

已经读了论坛里大量的关于winbugs的帖子，仍然想不通问题在哪里。。真心求助！

不知道能不能艾特出大神。。。好吧发现不能艾特。。算了。。

自己解决了该问题。。确实是原始数据的问题

空sutra 发表于 2015-3-16 20:15
自己解决了该问题。。确实是原始数据的问题

你好！我现在也遇到这样的问题了，请问您说的原始数据问题是什么意思呢？格式还是什么，如果是格式需要什么格式呢

撒哈拉lu 发表于 2017-12-5 20:53
你好！我现在也遇到这样的问题了，请问您说的原始数据问题是什么意思呢？格式还是什么，如果是格式需要什 ...

同问啊，请问是什么问题呢