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.