X=read. CSV("F: public_ html input. CSV head=F) Y=read. CSV (F: \public_htmloutputcS v",head=F)
◼ X=read.csv("F:\\public_html\\input.csv" ,head=F) ◼ Y=read.csv("F:\\public_html\\output.cs v",head=F)
■X=C y=ct (r) par(mfrow=c(1, 2)) ts plote(×) ts plot(y)
◼ x=c(t(X)) ◼ y=c(t(Y)) ◼ par(mfrow=c(1,2)) ◼ ts.plot(x) ◼ ts.plot(y)
par(mfrow=c(1, 2)) a acf(x) pacf(x) a par(mfrow=c(1, 2)) a acf(y) pacf(y)
◼ par(mfrow=c(1,2)) ◼ acf(x) ◼ pacf(x) ◼ par(mfrow=c(1,2)) ◼ acf(y) ◼ pacf(y)
一元分析 拟合输入序列 x.=-0.1228+ 1.97607B+1.37499B2-0.34336B3 拟合输出序列 y=53.90176+ 1-3.10703B+1.34005B2-021274B4
一元分析 ◼ 拟合输入序列 ◼ 拟合输出序列 2 3 1 1.97607 1.37499 0.34336 0.1228 B B B a x t t − + − = − + 2 4 1 3.10703 1.34005 0.21274 53.90176 B B B a y t t − + − = +
xar=arima(x, order=c(, 0,0)) yar=arima(y, order=c( 4,0, 0)) a xr=resid(xar) a yr=resid(yar)
◼ xar=arima(x,order=c(3,0,0)) ◼ yar=arima(y,order=c(4,0,0)) ◼ xr=resid(xar) ◼ yr=resid(yar)