[R] REML option in gstat

[ripley@stats.ox.ac.uk]

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