R help - Mar 2003 - REML option in gstat

Hi, please help!!

I've been trying to fit variogram models using the REML method in the gstat 
package.  Every time the Windows GUI crashes.  For example

library(gstat)
data(meuse)
x <- variogram(zinc ~ 1, ~x + y, meuse)
v <- vgm(140000, "Sph", 800, nug = 10000)
plot(x, model = fit.variogram(x, model = v, fit.method=5))

Other fit methods are non problematic (eg. fit.method=7 for WLS or 
fit.method=1 for OLS)

I've tried the code on Windows XP, 98 and NT and all fail, how about 
non-Windows platforms?
Any suggestions?

Many thanks
Dave

That example segfaults on Linux.

On Fri, 7 Mar 2003, David wrote:
> I've been trying to fit variogram models using the REML method in the
gstat
> package.  Every time the Windows GUI crashes.  For example
> 
> library(gstat)
> data(meuse)
> x <- variogram(zinc ~ 1, ~x + y, meuse)
> v <- vgm(140000, "Sph", 800, nug = 10000)
> plot(x, model = fit.variogram(x, model = v, fit.method=5))
> 
> Other fit methods are non problematic (eg. fit.method=7 for WLS or 
> fit.method=1 for OLS)
> 
> I've tried the code on Windows XP, 98 and NT and all fail, how about 
> non-Windows platforms?
> Any suggestions?
See Q7.3 in the rw-FAQ on how to debug compiled code.  Since you have
already re-built from the sources, this should be easy for you to do.

Alternatively, talk to the author about this.

-- 
Brian D. Ripley,                  ripley at stats.ox.ac.uk
Professor of Applied Statistics,  http://www.stats.ox.ac.uk/~ripley/
University of Oxford,             Tel:  +44 1865 272861 (self)
1 South Parks Road,                     +44 1865 272866 (PA)
Oxford OX1 3TG, UK                Fax:  +44 1865 272595