winbugs 求助

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

大虾们，我现有一模型需要应用winbugs编程，难住小弟了，特此请求高人指点。模型以及小弟编的部分程序如附件。谢谢
model{
    c<-10
    for(i in 1:n){zeros<-0
    zeros~dpois(zeros.mean)
    zeros.mean<--log(l)+c  
    l<-alpha*beta*(pow(y,(beta-1)))*(exp(-alpha*((pow(y,beta)-(pow(t,beta))))))}
    for(j in 2:n){log(beta[j])<-log(beta[j-1])+(eta[j])
    eta[j]~dnorm(0,lamda)}     
    s<-sqrt(1/lamda)
    lamda~dgamma(1,0.5)
    alpha~dgamma(1,0.5)
    beta[1]~dlnorm(0,0.1)
}
该程序难在lamda应该有两个值，我这里的程序只给出了一个。具体描述如下：如果认为y［i］为变点，则z［i］=1，lamda=0.005*dgamma(1,0.5),否则为z=0,lamda=dgamma(1,0.5).Z服从贝努力分布，p又是参数beta的函数。这样，程序中就出现了循环和选择了，困住小弟了，请高手指点。该模型的详细描述见PDF文件。模型原理很简单，很容易看懂。谢谢

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：winbugs WINBUG BUGS bug Win 求助 winbugs