Statystyka od zera · Lekcja 31 z 31

Szeregi czasowe: trendy, sezony i fałszywe korelacje

W jednym zdaniu: Dane zbierane w czasie mają trendy, sezony i pamięć, więc porównujesz podobne z podobnym (ten styczeń z zeszłym styczniem) i nigdy nie liczysz korelacji dwóch rosnących szeregów bez spojrzenia najpierw na ich zmiany.

Pomysł

Większość tego kursu traktowała dane jak worek osobnych obserwacji. Szereg czasowy jest inny: to ta sama rzecz mierzona raz za razem w czasie, na przykład miesięczna sprzedaż, dzienne temperatury czy cena akcji. Szeregi czasowe mają trzy nawyki, które łamią zwykłe reguły:

Budka z lodami cioci Grażyny

Ciocia Grażyna prowadzi budkę z lodami nad jeziorem. W styczniu trzeciego roku sprzedała 50 lodów, mniej niż 55 w grudniu. Wujek Janusz, jej samozwańczy doradca finansowy, nazywa to „9% zapaścią” i proponuje zamknięcie.

Grażyna zauważa, że jest styczeń.

Oto jej miesięczna sprzedaż w trzecim roku, z prostego modelu budki: wzór sezonowy plus stały wzrost o 1% miesięcznie.

Po­rów­na­nieZmiana
Styczeń wobec po­przed­nie­go grudnia9,1% w dół
Wrze­sień wobec sierp­nia26,9% w dół
Styczeń wobec stycz­nia rok wcze­śniej11,1% w górę
Lipiec wobec lipca rok wcze­śniej10,0% w górę
Cały rok wobec po­przed­nie­go10,2% w górę (1448 wobec 1314)

W porównaniach miesiąc do miesiąca dominuje sezon. Każda jesień wygląda na krach, a każda wiosna na boom. Rozwiązanie to porównywać podobne z podobnym: ten styczeń z zeszłym styczniem (rok do roku) albo cały rok z całym rokiem. Urzędy statystyczne publikują dane „wyrównane sezonowo” z tego samego powodu.

Sprzedane lody miesięcznie, trzeci rok

Sty50
Lut56
Mar76
Kwi102
Maj141
Cze194
Lip221
Sie216
Wrz158
Paź106
Lis67
Gru61
Model: wzór sezonowy plus wzrost o 1% miesięcznie. Rok był o 10,2% lepszy od poprzedniego, choć jesień wygląda na krach.

Pamięć: autokorelacja

Dzisiejsza cena akcji to prawie dokładnie wczorajsza plus mała zmiana. Policz korelację szeregu cen z nim samym przesuniętym o jeden dzień, a w naszej symulacji wyjdzie około 0,98: niezwykle silny związek zwany autokorelacją. Policz ją dla dziennych zmian, a wyjdzie około 0,00. Na tym polega cała różnica między poziomem, który pamięta wszystko, a jego zmianą, która może nie pamiętać nic.

Ma to znaczenie, bo większość testów zakłada niezależne obserwacje. Rok dziennych cen to nie 250 niezależnych informacji. To bliżej jednej błądzącej linii.

Pułapka fałszywej korelacji

Oto najdroższy błąd w tym temacie. Wygenerowaliśmy 2000 par zupełnie niezwiązanych błądzeń losowych po 250 kroków (pomyśl o dwóch wymyślonych cenach akcji, które nie mają ze sobą nic wspólnego) i policzyliśmy korelację każdej pary:

Co ko­re­lo­wa­li­śmyPary z ko­re­la­cją powyżej 0,5 (dowolny znak)
Poziomy („ceny”)41%
Dzienne zmiany0%

Dwa szeregi, które oba błądzą, często przez jakiś czas błądzą w tę samą stronę, a ich korelacja wygląda na silną przez czysty przypadek. Typowa para miała korelację 0,42. To korelacja pozorna (albo regresja pozorna) i dlatego „te dwa wykresy idą razem” samo w sobie niczego nie dowodzi. Policz korelację zmian, co usuwa wspólne błądzenie, a fałszywy związek znika.

Trzy nawyki dla danych w czasie

Spróbuj sam

Blog pokazuje dwie linie rosnące razem od 2000 do 2020 roku: liczbę subskrypcji streamingu i średnią krajową cenę domów. Korelacja wynosi 0,95. Autor sugeruje, że streaming podnosi ceny domów.

1. Dlaczego 0,95 nie robi tu wrażenia? 2. Co sprawdziłbyś zamiast tego?

Pokaż odpowiedź

1. Oba szeregi mają silne trendy wzrostowe, a prawie każde dwa rosnące szeregi są mocno skorelowane. W naszej symulacji nawet niezwiązane błądzenia losowe przekraczają 0,5 w 41% przypadków; dwa szeregi, które oba równo rosną, wypadną jeszcze lepiej.

2. Policz korelację zmian rok do roku: czy ceny domów rosły szybciej w latach, gdy subskrypcje rosły szybciej? Potem zapytaj o wspólne przyczyny, takie jak rosnące dochody i niskie stopy procentowe, które pchają w górę jedno i drugie (lekcja 7).

Słowa, których używa wykładowca

Wy­kła­dow­ca mówiTo znaczy
Szereg czasowyTa sama wiel­kość mie­rzo­na wie­lo­krot­nie w czasie
TrendDłu­go­ter­mi­no­wy dryf w górę albo w dół
Se­zo­no­wośćWzór po­wta­rza­ją­cy się w stałych od­stę­pach
Wy­rów­na­nie se­zo­no­weUsu­nię­cie ty­po­we­go wzoru se­zo­no­we­go, żeby dało się po­rów­ny­wać mie­sią­ce
Au­to­ko­re­la­cjaKo­re­la­cja szeregu z jego wła­sny­mi wcze­śniej­szy­mi war­to­ścia­mi
Błą­dze­nie losoweSzereg, w którym każda wartość to po­przed­nia plus losowy krok
Sta­cjo­nar­nySzereg, którego średnia i rozrzut nie dryfują w czasie; zmiany często są, poziomy często nie
Ko­re­la­cja (re­gre­sja) pozornaSilny, ale bez­sen­sow­ny związek dwóch sze­re­gów z trendem
Róż­ni­co­wa­niePraca na zmia­nach z okresu na okres zamiast na po­zio­mach

Trzy rzeczy do zapamiętania

Część 3 na razie

Lekcje 20-31 dodają narzędzia, które spotkasz na prawie każdym kursie i w każdym raporcie z badań: rozkład dwumianowy i rozkład Poissona do liczenia, ANOVA do porównywania więcej niż dwóch grup, wielkość efektu do oceny, czy istotny wynik ma znaczenie, regresję do średniej do wyłapywania zmian, które nie potrzebują przyczyny, i różnicę między obserwacją a eksperymentem do rozstrzygania, co jest przyczyną czego. Potem praktyczny zestaw: histogramy do oglądania kształtu danych, regresję wieloraką do ustalania innych czynników, bootstrap do przedziałów bez wzorów, ilorazy szans do czytania nagłówków o zdrowiu, testy rangowe do danych, które nie chcą się zachowywać, i szeregi czasowe do wszystkiego, co mierzy się w czasie.

Część kursu Statystyka od zera, darmowego kursu dla początkujących od WrongStat. Nie zastąpi wykładowcy, ale mamy nadzieję, że będzie dobrym tłumaczem.