我正在做一个作业，在找到一个最多具有3个交互项的模型之后，我需要进行boxcox转换。数据集包含19个独立变量E1-E4和G1-G15。
我的代码：
model <- lm( Y ~(E1+E2+E3+E4+G1+G2+G3+G4+G5+G6+G7+G8+G9+G10+G11+G12+G13+G14+G15)^3)
bc <- boxcox(model)
我希望boxcox（）返回一个图，并且我希望模型的摘要能够为我提供标准误差，t值和p值。
该模型的摘要输出的前几行如下：
“残差：所有1000个残差均为0：没有残差自由度！
系数：（由于奇异，未定义160）”
该表然后给出标准误差，t值和p值的所有“ NA”。
boxcox()返回
“ plot.window（...）中的错误：需要有限的'ylim'值”