使用预引导数据的BCa置信区间

时间:2018-01-11 10:27:09

标签: r confidence-interval statistics-bootstrap

我想使用 boot 包中的 boot.ci 函数来计算BCa置信区间。我已经进行了bootstrap重复(使用程序包提供的 boot 函数计算)。 要计算BCa置信区间,将引导复制与观察到的统计信息一起传递给 boot.ci 函数(参数 t t0 < / em>的)。但是,从技术上讲, boot.ci 需要一个类&#34; boot &#34;成为函数调用中的第一个参数。似乎可以将一些任意的 boot 对象传递给该函数,只要重复的数量( R )是合适的。但是,计算实际上不需要的伪引导对象似乎效率有点低。另一种可能性是从头开始构造一个伪引导对象,如此处所讨论的(How can I use pre bootstrapped data to obtain a BCa confidence interval?)和此处(Bootstrapping confidence intervals in R: BCa method and prescribed resamples)。这避免了至少浪费运行时间,但仍然不是很优雅。我想知道是否有人知道更智能的解决方案(例如 boot.ci 的修改代码)。非常感谢!

0 个答案:

没有答案