Re: R-beta: more glm bugs

Thomas Lumley (thomas@biostat.washington.edu)
Mon, 26 Jan 1998 08:34:35 -0800 (PST)


Date: Mon, 26 Jan 1998 08:34:35 -0800 (PST)
From: Thomas Lumley <thomas@biostat.washington.edu>
To: Jim Lindsey <jlindsey@luc.ac.be>
Subject: Re: R-beta: more glm bugs
In-Reply-To: <9801260749.AA11322@alpha.luc.ac.be>

On Mon, 26 Jan 1998, Jim Lindsey wrote:

>   1. At least on my computer, fitting a saturated gamma model with glm
> gives a core dump:
> glm(rgamma(20,1,1)~gl(20,1,20), family=Gamma)

I can't reproduce this problem on 0.50-a4, 0.60 or 0.61.1.

>   2. Factor variables are not vectors!!
> is.vector(gl(20,2,40)) gives FALSE

Even stranger 
R> is.vector(factor(1:3))
[1] FALSE
R> is.factor(factor(1:3))
[1] TRUE
R> is.vector(as.vector(factor(1:3)))
[1] TRUE
R> is.factor(as.vector(factor(1:3)))
[1] TRUE
so that factor variables *can* be vectors, but aren't naturally.


Thomas Lumley
------------------------
Biostatistics		
Uni of Washington	
Box 357232		
Seattle WA 98195-7232	
------------------------

-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
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
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._