[問題] 倖存分析 in R軟體
假定考慮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)