Я пытаюсь выполнить расчет мощности для моей продольной модели Lmer в R, используя функцию пакета longpower lmmpower. У меня есть 2 вопроса. Мои данные измеряются в пяти временных точках («время» = 4, 8, 16, 26, 52), и моя моя модель:
L4 <-lmer (результат ~ время + X1 + X2 + наркотик + наркотик * время + X2 * время + время ^ 2 + (время | ID), данные) </p>
Цель моей модели - выяснить, имеет ли время, связанное с наркотиками, существенные взаимодействия (т.е. влияет ли наркотик на траекторию). Я запутался насчет «т» - или вектора времени в амплитуде. Правильно ли я определил это ниже?
lmmpower (L4, pct.change = 0,30, t = c (4, 8, 16, 26, 52), мощность = 0,80)
Когда я запускаю их, это говорит мне, что n = 19,999 ... Так что мне нужно 20 человек в каждой группе, чтобы обнаружить изменение результата на 30%, верно? Но я думаю, что я делаю это неправильно, так как я специально хочу знать, как бы я обнаружил значительную разницу для взаимодействия лекарства * время. Есть ли способ, которым я могу указать это?
Спасибо!