Re: [Rd] gaussian family change suggestion

2006-04-11 Thread Simon Wood
> >>> Currently the `gaussian' family's initialization code signals an error if > >>> any response data are zero or negative and a log link is used. Given that > >>> zero or negative response data are perfectly legitimate under the GLM > >>> fitted using `gaussian("log")', this seems a bit unsatisf

Re: [Rd] gaussian family change suggestion

2006-04-11 Thread Simon Wood
> > Currently the `gaussian' family's initialization code signals an error if > > any response data are zero or negative and a log link is used. Given that > > zero or negative response data are perfectly legitimate under the GLM > > fitted using `gaussian("log")', this seems a bit unsatisfactory.

Re: [Rd] gaussian family change suggestion

2006-04-11 Thread Prof Brian Ripley
On Tue, 11 Apr 2006, Simon Wood wrote: >>> Currently the `gaussian' family's initialization code signals an error if >>> any response data are zero or negative and a log link is used. Given that >>> zero or negative response data are perfectly legitimate under the GLM >>> fitted using `gaussian("l

Re: [Rd] gaussian family change suggestion

2006-04-11 Thread Thomas Lumley
On Tue, 11 Apr 2006, Prof Brian Ripley wrote: > On Tue, 11 Apr 2006, Simon Wood wrote: > >> Hi, >> >> Currently the `gaussian' family's initialization code signals an error if >> any response data are zero or negative and a log link is used. Given that >> zero or negative response data are perfect

Re: [Rd] gaussian family change suggestion

2006-04-11 Thread Prof Brian Ripley
On Tue, 11 Apr 2006, Simon Wood wrote: > Hi, > > Currently the `gaussian' family's initialization code signals an error if > any response data are zero or negative and a log link is used. Given that > zero or negative response data are perfectly legitimate under the GLM > fitted using `gaussian("l

[Rd] gaussian family change suggestion

2006-04-11 Thread Simon Wood
Hi, Currently the `gaussian' family's initialization code signals an error if any response data are zero or negative and a log link is used. Given that zero or negative response data are perfectly legitimate under the GLM fitted using `gaussian("log")', this seems a bit unsatisfactory. Might it be