Model regresji wielorakiej: y = a + Ä…1 1i 2i x + Ä…2 x + & Ä…n ni n x + µ i 0 df liczba stopni swobody (degree of freedom): liczba obserwacji minus liczba zmiennych objaÅ›niajÄ…cych minus jeden, df = n (k+1) ZaÅ‚ożenia klasycznej MNK (metody najmniejszych kwadratów): dot. zmiennych i specyfikacji modelu: " model jest liniowy w stosunku do parametrów " zmienne egzogeniczne (objaÅ›niajÄ…ce) nie sÄ… zmiennymi losowymi " zmienne egzogeniczne nie sÄ… współliniowe (żadna z nich nie jest kombinacjÄ… pozostaÅ‚ych) " zmienne egzogeniczne nie sÄ… skorelowane ze skÅ‚adnikiem losowym (skoro nie sÄ… losowe to nie powinny być skorelowane ze skÅ‚adnikiem losowym) " liczba obserwacji jest wiÄ™ksza od liczby parametrów strukturalnych dot. skÅ‚adnika losowego: " wartość oczekiwana skÅ‚adnika losowego jest równa zero Å›redni wpÅ‚yw nie uwzglÄ™dnionych w modelu zmiennych, dla każdej zmiennej X jest równy zero " wariancja skÅ‚adnika losowego jest skoÅ„czona i staÅ‚a zaÅ‚ożenie czÄ™sto nie speÅ‚nione dla danych przekrojowych " brak autokorelacji skÅ‚adnika losowego należy sprawdzać dla szeregów czasowych, w danych przekrojowych to zaÅ‚ożenie jest zawsze speÅ‚nione WÅ‚asnoÅ›ci estymatorów MNK. Na mocy twierdzenia Gaussa Markowa estymator KMNK jest estymatorem: " liniowym " zgodnym (zbieżnym stochastycznie do Ä…) " nieobciążonym: E(a) = Ä… " najefektywniejszym (o najmniejszej wariancji) w kl. linio. i nieobciążonych estymat. Weryfikacja statystyczna: " interpretacja współczynnika determinacji i/lub kryteriów informacyjnych " test stopnia współliniowoÅ›ci zmiennych objaÅ›niajÄ…cych " testy istotnoÅ›ci: (a) podzbioru zmiennych objaÅ›niajÄ…cych (b) poszczególnych zmiennych objaÅ›niajÄ…cych " testy wÅ‚asnoÅ›ci skÅ‚adnika losowego modelu: (a) autokorelacji (b) heteroskedastycznoÅ›ci (c) normalnoÅ›ci rozkÅ‚adu Weryfikacja hipotez. JeÅ›li wynik próby należy do obszaru krytycznego hipoteza zerowa H0 jest odrzucana. BÅ‚Ä…d I rodzaju: odrzucenie poprawnej hipotezy zerowej BÅ‚Ä…d II rodzaju: przyjÄ™cie hipotezy H0, która jest faÅ‚szywa Test istotnoÅ›ci zmiennej: H0 badana zmienna nieistotna statystycznie; H1 zmienna istotna statystycznie; empiryczny poziom istotnoÅ›ci wynikajÄ…cy z przyjÄ™tej konwencji p- value = 0,05 lub 0,01 WzglÄ™dne bÅ‚Ä™dy szacunku: bÅ‚Ä…d standardowy współczynnika / wartość współczynnika (dopuszczalne sÄ… bÅ‚Ä™dy wzglÄ™dne < 50%) Kryterium informacyjne Akaike'a (AIC) Schwarza (SIC) miara 'odlegÅ‚oÅ›ci' danego modelu od modelu idealnego (prawdziwego), im niższe tym lepiej Klasyczny współczynnik determinacji (R2) miara dopasowania obserwacji należących do próby, interpretowalny gdy: (a) zależnoÅ›ci miÄ™dzy objaÅ›nianÄ… a objaÅ›niajÄ…cymi jest liniowa (b) parametry oszacowane MNK z wyrazem wolnym Skorygowany R2 kara za dodawania zbÄ™dnych zmiennych, w maÅ‚ych modelach (df<30) bardziej wiarygodny niż wersja klasyczna Typowe bÅ‚Ä™dy specyfikacji: " bÅ‚Ä™dna postać funkcyjna (np. liniowa zamiast nieliniowej) " pominiÄ™cie kluczowych zmiennych objaÅ›niajÄ…cych (=> obciążone estymatory MNK) " uwzglÄ™dnienie zbÄ™dnych zmiennych (=> nieefektywne estymatory) " bÅ‚Ä™dnie dobrana struktura dynamiczna " bÅ‚Ä™dy pomiaru zmiennych Z dwojga zÅ‚ego lepiej mieć estymatory nieefektywne niż obciążone. Test RESET (Regression Specification Error Test). Ogólny test wykrywajÄ…cy wiele bÅ‚Ä™dów specyfikacji modelu m.in.: (a) pominiÄ™te zmienne objaÅ›niajÄ…ce (b) korelacjÄ™ miÄ™dzy zmiennÄ… objaÅ›niajÄ…cÄ… a skÅ‚adnikiem losowym (c) bÅ‚Ä™dnÄ… postać funkcyjnÄ… 1. utworzenie modelu pomocniczego: model podstawowy + kolejne potÄ™gi zmiennej 2 3 objaÅ›nianej (przeważnie y oraz y . Wyznaczenie R2 modelu pomocniczego. i i 2. Hipoteza zerowa: współczynniki przy dodatkowych zmiennych sÄ… równe zero, model liniowy, hipoteza alternatywa: min. 1 wsp. różny od zera model nieliniowy Test Davidsona MacKinnona. Sprawdzenie kompletnoÅ›ci modelu, czy konkurencyjny model nie tÅ‚umaczy lepiej zjawiska. Dwa blizniacze modele o odmiennych zbiorach zmiennych. Y z pierwszego modelu dodawane jest jako zmien. objaÅ›niajÄ…ca do 2 modelu (i vice versa), jeżeli w nowym modelu Y nieistotne model pierwotny kompletny Test pominiÄ™tych zmiennych. Sprawdzenie czy któraÅ› ze zmiennych objaÅ›niajÄ…cych nie powinna zostać usuniÄ™ta, weryfikacja statystycznej istotnoÅ›ci zbioru wÄ…tpliwych zmiennych, wskazówki, że kluczowa zmienna objaÅ›niajÄ…ca zostaÅ‚a pominiÄ™ta: 'podejrzane' znaki parametrów, sÅ‚abe dopasowanie modelu do danych empirycznych Metodologia od ogólnego do szczególnego. Zaczynamy od modelu z nadmiarem parametrów, w kolejnych etapach model upraszczany na podstawie testów istotnoÅ›ci. Zmienne jakoÅ›ciowe. Zastosowanie: (a) cechy jednostek (b) wystÄ…pienie wydarzeÅ„ nietypowych (c) sezonowość Postać: (a) binarna (b) uporzÄ…dkowana klasyfikacja (c) nieuporzÄ…dkowana klasyfikacja. Sytuacja komplikuje siÄ™ przy autokorelacji i heteroskedastycznoÅ›ci skÅ‚adnika losowego. Gdy wiÄ™cej zmien. jakoÅ›ciowych należy rozważyć uwzglÄ™dnienie ich Å‚Ä…cznego efektu na zmiennÄ… endogennÄ…. LMP liniowy model prawdopodobieÅ„stwa: zastosowanie MNK do modelu regresji liniowej. Zalety: Å‚atwość estymacji, bezpoÅ›rednia interpretacja. Wady: wartoÅ›ci teoretyczne mogÄ… wykraczać poza przedziaÅ‚ (0,1), bÅ‚Ä™dy prognoz mogÄ… być bardzo duże Zmienne ukryte. Yi* > 0 => y = 1 oraz zero w przeciwnym przypadku. SkÅ‚adnik losowy ma rozkÅ‚ad normalny model probitowy; skÅ‚adnik losowy ma rozkÅ‚ad logistyczny m. logitowy WÅ‚asnoÅ›ci logit/probit: jeżeli x należy (od minus nieskoÅ„czoność do plus nieskoÅ„czoność) to 0=< y =< 1. Dystrybuanty rozkÅ‚adu normalnego i logistycznego sÄ… bardzo do siebie podobne. Częściej stosowany model logitowy. Skomplikowana interpretacja. BezpoÅ›rednia interpretacja znaku: znak + - zwiÄ™kszenie odpowiadajÄ…cej mu zmiennej zwiÄ™ksza prawdopodobieÅ„stwo zjawiska opisanego przez zmiennÄ… objaÅ›nianÄ…. Wyniki dla modelu logit i probit sÄ… porównywalne po pomnożeniu parametrów logitu przez pierwiastek z trzech / pi. Aby porównać oceny parametrów LPM z ocenami parametrów probitu należy je pomnożyć przez 2,5 a od wyrazu wolnego odjąć 1,25. Ocena jakoÅ›ci modeli logit/probit: (a) współczynnik pseudo-R^2, (b) wartość statystyki testu ilorazu wiarygodnoÅ›ci (c) tablica trafnoÅ›ci o wymiarach 2x2 z liczbÄ… przypadków trafionych / nietrafionych Współliniowość liniowa zależność (korelacja) pomiÄ™dzy zmiennymi objaÅ›niajÄ…cymi. Cecha danych nie modelu. Skutki: (a) utrudniona wiarygodna interpretacja ocen poszczególnych parametrów (b) zawyżone oceny Å›rednich bÅ‚Ä™dów szacunku parametrów CIW (z ang. VIF variance inflation factor) porównanie sytuacji faktycznej z idealnÄ…, w 2 której wszystkie zmienne objaÅ›niajÄ…ce sÄ… ze sobÄ… nieskorelowane: VIF (²^ ) = 1 / (1 R ) j i 1. model pomocniczy: dana zmienna objaÅ›niana przez pozostaÅ‚e. 2. użyć otrzymany R^2. Brak współliniowoÅ›ci CIW = 1, CIW > 5 => R^2 > 80 silne zakłócenia współliniowoÅ›ciÄ… Radzenie sobie ze współliniowoÅ›ciÄ…: (a) regresja grzbietowa dodanie staÅ‚ej do wariancji zmiennych (b) regresja wzglÄ™dem głównych skÅ‚adowych przeksztaÅ‚cenie zmiennych objaÅ›niajÄ…cych na zbiór zmiennych nieskorelowanych (c) usuwanie zmiennych powodujÄ…cych wystÄ™powanie współliniowoÅ›ci (d) zastÄ™powanie zmiennymi zastÄ™pczymi o podobnej info. merytorycznej (ale sÅ‚abej skorelowanymi) (e) zwiÄ™kszenie próby Heteroskedastyczność duża rozbieżność pomiÄ™dzy najmniejszymi i najwiÄ™kszymi obserwowanymi wartoÅ›ciami, najczęściej dot. danych przekrojowych, duży rozrzut reszt w oszacowanym modelu, wariancja skÅ‚adnika losowego nie staÅ‚a dla wszystkich obserwacji. Skutki: (a) estymatory nieefektywne (ale i nieobciążone) (b) oceny wariancji estymatorów obciążone (c) odchylenia standardowe nie sÄ… wiarygodne. Test hetero może wyÅ‚apać bÅ‚Ä™dnÄ… postać funkcyjnÄ… lub pominiÄ™te zmienne objaÅ›niajÄ…ce Test White'a. Hipoteza zerowa: homoskedastyczność; alternatywna: heteroskedastyczność. Procedura: 1. Model MNK, wyznaczenie reszt 2. model pomocniczy: zmienna objaÅ›niana kwadrat MNK-reszt (ej2) objaÅ›niajÄ…ce: objaÅ›niajÄ…ce w modelu podstawowym + ich kwadraty + ich iloczyny. Na podstawie Chi^2 decyzja czy odrzucić H0 jeżeli poziom istotnoÅ›ci <0,05 odrzucamy. Eliminacja heteroskedastycznoÅ›ci: ważona MNK, użycie deflatorów, transformacja do postaci logarytmicznej, estymatory HAC uwzglÄ™dniajÄ…ce autokorelacjÄ™ i heteroskedstyczn. Normalność rozkÅ‚adu skÅ‚adnika losowego: (a) nie jest niezbÄ™dna do uzyskania estymatorów o pożądanych wÅ‚asnoÅ›ciach (b) rozkÅ‚ad normalny skÅ‚adnika losowego ma duże znaczenie przy weryfikacji modelu ekonometrycznego. Test Jarque Bery, H-zero: rozkÅ‚. normalny Autokorelacja: korelacja miÄ™dzy skÅ‚adnikami losowymi modelu; autokorelacja miÄ™dzy et i et-k jest rzÄ™du k i oznaczana jest przez pk ; stanowi zÅ‚amanie zaÅ‚ożenia o sferycznoÅ›ci skÅ‚adnika losowego podstawy twierdzenia Gaussa Markowa. Przyczyny autokorelacji: (a) natura procesów gospodarczych: np. dÅ‚ugofalowe skutki decyzji gospodarczych, inercja procesów gosp. (b) bÅ‚Ä™dy specyfikacji modelu: (b1) niepoprawna postać analityczna (b2) niepeÅ‚ny zestaw zmiennych objaÅ›niajÄ…cych (b3) niewÅ‚aÅ›ciwa struktura dynamiczna. Skutki autokorelacji: (a) w przypadku braku opóznieÅ„: nieefektywność estymatorów (b) przy opóznieniach zmiennej objaÅ›nianej: estymatory MNK nie sÄ… zgodne Test DW (Durbina Watsona). Wady: (a) tylko modele z wyrazem wolnym, bez opóznieÅ„ zmiennej objaÅ›nianej, z normalnym rozkÅ‚adem skÅ‚adnika losowego (b) nie pozwala wykryć autokorelacji rzÄ™du wyższego niż 1 (c) nie zawsze prowadzi do jednoznacznego wyniku. H-zero: p = 0 (brak autokorelacji) H-1: p<0 auto ujemna v p>0 auto dodatnia (wartość DW decyzja) 0 p = 1; (0, DL) autokorelacja dodatnia; (DL; DU) test nie dziaÅ‚a; (DU, 2) brak autokorelacji; 2 p = 0; (2, 4-DU) brak autokorelacji; (DU, DL) nie dziaÅ‚a; (4- DL, 4) autokorelacja ujemna; 4 p = -1. Test LM (Lagrange multiplier), mnożnika Lagrange'a. Brak ograniczeÅ„ z testu DW. Hipoteza zerowa: brak autokorelacji, alternatywna wystÄ™powanie. Gdy Chi^2 < poziom istotnoÅ›ci odrzucamy hipotezÄ™ zerowÄ… => model z autokorelacjÄ…. Przeprowadzanie: 1. wyjÅ›ciowe równanie przy pomocy MNK + wyznaczenie jego reszt 2. model pomocniczy: model z pkt. 1 + p dodatkowych zmiennych bÄ™dÄ…cych opóznieniami oszacowanych reszt. Niemożliwość zastosowania w maÅ‚ych próbach (n < 30). Eliminacja autokorelacji: (a) korekta metody estymacji parametrów (b) zmiana postaci analitycznej (c) HAC: estymatory bÅ‚Ä™dy standardowego uwzgl. autokorelacjÄ™ i heterosk. Szeregi czasowe. Wyróżnia siÄ™ nastÄ™pujÄ…ce skÅ‚adowe: przeciÄ™tny poziom (M), dÅ‚ugookresowy trend (T), wahania sezonowe (S), wahania cykliczne (C), zm. nieregularne (I) yt = M + T + S + C + I Przyczyny opóznieÅ„: (a) psychologiczne: oczekiwania, przyzwyczajenia, inercja instytucji publicznych (b) technologiczne: koszty dostosowaÅ„ (c) instyt. - prawne: zobow., umowy Przyczynowość w sensie Grangera. X jest przyczynÄ… dla Y jeżeli X pomaga prognozować Y. Nie oznacza to, że Y jest skutkiem lub efektem X, a jedynie: X poprzedza Y. Stacjonarność Å›rednia i wariancja staÅ‚e i niezależne od upÅ‚ywu czasu Niestacjonoarność w miarÄ™ systematyczna zmiana wraz z upÅ‚ywem czasu np. trend rosnÄ…cy. Skutki regresja pozorna: (a) zawyżanie współczynnika determinacji (b) zawyżanie wartoÅ›ci statystyk k-Studenta i obciążenie innych statystyk wyznaczanych na podstawie odchyleÅ„ standardowych. WiÄ™kszość szeregów czasowych jest niestacjonarna. Test DF (Dickeya Fullera) z wyrazem wolnym lub bez / z trendem lub bez. "yt = Ä… + ´yt-1 + µt H-zero: ´ = 0 szereg niestacjonarny (z pierwiastkiem jednostkowym); H-jeden: ´ < 0 szereg stacjonarny Statystyka testu pierwiastka jednostkowego DF = ´ / bÅ‚Ä…d std. ´. RozstrzygniÄ™cie testu przy pomocy wartoÅ›ci krytycznych najczęściej w wersji ADF, uwzglÄ™dniajÄ…cej autokorelacjÄ™ skÅ‚adnika losowego w modelu testowym. StopieÅ„ integracji: ilość iteracji po której z szeregu niestacjonarnego otrzymamy szereg stacjonarny. Integracja sezonowa: obliczanie różnic dotyczy obserwacji oddalonych o s okresów, gdzie s jest dÅ‚ugoÅ›ciÄ… cyklu (s = 4 kwartaÅ‚y, s = 12 miesiÄ…ce) wówczas: "y = Ä… + ´y + µ Oba t t-s t rodzaje integracji mogÄ… wystÄ™pować jednoczeÅ›nie. Kointegracja dwa lub wiÄ™cej szeregi sÄ… niestacjonarne i zintegrowane w tym samym stopniu, ale ich liniowa kombinacja jest stacjonarna. WystÄ…pienie relacji kointegrujÄ…cej sugeruje obecność dÅ‚ugookresowej relacji równowagi (braku tendencji do zmiany) np. pÅ‚ace i ceny, konsumpcja i oszczÄ™dnoÅ›ci, deficyt budżetowy i inflacja, ceny akcji i dywiden. Test CRDW / kointegracji. Test niestacjonarnoÅ›ci reszt potencjalnej regresji kointegrujÄ…cej: Yt = ²0 + ²1 Xt + µt H-zero: pierwiastek jednostkowy w resztach niestacjonraność reszt; odrzucane gdy CRDW wiÄ™ksza niż wartość krytyczna (ok. 0,5) W przypadku podstaw do odrzucenia hipotezy zerowej: istnieje stacjonarna kombinacja niestacjonarnych zmiennych odpowiednia strategia modelowania: model korekty bÅ‚Ä™dem. Metoda Engle'a Grangera. Prostota zastosowania. Wady: (a) niska moc testów pierwiastka jednostkowego w maÅ‚ych próbach (b) asymetryczne traktowanie zmiennych (b) brak możliwoÅ›ci weryfikacji hipotez na temat relacji kointegrujÄ…cej. NOTATKI NIE UWZGLDNIAJ: testu istotnoÅ›ci zmiennej (nie chodzi tu o interpretacjÄ™ p- value t-Studenta), dokÅ‚adnego omówienia wykorzystania HAC estymatorów nieobciążonych (generalnie chodzi o to, że sÄ… stosowane gdy autokorelacja lub heteroskedastycznoÅ›ci wynika z samej natury procesów gospodarczych: tzn brak współliniowoÅ›ci, RESET na plus itp.), testu Beuschera (coÅ› tam) na autokorelacjÄ™; oraz pewnie kilku innych mniej ważnych aspektów H-zero testu Walda mówi, że żadna zmienna nie jest istotna: zaprzeczenie przynajmniej jedna jest (a nie, że wszystkie sÄ…) negacja zaÅ‚ożenia z kwantyfikatorem ogólnym jest ekwiwalentna do kwantyfikatora egzystenacalnego + negacji zaÅ‚ożenia tzn. wystarczy znalezć jeden kontrprzykÅ‚ad by caÅ‚e zaÅ‚ożenie byÅ‚o nieprawdziwe. POWODZENIA na egz.!!! pani Tomczyk jest naprawdÄ™ w porzÄ…dku osobÄ…!!! Nazwa testu Co testuje Hipoteza zerowa / Sposób Uwagi alternatywna liczenia t-Studenta Istotność Nieistotna ?? BÅ‚Ä…d I rodzaju: odrzucenie hipotezy danej (współczynnik = prawdziwej; zmiennej 0) / istotna (różny bÅ‚Ä…d II rodzaju: nieodrzucenie od zera) faÅ‚szywej; <0,05 zmienna istotna Walda Istotność Å»adna zmienna (R^2)/k : (1- F(k, n (k+1)) n ilość obserwacji, statystyka F modelu nie jest istotna / R^2)(n-k-1) k ilość zmiennych objaÅ›niajÄ…cych chociaż 1 jest <0,05 któraÅ› ze zmiennych istotna Chi^2 RozkÅ‚ad Normalny ?? <0,05 brak rozkÅ‚adu normalnego reszt rozkÅ‚ad reszt / nie normalny rozkÅ‚ad reszt RESET Adekwatn adekwatność / Dodanie Jeżeli współczynnik przy y różny od t ość postaci brak kwadratu zera zÅ‚a postać funkcyjna; testuje funkcyjnej, adekwatnoÅ›ci oraz również pominiÄ™te zmienne objaÅ›., szeÅ›cianu y t korelacjÄ™ miÄ™dzy zm. o-cÄ… a skÅ‚adnikiem losowym <0,05 zle DW Durbina Autokorela P = 0 brak ?? Interpretowany na podstawie tablic. Watsona cjÄ™ autokorelacji / Wady: tylko modele z wyr. wolnym, p < 0 ujemna, bez opóznieÅ„, z norm. rozkÅ‚adem skÅ‚. p > 0 dodatnia losowego, nie wykrywanie a-k rzÄ™du > 1, nie zawsze rozstrzyga DF (Dickeya niestacjona niestatcjonarny z DF = ´ / bÅ‚Ä…d. stand. ´ "yt = Ä… + Fullera) rność pierw jedn ´ = 0 / ´yt-1 + µt rozstrzygajÄ… wartoÅ›ci krytyczne stacjonarny ´ < 0 najczęściej w wersji ADF LM Autokorela Brak Niemożność stosowania w przypadku maÅ‚ych prób. cjÄ™ autokorelacji / 1. model MNK, wyznaczenie reszt 2. dodanie p wystÄ™powanie dodatkowych zmiennych bÄ™dÄ…cych opóznieniami autokorelacji oszacowanych reszt. Rozstrzyga Chi^2 (<0,05 - autokorelacja) VIF (en) Współlinio CIW = 1 brak VIF (²^ ) = 1. dla każdej objaÅ›niajÄ…cej model, j 2 CIW (pl) wość CIW > 5 1 / (1 R ) w którym jest objaÅ›nianÄ…. 2. z i poważnie dla każdej otrzymanego współczynnika zakłócajÄ…ca zmiennej obliczyć CIW interpretacjÄ™ White'a Heteroske Homoskedastycz 1. Model MNK, wyznaczenie reszt 2. zmienna dastycznoÅ› ność / objaÅ›niana kwadrat reszt z MNK wobec zmiennych ć heteroskedastycz objaÅ›niajÄ…cych, ich kwadratów oraz iloczynów. ność Decyzja na podstawie Chi^2 (<0,05 heteroskedastyczność) CRDW / Kointegrac Reszty Yt = ²0 + ²1 Gdy CRDW > od ok. 0,5 istnieje kointegracji ji niestacjonarne / stacjonarna kombinacja dwóch Xt + µt kombinacja zmiennych niestacjonarnych dwóch zmiennych stacjonarna