Dlaczego centrowanie zmiennych niezależnych może zmieniać główne efekty z umiarem?


28

Mam pytanie związane z regresją wielokrotną i interakcją, zainspirowane tym wątkiem CV: Pojęcie interakcji za pomocą analizy hierarchicznej regresji zmiennych centrowanych? Jakie zmienne powinniśmy wyśrodkować?

Podczas sprawdzania efektu moderacji centruję zmienne niezależne i mnożę zmienne wyśrodkowane, aby obliczyć termin interakcji. Następnie przeprowadzam analizę regresji i sprawdzam efekty główne i interakcyjne, które mogą wykazywać moderację.

Powtórzyć, jeśli analizę bez centrowania najwyraźniej współczynnik korelacji ( R2 ) nie ulegają zmianie, ale współczynniki regresji ( β e) do. To wydaje się jasne i logiczne.

Czego nie rozumiem: Wartości p głównych efektów zmieniają się zasadniczo wraz z centrowaniem, chociaż interakcja nie (co jest słuszne). Tak więc moja interpretacja głównych efektów może się diametralnie zmienić - po prostu określona przez centrowanie, czy nie. (To wciąż te same dane, w obu analizach!)

Czy ktoś może to wyjaśnić? - Ponieważ oznaczałoby to, że opcja wyśrodkowania moich zmiennych byłaby obowiązkowa i każdy powinien to zrobić, aby uzyskać te same wyniki z tymi samymi danymi.


Bardzo dziękuję za wyjaśnienie tego problemu i wyczerpujące wyjaśnienia. Zapewniamy, że twoja pomoc jest bardzo doceniana!

Dla mnie największą zaletą centrowania jest unikanie wielokoliniowości. Nadal dość mylące jest ustalanie reguły, czy centrować czy nie. Mam wrażenie, że większość zasobów sugeruje wyśrodkowanie, chociaż przy tym jest pewne „ryzyko”. Ponownie chcę podkreślić fakt, że 2 badaczy zajmujących się tym samym materiałem i danymi może wyciągać różne wyniki, ponieważ jeden centruje, a drugi nie. Właśnie przeczytałem część książki Bortza (był profesorem i swego rodzaju gwiazdą statystyki w Niemczech i Europie) i nawet nie wspomina o tej technice; wskazuje tylko, że należy zachować ostrożność przy interpretacji głównych efektów zmiennych, gdy są one zaangażowane w interakcje.

W końcu, jeśli przeprowadzasz regresję z jednym IV, jednym moderatorem (lub drugim IV) i DV, czy zaleciłbyś wyśrodkowanie, czy nie?


5
Prawie nigdy nie używam centrowania, uważając to za całkowicie niepotrzebne i mylące.
Frank Harrell,

3
Przeczytaj uważnie odpowiedzi ponownie. Twoje wnioski nie zmieniają się, gdy centrujesz zmienne niezależne lub nie stosujesz transformacji liniowej - jeśli są one poprawnie narysowane. Wielokoliniowość ze względu na brak centrowania jest czysto numerycznym problemem i jest rozwiązywana automatycznie przez każde przyzwoite oprogramowanie.
Scortchi - Przywróć Monikę

1
Zjawisko (p-wartościach zmieniających) mogą być rozumiane jako konsekwencja kwadratowego rodzaju interakcje, jak wyjaśniono w stats.stackexchange.com/questions/28730/... .
whuber

Odpowiedzi:


23

W modelach bez terminów interakcji (to znaczy bez terminów konstruowanych jako iloczyn innych terminów) współczynnik regresji każdej zmiennej jest nachyleniem powierzchni regresji w kierunku tej zmiennej. Jest stały, niezależnie od wartości zmiennych, i dlatego można powiedzieć, że mierzy ogólny efekt tej zmiennej.

W modelach z interakcjami tej interpretacji można dokonać bez dalszej kwalifikacji tylko dla tych zmiennych, które nie są zaangażowane w żadne interakcje. Dla zmiennej, która bierze udział w interakcjach, współczynnik regresji „głównego efektu” - to znaczy współczynnik regresji samej zmiennej - jest nachyleniem powierzchni regresji w kierunku tej zmiennej, gdy wszystkie inne zmienne, które oddziaływanie z tą zmienną ma wartości zerowe , a test istotności współczynnika odnosi się do nachylenia powierzchni regresji tylko w tym obszarze przestrzeni predyktora. Ponieważ nie ma wymogu, aby faktycznie istniały dane w tym obszarze przestrzeni, współczynnik efektu głównego może w niewielkim stopniu przypominać nachylenie powierzchni regresji w obszarze przestrzeni predyktora, w którym faktycznie zaobserwowano dane.

Innymi słowy, współczynnik efektu głównego jest analogiczny do prostego efektu głównego, a nie ogólnego efektu głównego. Co więcej, może odnosić się do tego, co w projekcie anova byłoby pustymi komórkami, w których dane były dostarczane przez ekstrapolację z komórek z danymi.

