贝叶斯地理统计模型R-INLA-2

Validation

上一期我们介绍了，如何利用100空间位置，来估计研究区域内的降雨量。
但是并没有做模型的validation
虽然我们已经将数据集分成test与train两个部分；接下来我们将介绍如何测试模型的好坏及与glm模型对比。

validation

首先我们将367个点绘制出来，看一下效果。

image.png

然后根据前述的SPDE函数，将367个空间效应给提取出来 然后整合放在stack.train里面，提示这里的y为NA但是X变量还是来源于train数据里面

# plot train
dim(train)
# train plot
ggplot() + 
  geom_point(data=train, aes(y = lat, x = lon,size=rainfall))