winbugs求助：winBUGS中怎么表示逆高斯分布及逆wishart分布？

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

winbugs求助：winBUGS中怎么表示逆高斯分布及逆wishart分布？如A~IWishart(R,2)，B~IG(0.5,0.5)

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：winbugs WINBUG BUGS wish Hart 高斯

我给楼主说下逆高斯怎么弄，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【i】，有什么问题还可以问：1530432346

多谢指教

强！！！

WinBUGS中可以实现逆伽马分布么？

libaijun19 发表于 2011-10-8 11:11
我给楼主说下逆高斯怎么弄，winbugs里，有些分布没有预设，需要我们自己另外定义。
# ---------------  IN ...

我有一个问题，想请教一下，如果一个参数x服从高斯分布，直接x~dnorm(mu,tau);
如果参数x服从逆高斯分布，该怎么表示呢？
我这里的x是一个参数，而不是似然函数服从逆高斯分布，谢谢！

使用0-1技巧

libaijun19 发表于 2011-10-8 11:11
我给楼主说下逆高斯怎么弄，winbugs里，有些分布没有预设，需要我们自己另外定义。
# ---------------  IN ...

dbern()是什么意思？