p值不存储在survdiff类中，因此必须在输出时动态计算。要重现p值，可以使用chisq分布函数：“pchisq”

diff = survdiff(Surv(Time, Censored) ~ Treatment+Gender, data = dat)
pchisq(diff$chisq, length(diff$n)-1, lower.tail = FALSE)

函数print.survdiff中显示p值的代码为：

cat("\n Chisq=", format(round(x$chisq, 1)), " on", df, 
            "degrees of freedom, p=", format(signif(1 - pchisq(x$chisq, 
                df), digits)), "\n")

引导它的代码：

if (is.matrix(x$obs)) {
            otmp <- apply(x$obs, 1, sum)
            etmp <- apply(x$exp, 1, sum)
        }         else {
            otmp <- x$obs
            etmp <- x$exp
        }
        df <- (sum(1 * (etmp > 0))) - 1

并且参数列表中的“digits”设置为3，因此使用surv.diff帮助页面上的示例：

x <- survdiff(Surv(time, status) ~ pat.karno + strata(inst), data=lung) 
cat( "p=", format(signif(1 - pchisq(x$chisq, 
                 df), digits)) )
#p= 0.00326

回应评论：在该示例中，第二码块简化为：

df <- with(x,    (sum(1 * (apply(x$exp, 1, sum) > 0))) - 1 )
> df
[1] 7

使用broom包中的glance()函数，可以非常容易地获得p.value。

diff = survdiff(Surv(Time, Censored) ~ Treatment+Gender, data = dat)

broom::glance(diff)$p.value