[R-sig-Geo] local Universal Kriging with krige()

[Sven Burbeck]

Dear all,
I have a problem when trying to perform local Universal kriging with krige().
The background:
For model definition I have to use a larger dataset since this is the only
one available with sufficient information. However, the interpolation
should only be performed for a smaller subset and as local UK (although I
use the global datset for model definition, I know that the data has local
differences).

MyExample:
myResult <- krige(log(myValue) ~ GEO_UNIT - 1,

locations = myGlobalDataset,

newdata = myPredictionGrid_1000_utm,

model = myVgmf,

nmax = Inf,

nmin = 0,

maxdist = 18000,

indicators = FALSE,

na.action = na.pass,

debug.level = -1)

I get the following error:
"solve.c", line 87: singular matrix in function Usolve()
"lufactor.c", line 207: singular matrix in function m_inverse()
Fehler in predict.gstat(g, newdata = newdata, block = block, nsim = nsim,  :
  m_inverse

When "maxdist = Inf" everything works well.

Any hint where the problem could be would be welcome.

Sven