R语言代码是:model1 <- auto.arima(quet.ts, xreg=cbind(step,ramp), max.d=1, max.D=1, stepwise=FALSE, trace=TRUE)
我初步写的stata代码是:autoarima y step ramp,noseas 但是发现时间序列y需要差分,根据ac,pac图定阶为1阶,所以修正为 arima y step ramp,arima(1,1,0) 但是这样就相当于把y 和step ramp三个变量都差分了,而不是R语言代码中的xreg表达的意思,那stata代码该怎么写呢球球大神们,要毕不了业了