![](/img/trans.png)
[英]How to Summarize Univariate Results from coxph model summary in R - multilevel covariates?
[英]how to extract formula from coxph model summary in R?
我能夠按照以下教程中的概述復制coxph模型和提取組件: http : //www.sthda.com/english/wiki/cox-proportional-hazards-model 。 但是,我正在努力從模型中提取公式對象。
library("survival")
library("survminer")
data("lung")
head(lung)
covariates <- c("age", "sex", "ph.karno", "ph.ecog", "wt.loss")
univ_formulas <- sapply(covariates,
function(x) as.formula(paste('Surv(time, status)~', x)))
univ_models <- lapply( univ_formulas, function(x){coxph(x, data = lung)})
# Extract data
univ_results <- lapply(univ_models,
function(x){
x <- summary(x)
p.value<-signif(x$wald["pvalue"], digits=2)
wald.test<-signif(x$wald["test"], digits=2)
beta<-signif(x$coef[1], digits=2);#coeficient beta
HR <-signif(x$coef[2], digits=2);#exp(beta)
HR.confint.lower <- signif(x$conf.int[,"lower .95"], 2)
HR.confint.upper <- signif(x$conf.int[,"upper .95"],2)
HR <- paste0(HR, " (",
HR.confint.lower, "-", HR.confint.upper, ")")
res<-c(beta, HR, wald.test, p.value)
names(res)<-c("beta", "HR (95% CI for HR)", "wald.test",
"p.value")
return(res)
#return(exp(cbind(coef(x),confint(x))))
})
res <- t(as.data.frame(univ_results, check.names = FALSE))
as.data.frame(res)
我嘗試將以下語句添加到function(x)而不成功:
formula_extract <- paste( c(x$formula[[2]],
x$formula[[1]],
x$formula[[3]]), collapse='')
添加語句看起來像這樣,結果數據幀“ res”包含該列,但全為空(甚至不適用)。
# Extract data
univ_results <- lapply(univ_models,
function(x){
x <- summary(x)
p.value<-signif(x$wald["pvalue"], digits=2)
wald.test<-signif(x$wald["test"], digits=2)
beta<-signif(x$coef[1], digits=2);#coeficient beta
HR <-signif(x$coef[2], digits=2);#exp(beta)
HR.confint.lower <- signif(x$conf.int[,"lower .95"], 2)
HR.confint.upper <- signif(x$conf.int[,"upper .95"],2)
HR <- paste0(HR, " (",
HR.confint.lower, "-", HR.confint.upper, ")")
formula_extract <- paste( c(x$formula[[2]],
x$formula[[1]],
x$formula[[3]]), collapse='')
res<-c(beta, HR, wald.test, p.value, formula_extract)
names(res)<-c("beta", "HR (95% CI for HR)", "wald.test",
"p.value", "formula_extract")
return(res)
#return(exp(cbind(coef(x),confint(x))))
})
res <- t(as.data.frame(univ_results, check.names = FALSE))
as.data.frame(res)
但是,它本身可以起作用:
paste( c(univ_models$age$formula[[2]],
univ_models$age$formula[[1]],
univ_models$age$formula[[3]]), collapse='')
任何幫助,將不勝感激!
您需要做的就是使用其他變量:即
univ_results <- lapply(univ_models,
function(y){
x <- summary(y)
:
:
formula_extract <- deparse(y$formula)
:
}
as.data.frame(res)
beta HR (95% CI for HR) wald.test p.value formula_extract
age 0.019 1 (1-1) 4.1 0.042 Surv(time, status) ~ age
sex -0.53 0.59 (0.42-0.82) 10 0.0015 Surv(time, status) ~ sex
ph.karno -0.016 0.98 (0.97-1) 7.9 0.005 Surv(time, status) ~ ph.karno
ph.ecog 0.48 1.6 (1.3-2) 18 2.7e-05 Surv(time, status) ~ ph.ecog
wt.loss 0.0013 1 (0.99-1) 0.05 0.83 Surv(time, status) ~ wt.loss
聲明:本站的技術帖子網頁,遵循CC BY-SA 4.0協議,如果您需要轉載,請注明本站網址或者原文地址。任何問題請咨詢:yoyou2525@163.com.