r - 如何从线性模型在国家地图上进行克里金法

Question

我想在网格上使用普通克里金法获取模型预测。在我下面的代码中，GPS 是“SpatialPointsDataFrame”，而 EU 是“SpatialPolygonDataFrame”。用于估计模型的所有协变量都包含在“All_data”中。这些协变量仅涉及 GPS 中的点数据。

# Merge points with data 
GPS = merge(GPS, All_data, by.x = "nuts", by.y = "geo")
GPS <- GPS[!is.na(GPS@data$Value),]
proj4string(GPS) = CRS("+init=epsg:28992")

# Merge polygons with data 
EU = merge(EU, All_data, by.x = "nuts", by.y = "geo")
EU = subset(EU, EU@data$nuts %in% GPS@data$nuts)
proj4string(EU) = proj4string(GPS)

# Create a grid     
ext <- extent(EU)
r <- raster(ext, res=0.1)  
r <- rasterize(EU, r, field=1)
EU_grid = as(r, 'SpatialPixels')

# Estimate linear model
OLS_model = lm(Value ~ pop + unemp, data = GPS_coords) 

有谁知道现在如何进行以获得带有线性模型估计的普通欧洲 kridged 地图？