我正在尝试运行 GWR 分析(spgwr 包)来预测哪些生态变量(即森林破碎化、NPP、NDVI 等)可以预测精确变量(生态系统服务供应单位)的地理位置。我做了一个非常简单的代码示例来测试模型。我可以毫无错误地计算 GWRbandwidth,但是当我使用 gwr 函数运行模型时,我收到以下错误消息:
“新数据矩阵行不匹配”
我该如何解决这个错误?我无法弄清楚问题是出在数据中还是出在代码中。谁能帮我解决它?
我在这里复制代码示例,以便您了解更多详细信息。
#library(spgwr)
#library(stars)
#set working directory
setwd("C:/Users/migno/OneDrive/Dottorato/++ELABORAZIONI_NUOVE/INSIEME_Muga_Fluvia/PPGIS_SRC")
raster1 <- readRDS("rasterobj.RDS")
coordinates <- as.matrix(st_coordinates(st_centroid(raster1))[,1:2])
GWRbandwidth <- gwr.sel(SPUH_MF.tif~forest_fragmentation.tiff, data = raster1, adapt = T, coords = coordinates)
GWRbandwidth
gwr.model = gwr(SPUH_MF.tif~forest_fragmentation.tiff,
data = raster1,
coords = coordinates,
adapt = GWRbandwidth,
hatmatrix = TRUE,
se.fit = TRUE)```
I obtain this errror: Error in gwr(SPUH_MF.tif ~ forest_fragmentation.tiff, data = raster1, : new data matrix rows mismatch
Thanks in advance for your help!
Enrica