I have never really understood of what use is the null deviance
printed out by glm or how it could be interpreted. No such idea exists
in GLIM. Here is exactly the same model fitted 4 time with 4 different
answers.

y <- rbinom(20,1,0.5)
ind <- rep(1,20)
glm(y~1,family=binomial)
glm(y~ind-1,family=binomial)
glm(y~ind-1,family=binomial(link=log))
glm(y~ind-1,family=binomial(link=cloglog))

Jim
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !)  To: r-help-request@stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._

