Closed EconProf closed 3 years ago
That's odd, I had thought that specifying X_13ARIMA_SEATS(y~seats())
would give the same results as X_13ARIMA_SEATS(y)
, as seats()
is the default. Investigating.
Looks like specifying seats = ""
overwrites the seats.noadmiss = "yes"
default option.
library(seasonal)
library(tsibble)
#>
#> Attaching package: 'tsibble'
#> The following objects are masked from 'package:base':
#>
#> intersect, setdiff, union
library(magrittr)
as.ts(tsibbledata::aus_production)[,"Beer"] %>%
seas(seats = "") %>%
final()
#> Error: X-13 run failed
#>
#> Errors:
#> - Available options for noadmiss are yes or no.
#>
#> Notes:
#> - Correct input errors in the order they are detected since the
#> first one or two may be responsible for the others (especially
#> if there are errors in the SERIES or COMPOSITE spec).
Created on 2021-03-26 by the reprex package (v1.0.0)
Oh, it actually is partial argument matching seats=""
to seats.noadmiss=""
in the above reprex. :vomiting_face:
Fixed, thanks as always.
I think this is more of a bug in the {seasonal}
package, but it can be fixed in {feasts}
.
The issue was that the default seas.noadmiss= "yes"
argument was dropped when list = list(seas = "")
is provided to seasonal::seas()
. If we also pass seas.noadmiss
into the list()
, it works without issue.
I followed the example in Section 3.5 of fpp3 to specify seats as a special, but using aus_production data for beer instead of us_retail_employment. It did not remove the seasonality (QS is very high) and reports that it was not a decomposition. Leaving out the seas special it worked, but with a note I do not understand.
Created on 2021-03-25 by the reprex package (v1.0.0)