记录一个在做逐步回归方程时的报错和可能的原因,如有错误还请批评指正。
报错:Error in step(lm(LW ~ ., data = data), direction = "both", trace = 0) : 这个模型AIC的值是负无限的,因此'step'不能继续进行
data <- data1[,c(2,11:20)];data
model2 = step(lm(LW~., data = data),direction = 'both' ,trace = 0); model2
# Error in step(lm(LW ~ ., data = data), direction = "both", trace = 0) :
# 这个模型AIC的值是负无限的,因此'step'不能继续进行
可能的原因:由于用于回归的数据条数只有11条,但回归方程中除了 "LW",剩余“预测变量”的数量为10个,经过多组数据测试分析发现,预测变量的个数不能大于数据条数(或数据行数)- 2,这里对应的为9个预测变量个数。
减少预测变量个数即不会出现报错
为了保证回归方程的可靠性,还是尽量保证数据量吧,这里我是计算了11组数据的平均值进行计算,所以本来有将近20个预测变量但无法进行逐步回归,最后选择了使用原始数据(122条/行)进行计算得到的预测结果和预期还挺接近的。