我建立了一个广义的加性混合效应模型,它具有固定效应和随机拦截效应(即一个分类变量).运行模型后,我可以使用每个类别提取随机截取ranef(m1$lme)$x[[1]].但是,当我尝试使用时提取随机效果的标准误差时se.ranef(m1$lme),该功能不起作用.其他尝试使用se.ranef(m1)和se.ranef(m1$gam)不工作.我不知道这是否因为这些函数只适用于lmer类的模型？
任何人都可以帮助我,以便我可以从"gamm"课程中抽出我随机拦截的标准错误吗？我想使用随机截距和标准误差来绘制我的gamm模型的最佳线性无偏预测器.
我最初的模型是:gamm(y ~ s(z), random = list(x = ~1), data = dat).
library(mgcv)
library(arm)
example <- gamm(mag ~ s(depth), random = list(stations = ~1), data = quakes)
summary(example$gam)
#Family: gaussian
#Link function: identity
#Formula:
# mag ~ s(depth)
#Parametric coefficients:
# Estimate Std. Error t value Pr(>|t|)
#(Intercept) 5.02300 0.04608 109 <2e-16 ***
# ---
# Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ …