R boxcox 线性模型的 Box-Cox 变换

说明

计算并可选择绘制 Box-Cox 幂变换参数的剖面对数似然。

用法

boxcox(object, ...)

## Default S3 method:
boxcox(object, lambda = seq(-2, 2, 1/10), plotit = TRUE,
       interp, eps = 1/50, xlab = expression(lambda),
       ylab = "log-Likelihood", ...)

## S3 method for class 'formula'
boxcox(object, lambda = seq(-2, 2, 1/10), plotit = TRUE,
       interp, eps = 1/50, xlab = expression(lambda),
       ylab = "log-Likelihood", ...)

## S3 method for class 'lm'
boxcox(object, lambda = seq(-2, 2, 1/10), plotit = TRUE,
       interp, eps = 1/50, xlab = expression(lambda),
       ylab = "log-Likelihood", ...)

参数

object	公式或拟合模型对象。目前仅处理lm 和aov 对象。
lambda	lambda 值向量 - 默认 (-2, 2)，步长为 0.1。
plotit	逻辑控制是否应绘制结果。
interp	控制是否使用样条插值的逻辑。如果使用长度小于 100 的 lambda 进行绘图，则默认为 TRUE。
eps	lambda = 0 的容差；默认为 0.02。
xlab	默认为 "lambda" 。
ylab	默认为 "log-Likelihood" 。
...	模型拟合中使用的附加参数。

值

lambda 向量和计算的轮廓对数似然向量的列表，如果绘制结果，则不可见。

副作用

如果 plotit = TRUE 绘制对数似然与 lambda 的关系，并指示关于 lambda 最大观测值的 95% 置信区间。如果 interp = TRUE ，则使用样条插值来给出更平滑的绘图。

例子

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

boxcox(Days+1 ~ Eth*Sex*Age*Lrn, data = quine,
       lambda = seq(-0.05, 0.45, length.out = 20))

参考

Box, G. E. P. and Cox, D. R. (1964) An analysis of transformations (with discussion). Journal of the Royal Statistical Society B, 26, 211-252.

Venables, W. N. and Ripley, B. D. (2002) Modern Applied Statistics with S. Fourth edition. Springer.