... związek jest nieliniowy, ale istnieje wyraźny związek między xiy, jak mogę przetestować powiązanie i oznaczyć jego naturę?
Jednym sposobem osiągnięcia tego byłoby dopasować jako częściowo parametrycznie przybliżony funkcji z wykorzystaniem, na przykład, uogólnione dodatków modelu i testowania, czy nie, że oszacowanie funkcjonalna jest stała, które nie wskazują na stosunek między i . Takie podejście uwalnia cię od konieczności wykonywania regresji wielomianowej i podejmowania czasami arbitralnych decyzji dotyczących kolejności wielomianu itp.x y xyxyx
W szczególności, jeśli masz obserwacje , możesz dopasować model:(Yi,Xi)
E(Yi|Xi)=α+f(Xi)+εi
i przetestuj hipotezę . Za pomocą tej funkcji można to zrobić . Jeśli Twój wynik jest predyktorem, możesz wpisać:H0:f(x)=0, ∀xR
gam()
y
x
library(mgcv)
g <- gam(y ~ s(x))
Wpisanie summary(g)
daje wynik testu hipotezy powyżej. Jeśli chodzi o charakter relacji, najlepiej byłoby to zrobić za pomocą fabuły. Jednym ze sposobów na to R
(przy założeniu, że powyższy kod został już wprowadzony)
plot(g,scheme=2)
Jeśli twoja zmienna odpowiedzi jest dyskretna (np. Binarna), możesz to uwzględnić w tym frameworku, dopasowując logistyczną GAM (w R
, dodajesz family=binomial
do swojego wywołania gam
). Ponadto, jeśli masz wiele predyktorów, możesz dołączyć wiele dodatków (lub zwykłe terminy liniowe) lub dopasować funkcje wielu zmiennych, np. jeśli masz predyktory . Złożoność relacji jest automatycznie wybierana przez krzyżową weryfikację, jeśli korzystasz z metod domyślnych, chociaż jest tutaj duża elastyczność - w razie zainteresowania przejrzyj plik pomocy .f(x,z)x, z
gam