Cox-Regression in MATLAB
-
28-09-2019 - |
Frage
Ich weiß, es ist COXPHFIT Funktion in MATLAB Cox-Regression zu tun, aber ich habe Probleme zu verstehen, wie sie anzuwenden.
1) Wie zwei Gruppen von Proben mit Überlebensdaten in Tagen (survdays
) zu vergleichen, zu zensieren (cens
) und einiger Prädiktor-Wert (x
)? Die Gruppen, die durch groups
logische Variable definiert. Gruppen haben unterschiedliche Anzahl von Proben.
2) Welches ist das Basisparameter in coxphfit? Ich habe die Dokumentation lesen, aber wie soll ich die Basislinie richtig wählen?
Es wäre toll, wenn Sie eine Seite mit guten detaillierten Beispielen für medizinische Überlebensdaten kennen. Ich fand nur die Mathworks Demo , die nicht einmal coxphfit nicht erwähnt.
Möchten Sie wissen, kann eine andere 3rd-Party-Funktion für die Cox-Regression sein?
UPDATE . Der r
Tag seit der Antwort hinzugefügt, was ich habe ist für R
Lösung
Mit Überlebensanalyse, die Hazard-Funktion ist die momentane Todesrate.
Bei diesen Analysen Sie messen typischerweise was Wirkung etwas auf dieser Gefahrenfunktion. Zum Beispiel können Sie „nicht schlucken Arsen erhöhen die Rate, mit der die Menschen sterben?“ Fragen. Ein Hintergrund Gefahr ist die Ebene, auf der die Menschen sowieso sterben würde (ohne Arsen beim Schlucken, in diesem Fall).
Wenn Sie die Dokumentation für coxphfit
sorgfältig zu lesen, werden Sie feststellen, dass diese Funktion versucht, die Basis-Hazard zu berechnen; es ist nicht etwas, dass Sie eingeben.
Baseline Der X-Wert, bei denen zu Berechnen der Basis-Hazard.
EDIT: MATLAB coxphfit
Funktion offensichtlich nicht mit gruppierten Daten arbeiten. Wenn Sie gerne Schalter auf R sind, dann ist der anaylsis ein Einzeiler.
library(survival)
#Create some data
n <- 20;
dfr <- data.frame(
survdays = runif(n, 5, 15),
cens = runif(n) < .3,
x = rlnorm(n),
groups = rep(c("first", "second"), each = n / 2)
)
#The Cox ph analysis
summary(coxph(Surv(survdays, cens) ~ x / groups, dfr))
ANOTHER EDIT: Der baseline
Parameter auf MATLAB coxphfit
erscheint eine Normierungskonstante zu sein. R coxph
Funktion hat keine gleichwertigen Parameter. Ich schaute in Statistical Computing von Michael Crawley und es scheint darauf hinzudeuten, dass der Basis-Hazard isn ‚t wichtig, da es aufhebt, wenn Sie die Wahrscheinlichkeit, dass Ihre individuellen Sterben berechnen. Siehe Kapitel 33 und p615-616 insbesondere. Mein Wissen darüber, wie das Modell funktioniert, ist nicht tief genug, um die Diskrepanz in den MATLAB und R-Implementierungen zu erklären; vielleicht könnten Sie auf die fragen Stapel Exchange-Statistik Analyse Website.