Czy w modelu logit jest lepszy sposób na określenie wpływu niezależnej zmiennej porządkowej niż użycie zmiennych zastępczych dla każdego poziomu?
logistic
logit
ordinal-data
fgregg
źródło
źródło
Odpowiedzi:
Aby dodać do odpowiedzi @ dmk38, „dowolny zestaw wyników daje prawidłowy test, pod warunkiem, że są one konstruowane bez konsultacji z wynikami eksperymentu. Jeśli zestaw wyników jest słaby, to źle zniekształca skalę numeryczną, która naprawdę leży u podstaw uporządkowana klasyfikacja, test nie będzie wrażliwy. Dlatego wyniki powinny zawierać najlepszy możliwy wgląd w sposób, w jaki klasyfikacja została zbudowana i zastosowana. ” (Cochran, 1954, cytowany przez Agresti, 2002, s. 88–89). Innymi słowy, traktowanie uporządkowanego czynnika jako zmiennej punktowanej jest jedynie kwestią modelowania. Pod warunkiem, że ma to sens, wpłynie to tylko na sposób interpretacji wyniku i nie ma ostatecznej reguły wyboru najlepszego reprezentacji zmiennej porządkowej.
Rozważ następujący przykład dotyczący spożycia alkoholu przez matkę oraz obecności wrodzonych wad rozwojowych lub ich braku (Agresti, analiza danych kategorycznych , tabela 3.7, str. 89):
W tym konkretnym przypadku możemy modelować wynik za pomocą regresji logistycznej lub prostej tabeli asocjacji. Zróbmy to w R:
Zwykła statystyka (12,08, p = 0,016751) lub LR (6,20, p = 0,184562) (z 4 df) nie uwzględnia uporządkowanych poziomów spożycia alkoholu.χ2)
Traktując obie zmienne jako porządkowe z równomiernie rozmieszczonymi wynikami (nie ma to wpływu na zmienne binarne, takie jak zniekształcenie, i wybieramy linię bazową jako 0 = nieobecna), moglibyśmy przetestować liniowy przez asocjację liniową. Najpierw stwórzmy rozłożoną wersję tej tabeli awaryjnej:
Następnie możemy przetestować połączenie liniowe za pomocą
Jak widać, niewiele jest dowodów na wyraźne powiązanie między tymi dwiema zmiennymi. Jak uczyniło Agresti, jeśli zdecydujemy się przekodować poziomy alkoholu na {0,0,5,1,5,4,7}, czyli wykorzystujemy wartości średnie dla hipotetycznej skali ciągłej, przy czym ostatni wynik jest nieco arbitralny, to doszlibyśmy do wniosku na większy wpływ spożywania alkoholu przez matkę na rozwój wrodzonych wad rozwojowych:
daje statystykę testową 6,57 z powiązaną wartością p 0,01037.
Stosując podejście GLM, postępowalibyśmy w następujący sposób. Ale najpierw sprawdź, jak alkohol jest kodowany w R:
Jest to prosty nieuporządkowany czynnik (
"factor"
), stąd nominalny predyktor. Oto trzy modele, w których uważamy alkohol za predyktor nominalny, porządkowy lub ciągły.Jednak w dwóch innych przypadkach (
mod1
imod2
) otrzymujemy różne dane wyjściowe, ponieważ macierz projektowa zastosowana do modelowania predyktora jest różna, co można sprawdzić za pomocą:mod1
mod2
mod1
mod2
Możesz spróbować przypisać nowe wyniki do Alkoholu i zobaczyć, jak wpłynie to na przewidywane prawdopodobieństwo wad rozwojowych.
źródło
doskonale jest używać predyktora jakościowego w modelu regresji logit (lub OLS), jeśli poziomy są porządkowe. Ale jeśli masz powód, aby traktować każdy poziom jako dyskretny (lub jeśli w rzeczywistości twoja zmienna kategorialna jest nominalna, a nie porządkowa), to jako alternatywę dla kodowania pozornego możesz również użyć ortogonalnego kodowania kontrastowego. Aby uzyskać bardzo kompletną i dostępną dyskusję, patrz Judd, CM, McClelland, GH & Ryan, CS Analiza danych: metoda porównywania modeli, Edn. 2. miejsce (Routledge / Taylor and Francis, New York, NY; 2008), lub po prostu Google „kodowanie kontrastowe”
źródło