我有一个关于JoséA.Sánchez-Espigares和JordiOcaña的R演示文稿的问题,题为“混合模型的自助程序的R实现”。
请参阅:http://www.r-project.org/conferences/useR-2009/slides/SanchezEspigares+Ocana.pdf
在示例的幻灯片22上,他们使用函数bootstrap()
(在幻灯片中:sleep.boot=bootstrap(model,B=1000)
)。
他们引用的唯一包是lme4
,但该包不包含bootstrap()
函数,我得到:
Error: could not find function "bootstrap"
有人知道他们在这里使用的套餐吗?
答案 0 :(得分:0)
我99%确定lme4
从未拥有bootstrap()
函数,并且它实际上来自演示者编写的代码(我之前已经在线查找过该代码但从未找到了)。内置的bootMer
函数尝试执行此演示文稿所描述的大部分内容,但不会执行所有操作(例如,它不包含Wild bootstrap的实现。)
如果您想在不指定函数bootMer
的情况下使用FUN
,可以查看confint(.,method="boot")
...它看起来像作者使用的摘要函数可能是
function(x) c(fixef(x),getME(x,"theta"),sigmaREML=sigma(x))