 

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]]


