2
I have the following matrix:
DIV RIQ EXC FX_FREQ AMPL FQMIN FQMAX FQDOM DUR
BM 1.16 14 1 6 9950 250 10200 3254.4 173.77
CP 1.23 12 1 6 9840 260 10100 3417.5 84.06
MT 2.03 14 5 9 9955 160 10115 1359.8 60.70
From there I want to run a GLM. My global model is as follows (I am not using either RIQ or FX_FREQ):
M1<-glm(ADI ~ EXC+AMPL+FQMIN+FQMAX+FQDOM+DUR, family=poisson, data=env)
summary(M1)
When I call the Summary(M1) command, I get: Call:
glm(formula = ADI ~ EXC + AMPL + FQMIN + FQMAX + FQDOM + DUR,
family = poisson, data = env)
Deviance Residuals:
[1] 0 0 0
Coefficients: (4 not defined because of singularities)
Estimate Std. Error z value Pr(>|z|)
(Intercept) 5.308e+00 1.162e+02 0.046 0.964
EXC 1.406e-01 2.996e-01 0.469 0.639
AMPL -5.327e-04 1.177e-02 -0.045 0.964
FQMIN NA NA NA NA
FQMAX NA NA NA NA
FQDOM NA NA NA NA
DUR NA NA NA NA
(Dispersion parameter for poisson family taken to be 1)
Null deviance: 3.0247e-01 on 2 degrees of freedom
Residual deviance: -3.2641e-16 on 0 degrees of freedom
AIC: Inf
Number of Fisher Scoring iterations:
What bothers me is this "NA". I try to follow the analysis by creating new models, the criterion I use is: always in the new model, I do not insert the variable with the highest p value obtained in the previous model. Ex: my M2 is:
M2<-glm(ADI ~ EXC+FQMIN+FQMAX+FQDOM+DUR, family=poisson, data=env)
because in Summary(M1), the variable AMPL had higher p. I proceed until I have 6 models.
At the end I want to select based on Akaike (Aicc). Then I apply the command:
AICc <- ICtab(M1,M2,M3,M4,M5,M6,type = c("AICc"),base=T, weights = TRUE,
delta = TRUE, sort = TRUE, nobs = 4)
The biggest problem occurs when I call Aicc. What appears is:
AICc dAICc df weight
M1 Inf NaN 3 NA
M2 Inf NaN 3 NA
M3 Inf NaN 3 NA
M4 Inf NaN 3 NA
M5 Inf NaN 3 NA
M6 Inf NaN 2 NA
I need to fix this. Could you help me?