Bez modelowania wielopoziomowego, jak poradzić sobie z replikacją w ramach badania w metaanalizie, gdzie badanie jest jednostką replikacji?

13

Opis badania:

Zauważyłem częsty błąd w metaanalizach w odniesieniu do obsługi replikacji w ramach badania. Nie jest dla mnie jasne, czy błąd unieważnia badania po sformułowaniu założeń. Jednak, jak rozumiem, założenia te naruszają podstawową przesłankę statystyki.

Na przykład, badania Badania efektów chemicznej w reakcji .YXY

Analiza jest przeprowadzana na logarytmicznym stosunku odpowiedzi: stosunek leczenia (w obecności ) do kontroli (bez ): X Y 0 XY+XXY0X

R=ln(Y+XY0)

Niektóre z tych badań zawarte w meta-analizy zawierają wiele zabiegów, na przykład różne poziomy lub chemiczne formy . Dla każdego leczenia istnieje inna wartość , chociaż zawsze używa tej samej wartości .R R Y 0XRRY0

Metody określają:

odpowiedzi na różne terapie (poziomy i formy ) w ramach jednego badania uznano za niezależne obserwacje.X

Pytania:

  • Czy to nie jest pseudoreplikacja?
  • Czy jest to niewłaściwe, nawet jeśli w metodach stwierdzono naruszenie niezależności?
  • Jaki byłby łatwy sposób (np. W ramach możliwości prostego pakietu oprogramowania do metaanalizy) obsługi replikacji badania?

Wstępne przemyślenia:

  • Podsumuj wyniki każdego badania, np. Biorąc średnią odpowiedź
  • Wybrać tylko jeden zabieg z każdego badania w oparciu o kryteria a priori (np. Najwyższa dawka, pierwszy pomiar)?

Czy są jakieś inne rozwiązania?

David LeBauer
źródło
To tylko szybkie przypuszczenie, ale możesz chcieć sprawdzić Kim / Becker 2010: Stopień zależności między wielkościami efektów wielokrotnego leczenia ; Nie przeczytałem tego artykułu, ale może być związany z twoim pytaniem.
Bernd Weiss,
Czy metaanaliza naprawdę po prostu uśrednia wszystkie wartości różnic R? Wydaje się to dość dziwne w porównaniu do np. Próby meta-regresji - w takim przypadku różnice między R na różnych poziomach X mogą być tym, co chcesz połączyć w różnych badaniach.
gość
@ gość tak, naprawdę są; Interesujące byłoby, w jaki sposób różne poziomy X wpływają na R, ale pytanie brzmi po prostu: „czy jest efekt X”? W tym kontekście może istnieć ograniczona moc do testowania wpływu X na R (reakcja ekosystemu na dodatek składników odżywczych), ze względu na różnorodność metod i warunków badań.
David LeBauer
1
Masz rację, to jest problem. Nie tyle z oszacowaniami punktowymi, ale miary precyzji (tj. Błędy standardowe) będą zbyt małe; ignoruje wielokrotne wykorzystanie danych grupy kontrolnej. Nie powinno to jednak stanowić wiadomości dla nikogo w metaanalizie. Powyższy artykuł Kim / Becker jest w zasadzie powtórnym oświadczeniem Glesera i Olkina (1994). Rozmiary efektów zależne od stochastycznie. W Cooper & Hedges (Eds), Podręcznik syntezy badań (s. 339–355). Ta książka jest standardowym tekstem w tej dziedzinie, wierzę teraz w drugie wydanie.
gość

Odpowiedzi:

3

