Функция
lmer пакета lme4 дает интервал вероятности профиля параметров помех.Например,
fm1 <- lmer(Reaction ~ Days + (Days | Subject), sleepstudy)
pp <- profile(fit,"theta_")#profiling of the deviance
ci <- confint(pp)
Функция выдает статистику Вальда оценок по следующим кодам:
res <- summary(fit)
se <- coefficients(res)[,2]
Можно ли также получить интервал вероятности профиля коэффициента регрессии?