Kontekst
Próbuję zrozumieć, w jaki sposób Rxx () przyjmuje i obsługuje powtarzające się wpisy dla podmiotów (lub pacjenta / klienta, jeśli wolisz). Niektórzy nazywają to długim formatem, inni nazywają to „powtarzanymi pomiarami”.
Zobacz na przykład zestaw danych zawierający kolumnę identyfikatora w sekcji Odpowiedzi na:
Najlepsze pakiety dla modeli Coxa z zmiennymi zmiennymi w czasie
Załóżmy również, że zmienne towarzyszące są zmienne w czasie i istnieje dokładnie jedna zmienna cenzora (tj. Zdarzenie), która jest binarna.
pytania
1) Jeśli w powyższym linku odpowiedź nie zostanie podana jako parametr w wywołaniu coxph (), to czy wyniki powinny być takie same, jak uwzględnienie klastra (ID) jako parametru w coxph ()?
Próbowałem wyszukać dokumentację, ale następujące elementy nie wydają się wyraźnie adresować (1): https://stat.ethz.ch/pipermail/r-help//2013-July/357466.html
2) Jeśli odpowiedź na (1) brzmi „nie”, to (matematycznie) dlaczego? Wydaje się, że klaster () w coxph () szuka korelacji między podmiotami zgodnie z podrozdziałem „klaster” na stronie. 20 at
https://cran.r-project.org/web/packages/survival/survival.pdf
3) Niejasne pytanie: w jaki sposób coxph () z powtarzanymi pomiarami wypada w porównaniu z metodami regresji kruchego pakietu R.
Dodatki
Następujące wskazówki dotyczące korzystania z klastra (ID):
Czy istnieje wersja testu logrank uwzględniająca powtarzane pomiary?
podobnie jak:
https://stat.ethz.ch/pipermail/r-help//2013-July/357466.html
Podejście GEE: dodaj „+ klaster (podmiot)” do instrukcji modelu w Coxph Podejście modeli mieszanych: dodaj „+ (1 | podmiot)” do statystyki modelu w Coxme.
Z góry dziękuję!
Oto odpowiedź z
survival
winiety paczki, która okazała się pomocna - jest powiązana z pierwszą odpowiedzią na pierwsze pytanie, do którego dałeś link:Najlepsze pakiety dla modeli Coxa z zmiennymi zmiennymi w czasie
Odnoszą się do konfiguracji danych w formie długiej lub danych z powtarzającymi się wpisami dla tematów.
Podają przykład
sugerując, że jeśli podasz dwa razy (początek i koniec okresu)
Surv
zamiast zamiast jednego,coxph()
rozliczy resztę.źródło