Re: [R-sig-Geo] Bug in surf.ls (spatial package)?

2007-04-10 Thread ONKELINX, Thierry
Onderwerp: Re: [R-sig-Geo] Bug in surf.ls (spatial package)? On Tue, 10 Apr 2007, ONKELINX, Thierry wrote: > > -Oorspronkelijk bericht- > > Van: Prof Brian Ripley [mailto:[EMAIL PROTECTED] > > Verzonden: dinsdag 10 april 2007 16:20 > > Aan: ONKELINX, T

Re: [R-sig-Geo] Bug in surf.ls (spatial package)?

2007-04-10 Thread Roger Bivand
On Tue, 10 Apr 2007, ONKELINX, Thierry wrote: > > -Oorspronkelijk bericht- > > Van: Prof Brian Ripley [mailto:[EMAIL PROTECTED] > > Verzonden: dinsdag 10 april 2007 16:20 > > Aan: ONKELINX, Thierry > > CC: r-sig-geo@stat.math.ethz.ch > > Onderwerp: Re: Bug in surf.ls (spatial package)? >

Re: [R-sig-Geo] Bug in surf.ls (spatial package)?

2007-04-10 Thread ONKELINX, Thierry
> -Oorspronkelijk bericht- > Van: Prof Brian Ripley [mailto:[EMAIL PROTECTED] > Verzonden: dinsdag 10 april 2007 16:20 > Aan: ONKELINX, Thierry > CC: r-sig-geo@stat.math.ethz.ch > Onderwerp: Re: Bug in surf.ls (spatial package)? > > On Tue, 10 Apr 2007, ONKELINX, Thierry wrote: > > > Dea

Re: [R-sig-Geo] Bug in surf.ls (spatial package)?

2007-04-10 Thread Prof Brian Ripley
On Tue, 10 Apr 2007, ONKELINX, Thierry wrote: > Dear listers, > > When fitting a trend surface I noticed that surf.ls() generates other > values for the model parameters than an ordinary lm(). Since both fit, > in my opinion, the same model to the same data. So I would suppect that > both would ge

Re: [R-sig-Geo] Bug in surf.ls (spatial package)?

2007-04-10 Thread Gabor Grothendieck
They are using different coefficient parameterizations but the models are otherwise the same. For example, they give the same fitted values: > all.equal(fitted(surface.surf), as.vector(fitted(surface.lm))) [1] TRUE On 4/10/07, ONKELINX, Thierry <[EMAIL PROTECTED]> wrote: > Dear listers, > > Whe

[R-sig-Geo] Bug in surf.ls (spatial package)?

2007-04-10 Thread ONKELINX, Thierry
Dear listers, When fitting a trend surface I noticed that surf.ls() generates other values for the model parameters than an ordinary lm(). Since both fit, in my opinion, the same model to the same data. So I would suppect that both would generate the same parameter values. Or am I missing somethin