如何在R中运行步骤（顺序）回归？

Question

我想运行顺序回归

Step1: y ~ a + b + c

Step2: y ~ a + b + c + a:b + a:c + b:c

Step3: y ~ a + b + c + a:b + a:c + b:c + a:b:c

我试过

firstModel  = lm(y ~ a + b + c, data=data1)

secondModel = lm(y ~ a + b + c + a:b + a:c + b:c, data=data1) 

fullModel  =  lm(y ~ a + b + c + a:b + a:c + b:c + a:b:c, data=data1) 

scope1 = list(firstModel, secondModel, fullModel)

step1 = step(firstModel, scope = scope1, data=data1, direction="forward")

我得到的结果如下：

Call:
lm(formula = y ~ a + b + c + a:c, data = data1)

在这个结果中，我没有看到我指定的3个模型。这是为什么？ 有没有办法“强迫”R来计算3个模型的回归？