Мне нужно использовать пакет metafor
для анализа некоторых вложенных данных и представления результатов довольно наивной статистике, поэтому я надеюсь получить некоторую обратную связь по модели и способам суммирования результатов пакета.Данные получены из 50 единиц в 15 палатах в 5 больницах, и у меня есть только данные, соответствующие среднемесячным значениям каждой переменной интереса.Вот текущая модель:
rma.mv(datES, datESvar, mods = cbind(avgLoS, avgAge, avgEdu, avgMedComp), random = ~ 1 | Hospital/Ward/Unit, struct = "AR", data = dat)
Вкратце, «avgLoS» соответствует средней продолжительности пребывания пациентов в каждом отделении, «avgAge» и «avgEdu» относятся к среднему возрасту.и образование каждого подразделения, и «avgMedComp» относится к средней приверженности к лечению пациентов в каждом отделении.Я буду использовать функцию robust()
, чтобы учесть зависимости в переменных измерениях.
У меня есть вопрос относительно моделирования продольных эффектов с помощью модераторов, так как переменная для средней продолжительности пребывания также собираетсябыть связанным с последовательными ежемесячными измерениями.metafor
Рассчитывает ли автокорреляцию между фиксированными эффектами, а также случайными эффектами?
Я впервые использую пакет metafor
, поэтому я несколько незнаком с некоторыми функциями, включенными впакет для обобщения и построения модели.Есть ли у кого-нибудь рекомендации по представлению результатов таких моделей, которые делают их более интуитивно понятными или более удобными для восприятия, особенно для тех, у кого ограниченный статистический фон?