Re: R-alpha: glm and substitute bugs

Kurt Hornik (Kurt.Hornik@ci.tuwien.ac.at)
Wed, 15 Jan 1997 08:35:37 +0100


Date: Wed, 15 Jan 1997 08:35:37 +0100
Message-Id: <199701150735.IAA31557@aragorn.ci.tuwien.ac.at>
From: Kurt Hornik <Kurt.Hornik@ci.tuwien.ac.at>
To: Thomas Lumley <thomas@biostat.washington.edu>
Subject: Re: R-alpha: glm and substitute bugs
In-Reply-To: <Pine.SUN.3.91.970114182458.4333A-100000@chekov>
	<Pine.SUN.3.91.970114182458.4333A-100000@chekov>

>>>>> On Tue, 14 Jan 1997 18:32:05 -0800 (PST),
>>>>> Thomas Lumley <thomas@biostat.washington.edu> said:

> On Mon, 13 Jan 1997, Jim Lindsey wrote:
>> 1. Kurt Hornik has asked me to look at a  data set where glm with
>> binomial gives all zero coefficients and S-Plus does not. It turns out
>> that R gives all zero coefficients (except the intercept) for all
>> binomial data sets! 

> I'm not getting this problem:

> shuksan% uname -a
> SunOS shuksan 5.5.1 Generic_103640-03 sun4m sparc SUNW,SPARCstation-20
> shuksan% ~/R15
> R : Copyright 1996, Robert Gentleman and Ross Ihaka
> Version 0.15  (Released: December 25, 1996)

> R is free software and comes with ABSOLUTELY NO WARRANTY.
> You are welcome to redistribute it under certain conditions.
> Type "license()" for details.

R> y<-rep(0:1,20)
R> x<-rnorm(40)
R> glm(y~x,family=binomial())

> Call:  glm(y ~ x, family = binomial())

> ...

Guess Jim meant something like `for all binomial data sets close to the
one Kurt sent'.  It definitely did not work on that one, and now does
following Ross's hint about commenting out the lines in coerce.c.

I only had the problems in cases where y is not a binary vector, but a
matrix of successes and errors ...

Jim?

-k
=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-
r-testers mailing list -- For info or help, send "info" or "help",
To [un]subscribe, send "[un]subscribe"
(in the "body", not the subject !)  To: r-testers-request@stat.math.ethz.ch
=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-