r语言时间序列预测

1531

收藏 2016-10-07

> w[,1]
Time Series:
Start = 1
End = 288
Frequency = 1
  [1] 0.181 0.180 0.180 0.175 0.173 0.170 0.170 0.166 0.209 0.202
[11] 0.192 0.174 0.163 0.189 0.186 0.164 0.185 0.222 0.235 0.246
[21] 0.254 0.224 0.203 0.194 0.203 0.214 0.213 0.215 0.234 0.274
[31] 0.281 0.264 0.249 0.234 0.237 0.225 0.213 0.240 0.217 0.225
[41] 0.283 0.285 0.289 0.283 0.269 0.245 0.219 0.221 0.233 0.237
[51] 0.234 0.256 0.261 0.284 0.279 0.278 0.251 0.242 0.274 0.250
[61] 0.230 0.262 0.231 0.266 0.295 0.298 0.292 0.270 0.250 0.243
[71] 0.241 0.224 0.219 0.215 0.217 0.230 0.220 0.232 0.255 0.262
[81] 0.257 0.259 0.255 0.238 0.223 0.209 0.210 0.205 0.199 0.197
[91] 0.213 0.253 0.221 0.215 0.216 0.220 0.224 0.217 0.213 0.211
[101] 0.211 0.206 0.205 0.203 0.202 0.204 0.203 0.201 0.202 0.200
[111] 0.205 0.208 0.206 0.209 0.210 0.205 0.207 0.210 0.213 0.193
[121] 0.183 0.215 0.221 0.219 0.213 0.219 0.213 0.209 0.219 0.220
[131] 0.219 0.213 0.222 0.214 0.200 0.205 0.210 0.208 0.204 0.195
[141] 0.195 0.193 0.192 0.188 0.180 0.181 0.182 0.180 0.178 0.178
[151] 0.183 0.185 0.181 0.182 0.185 0.181 0.182 0.180 0.178 0.178
[161] 0.180 0.178 0.178 0.180 0.174 0.172 0.169 0.168 0.167 0.168
[171] 0.169 0.167 0.168 0.167 0.167 0.164 0.163 0.169 0.164 0.165
[181] 0.165 0.160 0.162 0.159 0.164 0.159 0.155 0.155 0.154 0.153
[191] 0.150 0.145 0.135 0.171 0.214 0.226 0.220 0.232 0.229 0.215
[201] 0.191 0.189 0.171 0.160 0.165 0.176 0.181 0.192 0.221 0.221
[211] 0.185 0.215 0.241 0.187 0.167 0.162 0.190 0.160 0.173 0.216
[221] 0.235 0.248 0.253 0.244 0.226 0.220 0.205 0.205 0.210 0.229
[231] 0.243 0.262 0.272 0.250 0.252 0.216 0.198 0.192 0.190 0.187
[241] 0.183 0.176 0.174 0.176 0.179 0.180 0.182 0.180 0.175 0.173
[251] 0.170 0.170 0.172 0.170 0.172 0.211 0.195 0.212 0.205 0.192
[261] 0.201 0.270 0.238 0.240 0.259 0.250 0.218 0.208 0.204 0.205
[271] 0.203 0.203 0.205 0.210 0.205 0.198 0.198 0.209 0.210 0.221
[281] 0.218 0.216 0.248 0.250 0.239 0.255 0.235 0.234
> plot(diff(w[,1],difference=2),type="l")
> abline(h=0)

> acf(diff(w[,1],difference=2))
> pacf(diff(w[,1],difference=2))
> eacf(diff(w[,1],difference=2))
AR/MA
  0 1 2 3 4 5 6 7 8 9 10 11 12 13
0 x x x o o o o o o o o  o  o  o
1 x x x o o o o o o o o  o  o  o
2 x o x o o o o o o o o  o  o  o
3 x o o o x o o o o o o  o  o  o
4 x o o o o o o o o o o  o  o  o
5 x o o o x o o o o o o  o  o  o
6 x o x x o o x o o o o  o  o  o
7 x x o x o o x x o o o  o  o  o

> m1=arima((w[,1]),order=c(0,2,3))

> plot(m1,n.ahead=50,type="b")