[R] REML option in gstat

David d.r.j.pleydell at pgr.salford.ac.uk
Fri Mar 7 20:10:19 CET 2003

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

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

More information about the R-help mailing list