[問題] 急問SAS的zero-inflalted negative binomial r …
我要跑有無某疾病對於住院次數 跟住院費用的相關性
因為住院次數跟住院費用都不是常態分佈
加上研究對象裡面很多人都沒有住院
所以住院費用跟住院次數有很多零
因此要使用到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)