如何在R中提取boxcox的95％CI

Question

boxcox(Volume ~ log(Height) + log(Girth), data = trees,
   lambda = seq(-0.25, 0.25, length = 10))

我的问题是，如果不绘制结果，如何提取λ的95％置信区间（CI）？ 我想检查1是否在CI内。如果是这样，我将不会转换数据，否则，我将使用具有最高可能性的lambda数据进行转换。

上面的代码只是R中的boxcox的一个例子。

谢谢！

Answer 1

我想我明白了。

bc = boxcox(Volume ~ log(Height) + log(Girth), data = trees,lambda = seq(-0.25, 0.25, length = 10))

bc$x[bc$y > max(bc$y) - 1/2 * qchisq(.95,1)]

会起作用。

参考：线性模型与R by Julian Faraway第8.1节。