lmer软件包的 lme4函数给出了干扰参数的轮廓似然间隔。例如，
fm1 <- lmer(Reaction ~ Days + (Days | Subject), sleepstudy)
pp <- profile(fit,"theta_")#profiling of the deviance
ci <- confint(pp)
该函数通过以下代码提供估算值的wald统计：
res <- summary(fit)
se <- coefficients(res)[,2]
是否也可以获得回归系数的轮廓似然区间？