我使用. lmer_ 使用我可以获得模型中包含的固定效果的平方和。不幸的是，没有提供总平方和 (SSTO) 和误差平方和 (SSE)。有没有办法获得这些值？我在下面包含了一个人为的数据集，用于检查解决方案。lme4RlmerTest::anova
# libraries
library(lme4)
library(lmerTest)
# factors
fixed <- c(rep('a', 4), rep('b', 4))
random <- c('r1', 'r1', 'r2', 'r2', 'r3', '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