Я запускаю модель смешанных эффектов в R с использованием пакета 'nlme' и пытаюсь визуализировать предсказанные результаты с помощью пакета 'visreg'. В моей модели SEA c в качестве переменной отклика объясняется взаимодействием между «точкой» (числовой) и «местом» (факториал) и ковариатой «богатство» (числовой) с «рекой» в качестве случайного эффекта, а именно:
mod1 <- lme(SEAc~point*site+richness, random =~1|river, na.action=na.omit, data = alldata)
Я хотел построить регрессию точки SEAc ~ для каждого сайта на основе прогнозов интерактивной модели, используя следующий код visreg:
visreg(fit = SEAc, xvar = "point", by = "site", re.form = ~ (1|river),
cond = list(richness = mean(richness),na.rm = TRUE)), scale = "response")
Однако я продолжаю получать это сообщение об ошибке :
Error in UseMethod("family"):no applicable method for 'family' applied to an object of class "lme"
Я не понимаю, что этот код раньше работал отлично и до сих пор работает на других машинах. Возможно ли, что обновление «visreg», «nlme» или «R» помешало работе функции? Я попытался использовать предыдущие версии всех трех безрезультатно.
Любые советы относительно того, что происходит, или почему можно получить эту ошибку, используя visreg с nlme, были бы очень признательны!