如何在R中提取boxcox的95%置信区间

如何在R中提取boxcox的95%置信区间,r,R,我的问题是,在不绘制结果的情况下,如何提取lambda的95%置信区间(CI)? 我想检查1是否在CI中。如果是这样,我将不转换数据,否则,我将使用具有最大可能性的lambda的数据进行转换 上面的代码只是R中boxcox的一个示例 谢谢大家! 我想我明白了 boxcox(Volume ~ log(Height) + log(Girth), data = trees, lambda = seq(-0.25, 0.25, length = 10)) 会有用的 参考:Julian Faraw

我的问题是,在不绘制结果的情况下,如何提取lambda的95%置信区间(CI)? 我想检查1是否在CI中。如果是这样,我将不转换数据,否则,我将使用具有最大可能性的lambda的数据进行转换

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

谢谢大家!

我想我明白了

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

参考:Julian Faraway第8.1节中的线性模型和R

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)]