[R] Spatial Statistics: surf.gls

[Stegemann, Berthold]

 In an experimental setup we obtain z-data samples at equidistant grid points.

The surf.gls (Kriging) algorithm produces an error under this circumstance when performing the Choleski decomposition.

A workaround is to dither the grid coordinates using  (x <- rnorm(length(x)) ; y<- rnowm(length(y))).

Question: Is this an expected behaviour  of the  surf.gls function ?

Regards,
Berthold 

---------------------------------
Dr. Berthold Stegemann 
Senior Scientist
Bakken Research Center
Heart Failure Management
Endepolsdomein 5
6229 GW Maastricht
The Netherlands
Office: +31.43.3566.613
Mobile: +31.6.5335.1669
Fax: +31.43.3566.509