如何解决如何获得lme4中的所有平方和
我使用 lmer
中 lme4
包的 R
拟合混合模型。使用 lmerTest::anova
我可以获得模型中包含的固定效应的平方和。不幸的是,没有提供总平方和 (SSTO) 和误差平方和 (SSE)。有没有办法获得这些值?我在下面包含了一个人为的数据集,用于检查解决方案。
# libraries
library(lme4)
library(lmerTest)
# factors
fixed <- c(rep('a',4),rep('b',4))
random <- c('r1','r1','r2','r3','r4','r4')
# constructing response variable
fixed_effect <- ifelse(fixed == 'a',1,2)
random_effect <- as.numeric(gsub('r','',random))
error <- rnorm(length(fixed))
y <- fixed_effect + random_effect + error
# modeling
mod <- lmer(y ~ fixed + (1|random))
anova(mod) # only provides sums of squares for fixed effects
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。