Я пытаюсь проанализировать некоторые данные, используя модель смешанного эффекта. Собранные мной данные показывают вес некоторых молодых животных с различным генотипом с течением времени.
Я использую предлагаемый здесь подход: https://gribblelab.wordpress.com/2009/03/09/repeated-measures-anova-using-r/
В частности я использую решение № 2
Так у меня что то типа
require(nlme)
model <- lme(weight ~ time * Genotype, random = ~1|Animal/time,
data=weights)
av <- anova(model)
Теперь я хотел бы провести несколько сравнений. Используя multcomp
я могу сделать:
require(multcomp)
comp.geno <- glht(model, linfct=mcp(Genotype="Tukey"))
print(summary(comp.geno))
И, конечно, я мог сделать то же самое со временем.
У меня есть два вопроса:
- Как мне использовать,
mcp
чтобы увидеть взаимодействие между временем и генотипом? Когда я бегу,
glht
я получаю это предупреждение:covariate interactions found -- default contrast might be inappropriate
Что это означает? Могу ли я безопасно проигнорировать это? Или что я должен сделать, чтобы избежать этого?
РЕДАКТИРОВАТЬ: я нашел этот PDF, который говорит:
Поскольку в этом случае невозможно автоматически определить интересующие параметры, mcp () в multcomp по умолчанию будет генерировать сравнения только для основных эффектов, игнорируя ковариаты и взаимодействия . Начиная с версии 1.1-2, можно указать усреднение по терминам взаимодействия и ковариатам с использованием аргументов взаимодействия_верация = ИСТИНА и covariate_average = ИСТИНА соответственно, тогда как версии старше 1.0-0 автоматически усредняются по условиям взаимодействия. Тем не менее, мы предлагаем пользователям вручную написать набор желаемых контрастов.Это следует делать всякий раз, когда возникает сомнение в том, что измеряют контрасты по умолчанию, что обычно происходит в моделях с терминами взаимодействия более высокого порядка. Мы обращаемся к Hsu (1996), глава ~ 7, и Searle (1971), глава ~ 7.3, для дальнейшего обсуждения и примеров по этому вопросу.
У меня нет доступа к этим книгам, но, может быть, кто-то здесь есть?