Displaying 1 result from an estimated 1 matches for "revtotal".
2002 Nov 18
1
Prediction from arima() object (library ts) (PR#2305)
...ata.
Example:
total <- c(750, 775, 971, 1099, 1344, 610, 910, 1056, 1589, 1006, 1469, 1598,
876, 1104, 1197, 1062, 1783, 1554, 1400, 1340, 1013, 1031, 1030, 860, 943, 1542,
1400, 1301, 1443, 1057, 1296, 710, 1038, 1447, 694, 889, 1459, 957, 1284, 1520,
1172, 1826, 1751, 1968, 933, 1409, 1889)
revtotal <- numeric(length(total))
for(i in 1:length(total))
revtotal[i] <- total[length(total)-i+1]
# first model
m1 <- arima(revtotal, c(1,1,1))
pm1 <- predict(m1, n.ahead=15)
plot(1:47, y=revtotal, type="l", xlim=c(0,65))
lines(48:62, y=pm1$pred, col="red")
# second mod...