Agriculture Reference
In-Depth Information
Fig. 11.3 Comparison of the variance of the direct estimator with the MSE estimate of the
SEBLUP estimator
loglike AIC BIC
-28.11950 62.23900 62.83067
The corresponding estimates of the MSE of the SEBLUP estimates are as
follows.
> domSFH < - mseSFH(yobs ~ auxFH - 1, vardir ¼ var, proxmat ¼ cmatr,
+
data¼datFH)
> domSFH$mse
[1] 20.0397833 15.3639483 23.4089002 18.7782084 0.6928073 16.3429502
15.0361337 23.3314683 18.2791517
We can compare the variance estimates of the direct estimator with the MSE
estimates of the SEBLUP (see Fig. 11.3 ) using the following code.
> plot(datFH$var,domSFH$mse,axes ¼ T,cex ¼ 0.5,pch ¼ 19,
+ xlab ¼ "Variance Direct",ylab ¼ "MSE Spatial EBLUP")
> abline(a ¼ 0, b ¼ 1)
Note that of the exercise described in this section is only used for illustrative
purposes. To evaluate the capacity of the spatial SAE, we must define a more complex
structure for the neighborhood system of the SAs. In fact, for a 3 3system
(as considered here), the presence of a strong edge effect may influence the results.
Search WWH ::




Custom Search