Box-Cox変換

なかなか正規化しない分布を無理やり正規化する方法。
λを用いてこんな変換をする

T(x)=\frac{x^{\lambda}-1}{\lambda}

MASSライブラリをロードし、

library(MASS)
box_par <- boxcox(y~x,lambda = seq(-0.25, 1, length = 100))

などとするとλの尤度が計算される。

p <- box_par$x[which.max(box_par$y)]

で、最尤推定量を取り出す。