Dear All,
I work on four-year spatial yield monitor data from the same field. I use
GSTAT library of  R to fit variogram models (spherical, Gaussian,
exponential, and linear) for each year and I identified the exponential
function to be optimal in explaining the spatial yield pattern in each year.

I intend to carry out the combine analysis of the data over years to assess
the stability of spatial pattern. I have the attention of using the kriging
approach to obtain predicted yield value from each year and possibly
looking at their correlations  to assess the stability of yield pattern
across the yield. My worry is that the predicted value may not likely be
from the same locations over the years for correlation analysis. How do I
achieve this?

Please could anyone guide me on how to use GSTAT library to achieve this?
Or is there any other approach of assessing the stability of spatial yield
pattern over years with the same field?

Thanks while looking forward to your response and advice.

	[[alternative HTML version deleted]]

