У меня есть вопрос, касающийся того, как я могу получить остатки из линейной модели смешанного эффекта (например, используя lme из nlme в R), чтобы построить данные после того, как модель lme удалила определенные переменные, которые не представляют интереса. В приведенном ниже коде я показал модель, которую я использую. Переменная, о которой я забочусь, это grp, а переменные, которые я пытаюсь вычеркнуть (потому что они являются неприятными переменными), - это такие вещи, как пол и возраст с фиксированным и случайным эффектом сайта. Я хотел бы по существу построить график y_var относительно переменной grp, но после того, как модель учитывает пол, возраст и сайт. Может кто-нибудь предложить, как я мог бы реализовать это в R?
lme_model = lme(fixed = y_var ~ grp + sex + age, random = ~1|site, data = my_df)