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