Я бы подошел к этому, сначала посмотрев, как пакет survival
создает таблицу, которую он печатает по умолчанию.
Чтобы найти функцию, выполняющую эту печать, изучите класс подходящего объекта.и затем найдите метод печати для этого класса:
class(fit.pbc)
# [1] "coxph"
grep("coxph", methods("print"), value=TRUE)
# [1] "print.coxph" "print.coxph.null"
# [3] "print.coxph.penal" "print.summary.coxph"
После просмотра print.coxph
вот что я придумал:
cox <- fit.pbc
# Prepare the columns
beta <- coef(cox)
se <- sqrt(diag(cox$var))
p <- 1 - pchisq((beta/se)^2, 1)
CI <- round(confint(cox), 3)
# Bind columns together, and select desired rows
res <- cbind(beta, se = exp(beta), CI, p)
res <- res[c("age", "log(protime)"),]
# Print results in a LaTeX-ready form
xtable(res)