Tak, jest to problem, ponieważ w odpowiedziach należy wziąć pod uwagę zależność próbkowania (chociaż czasami efekt może być nieistotny i cały czas naruszamy założenie podczas przeprowadzania analiz statystycznych). Istnieją metody, aby sobie z tym poradzić, jednym podejściem jest uwzględnienie kowariancji między powiązanymi eksperymentami (bloki nie przekątne) w macierzy wariancji błędu-kowariancji (patrz np. Hedges i in., 2010). Na szczęście przy stosunkach logów jest to dość łatwe. Można uzyskać przybliżone kowariancje między eksperymentami, ponieważ wariancja (var) log R wynosi (jeśli Yx i Y0 są niezależnymi grupami): log Yx - log Y0, aby postępować zgodnie z zapisem w pytaniu, Yx odnosi się do grupy eksperymentalnej, a Y0 Grupa kontrolna. Kowariancja (cov) między dwiema wartościami (np. Leczenie 1 i leczenie 2) dla log R wynosi cov (loge Yx_1 - log Y0, log Yx_2 - log Y0), który jest równy var (log Y0), i jest obliczany jako SD_Y0 / (n_Y0 * Y0), gdzie SD_Y0 jest standardowym odchyleniem Y0, n_Y0 jest rozmiarem próbki w traktowaniu kontrolnym, a Y0 wynosi wartość w leczeniu kontrolnym. Teraz możemy podłączyć całą macierz wariancji-kowariancji do naszego modelu, zamiast używać tylko wariancji (ei), co jest klasycznym sposobem przeprowadzania metaanalizy. Przykład tego można znaleźć wLimpens i in. 2011 przy użyciu pakietu metahdep w R (na bioprzewodniku) lub Stevens i Taylor 2009 dla Hedge's D.

Jeśli chcesz to uprościć, pokusiłbym się zignorować problem i spróbować ocenić efekt zależności próbkowania (np. Ile zabiegów jest w ramach badań? Jak zmieniają się wyniki, jeśli używam tylko jednego leczenia? Itp.) .

GGeco
źródło
2

Tak, to jest problem.

Tak, jest nieodpowiedni, chociaż przynajmniej jest przejrzysty na temat tego, co robi (dostaje punkty za przejrzystość, ale nadal nie jest zadowalający).

Wątpię, czy istnieje „łatwy sposób”, aby to naprawić. Nie wiem wiele na temat podejść do metaanalizy, ale jeśli istnieje specjalne oprogramowanie do metaanalizy, a badania tego typu są tworzone przy jego użyciu i publikowane, może to być wspólne podejście. Każda z proponowanych odpowiedzi traci pewną ziarnistość informacji z każdego badania (tj. Odwrotny problem tego, co zrobili wydawcy).

Oczywistym rozwiązaniem jest model z efektami mieszanymi (tj. Wielopoziomowy) z badaniem jako czynnikiem losowym. Sugerowałbym użycie do tego specjalistycznego pakietu statystycznego, jeśli oprogramowanie do metaanalizy nie może tego zrobić. Nadal możesz używać oprogramowania do metaanalizy do przechowywania i przetwarzania danych, a także eksportować dane do R, Stata lub SAS w celu przeprowadzenia analizy.

Peter Ellis
źródło
Zastanawiałem się nad próbami klinicznymi i zastanawiałem się, czy nie jest to właściwe w sytuacji, w której krzywa dawka-odpowiedź jest wynikiem, ponieważ wtedy można porównać funkcje krzywej. Czy to jest możliwe?
Michelle,
Nie sądzę, żeby miało to duży wpływ na problem polegający na tym, że kilka wyników jednego badania będzie w jakiś sposób skorelowanych, a zatem nie będą to „nowe” informacje. Ale metaanaliza funkcji krzywych z pewnością byłaby możliwa, o ile w jakiś sposób kontrolujesz korelację między różnymi oszacowaniami tych krzywych. Jeśli wszystkie mają tę samą formę i jest to tylko kwestia oszacowania parametrów, powinno to być możliwe.
Peter Ellis,
@Michelle Zgadzam się z Peterem: jeśli podsumowujesz parametry krzywej, to otrzymujesz oszacowanie jednego parametru z każdej krzywej, i powinno być dobrze.
Abe