[R] ANCOVA/glm missing/ignored interaction combinations

lara harrup (IAH-P) lara.harrup at bbsrc.ac.uk
Wed Sep 3 10:46:54 CEST 2008


Hi

I am using R version 2.7.2. on a windows XP OS and have a question
concerning an analysis of covariance with count data I am trying to do,
I will give details of a scaled down version of the analysis (as I have
more covariates and need to take account of over-dispersion etc etc) but
as I am sure it is only a simple problem but I just can't see how to fix
it.

I have a data set with count data as the response (total) and two
continuous covariates and one categorical explanatory variable (semio).
When I run the following lines, after loading the data and assigning
'semio' as a factor, taking into account that I want to consider two way
interactions:

> model.a<-glm(total~(temperature+humidity+semio)^2,family=poisson)
> summary(model.a)

I get the output below. But not all interactions are listed there are 4
semio categories, 1,2,3 and 4 but only 2,3,and 4 are listed in the
summary (semio2,semio3,semio4). And I cant for the life of me work out
why category one (semio1) is being ignored, missing etc.

Any help or suggestions would be most appreciated. Thanks in advance

Lara
lara.harrup at bbsrc.ac.uk

Call:
glm(formula = total ~ (temperature + humidity + semio)^2, family =
poisson)

Deviance Residuals: 
    Min       1Q   Median       3Q      Max  
-22.212   -5.132   -2.484    3.200   18.793  

Coefficients:
                      Estimate Std. Error z value Pr(>|z|)    
(Intercept)          23.848754   2.621291   9.098  < 2e-16 ***
temperature          -1.038284   0.150465  -6.901 5.18e-12 ***
humidity             -0.264912   0.033928  -7.808 5.81e-15 ***
semio2               22.852664   1.291806  17.690  < 2e-16 ***
semio3                3.699901   1.349007   2.743   0.0061 ** 
semio4               -1.851163   1.585997  -1.167   0.2431    
temperature:humidity  0.012983   0.001983   6.545 5.94e-11 ***
temperature:semio2   -0.870430   0.037602 -23.149  < 2e-16 ***
temperature:semio3   -0.060846   0.038677  -1.573   0.1157    
temperature:semio4    0.055288   0.046137   1.198   0.2308    
humidity:semio2      -0.114718   0.013369  -8.581  < 2e-16 ***
humidity:semio3      -0.031692   0.013794  -2.298   0.0216 *  
humidity:semio4       0.008425   0.016020   0.526   0.5990    
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 

(Dispersion parameter for poisson family taken to be 1)

    Null deviance: 10423.5  on 47  degrees of freedom Residual deviance:
2902.2  on 35  degrees of freedom
AIC: 3086.4

Number of Fisher Scoring iterations: 7



More information about the R-help mailing list