作者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/m.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