作者saintchu ()
看板Statistics
标题[问题] 幸存分析in R软体
时间Mon Jan 29 20:10:21 2007
假定考虑Cox提出的discrete logistic hazard model
也就是：
h(t,x)/1-h(t,x) = h0(t)/1-h0(t) *exp(bx)
t表时间，x为解释变数，h0(t)代表baseline hazard。
用partial likelihood的方法估计b，R的function为：
library(survival) # 需要先download这个package
cox<-coxph(Surv(time,default)~LTV, method="exact")
method="exact"....这是discrete time时的估计方法
问题：baseline hazard要怎麽估计?
是用"basehaz(cox)"这个function吗?
我试了几笔模拟出来的资料，都假设baseline hazard为t的二次式
h0(t)=0.001*(t-0.06*t^2)
但是估计出的baseline hazard却是monotone increasing
所以到底discrete time中的baseline hazard要用哪个function去估计?
另问： R有没有讨论区?
R网页中的search可搜寻过去别人问过的问题
那要怎麽发问?去哪发问?
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 61.57.133.136
※ 编辑: saintchu 来自: 61.57.133.136 (01/29 20:11)