如何在 R 中獲取僅攔截 cox ph 模型，這在生存包中意味着什么？

Question

一個人如何設置參數以使其僅在生存包或其他情況下攔截，以及關於 cox 比例風險模型，攔截僅真正意味着什么。

Answer 1

R 回歸公式中制作我稱之為“空模型”而不是“僅截取”模型的標准方法是使用~1作為唯一術語。 （事實證明，您也可以使用~0作為唯一的術語。）在生存模型中，“截距”是基線風險（所有協變量估計值都參考的估計值），可以在survival包中找到對於包含使用basehaz函數的協變量的模型。 因此，如果我理解您的意圖，“僅攔截”Cox 模型將只是未經調整的 Kaplan-Meier 估計。

 test1 <- list(time=c(4,3,1,1,2,2,3),   # example code ?coxph
               status=c(1,1,1,0,1,1,0), 
               x=c(0,2,1,1,1,0,0), 
               sex=c(0,0,0,0,1,1,1)) 
 # Fit a stratified model 
 nullfit <- coxph(Surv(time, status) ~ 0, test1) # Null model
 survfit( nullfit)
# ----
Call: survfit(formula = nullfit)

     n events median 0.95LCL 0.95UCL
[1,] 7      5      3       2      NA
# ----
 png(); plot( survfit( nullfit) ); dev.off()

Answer 2

沒有任何其他上下文，“僅攔截”模型只是平均值。 基本上，如果您在 x 上對 y 進行回歸，則僅截距模型是 y 平均值處的一條平線（在 x 處的所有值處）。