我正在尝试使用 survminer 包中的 ggsurvplot 函数中的趋势测试选项对有序分类变量进行一些生存分析。在计算趋势的对数秩检验时，如果我的变量是数字 (1,2,3,4) 或有序因子 ("1","2","3","4")，则获得的 p 值不同。我应该使用哪种格式，为什么？
我还使用了 survMisc 包中的 comp 函数，因为 ggsurvplot 似乎是基于这个包并获得了预期的相同结果。
非常感谢您的帮助。
这是重现该问题的代码：
library(survminer)
data("larynx", package="KMsurv")
larynx.table=larynx
larynx.table$stage.cat=as.character(larynx.table$stage)
larynx.table$stage.fact=factor(larynx.table$stage.cat,levels=c("1","2","3","4"))
larynx.table$stage.fact.inv=factor(larynx.table$stage.cat,levels=c("4","3","2","1"))
fit1<- survfit(Surv(time, delta) ~ stage, data = larynx.table)
fit2<- survfit(Surv(time, delta) ~ stage.cat, data = larynx.table)
fit3<- survfit(Surv(time, delta) ~ stage.fact, data = larynx.table)
fit4<- survfit(Surv(time, delta) ~ stage.fact.inv, data = larynx.table)
SA=ggsurvplot(fit1, data = larynx.table, pval = TRUE,test.for.trend = TRUE)
SB=ggsurvplot(fit2, data = larynx.table, pval = TRUE,test.for.trend = TRUE)
SD=ggsurvplot(fit3, data = larynx.table, pval = TRUE,test.for.trend = TRUE)
SE=ggsurvplot(fit4, data = larynx.table, pval = TRUE,test.for.trend = TRUE)
arrange_ggsurvplots(list(SA,SB,SD,SE),ncol=2,nrow=2)
