小编Vic*_*Vic的帖子

我想使用新lme4包的新bootMer()功能(当前的开发人员版本).我是R的新手,不知道我应该为其FUN参数编写哪个函数.它说它需要一个数字向量,但我不知道该函数将执行什么.所以我有一个混合模型公式,它被转换为bootMer(),并且有许多重复.所以我不知道那个外部函数是做什么的？它应该是引导方法的模板吗？他的bootMer中是否已经实现了自举方法？那么为什么他们需要外部的"利益统计"呢？我应该使用哪种统计数据？

以下语法是否适用？R继续产生错误,产生FUN必须是数字向量.我不知道如何将估计与"适合"分开,甚至我应该首先做到这一点？我可以说我迷失了那种"有趣"的说法.另外我不知道我应该使用变量"Mixed5"传递混合模型glmer()公式还是应该传递一些指针和引用？我在示例中看到X(bootMer()的第一个参数是*lmer()对象.我想编写*Mixed5但是它呈现错误.

非常感谢.

我的代码是:

library(lme4)
library(boot)

(mixed5 <- glmer(DV ~ (Demo1 +Demo2 +Demo3 +Demo4 +Trt)^2 
                 + (1 | PatientID) + (0 + Trt | PatientID)
                 , family=binomial(logit), MixedModelData4))


FUN <- function(formula) {
  fit <- glmer(DV ~ (Demo1 +Demo2 +Demo3 +Demo4 +Trt)^2 
               + (1 | PatientID) + (0 + Trt | PatientID)
               , family=binomial(logit), MixedModelData4)
  return(coef(fit))
}

result <- bootMer(mixed5, FUN, nsim = 3, seed = NULL, use.u = FALSE,
        type = c("parametric"),
        verbose = T, .progress = "none", PBargs = …