From thomas@biostat.washington.edu  Tue Mar  7 19:17:30 2000
From: thomas@biostat.washington.edu (Thomas Lumley)
Date: Tue, 7 Mar 2000 11:17:30 -0800 (PST)
Subject: [Rd] lm(rnorm(1000)~rnorm(1000)) kills rw1000 (PR#476)
In-Reply-To: <Pine.WNT.4.05.10003071851000.98-100000@tern.stats>
Message-ID: <Pine.GSO.4.21.0003071108270.27319-100000@mail.biostat.washington.edu>

On Tue, 7 Mar 2000, Prof Brian D Ripley wrote:

> On Tue, 7 Mar 2000 cberry@tajo.ucsd.edu wrote:
> 
> > It does not happen on rw0901 or rw0651 (? version number ?) on my HOME
> > PC. Nor on R 0.90.1 under Solaris
> 
> It does happen on R 1.0.0 on Solaris and on rw1000 on a real PC.

Something strange is happening where we ensure that the response does not
appear on the rhs of the formula.  

At lines 1643ff of model.c we count how many columns of the design matrix
are needed for each model term.  This count is not done if the term is the
response, which leaves the count uninitialised.

We then use this uninitialised count at 1673ff to construct the assign
attribute.

	-thomas



-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-devel 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-devel-request@stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._


[Rd] lm(rnorm(1000)~rnorm(1000)) kills rw1000 (PR#476)
Thomas Lumley