[Gstat-info] krige with KED

[G. Allegri]

Hi.
I need to run SGS on a variable with a linear external drift:

krige(EceZ~dist_costa,t,newdata=grdPIX,model=Vgm,nsim=150,nmax=50,maxdist=600)

where both EceZ and dist_costa are inside t (a SpatialPointsDataFrame).
When I launch the above function I receive an error saying that the object
"dist_costa" isn't found. I've tried with t$EceZ~t$dist_costa, and in this
case the message is "NROW(locs) != NROW(X)".
What am I doing wrong?

Giovanni
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.geo.uu.nl/pipermail/gstat-info/attachments/20080314/e16382f7/attachment.html