Wie 95% CI von R coxph bekommen?

Ich habe die folgende Funktion in R coxph() verwendet, um ein Cox-Hazard-Modell anzupassen. Ich möchte die richtigen Statistiken melden; In der Ausgabe ist jedoch kein 95% CI enthalten.Wie 95% CI von R coxph bekommen?

Surv(days, censor) ~ gender + age + treatment, data_1)

ich die folgenden Spalten nur bekommen.

coef exp(coef) se(coef) z p

Quelle

2016-04-28 ejg

Eine einfache Möglichkeit, Vertrauensintervalle für den Hazard Ratios mit Prädiktorvariablen verbunden zu bekommen wäre, die „Zusammenfassung“ Funktion paßt auf Ihrem Modell zu verwenden. Wenn Sie Konfidenzintervalle für die Koeffizientenschätzungen selbst benötigen, können Sie die Funktion "confint" verwenden. Die Exponentiation der Ergebnisse von confint kann auch verwendet werden, um die Hazard-Ratio-Konfidenzintervalle zu erhalten.

fit <- coxph(Surv(t,y) ~ x) 
summary(fit) #output provides HR CIs 
confint(fit) #coefficient CIs 
exp(confint(fit)) #Also HR CIs

Quelle

2016-04-28 23:40:15

Wenn Sie dies für die weitere Verarbeitung benötigen, sollten Sie eine tidy Lösung mit broom:

mod <- coxph(Surv(time, status) ~ as.character(sex) + age, data = lung) 

mod %>% 
    tidy %>% 
    mutate(
     estimate=exp(estimate), 
     conf.low=exp(conf.low), 
     conf.high=exp(conf.high) 
    ) %>% 
    select(term, estimate, starts_with("conf")) 

## term estimate conf.low conf.high 
## 1 sex2 0.598566 0.4310936 0.8310985 
## 2 age 1.017191 0.9989686 1.0357467

Weitere Details unter this reference

Quelle

2017-10-14 22:57:42 Paul

Sie einfach CI in CoxPH Modell Zusammenfassung in R finden (in Fall mit survival Paket.Ist unter lower .95 und upper .95

Quelle

2017-10-18 18:29:50 Abdal

Wie 95% CI von R coxph bekommen?

Antwort

Verwandte Themen