[软体程式类别]: openBUGS [程式问题]: openBUGS要进行update却出现 "***Sorry something went wrong in procedure Ln in module Math***" [软体熟悉度]: 新手 [问题叙述]: 我想用openBUGS来处理含有 right censored和interval censored资料型态 的model里的参数估计问题 假设资料来自log-logistic(α,β)分配 资料的型态为(x,t1,t2,d),x为解释变数, t1为区间开始时间,t2为区间结束时间,d为指标,right censored=0,interval censored=1 即right censored=(x,t1,inf,0),interval censored=(x,t1,t2,1) 为了输入data方便我把它改成t1=t,t2=t1+k=t+k,inf=0 right censored=(x,t,0,0),interval censored=(x,t,t+k,1) 令α=b0+b1*x+b2*x^2,π=(a2+(1-a2)*exp(-exp(a0-a1*x))) 总之,每笔资料对应到的概似机率为 L=(((π)+(1-π)*(1+(t/α)^β)^(-1)-k)^(1-d)) *((1-π)*((1+(t/α)^β)^(-1)-(1+((t+k)/α)^β)^(-1)))^d 所以t~dloglik(logL),要估β,a0,a1,a2,b0,b1,b2 打完code放到openBUGS执行 check model,load data,compile都可以执行 但要进行update却出现 "***Sorry something went wrong in procedure Ln in module Math***" 爬过文有试过用较简单的model减少要估计的参数的个数 还是出现同样问题 完全不清楚到底哪里有问题....求版上高手解救.... [程式范例]: 以下是我打的code http://pastie.org/10953126 但却无法updata,求高手解救,谢谢~ --

※ 发信站: 批踢踢实业坊(ptt.cc), 来自: 140.112.74.231 ※ 文章网址: https://webptt.com/cn.aspx?n=bbs/Statistics/M.1477907578.A.E57.html