On 01/09/2013 07:26 PM, halim10-fes wrote:
> Dear Ashton, Andy and the group,
>
> Thank you very much for your suggestions. Yes, specifying ‘data=data’ is
> throwing the error and
>
> prk.rich<-krige(Richness~Elev, data, newdata=grd, model=fit.sph1)
>
> is not giving that. But regarding your
Your prediction grid must have elevation data, that you predicted by OK.
Newdata file (prediction grid file) may be looks like this:
x
y
Elev
Data file: x y Richness Elev
Zia
Sent from my iPad
On Jan 9, 2013, at 2:33 PM, "halim10-fes" wrote:
> Dear Ashton, Andy and the group,
>
> Thank y
Dear Ashton, Andy and the group,
Thank you very much for your suggestions. Yes, specifying data=data is
throwing the error and
prk.rich<-krige(Richness~Elev, data, newdata=grd, model=fit.sph1)
is not giving that. But regarding your second suggestion to have elevation
covariate at all locati
Dear Abdul,
there are a couple of strange things about your code (the cross-variogram
section actually performs cross-validation), but there are two problems with
your code for the regression kriging in particular:
1. in my version of R at least, specifying 'data=data' throws the error you
desc
Dear Group,
Happy New Year to Everybody.
I am a novice user of gstat and geostatistics. I have a data set (please see
attached: data.csv) containing coordinates (East_x, North_y),Elev(Elevation),
and Richness (Plant Species) collected from 125 plots of 10m x 10m plots. I
want to interpolate pl