急求大神帮忙，小弟拜谢，winbugs提示有函数未定义，不明白

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

###### write the model file ######
require(R2Open BUGS)
ckls_jump<-function() {   
   y0<-alpha   
   ymean[1]<-alpha+beta*y0
   ytau[1]<-pow(y0,-2*gamma)   
   y[1]~dnorm(ymean[1],ytau[1])  

   for(t in 2:T){     
     j[t] ~dnorm(jmean,jtau)
     q[t]~dbern(k)     
     ymean[t] <-alpha+beta*y[t-1]+j[t]*q[t]
     ytau[t] <-pow(y[t-1],-2*gamma)     
     y[t] ~ dnorm(ymean[t],ytau[t])
   }   
   alpha ~ dunif(0,1)
   beta ~ dnorm(2,0.25)
   gamma ~ dunif(0.5,1.5)   
   jmean ~ dnorm(2,2.5)   
   jtau~dgamma(2.5,0.025)
   jsimga2 <- 1/jtau   
   k~dbeta(2,1)
}
write.model(ckls_jump,"ckls-j.txt")
###### prepare data   #######
data = list(T=length(r),y= r)
###### inits #####
inits = NULL
##### sample with Open BUGS
  
##########
cklsj <- bugs(data,
               inits,
               parameters=c("alpha","beta","gamma","jmean","jsigma2","k"),
               model.file="ckls-j.txt",            
               n.chains=2,
               n.iter=10000,
               debug=T
           )

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：winbugs WINBUG BUGS 不明白 bug require write file