Sparowana, ANOVA z powtarzanymi pomiarami czy model mieszany?

9

Poproszono mnie o przeanalizowanie niektórych danych z badania klinicznego dotyczącego dwóch metod pomiaru ciśnienia krwi. Mam dane od 50 osób, z których każda zawiera od 2 do 57 miar przy użyciu każdej metody.

Zastanawiam się, jak najlepiej postępować.

Oczywiście potrzebuję rozwiązania, które uwzględni fakt, że miara ciśnienia krwi jest sparowana (dwie metody mierzone jednocześnie), a także zmienną towarzyszącą w czasie (z różną liczbą obserwacji na pacjenta), a także uwzględniającą intra i inter- zmienność pacjenta.

Zastanawiałem się nad tym, jak przerobić to na ANOVA z powtarzanymi pomiarami, ale wydaje mi się, że może to być model mieszany.

Byłbym wdzięczny za wszelkie pomocne porady, które możesz zaoferować.

Jestem kompletnym nowicjuszem R, ale jestem bardzo podekscytowany rozwijaniem umiejętności i mam umiarkowane doświadczenie w Stacie, więc zawsze mogę na tym polegać.

Sam
źródło

Odpowiedzi:

11

Nie sądzę, że możesz łatwo zrobić to, co chcesz zrobić z RM-ANOVA, ponieważ liczba powtórzeń nie jest taka sama dla wszystkich przedmiotów. Uruchamianie modeli z efektami mieszanymi jest bardzo łatwe w R. W rzeczywistości, poświęcając trochę czasu na naukę podstaw i poleceń, otworzy przed Tobą wiele możliwości. Uważam też, że modelowanie mieszane jest o wiele prostsze w użyciu i bardziej elastyczne i prawie nigdy nie muszę bezpośrednio wykonywać RM-ANOVA. Na koniec należy wziąć pod uwagę, że w przypadku modelowania mieszanego można również uwzględnić strukturę kowariancji reszt (RM-ANOVA po prostu zakłada strukturę diagonalną), co może być ważne w wielu zastosowaniach.

Istnieją dwa główne pakiety liniowego modelowania mieszanego w R: nlmei lme4. Do lme4paczki jest bardziej nowoczesna, które jest idealne dla dużych zbiorów danych, a także dla przypadków zajmują się danymi klastrowych. Nlmejest starszym pakietem i jest w większości przestarzały na korzyść lme4. Jednak w przypadku projektów z powtarzanymi pomiarami jest to wciąż lepsze niż, lme4ponieważ nlmepozwala tylko modelować strukturę kowariancji reszt. Podstawowa składnia nlmejest bardzo prosta. Na przykład:

fit.1 <- lme(dv ~ x + t, random=~1|subject, cor=corCompSymm())

Tutaj modeluję związek między zmienną zależną dva czynnikiem xi zmienną zależną od czasu t. Subjectjest efektem losowym i zastosowałem złożoną strukturę symetrii dla kowariancji reszt. Teraz możesz łatwo uzyskać niesławne wartości p poprzez:

anova(fit.1)

Na koniec mogę zasugerować, aby przeczytać więcej o nlme, korzystając z jego ostatecznego przewodnika, modeli efektów mieszanych w S i S-Plus . Innym dobrym źródłem informacji dla początkujących jest liniowy model mieszany - praktyczny przewodnik z wykorzystaniem oprogramowania statystycznego, który kompiluje wiele przykładów różnych aplikacji modelowania mieszanego z kodem w R, SAS, SPSS itp.

AlefSin
źródło
Dzięki Alef - te dwie referencje są niesamowite - podobnie jak powyższe Wolfa. Zastanawiam się, czy mogę nieco rozszerzyć moje pytanie dotyczące struktury modelu. Nie mogę rozpoznać dv !! Mam dwa zestawy pomiaru BP (dwie metody), a także identyfikator pacjenta i czas obserwacji. Jak mogę modelować różnicę między dwoma pomiarami BP (analogicznie do testu t dla jednej próbki, że różnica = 0)? Przepraszam, że cię prześladuję - zabiorę się teraz do czytania!
Sam
Nie martw się wszyscy - myślę, że to rozgryzłem !!! Miałem swoje dane w niewłaściwym formacie. Kiedy w końcu wymyśliłem i zmanipulowałem go do długiego formatu, wszystkie te posty miały o wiele więcej sensu !! Jeszcze raz dziękuję.
Sam
Cieszę się, że to rozgryzłeś. Wydaje się, że co do zasady większość pakietów w R działa z danymi w długim formacie.
AlefSin
1

Jeśli szukasz RM-ANOVA z modelem mieszanym za pomocą R. Możesz to sprawdzić http://blog.gribblelab.org/2009/03/09/repeated-measures-anova-using-r/ Istnieją świetne przykłady pokazujące, jak stosować model mieszany do osiągnięcia RM-ANOVA.

Bazując na moim doświadczeniu, SAS jest lepszym narzędziem do radzenia sobie z modelem mieszanym. Jeśli używasz SAS, możesz sprawdzić pomoc SAS „Proc Mixed” dla RM-ANOVA.

Tu.2
źródło