Ich versuche, eine ähnliche Funktion zu verwenden, wenn nicht tatsächlich, Regsubsets im Paket springt in Programm R bei der Auswahl der oberen Cox Proportional Hazards-Modelle für meine Daten. Ist das möglich? und wenn ja, existiert bereits eine Funktion?Kann ich die Auswahl aller Teilmengen für ein Cox Proportional Hazards Model in R durchführen?
Ich vermute, Sie bereits vertraut sind mit dem folgenden ... Wenn Sie AIC als Kriterium für ‚Topmodel‘ mit du dann wäre dies ein vernünftiger Ausgangspunkt:
library(survival)
data(colon)
c1 <- coxph(Surv(time=time, event=status) ~
as.factor(extent) + age + sex, data=colon)
step(c1)
Vorsicht wenn Sie fehlende Werte haben (NA). Es kann natürlich ein besseres Modell geben, das mit dieser Methode nicht gefunden wird, aber mit einer kleinen Anzahl potenzieller Prädiktoren werden Sie es wahrscheinlich nicht verpassen. Vorbehalte wie oben (danke @DWin) über die Verwendung von numerischen Methoden, wo eine fundierte Meinung zuverlässiger sein kann.
Sie sollten zuerst ernsthaft in Frage stellen, ob das eine gute Idee ist. –
Welches Modellauswahlverfahren empfehlen Sie? Ist es nicht angemessen zu sehen, wie nah die konkurrierenden Modelle sind, bevor ich ein Modell als das "beste" akzeptiere? – marcellt
Sie sollten zuerst über die zugrunde liegende Wissenschaft auf dem Gebiet der Untersuchung nachdenken. Es ist nicht möglich, ein "bestes Verfahren" zu empfehlen, es sei denn, die Ziele sind klar. –