Niedawno dopasowałem 4 modele regresji wielokrotnej dla tych samych danych predykcyjnych / odpowiedzi. Dwa modele pasuję do regresji Poissona.
model.pois <- glm(Response ~ P1 + P2 +...+ P5, family=poisson(), ...)
model.pois.inter <- glm(Response ~ (P1 + P2 +...+ P5)^2, family=poisson(), ...)
Dwa modele pasuję do ujemnej regresji dwumianowej.
library(MASS)
model.nb <- glm.nb(Response ~ P1 + P2 +...+ P5, ...)
model.nb.inter <- glm.nb(Response ~ (P1 + P2 +...+ P5)^2, ...)
Czy istnieje test statystyczny, którego można użyć do porównania tych modeli? Użyłem AIC jako miary dopasowania, ale AFAIK nie reprezentuje rzeczywistego testu.
model.nb.inter
jest znacznie lepsze niż dopasowanie model.pois.inter
. Tak, AIC jest niższy, ale o ile niższy jest znacznie lepszy ?