使用bootMer进行自举的预测值？

Question

我正在尝试为lmer中开发的混合模型运行bootstrap模拟。该模型如下所示：

model1<-lmer(Hg~A+B+C+D+(1|K),data=mydata)

我写了这样的bootMer命令：

bootfit1<-bootMer(model1,FUN=function(x)predict(x,re.form=NA),nsim=1000)

有效！但是当我使用实际和模型预测值绘制bootfit1值时，bootfit1值比实际值和预测值更接近0。在我的数据中Hg在技术上是正态分布的，但实际上偏向0。

所以这是我的问题：有没有办法控制从实际数据中抽取引导程序中每个模拟数据集的数量以及纯模拟的数量？我想知道是否因为我的数据偏差，引导程序模拟没有获得远高于0的许多值。