Aby zmierzyć ogólny efekt zmiennej, który jest analogiczny do ogólnego efektu głównego w anova i nie dokonuje ekstrapolacji poza region, w którym zaobserwowano dane, musimy spojrzeć na średnie nachylenie powierzchni regresji w kierunku zmiennej , gdzie uśrednianie jest dla N zaobserwowanych przypadków. To średnie nachylenie może być wyrażone jako ważona suma współczynników regresji wszystkich terminów w modelu, które obejmują daną zmienną.

Ciężary są trudne do opisania, ale łatwe do zdobycia. Współczynnik głównego efektu zmiennej zawsze otrzymuje wagę 1. Dla każdego innego współczynnika terminu obejmującego tę zmienną, waga jest średnią iloczynu innych zmiennych tego terminu. Na przykład, jeśli mamy pięć „surowych” zmiennych x1, x2, x3, x4, x5, plus cztery interakcje dwukierunkowe (x1,x2), (x1,x3), (x2,x3), (x4,x5)i jedną interakcję trójstronną (x1,x2,x3), wówczas model jest

y = b0 + b1*x1 + b2*x2 + b3*x3 + b4*x4 + b5*x5 +
    b12*x1*x2 + b13*x1*x3 + b23*x2*x3 + b45*x4*x5 +
    b123*x1*x2*x3 + e

a ogólne główne efekty to

B1 = b1 + b12*M[x2] + b13*M[x3] + b123*M[x2*x3],

B2 = b2 + b12*M[x1] + b23*M[x3] + b123*M[x1*x3],

B3 = b3 + b13*M[x1] + b23*M[x2] + b123*M[x1*x2],

B4 = b4 + b45*M[x5],

B5 = b5 + b45*M[x4],

gdzie M [.] oznacza średnią próbki ilości w nawiasach. Wszystkie warunki produktu w nawiasach należą do tych, które zostały zbudowane w celu przeprowadzenia regresji, więc program regresji powinien już o nich wiedzieć i być w stanie wydrukować swoje środki na żądanie.

W modelach, które mają tylko główne efekty i interakcje dwukierunkowe, istnieje prostszy sposób na uzyskanie ogólnych efektów: wyśrodkuj [1] surowe zmienne na ich średnich. Należy to zrobić przed obliczeniem warunków produktu i nie można tego zrobić w odniesieniu do produktów. Wtedy wszystkie wyrażenia M [.] Staną się 0, a współczynniki regresji będą interpretowane jako efekty ogólne. Wartości b's zmienią się; wartości B nie będą. Tylko zmienne, które są zaangażowane w interakcje, muszą być wyśrodkowane, ale zwykle nie ma szkody w centrowaniu innych mierzonych zmiennych. Ogólnym efektem centrowania zmiennej jest to, że oprócz zmiany punktu przecięcia zmienia tylko współczynniki innych zmiennych, które oddziałują ze zmienną centrowaną. W szczególności, nie zmienia współczynników żadnych terminów, które dotyczą zmiennej centrowanej. W powyższym przykładzie centrowanie x1 zmieniłoby b0, b2, b3 i b23.

