Dear Subscribers,
I have some trouble in specifying the parameters trend.d and trend.l
correctly.
As global trend I assume row and column effects which translate into the
formula:
trend=~factor(Row)+factor(Column)
Consequently, I obtain the parameters of a the model via variofit() and want
to do OK-estimation
for each point of a grid. Therefore, I call krige.cov() n times, for my n
data points. Each time leaving
out the point which I want to estimate. Actually, I want to smooth the
surface using information
of neighboring locations.
parms <- krige.control( type.krige="ok", obj.model=my.model, aniso.pars=c(
ai.axis,ai.ratio),
trend.d =
.?, trend.l =
.? )
I always get the error: ..trend.spatial(trend = krige$trend.d, geodata =
geodata) :
trend elements not found.
The geodata object I pass to krige.conv() contains the covariates Row and
Column.
How can I specify the parameters trend.d and trend.l correctly in order
to do the OK-estimation
location by location?
I would appreciate to get a hint how to do that.
Thanks,
Andre
[[alternative HTML version deleted]]
_______________________________________________
R-sig-Geo mailing list
[email protected]
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
