原文鏈接:http:///?p=14683預(yù)期繪制生存曲線圖 實際object of type ‘symbol‘ is not subsettable 重現(xiàn)問題的步驟
library(survminer)
require("survival")
survie <- Surv(time, status) ~ sex fit<- survfit(survie, data = lung)
ggsurvplot(fit, data = lung)
ggsurvplot()嘗試從對象中捕獲生存公式Surv(time, status) ~ sex,如下所示: fit$call$formula ,。
通過分解R代碼,執(zhí)行fit $ call $ formula的結(jié)果是“ survie”而不是 Surv(time, status) ~ sex ,。這將導(dǎo)致錯誤消息,。 請擬合并可視化生存曲線,如下所示:
library(survival)
fit <- survfit(Surv(time, status) ~ sex, data = lung)
ggsurvplot(fit, data = lung) 如果您仍然想分解R代碼,,那么我建議使用surv_fit() [survminer軟件包中的] 函數(shù) ,,它是對R基本函數(shù)的擴展, survfit() 具有更多功能
您可以按以下方式使用它:
survie <- Surv(time, status) ~ sex
fit <- surv_fit(survie, data = lung)
ggsurvplot(fit, data = lung)
|