Then run imapa:
testModel <- imapa(ts(x[1:obsInsample],frequency=12),h=3)
The returned values in the frc.out are:
[,1] [,2] [,3]
AL1 0.7472482 0.7472482 0.7472482
AL2 NA NA NA
AL3 NA NA NA
AL4 NA NA NA
AL5 NA NA NA
AL6 NA NA NA
AL7 NA NA NA
AL8 NA NA NA
AL9 NA NA NA
AL10 NA NA NA
AL11 NA NA NA
AL12 NA NA NA
AL13 NA NA NA
AL14 NA NA NA
AL15 NA NA NA
AL16 NA NA NA
AL17 NA NA NA
AL18 NA NA NA
AL19 NA NA NA
AL20 NA NA NA
AL21 NA NA NA
AL22 NA NA NA
AL23 NA NA NA
AL24 NA NA NA
AL25 NA NA NA
AL26 NA NA NA
AL27 NA NA NA
AL28 NA NA NA
This is probably due to the huge number of zeroes in the series. But maybe make function more robust in such cases?
Use the following series:
x <- structure(c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 2, 18), .Tsp = c(1, 5.66666666666667, 12), class = "ts")
Then run imapa:
testModel <- imapa(ts(x[1:obsInsample],frequency=12),h=3)
The returned values in the
frc.out
are: