Jak odstranit sezónu z řady? (sezónní očištění pomocí X-13, kód R s sebou)

Stručně

Jak odstranit sezónu z měsíční nebo čtvrtletní řady? X-13 v prohlížeči neběží, proto Statable sezónní očištění sám nepočítá. Vytvoří řadu jako soubor CSV a skript R s sebou, spouštěný v R s balíčkem seasonal.

Jak to funguje ve Statable

1. V analýze časových řad otevřete oddíl „Sezónní očištění metodou X-13 (v R na vašem počítači)“.
2. Zvolte metodu: SEATS (výchozí) nebo X-11.
3. „Uložit řadu jako CSV“ a „Uložit skript R“.
4. Oba soubory dejte do jedné složky a spusťte skript v R. Zapíše druhý CSV s očištěnou řadou.

Není totéž co klasická dekompozice

Dekompozice v analýze časové řady (trend, sezóna, zbytek) není sezónní očištění metodou X-13 a tak se nenazývá. Zůstává, jak je.

Kdy se nabízí

• Jen pro měsíční nebo čtvrtletní hodnoty.
• Řada nesmí mít mezery. Statable nic nedoplňuje.
• Je potřeba alespoň 3 celých let: 36 měsíčních hodnot nebo 12 čtvrtletních hodnot.
• Pokud některá z těchto podmínek neplatí, objeví se místo tlačítek věta s důvodem.

Všechny kroky podrobně

1. Zvolte metodu: SEATS (výchozí, standard balíčku seasonal) nebo X-11.
2. „Uložit řadu jako CSV“ a „Uložit skript R“.
3. Oba soubory dejte do stejné složky a spusťte skript v R. R k tomu potřebuje balíček seasonal (s x13binary).
4. Skript zapíše druhý soubor CSV s očištěnou řadou. Lze ho ve Statable znovu otevřít.

Soubor CSV obsahuje časový okamžik a hodnotu, s tečkou jako desetinným znakem, v UTF-8. Pomocí „Zobrazit skript R“ lze skript před uložením přečíst.

Varování

• Statable to nepočítá. Výsledek vzniká v R na vašem vlastním počítači a Statable ho nezkontroloval.
• Kalendář: seasonal testuje standardní efekty pracovních dnů a Velikonoc. Rakouské státní svátky nejsou zahrnuty; směrnice ESS doporučují národní kalendáře. Kdo je potřebuje, přidá vlastní regresory přes xreg.
• Volba modelu v X-13 je automatická. Je třeba si přečíst summary(m) a diagnostiku.

Typické příklady

• Odstranit roční vzorec z měsíční míry nezaměstnanosti dříve, než se vývoj interpretuje.
• Očistit čtvrtletní tržby za několik let a porovnat je s původní řadou.

Pro odborníky

library(seasonal)
d <- read.csv("reihe.csv")
y <- ts(d$wert, start = c(1949, 1), frequency = 12)   # 12 bei Monaten, 4 bei Quartalen
m <- seas(y)                                           # bei X-11: seas(y, x11 = "")
summary(m)
plot(m)
bereinigt <- final(m)

Pro kontrolu: seas(AirPassengers) dává model (0 1 1)(0 1 1); první očištěné hodnoty jsou 122,713, 124,766 a 125,073.