[Rd] Apparant bug in binomial model in GLM (PR#13434)

soren.faurby at biology.au.dk soren.faurby at biology.au.dk
Tue Jan 6 12:25:04 CET 2009


Full_Name: Søren Faurby
Version: 2.4.1 and 2.7.2
OS: 
Submission from: (NULL) (192.38.46.92)


There appear to be a bug in the estimation of significance in the binomial model
in GLM. This bug apparently appears when the correlation between two variables
is to strong.

Such as this dummy example
c(0,0,0,0,0,1,1,1,1,1)->a
a->b
m1<-glm(a~b, binomial)
summary(m1)

It is sufficient that all 1's correspond to 1's such as this example

c(0,0,0,0,0,1,1,1,1,1)->a
c(0,0,0,0,1,1,1,1,1,1)->c
m1<-glm(a~c, binomial)
summary(m1)

I hope that this message is understandable. 

Kind regards, Søren



More information about the R-devel mailing list