楼主: ayyqyyq
1116 1

winbugs中如何实现if else [推广有奖]

  • 0关注
  • 0粉丝

等待验证会员

小学生

28%

还不是VIP/贵宾

-

威望
0
论坛币
7 个
通用积分
0
学术水平
0 点
热心指数
0 点
信用等级
0 点
经验
73 点
帖子
3
精华
0
在线时间
5 小时
注册时间
2019-5-21
最后登录
2019-7-10

楼主
ayyqyyq 发表于 2019-5-23 12:12:18 |AI写论文
3论坛币
程序如下,红色加粗部分是需要修改的,感谢百忙之中帮我解答疑惑,或您发现了其他问题可私聊。

model
{
#likelihood
   for(i in 1:n){
          x ~ dnorm( mu, tau0 )
         mu<- mu1*t*I[0]+(mu1*tau-mu2*tau+mu2*t[i,j+1]-mu1*t[i,j])*I[1]+mu2*t*I[2]
         if tau> t[i, j+1]
            I[0]=1
         if t[i,j]<tau< t[i, j+1]
            I[1]=1
         if tau< t[i, j]
            I[2]=1
         sigma <-1/sqrt(tau0)
#prior
   tau~dnorm(mutau,tautau)
   mutau~dorm(0,  1.0E-3)
   tautau~dgamma(0.01, 0.01)
  sigmatau<-1/sqrt(tautau)
   }

   mu1~dnorm(mumu1, tau1)
   mumu1~dorm(0, 1.0E-3)
   tau1~dgamma(0.01,0.01)
   sigmatau1<-1/sqrt(tau1)

   mu2~dnorm(mumu2, tau2)
   mumu2~dorm(0, 1.0E-3)
   tau2~dgamma(0.01,0.01)
   sigmatau2<-1/sqrt(tau2)

   tau0~dgamma(alpha,beta)
   sigmatau0<-1/tau0

}



关键词:winbugs WINBUG BUGS else 如何实现

沙发
ayyqyyq 发表于 2019-5-28 15:16:34
最后还是自己解决了

您需要登录后才可以回帖 登录 | 我要注册

本版微信群
加好友,备注jltj
拉您入交流群
GMT+8, 2026-1-7 13:51