1 / 39

Dzisiaj na wykładzie

Dzisiaj na wykładzie. Regresja wieloraka – podstawy i założenia Przykładowe oblicznia w Statistica Weryfikacja założeń w Statistica Zadanie. Regresja wieloraka. Bada związki między zmienną objaśnianą (zależną) i zmiennymi objaśniającymi (niezależnymi).

hayes-serrano
Download Presentation

Dzisiaj na wykładzie

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. Dzisiaj na wykładzie Regresja wieloraka – podstawy i założenia Przykładowe oblicznia w Statistica Weryfikacja założeń w Statistica Zadanie

  2. Regresja wieloraka Bada związki między zmienną objaśnianą (zależną) i zmiennymi objaśniającymi (niezależnymi)

  3. Jak przewidzieć odsetek upadków zwierząt przy transporcie? Jakie zmienne związane są z upadkami zwierząt? Co decyduje o upadkach zwierząt? Które z nich są najważniejsze i w jakim stopniu są powiązane ze stratami zwierząt? Typowe zadanie

  4. regresja pierwszego rzędu y = a + b1x1 + b2x2 + b3x3 + ... + e drugiego rzędu y = a + b1x1 + b2x2x2 + ... y = a + b1x1 + b2x2 +b3x2x3 + Dlaczego te regresje są liniowe?

  5. y = a + b1x1 + b2x2 + b3x3 + ... + e a to wyraz wolny. Średnia? b1, b2, b3 – cząstkowe współczynniki regresji to niezależne wkłady każdej ze zmiennych objaśniających Dlaczego cząstkowe?

  6. Współczynniki cząstkowe obrazują zależności po uwzględnieniu pozostałych czynników. przykład: długość włosów i wzrost Jak je oszacować?

  7. Jak wyznaczyć prostą regresji metoda najmniejszych kwadratów metoda najmniejszych kwadratów ważonych metoda najmniejszych reszt bezwględnych

  8. Metoda najmniejszych kwadratów

  9. R-kwadrat czyli współczynnik determinacji • Mówi o jakości przewidywania • R2=30% znaczy 30% wariancji opisane przez zmienne opisujące i 70% wciąż nie opisane (błędy)

  10. R-kwadrat R-kwadrat = 1 - SSE / SST SSE - suma kwadratów reszt SST - suma kwadratów obserwacji

  11. R-kwadrat skorygowane R-kwadrat zależy od liczby zmiennych objaśniających! Im więcej zmiennych tym większy. R-kwadrat skorygowane NIE zależy od liczby zmiennych objaśniających Porównując dwa modelu o różnej liczbie zmiennych patrz na R-kwadrat skorygowane

  12. R czyli korelacja R to pierwiastek z R-kwadrat. Tylko dodatni ( 0 - 1 ) !!! wskazuje na stopień powiązania zmiennych

  13. Zależność między zmiennymi jest liniowa • Trudne do sprawdzenia • Małe odstępstwa niegroźne • Liniowość oceniamy na oko – wykresy rozrzutu • Co jeżeli zależność nie jest liniowa? • transformacja danych • regresja nieliniowa

  14. Regresja mówi o współwystępowaniu zjawisk, a nie o przyczynach i skutkach! Regresja liczby kradzieży na liczbę policjantów jest dodatnia!

  15. Zwodnicza regresja Jeżeli w modelu umieścisz dużą liczbę zmiennych objaśniających część z nich na pewno będzie istotna. Im więcej danych tym mniej złudne są wyniki. Ile?

  16. Reszty mają rozkład normalny • Ważne przy testowaniu, nie przy szacowaniu • Stosujemy histogramy reszt i wykresy normalności reszt • Niewielkie odchylenia nie są groźne • Dobry model daje duży R-kwadrat i normalność reszt. Czy taki potrafimy znaleźć?

  17. Nadmiarowość danych • % upadków przy transporcie tak samo dobrze opisuje liczba przejechanych kilometrów jak i dystans do ubojni (to to samo) • Statistica daje sygnał o „złym uwarunkowaniu macierzy”, ale nie zawsze.

  18. Odstające obserwacje Znacznie przekłamują oszacowania. Najczęściej to błędy powstałe przy wpisywaniu danych. Najlepiej usunąć je przed analizą regresji.

  19. Przykładowe dane

  20. Analiza bez interakcji z wieloma zmiennymi objaśniającymi ciągłymi

  21. Który model jest najlepszy? • Ten z max. R-kwadrat? Im szerszy model tym większy wsp. determ. • Lepiej stosować statystykę FjSSEj / MSE

  22. Jak wyznaczyć najlepszy model? • Najlepiej rozpatrywać każdy model z osobna (Statistica tego nie ułatwia). • Jeżeli więcej zmiennych kandydujących to korzystamy z regresji krokowej • wstecznej • postępującej

  23. R. krokowa postępująca • Najprostszy model – tylko wyraz wolny • Testujemy każdy z osobna, i dodajemy do modelu zmienną, której F>Fwprow. • Kontynuujemy - wprowadzamy następne zmienne i usuwamy te, dla których F<Fusun.

  24. y = -0,27 + 0,04  dystans + 0,45  obsada BETA to współ., po standaryzacji wszystkich zmiennych na średnią 0 i odch. std=1. Wskazują relatywny wkład każdej zmiennej.

  25. Po dopasowaniu modelu zawsze analizujemy reszty. Powtarzamy analizę jeżeli mamy duże wartości odstające!

  26. Wykres reszt wg przypadków Jeśli jedna lub więcej reszt wychodzi poza granicę ±3s usuwamy dane i powtarzamy analizę.

  27. Odległości Mahalanobisa • Odległości przypadków przewidywanych od średniej przewidywań. • Sposób na analizę wartości ekstremalnych po stronie zmiennych objaśniających.

  28. Wykres odstających wg. przypadków

  29. Reszty usunięte ...to reszty jakie byśmy uzyskali, gdyby dany przypadek pominąć przy obliczeniach regresji. Jeśli reszta usunięta znacznie się różni od zwykłej reszty standaryzowanej, to dany przypadek przekłamuje całą analizę!

  30. Reszty wzg. usuniętych reszt

  31. Problem na ćwiczenia udział tłuszczu w ciele Ważny dla zdrowia Trudny w pomiarze - wymaga ważenia ciała w wodzie. Czy można go przewidzieć na podstawie łatwych pomiarów.

  32. Density determined from underwater weighing Percent body fat from Siri's (1956) equation Age (years) Weight (lbs) Height (inches) Neck circumference (cm) Chest circumference (cm) Abdomen 2 circumference (cm) Hip circumference (cm) Thigh circumference (cm) Knee circumference (cm) Ankle circumference (cm) Biceps (extended) circumference (cm) Forearm circumference (cm) Wrist circumference (cm) dane BODYFAT http://lib.stat.cmu.edu/datasets/

  33. Zadania na ćwiczenia Skonstruuj dobry model predykcji udziału tłuszczu w ciele Które zmienne są najlepiej objaśniają udział tłuszczu w ciele człowieka? Dokonaj analizy reszt. Usuń przypadki zniekształcające przewidywanie i popraw model.

  34. zadanie dla chętnychZbrodnie Detroit http://lib.stat.cmu.edu/datasets/detroit The data are on the homicide rate in Detroit for the years 1961-1973. FTP - Full-time police per 100,000 population UEMP - % unemployed in the population MAN - number of manufacturing workers in thousands LIC - Number of handgun licences per 100,000 population GR - Number of handgun registrations per 100,000 population CLEAR - % homicides cleared by arrests WM - Number of white males in the population NMAN - Number of non-manufacturing workers in thousands GOV - Number of government workers in thousands HE - Average hourly earnings WE - Average weekly earnings HOM - Number of homicides per 100,000 of population ACC - Death rate in accidents per 100,000 population ASR - Number of assaults per 100,000 population skonstruuj model predykcji liczby zabójstw

More Related