作者musaki (我也需要充电...)
看板Statistics
标题急问SAS的zero-inflalted negative binomial regression撰写方法
时间Sat Jun 21 18:00:27 2008
我要跑有无某疾病对於住院次数 跟住院费用的相关性
因为住院次数跟住院费用都不是常态分布
加上研究对象里面很多人都没有住院
所以住院费用跟住院次数有很多零
因此要使用到zero-inflated negative binomial的回归模式
我这边有找到一个方法
但是我不知道指令的意义
因此无法应用 不知道有没有人可以替我解答
*** ZINB Model
proc nlmixed data=nhanes.data;
/* the starting values */
parms b0 0 b1 0 bp0 0 bp1 0 k 1;
/* the logit model. Agedec = age/10 */
eta_p = bp0 + bp1*agedec;
q0 = exp(eta_p)/(1+exp(eta_p));
/* the negative binomial model */
eta = b0 + b1*agedec;
mu = exp(eta);
/* the ZINB loglikelihood function that combines the two parts of the mdoel */
/* q0 = 1-p = the probability of a case
having the negative binomial distribution */
if numhosp=0 then loglike = log(1-q0 + q0*exp(-mu));
else loglike = log(q0) + lgamma(numhosp + (1/k) - lgamma(numhosp+1) -
lgamma(1/k) + numhosp*log(k*mu) - (numhosp+(1/k))*log(1+k*mu));
model numhosp ~ general(loglike);
run;
感谢！！
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 203.64.48.9
※ 编辑: musaki 来自: 203.64.48.9 (06/21 18:43)
※ 编辑: musaki 来自: 203.64.48.9 (06/22 00:12)