作者phil5566 (5566)
看板Statistics
标题[程式] openbugs怎麽建立log logistic的分配?
时间Sat Oct 29 21:01:40 2016
[软体程式类别]:
openBUGS
[程式问题]:
小弟想建立一个model,里面有变数T,而T~log-logistic分配
但翻了openBUGS的分配指令并没有属於log-logistic的指令
想请教版上各位高手该如何处理这个问题?
[软体熟悉度]:
新手
[问题叙述]:
以下是小弟举的一个例子
model {
for(i in 1:N) {
T[i] ~dloglogistic(shape, scale)
}
shape ~ dexp(1)
scale ~ dgamma(0.1,1.0)
}
list(t = c(94.3,15.7,62.9,126,5.24,31.4,1.05,1.05,2.1,10.5), N=10)
其中的T[i] ~dloglogistic(shape, scale)这行并没有
dloglogistic(shape, scale)这个指令
想请教该如何解决?
谢谢
[程式范例]:
http://pastebin.com/037N6k3j
请各位大大指点迷津
谢谢
--
※ 发信站: 批踢踢实业坊(ptt.cc), 来自: 140.112.74.231
※ 文章网址: https://webptt.com/cn.aspx?n=bbs/Statistics/M.1477746103.A.40A.html
1F：推 recorriendo: 可参考random numer generation的书 有可以给出任意 10/31 01:02
2F：→ recorriendo: 分布的方法e.g inverse method,acceptance-rejection 10/31 01:07
3F：→ clickhere: T[i] <- exp(Y[i]) 10/31 10:01
4F：→ clickhere: X[i] ~ dbinom(Y[i], N[i]) 10/31 10:02
5F：→ phil5566: 看不懂C大的解法?为何不是Y[i]<-log[T[i]] 10/31 15:30
6F：→ phil5566: Y[i]~dlogis(mu, tau)) 10/31 15:30
7F：→ phil5566: 虽然这样估出来是mu和tau 10/31 15:31
8F：→ phil5566: 而不是dloglogistic的shape和scale,但不知可否透过 10/31 15:33
9F：→ phil5566: mu=ln(scale),tau=1/shape的关系得到scale和shape的估计 10/31 15:37