[R] REML option in gstat

[David]

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