We can clearly see that there are some spatial trend that is not well captured by the second model.Īsreml-R fits a model by estimating the variance components using REML approach with the average information (AI) algorithm (Johnson and Thompson, 1995). Theme_bw(base_size=16) + xlab('Unit Lag') + ylab('Empiricial Variogram for Stud. m1 % ggplot(aes(x, gamma)) + geom_point() + We want to the following model to the yield \(\boldsymbol\).
