fit=auto.arima(as.numeric(datanew[,j]),# 第j列数据转数值。coef=matrix(0,16,12)# 新建16行12列的0矩阵，存12个序列的AR系数。max.q=0,d=1,# 强制q=0，不做MA。xlab='',ylab='',zlab='')# 不显示轴标签。# 核心：拟合 ARIMA(p,1,0) 模型 = 一阶差分后的 AR(p)coef=co