Как согласовать регрессионную модель повторных измерений со случайным перехватом - PullRequest
Новая
       3

Как согласовать регрессионную модель повторных измерений со случайным перехватом

0 голосов
/ 12 февраля 2019

Основная проблема, с которой я столкнулся, заключается в том, чтобы убедиться, что я правильно указываю модель повторных измерений со случайным эффектом для группы.Там нет группы лечения.Скорее,

У меня есть набор данных, в котором участники отвечали на опрос в течение одного года и снова в последующем опросе в следующем году.Участники принадлежат к разным группам, каждая из которых вела программу повышения качества между двумя исследованиями.Заинтересованные стороны хотят знать, произошло ли изменение между 1 годом и 2 годом. Участники вкладываются в группы.Группу считают случайной, потому что меня интересует, прежде всего, изменение во времени (группы могли управлять программами по-разному).(Одна из проблем, с которыми я сталкиваюсь, заключается в том, что раньше у меня был SAS, и я был более знаком с модельными процедурами в этом программном обеспечении.)

Данные выглядят так (всего насчитывается около 108 подходящих пар). Здесьнебольшой набор данных: 

mydata <- data.frame(
  Participant = c(21, 22, 23, 24, 25, 26, 27), 
  Year = c(rep(c("Year1", "Year2"), each=7)), 
  Group = c(rep(c("Group1", "Group2", "Group3", "Group4", "Group5", "Group6", "Group7"), 
                  each = 1)), 
  Outcome = c(1.2, 1.45, 2.1, 3.0, 3, 3.5, 2.8, 
                1.5, 1.9, 2.5, 3.6, 3.8, 4.1, 3.8 
            ))

Я пытался подобрать модель таким образом: 

 model <- lmer(Outcome ~ Year  +(1|Group:Year), mydata)

Однако мне не ясно, подходит ли мне подходящая модельдля моей настройки, потому что я не уверен, что понимаю термин (1 | Group: Year).Я думаю, это означает, что каждая группа имеет случайный перехват для каждого значения года, но я не уверен, что это более уместно, чем просто оценка одного случайного эффекта, который не является вложенным во времени.В SAS я смог увидеть оценки конкретных параметров для каждого уровня случайного эффекта, но coef(model) производит только оценки самих перехватов без вывода, что указывает, отличалась ли конкретная группа от смоделированного среднего.Кроме того, я не вижу, как указать, что для каждого значения Участника есть несколько измерений.

Буду признателен за советы о том, как правильно использовать lmer - это мощный, но новый для меня.

...