第4行出现错误，说“ Quarter [-1]”的可变长度不同。而且不太确定如何解决它。如果有人可以帮助，那就太好了。谢谢
co2.ts=ts(co2,start=c(2000,1),end=c(2016,4),frequency = 4)
Time = 1:68
Quarter=factor(c((1:4),rep(1:4,6)))
co2.fit1=lm(co2.ts[-1]~Time[-1]+Quarter[-1]+co2.ts[-68])
plot.ts(residuals(co2.fit1),main="Residual Series")
summary(co2.fit1)