[R] REML option in gstat
ripley at stats.ox.ac.uk
Fri Mar 7 20:37:28 CET 2003
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
> 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
More information about the R-help