楼主: yll95
1217 2

[问答] 【急求解答】请教一个winbugs问题 [推广有奖]

  • 1关注
  • 1粉丝

已卖:68份资源

大专生

86%

还不是VIP/贵宾

-

威望
0
论坛币
498 个
通用积分
37.5846
学术水平
5 点
热心指数
5 点
信用等级
5 点
经验
2362 点
帖子
52
精华
0
在线时间
61 小时
注册时间
2017-9-22
最后登录
2021-12-28

楼主
yll95 发表于 2019-2-16 21:15:59 |AI写论文
20论坛币
如图所示,想在原来的logistic函数中增加一个参数lambda,因为极大似然估计太复杂所以考虑用winbugs进行估计,请问我的model写的正确吗?感恩model
   {
    for (i in 1:N) {
    y~dbin(p,1)
   p<-1/pow((exp(-beta0-beta1*x)+1),lambda)
       }
beta1~dnorm(0,1)
beta0~dnorm(0,1)
lambda~dgamma(2,1)
}

list(y=c(1,1,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0),
x=c(1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23),
N=23)

#Inits

list(beta1=1,beta0=0,lambda=1) 微信图片_20190216210730.png


关键词:winbugs WINBUG BUGS bug Win winbugs 贝叶斯 贝叶斯估计 Gibbs Gibbs抽样

沙发
啊啊啊啊啊吖 发表于 2019-3-1 14:49:27
帮楼主顶一下帖子
~~

藤椅
gerryli 发表于 2021-10-17 21:29:21
This should work.

model
   {
    for (i in 1:N) {
    y[i]~dbin(p[i],1)
   p[i]<-1/pow((exp(-beta0-beta1*x[i])+1),lambda)
       }
beta1~dnorm(0,1)
beta0~dnorm(0,1)
lambda~dgamma(2,1)
}

list(y=c(1,1,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0),
x=c(1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23),
N=23)

#Inits

list(beta1=1,beta0=0,lambda=1)

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

本版微信群
加好友,备注cda
拉您进交流群
GMT+8, 2025-12-29 05:14