如何引导混合模型的 R 平方?

Question

我正在尝试为混合效应模型获取引导R^2。由于已经只有一种解决方法来获取条件和边际 R^2，因此我尝试根据 statmethods 给出的示例来引导这些统计数据，用于引导单个统计数据。代码有效，但偏差和标准误差始终为零。

library(lme4)
library(boot)
data(Dyestuff, package = "lme4")
model <- lmer(Yield ~ 1|Batch, data=Dyestuff)
summary(model)
r.squaredGLMM(model)
rsq <- function(formula, data, indices) {
d <- data[indices,] 
model.fit <- lmer(Yield ~ 1|Batch, data=Dyestuff)
fit.r.squared <- r.squaredGLMM(model.fit)
return(summary(fit.r.squared[,2]))
}
set.seed(101)
results <- boot(data=Dyestuff, statistic=rsq,
R=1000, formula=Yield ~ 1|Batch)
results

Bootstrap Statistics :
original  bias    std. error
t1* 0.4184874       0           0
t2* 0.4184874       0           0
t3* 0.4184874       0           0
t4* 0.4184874       0           0
t5* 0.4184874       0           0
t6* 0.4184874       0           0

当我引导模型时，条件和边际 R^2 不应该也改变吗？有没有其他方法可以获得自举条件和边际 R^2？

Answer 1

rsq <- function(formula, data, indices) {
d <- data[indices,] 
model.fit <- lmer(Yield ~ 1|Batch, data = d)
fit.r.squared <- r.squaredGLMM(model.fit)
return(fit.r.squared[,2])
}