我有一个大型栅格数据集，已将其转换为 SpatialPointsDataFrame 并从那里转换为常规数据框。我需要使用 spgwr 包运行地理加权回归 (GWR)。我已经使用 gwr.sel() 函数成功计算了 GWR 带宽，但是当我尝试使用 gwr() 运行 GWR 时，我收到一条错误消息，指出 R 无法分配大小为 25 Gb 的向量。我在想有可能只对整个数据框的随机子样本运行局部 GWR 回归，同时将所有其他点保留为回归的邻域值。
这是我到目前为止所拥有的，但我不确定如何在子样本上运行 GWR 而不会摆脱其他所有点。
##Calculate kernel bandwidth
gwr.bw <- gwr.sel(merge.df$nsEVI_YAproj ~ merge.df$aDEM_YAproj, data = merge.df, coords = cbind(merge.df$x, merge.df$y), adapt = TRUE)
##Run GWR model
gwr <- gwr(merge$nsEVI_YAproj ~ merge$aDEM_YAproj, data = merge.df, coords = cbind(merge.df$x, merge.df$y), adapt = gwr.bw, hatmatrix = TRUE, se.fit = TRUE)