r - Box-Cox 参数估计是否达到 lambda=0？

Question

Box-Cox 参数估计是否达到 lambda=0？这样推荐的转换是对数函数？还是这种可能性只是使函数连续的一种数学精确性，但很少在实践中使用？

我问是因为我一直在尝试在 R 中测试 Box Cox 变换，看看它是否在适当的时候估计 lambda==0:

> require(car)
> x <- exp(qnorm(runif(1000))) # 1000 numbers from a normal distribution then exp()
> p = powerTransform(x)
> p
Estimated transformation parameters
     x
-0.1098415
# clearly this does not equal zero

Answer 1

由于数据是随机的，因此您不会得到恰好为零。

> p <- powerTransform(rlnorm(100))
> coef(p)
rlnorm(100) 
-0.05007203 
> coef(p, round=TRUE) # Recommended transform
rlnorm(100) 
          0 

如果去除随机性，该值会更接近于零。

> p <- powerTransform(qlnorm(ppoints(100)))
> p
Estimated transformation parameters 
qlnorm(ppoints(100)) 
       -2.635191e-12