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(β′(xi−xj))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 .
survreg
h0(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
Aby zacytować dokumentację dotyczącą metody drukowania dla obiektu Coxph, uzyskaną w R, wpisując
?survival::print.coxph
: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
coef
kolumna daje powyższecoefficients
, aexp(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.źródło