求助winbugs里怎么表示 IG 逆高斯分布？

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

1：不知道winbugs里是否有IG（逆高斯分布），只是我不知道像   dnorm 这样来表示
2：如果程序里没有，那么我该怎么把逆高斯分布用到里面去，也即如何定义一个库中没有的分布函数


有这方面学习交流  加QQ 1530432346.说明是学习winbugs。

扫码加我 拉你入群

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

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

分享0 收藏1 回帖

关键词：winbugs WINBUG BUGS Win bug 分布函数 我不知道 高斯 程序 如何

我给楼主说下逆高斯怎么弄，winbugs里，有些分布没有预设，需要我们自己另外定义。
# ---------------  INVERSE GAUSSIAN MODEL ------------------------------------
# ---------------  ZEROS TRICK ------------------------------------
model{


      C <- 10000
      for (i in 1:n) {
          zeros <- 0
          zeros ~ dpois( zeros.mean)
          zeros.mean <- -l + C
          l <-   0.5*( log( lambda ) - log(2*3.14) - 3*log(y) ) - 0.5* lambda * pow( (y-mu)/mu, 2 )/y
          log(mu) <- beta[1] + beta[2]*x1 + beta[3]*x2 + beta[4]*x3 + beta[5]*x4
       }  
       # priors
       for (j in 1:5){ beta[j] ~ dnorm( 0.0 ,  0.001 ) }
       lambda ~ dgamma( 0.01, 0.01)
       s2 <- 1/lambda
}
INITS
list( beta=c(0,0,0,0,0), lambda=1 )


# ---------------  ONES TRICK ------------------------------------
model{
      C <- 10
      for (i in 1:n) {
          ones <- 1
          ones ~ dbern( ones.p)
          ones.p <- exp( l - C )
          # write here the expression of the log-likelihood for i observation
          l <-   0.5*( log( lambda ) - log(2*3.14) - 3*log(y) ) - 0.5* lambda * pow( (y-mu)/mu, 2 )/y
          log(mu) <- beta[1] + beta[2]*x1 + beta[3]*x2 + beta[4]*x3 + beta[5]*x4
       }  
       # priors
       for (j in 1:5){ beta[j] ~ dnorm( 0.0 ,  0.001 ) }
       lambda ~ dgamma( 0.01, 0.01)
       s2 <- 1/lambda
}




在自己的程序里这样运用如下：
#model
{
  
C <- 10
      for (i in 1:n) {
          ones <- 1
          ones ~ dbern( ones.p)
          ones.p <- exp( l - C )
          # write here the expression of the log-likelihood for i observation
          l <-   0.5*( log( lambda ) - log(2*3.14) - 3*log(y) ) - 0.5* lambda * pow( (y-mu)/mu, 2 )/y
           lambda<-2
          mu<-4
         


    delta2<-y


y,就是我们要用的逆高斯分布，我这里期望值：4，delta的平方：2，lambda与delta平方实际室友关系的，但是，转化后lambda还是为2，转化没有什么意义，相当于没变，在用逆高斯分布时，你只需要改这两个值，然后引用y