Złożyłem artykuł, który został odrzucony z powodu niewłaściwego sposobu przeprowadzenia analizy przeżycia. Sędzia nie pozostawił żadnych innych szczegółów ani wyjaśnień poza: „analiza przetrwania trendów czasowych wymaga bardziej wyrafinowanych sposobów cenzury”.
Pytanie:
Czy w ostatnich dziesięcioleciach zmniejszono nadmierne ryzyko śmierci wśród palaczy?
Dane:
25.000 palaczy w Niemczech. Zostali oni zapisani do kohorty w dowolnym momencie w latach 1995–2014. Każdy palacz został dopasowany (w momencie rejestracji) do płci i wieku dopasowanych kontroli z populacji ogólnej (która nie paliła). Mam dokładny czas śmierci dla wszystkich, którzy zmarli podczas całego okresu studiów. Ci, którzy nie umarli podczas obserwacji, zostaną ocenzurowani. Badanie ma na celu zbadanie nadmiernego ryzyka śmierci wśród palaczy każdego roku od 1995 do 2014 roku.
Celem jest obliczenie:
- Wskaźniki śmiertelności wśród palaczy i osób niepalących co roku i zbadanie tych trendów
- nadmierne ryzyko śmierci wśród palaczy każdego roku (lub okresu kilku kolejnych lat).
Jak należy analizować dane? Przypomnij sobie, że ktoś, kto jest uwzględniony w 1998 r., Może umrzeć w 2015 r. Czy właściwe podejście do korzystania z formatu procesu liczenia z uruchomieniem i zatrzymaniem jest aktualizowane każdego roku?
Takie podejście nie podobało się sędziemu:
Wskaźniki zapadalności obliczono za pomocą regresji Poissona. Uwzględniliśmy czas obserwacji jako przesunięcie w modelu oraz uwzględniliśmy wiek, płeć, status palenia i okres kalendarzowy (łącząc dwa kolejne lata) jako predyktory w modelu. Następnie obliczono stawki za 1000 osobolat za pomocą funkcji przewidywania () R. Przesunięcie (czas obserwacji) to całkowity czas obserwacji (dni) osób od rejestracji.
Zastosowano model Coxa do oszacowania względnego ryzyka dla palaczy w każdym okresie od początku do końca badania. Dla uproszczenia porównaliśmy współczynnik ryzyka w pierwszym okresie ze współczynnikiem ryzyka w okresie końcowym.
Problemy: - osoba (wraz z jej kontrolą) może zostać objęta w 1998 r., A zatem należy do tej grupy kalendarza, ale doznaje zdarzenia w 2006 r. - Jak należy rozłożyć dane do analizy regresji Poissona i Coxa? Proces liczenia dla Coxa? Jaki jest czas rozpoczęcia i zakończenia? - Jak można oceniać trendy w tej sytuacji?
Kilka wyjaśnień: załóżmy, że pacjenta zaobserwowano po raz pierwszy 15 czerwca 1998 r. I przeżył zdarzenie 31 grudnia 1998 r., Wartość naszej zmiennej czasowej dla tego pacjenta wynosi 182,5 z 730 możliwych dni, ponieważ okres składa się z 2 kolejnych lat. Maksymalny czas obserwacji w każdym okresie wynosi 730 dni.
Kiedy pacjent jest obserwowany w jednym przedziale czasowym, ale ocenzurowany (tj. Albo doświadczony i zdarzenie, albo porzucony) w innym przedziale czasowym, czy liczbę dni obserwowanych należy dodać do następnego okresu, czy co?
Zatem głównym problemem jest obsługa czasu obserwacji i roku kalendarzowego (który jest wykorzystywany jako zmienna kategoryczna, składająca się z dwóch kolejnych lat).
Odpowiedzi:
Z powyższego wynika kilka możliwości modelu Coxa:
OBRÓBUJ WZGLĘDNE ZMIANY W ZAGROŻENIU U PALĄCYCH I OSÓB NIEPALĄCYCH ODDZIELNIE : jedna obserwacja na osobę; oblicz czas obserwacji (niezależnie od tego, kiedy nastąpi cenzura / zdarzenie), a następnie wykorzystaj wszystkich pacjentów (od 1995 do 2014 r.) w modelu, użyj okresu czasu jako zmiennej kategorialnej i ustaw jeden z okresów jako wartość odniesienia.
źródło
Mimo, że przesadzanie z tajemniczymi komentarzami recenzenta jest niebezpieczne, zgaduję, że sprzeciw dotyczy tego, czy cenzura ma charakter informacyjny.
Interpretacja modeli przeżycia zazwyczaj opiera się na założeniu, że osoba ocenzurowana w czasie jest reprezentatywna dla wszystkich podmiotów, które przeżyły do czasu po wejściu do badania. (Brzmienie dostosowane od tego wstępu do analizy przeżycia). W takim razie cenzura nie ma charakteru informacyjnego.T T
W twojej analizie cenzurowani byli jednak ci, którzy przeżyli do 2014 r. Jeśli uważasz, że w ciągu ostatnich 20 lat nastąpiła zmiana nadmiernego ryzyka śmierci z powodu palenia (lub nawet jeśli nastąpiły równoległe zmiany w śmiertelności dla obu grup), wówczas osoby cenzurowane mogą nie być reprezentatywne dla tych, którzy przeżyli w tym samym czasie, ale weszli do badania wcześniej. Według twojej hipotezy cenzura może mieć charakter informacyjny.
Możliwe, że szczegóły projektu twojej analizy uniknęły tego problemu, ale nie było to jasne w recenzowanym manuskrypcie. A może recenzentowi nie spodobało się badanie z kilku dodatkowych powodów i stwierdził, że jest to sposób na odrzucenie go, którego redaktor by nie kwestionował. Niemniej jednak wydaje się, że jest to potencjalny sprzeciw wobec sposobu, w jaki analizowałeś te dane i powinieneś upewnić się, że są one właściwie przetwarzane. (To jest poza moją osobistą wiedzą; inni na tej stronie mogą mieć wskazówki, jak postępować. Bardziej precyzyjny tytuł tego pytania, z bardziej szczegółowymi informacjami na temat projektu i analizy, może uzyskać bardziej pomocne odpowiedzi).
Z twojego pytania i wyjaśniającego komentarza nie wynika jasno, że analizy Coxa dodają coś użytecznego do prostego modelowania wskaźników śmiertelności rocznie (lub w odstępach 2-letnich). Ponadto, twoja hipoteza wydaje się sugerować, że zagrożenia nie są proporcjonalne w czasie między osobami niepalącymi i palącymi, co stanowi podstawę standardowych analiz Coxa. Jeśli interesuje Cię różnica śmiertelności między palaczami a osobami niepalącymi w zależności od roku kalendarzowego, jest to najprostsza miara do modelowania (chociaż być może będziesz musiał wziąć pod uwagę zakładane wzbogacenie osób niepalących w badanej próbie gdy ich pasujący palacze zginą).
źródło