Jaka jest różnica między wyjściem „coef” i „(exp) coef” coxpha w R?

14

Próbowałem rozeznać, co dokładnie oznaczają „coef” i „(exp) coef” coxph. Wygląda na to, że „(exp) coef” to porównania pierwszej zmiennej w modelu zgodnie z grupą przypisaną w poleceniu.

W jaki sposób funkcja Coxph osiąga wartości dla „coef” i „(exp) coef”?

Dodatkowo, w jaki sposób coxph określa te wartości, gdy występuje cenzura?

annemphillip
źródło

Odpowiedzi:

21

Jeśli masz jedną zmienną objaśniającą, powiedzmy grupę leczenia, dopasowany jest model regresji Coxa coxph(); współczynnik ( coef) odczytuje się jako współczynnik regresji (w kontekście opisanego poniżej modelu Coxa), a jego wykładniczy poziom daje ryzyko w grupie leczonej (w porównaniu z grupą kontrolną lub placebo). Na przykład, jeśli p = - 1,80 , to ryzyko jest exp ( - 1,80 ) = 0,165 , to jest 16,5%.β^=1.80exp(1.80)=0.165

Jak zapewne wiesz, funkcja hazardu jest modelowana jako

h(t)=h0(t)exp(βx)

h0(t)exp(β1)x1ijexp(β(xixj))i nazywa się współczynnikiem ryzyka (lub współczynnikiem zapadalności). Przyjmuje się, że stosunek ten jest stały w czasie, stąd nazwa zagrożenia proporcjonalnego .

survregh0(t)

Parametry regresji są szacowane przez maksymalizację częściowego prawdopodobieństwa logarytmu określonego przez

=flog(exp(βxf)r(f)exp(βxr))

fr(f)βh0(t)h0(t)βh0(t)β

Jeśli chodzi o cenzurę, nie jest jasne, czy odwołujesz się do lewej cenzury (jak może się zdarzyć, jeśli weźmiemy pod uwagę początek skali czasowej, która jest wcześniejsza niż czas rozpoczęcia obserwacji, zwany również opóźnionym wejściem ), czy cenzurę prawą. W każdym razie więcej szczegółów na temat obliczania współczynników regresji i sposobu, w jaki pakiet przetrwania obsługuje cenzurę, można znaleźć w Therneau i Grambsch, Modeling Survival Data (Springer, 2000). Terry Therneau jest autorem poprzedniego pakietu S. Dostępny jest samouczek online .

χ2

chl
źródło
4
+1, to naprawdę pouczająca odpowiedź na dość proste i proste pytanie. Szczególnie doceniam linki do zasobów, aby uzyskać więcej informacji. NB link do samouczka online Terry'ego Therneau jest martwy; Pakiet do analizy przeżycia w S może być żywym łączem do tego samego materiału.
gung - Przywróć Monikę
5

Aby zacytować dokumentację dotyczącą metody drukowania dla obiektu Coxph, uzyskaną w R, wpisując ?survival::print.coxph:

coefficients współczynniki predyktora liniowego, które mnożą kolumny macierzy modelu.

To cała dokumentacja dostarczona przez autora pakietu. Opakowanie nie zawiera instrukcji obsługi ani winiety opakowania. R nie jest zaprojektowany jako przyjazny dla użytkownika, a dokumentacja zakłada, że ​​znasz już stosowane metody statystyczne.

Zakładam, że coefkolumna daje powyższe coefficients, a exp(coef)kolumna jest wykładnicza z nich. Ponieważ regresja Coxa obejmuje funkcję logarytmiczną, współczynniki są ilorazami ryzyka logarytmicznego . Wykładanie ich daje zatem współczynniki ryzyka pleców.

jeden przystanek
źródło