В соответствии с MuMIn
R
21
#load packages
library(lme4)
library(MuMIn)
#Fit Model
m <- lmer(mpg ~ gear + disp + (1|cyl), data = mtcars)
#Determine R2:
r.squaredGLMM(m)
R2m R2c
0.5476160 0.7150239
Выход для функции r.squaredGLMM
обеспечивает:
R2m : предельное значение R в квадрате, связанное с фиксированными эффектами
Условное значение R2c R2, связанное с фиксированными эффектами плюс случайные эффекты.
Примечание: комментарий к связанному сообщению в блоге предполагает, что альтернативный подход, вдохновленный Nakagawa & Schielzeth, разработанный Джоном Лефчеком (с использованием sem.model.fits
функции в piecewiseSEM
пакете), дал идентичные результаты. Итак, у вас есть варианты: р.
Я не тестировал эту последнюю функцию, но я тестировал r.squaredGLMM()
функцию в MuMIn
пакете и поэтому могу подтвердить, что она все еще функционирует сегодня (2018).
Что касается обоснованности этого подхода, я оставляю чтение Nakagawa & Schielzeth (2013) (и последующую статью Johnson 20142) вам решать.
1: Накагава, S., и Schielzeth, H. 2013. Общий и простой метод получения R2 из обобщенных линейных моделей смешанных эффектов. Методы в экологии и эволюции 4 (2): 133-142.
2: Джонсон, PCD 2014 Расширение модели R2GLMM компании Nakagawa & Schielzeth для моделей со случайными уклонами. Методы в экологии и эволюции 5: 44–946.