[Gstat-info] krige with KED

[Jon Olav Skoien]

Hi,

Do you have dist_costa as a column in grdPIX as well? For universal 
kriging you need to know the value of your secondary variable at all 
prediction locations in addition to your observation locations.

Cheers,
Jon

G. Allegri wrote:
> 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
> ------------------------------------------------------------------------
>
> _______________________________________________
> Gstat-info mailing list
> Gstat-info at geo.uu.nl
> http://mailman.geo.uu.nl/mailman/listinfo/gstat-info
>