Jeśli dopasuję moje dane do czegoś podobnego lm(y~a*b)
, w składni R, gdzie a
jest zmienną binarną i b
zmienną numeryczną, to a:b
termin interakcji jest różnicą między nachyleniem y~b
at a
= 0 i at a
= 1.
Teraz załóżmy, że relacje między y
i b
jest krzywoliniowy. Jeśli teraz pasuję lm(y~a*poly(b,2))
, to a:poly(b,2)1
jest zmiana zmiany y~b
warunkowej na poziomie a
jak wyżej, i a:poly(b,2)2
zmiana y~b^2
warunkowa na poziomie a
. Wymaga to trochę falowania ręcznego, ale jeśli którykolwiek z tych współczynników interakcji różni się znacznie od zera, mogę argumentować, że oznacza to a
wpływ nie tylko na przesunięcie w pionie, y
ale także na lokalizację piku i stromość podejścia do szczytu y~b+b^2
krzywej.
Co jeśli pasuję lm(y~a*bs(b,df=3))
? Jak mam interpretować a:bs(b,df=3)1
, a:bs(b,df=3)2
oraz a:bs(b,df=3)3
warunki? Czy są to pionowe przesunięcia y
splajnu, które można przypisać a
każdemu z trzech segmentów?