winbugs编程问题

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

对winbugs不是很懂 随便瞎编的，主要是想在已得到的最大似然函数去做参数估计，要怎么做啊，有没有哪位高手能指导一下下，万分感谢啊！

model
{
   for(i in 1:K)
{
                 mu1[i]<-b1+a1*exp(c1*T[i])
                 mu2[i]<-b2+a2*exp(c2*T[i])   
              }     

      for(l in 1:n)
         {
           for(h in 1:c)
             {   
            
                 U11[l,h]<-(y11[l,h]-mu1[1]/(sigma1)
                 U12[l,h]<-(y12[l,h]-mu1[2])/(sigma1)
                 U13[l,h]<-(y13[l,h]-mu1[3]/(sigma1)
                 U14[l,h]<-(y14[l,h]-mu1[4]/(sigma1)
                 U21[l,h]<-(y21[l,h]-mu2[1]/(sigma2)
                 U22[l,h]<-(y22[l,h]-mu2[2]/(sigma2)
                 U23[l,h]<-(y23[l,h]-mu2[3]/(sigma2)
                 U24[l,h]<-(y24[l,h]-mu2[4]/(sigma2)


                 U11[l,h]~dnorm(0,1)
                 U12[l,h]~dnorm(0,1)
                 U13[l,h]~dnorm(0,1)
                 U14[l,h]~dnorm(0,1)
                 U21[l,h]~dnorm(0,1)
                 U22[l,h]~dnorm(0,1)
                 U23[l,h]~dnorm(0,1)
                 U24[l,h]~dnorm(0,1)
                 
                u11[l,h]<-phi(U11[l,h],0,1)
                u12[l,h]<-phi(U12[i,h],0,1)
                u13[l,h]<-phi(U13[l,h],0,1)
                u14[l,h]<-phi(U14[l,h],0,1)
                u21[l,h]<-phi(U21[l,h],0,1)
                u22[l,h]<-phi(U22[l,h],0,1)
                u23[l,h]<-phi(U23[l,h],0,1)
                u24[l,h]<-phi(U24[l,h],0,1)

                  
                c1[l,h]<(-alpha*(1-exp(-alpha))*(exp(-alpha*(u11[l,h]+u21[l,h]))))/(((1-exp(-alhpa))-(1-exp(-alpha*(u11[l,h]))*(1-exp(-alpha*(u21[l,h])))^2)
               c2[l,h]<(-alpha*(1-exp(-alpha))*(exp(-alpha*(u12[l,h]+u22[l,h]))))/(((1-exp(-alhpa))-(1-exp(-alpha*(u12[l,h]))*(1-exp(-alpha*(u22[l,h])))^2)
               c3[l,h]<(-alpha*(1-exp(-alpha))*(exp(-alpha*(u13[l,h]+u23[l,h]))))/(((1-exp(-alhpa))-(1-exp(-alpha*(u13[l,h]))*(1-exp(-alpha*(u23[l,h])))^2)
               c4[l,h]<(-alpha*(1-exp(-alpha))*(exp(-alpha*(u14[l,h]+u24[l,h]))))/(((1-exp(-alhpa))-(1-exp(-alpha*(u14[l,h]))*(1-exp(-alpha*(u24[l,h])))^2)
               

            }
         }
         lnL.theta<-sum(log(c1[ ]))+sum(log(c2[ ]))+sum(log(c3[ ]))+sum(log(c4[ ]))-(K+n+c)*log(sqrt(2*pi))*(sigma1+sigma2)-1/2*(sum((U11[ ])^2)+sum((U12[ ])^2)+sum((U13[ ])^2)+sum((U14[ ])^2))-1/2*(sum((U21[ ])^2)+sum((U22[ ])^2)+sum((U23[ ])^2)+sum((U24[ ])^2))
   

#priors for regression parameters
                 a1~dnorm(0,0.0001)
                 b1~dnorm(0,0.0001)
                 c1~dnorm(0,0.0001)
                 a2~dnorm(0,0.0001)
                 b2~dnorm(0,0.0001)
                 c2~dnorm(0,0.0001)
                 alpha~dnorm(0,0.0001)
            

}

扫码加我 拉你入群

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

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

分享0 收藏1 回帖

关键词：winbugs WINBUG BUGS bug Win 编程