貌似不能直接做~~网上找到个R的，boxcox(线性模型,lambda=seq(起点,终点,by=间距))
或者boxcox(线性模型，lambda=seq(起点,终点,len=点的数目))
boxcox函数是计算若干个离散lambda点的值，给出boxcox变换后的线性模型的残差为正态分布N(0,sigma^2)的置信率。
例子：
boxcox(Y~X,lambda=seq(0,1,len=20))