2017-06-28 2 views
1

ich zum ersten Mal ein cox Modell in R etablieren:Verwendung cox Modell zur Schätzung Überleben

test1<- test[1:20,] 
model.1 <- coxph(Surv(test1$days,test1$status==1) ~ test1$MTT+test1$ADC,data=test1) 

und als ich versuchte nächsten Patienten das Überleben wie dies vorherzusagen:

covs1 <- data.frame(test[21,]$MTT,test[21,]$ADC) 
summary(survfit(model.1, newdata= covs1, type ="aalen")) 

es gab mir zu viele Überleben Ergebnisse und die Warnung ist "'newdata' hatte 1 Zeile, aber Variablen gefunden haben 20 Zeilen" fyi, gibt es 20 Ereignisse und die Ergebnisse enthalten 20 Überlebensergebnisse.

Antwort

1

Die Namen der Spalten im Datenrahmen, die als Grundlage für eine Vorhersage verwendet werden, müssen die gleichen Spaltennamen haben wie die RHS der Modellformel. Ich glaube nicht, Ihr wird qualifiy, wenn Sie so etwas tun:

test1<- test[1:20,] 
model.1 <- coxph(Surv(days, status==1) ~ MTT + ADC, data=test1) 

covs1 <- test[21, c("MTT", "ADC")] 
# then do your prediction 

Sie nicht $ verwenden sollten Argumente Surv zu liefern. Es ist wichtig, dass das Modell in der Umgebung des Datenrahmens konstruiert wird.

+0

es hat funktioniert, danke! – ineffable