我有一个包含 13 个自变量（其中 2 个是分类变量）和 678 个观察值的模型。所有变量都很重要。我打算检查系数的稳定性。我的方法是从具有一个变量的模型开始并找到该系数。然后我将添加另一个变量；但是，我也保留了以前的变量。我将继续，直到所有变量都在最终模型中。我正在使用以下 R 代码：

add.APV<-lm(HPV~log(APV),data=Mydata)
add1(add.APV,scope =~ AWD + Variety.bodyandchassis + Flexibility + Model.types + Year + log(CAC)+ Outsourcing + Platformstrategy + Hourly.total + Newownership + Vehicle.Launch + Segment.MD., test="F", k=2, trace=TRUE)

我收到以下错误：

Error in factor.scope(attr(terms1, "factors"), list(add = attr(terms2,  : 
  upper scope has term ‘NA’ not included in model

我想知道问题是什么。对于任何反馈，我们都表示感谢。