[1 - „Centrowanie” jest używane przez różnych ludzi w sposób, który różni się tylko na tyle, aby spowodować zamieszanie. W użytym tutaj znaczeniu „wyśrodkowanie zmiennej na #” oznacza odjęcie # od wszystkich wyników na zmiennej, przekształcając oryginalne wyniki na odchylenia od #.]

Dlaczego więc nie zawsze rutynowo koncentrować się na środkach? Trzy powody. Po pierwsze, interesujące mogą być same współczynniki efektu głównego zmiennych niecentrowanych. Centrowanie w takich przypadkach przyniosłoby efekt przeciwny do zamierzonego, ponieważ zmienia współczynniki efektu głównego innych zmiennych.

Po drugie, centrowanie spowoduje, że wszystkie wyrażenia M [.] 0, a tym samym przekształci proste efekty w ogólne, tylko w modelach bez interakcji trójstronnych lub wyższych . Jeśli model zawiera takie interakcje, należy wykonać obliczenia b -> B, nawet jeśli wszystkie zmienne są wyśrodkowane na ich średnich wartościach.

Po trzecie, centrowanie na wartości takiej jak średnia, która jest zdefiniowana przez rozkład predyktorów w przeciwieństwie do racjonalnego wyboru, oznacza, że ​​wszystkie współczynniki, na które ma wpływ centrowanie, będą specyficzne dla konkretnej próbki. Jeśli koncentrujesz się na średniej, to ktoś, kto próbuje powtórzyć twoje badanie, musi skoncentrować się na średniej, a nie na własnej średniej, jeśli chce uzyskać te same współczynniki, które masz. Rozwiązaniem tego problemu jest wyśrodkowanie każdej zmiennej na racjonalnie wybranej wartości centralnej tej zmiennej, która zależy od znaczenia wyników i nie zależy od rozkładu wyników. Jednak obliczenia b -> B nadal są konieczne.

Istotność ogólnych efektów można sprawdzić za pomocą zwykłych procedur testowania liniowych kombinacji współczynników regresji. Jednak wyniki należy interpretować ostrożnie, ponieważ ogólne efekty nie są parametrami strukturalnymi, ale zależą od projektu. Parametry strukturalne - współczynniki regresji (niecentryczne lub z racjonalnym centrowaniem) i wariancja błędu - mogą pozostać niezmienne przy zmianach rozkładu predyktorów, ale ogólne efekty na ogół się zmienią. Ogólne efekty są specyficzne dla konkretnej próbki i nie należy oczekiwać, że zostaną przeniesione na inne próbki o różnych rozkładach predyktorów. Jeśli ogólny efekt jest znaczący w jednym badaniu, a nie w innym, może odzwierciedlać jedynie różnicę w rozkładzie predyktorów.


10

β

y=β1x1+β2x2+β3x1x2+ϵβ1x1β3x1x2x1x1x2β

ββ1yx1 x2=0x1yx2β1x2

βx1yx2yx1x2


-1

Oszalałem z tego samego pytania, ale w końcu znalazłem rozwiązanie twojego i mojego problemu. TO WSZYSTKO, JAK OBLICZASZ SWOJE CENTRALNE ZMIENNE. Dostępne są dwie opcje:
1. ZNACZENIE - ZMIENNE INDYWIDUALNE 2. ZMIENNE INDYWIDUALNE - ZNACZONE
Prawdopodobnie zmienne wyśrodkowane obliczono jako (zmienna indywidualna - wartość średnia) , dlatego te o niskich wartościach uzyskałyby wyniki ujemne, a te o wysokich wartościach - dodatnie wyniki.
Wyjaśnię na przykładzie, aby ułatwić zrozumienie. Chcę zobaczyć, jak siła mięśni wpływa na masę kości i chcę wziąć pod uwagę płeć, aby zobaczyć, czy wpływa to inaczej u dziewcząt i chłopców. Chodzi o to, że im wyższa siła mięśni, tym wyższa masa kostna. Dlatego mam:

Zmienna zależna: masa kości Zmienne niezależne: płeć, siła mięśni, interakcja_SEKST_MUSCYLE siła.

Jak znalazłem wielokoliniowość (zazwyczaj dzieje się tak, gdy masz termin interakcji), wyśrodkowałem siłę mięśni (ŚREDNIA - ZMIENNA INDYWIDUALNA) i utworzyłem nowy termin interakcji z nową zmienną centrowaną. Moje współczynniki były

Stała: 0.902
Płeć: -0.010(Chłopcy = 0; Dziewczyny = 1)
Wyśrodkowany mięsień: -0.023
Interakcja: 0.0002
Dlatego jeśli chcesz oszacować masę kości chłopca, będziesz mieć następujące równanie:
Masa kości =0.902(00.010)(0.023musclecentredvalue)+(Interaction0.0002)

Patrząc na to, możesz pomyśleć, że mięsień wpływa negatywnie na kość, ale musisz pomyśleć o swoich środkowych zmiennych, a nie o pierwotnych zmiennych. Powiedzmy, że średnia siła grupy była 30 KG. I chcesz oszacować masę kości chłopca (WEAKBOY), który wykonał, 20 KGa drugiego, który wykonał 40KG(STRONGBOY). Wyśrodkowane wartości WEAKBOY będą wynosić (WARTOŚĆ GRUPY MEAN - WARTOŚĆ INDYWIDUALNA; 30-20 = 10), a dla STRONGBOY będzie wynosić -10. Zastosowanie tych wartości do równania:

WEAKBOY Masa kości = 0,902 - 0 - (0,023 * 10) + .... = 0,672

STRONGBOY Masa kości = 0,902 - (0,023 * (- 10)) + ... = 1,132

Jak widać STRONGBOY rzeczywiście miał silniejszą kość. Jeśli wyśrodkowałeś swoje zmienne na odwrót: (INDYWIDUALNY - ŚREDNI), wszystkie współczynniki będą takie same, ale symbole będą różne. Wynika to z tego, że po zastosowaniu zmiennej centralnej WEAKBOY będzie (-10), a STRONGBOY będzie (+10). Dlatego końcowe wyniki będą dokładnie takie same.

To wszystko ma sens, kiedy to zrozumiesz.

Mam nadzieję, że przykład jest wystarczająco jasny.


Ten błąd nie wyjaśniłby zmian wartości p. BTW, twoja opcja (1) nie jest wyśrodkowana, ponieważ obejmuje także pomnożenie wartości przez stałą. (Stała wynosi -1.)
whuber
Korzystając z naszej strony potwierdzasz, że przeczytałeś(-aś) i rozumiesz nasze zasady używania plików cookie i zasady ochrony prywatności.
Licensed under cc by-sa 3.0 with attribution required.