问题标签 [gwr]

我是 R 新手，我正在尝试使用 Roger Bivand 为 R 开发的一些代码来计算地理加权回归。有关代码的讨论，请参阅 help(gwr)。在我的教授的帮助下，我能够让它在自适应带宽下工作：

图书馆(GWmodel)

DM<-gw.dist(dp.locat=coords)

Adaptive.bs1 <- gwr.sel（OnTer_tory ~ CL + cinc + Allie_thUS + Island + Landlocked，adapt = TRUE，method = "cv"，gweight = gwr.bisquare，coords = coords）

gwr1 <- gwr（OnTer_tory ~ CL + cinc + Allie_thUS + Island + Landlocked，适应 = Adaptive.bs1，gweight = gwr.bisquare，coords = coords，hatmatrix = TRUE）

但是，当我尝试使用固定带宽时：fixed.bs1 <- gwr.sel(OnTer_tory ~ CL + cinc + Allie_thUS + Island + Landlocked, adapt = FALSE, method = "cv", gweight = gwr.bisquare, coords = coords )

gwr2 <- gwr（OnTer_tory ~ CL + cinc + Allie_thUS + Island + Landlocked，adapt = fixed.bs1，gweight = gwr.bisquare，coords = coords，hatmatrix = FALSE）

我收到以下错误：gwr2 <- gwr(OnTer_tory ~ CL + cinc + Allie_thUS + Island + Landlocked, + adapt = fixed.bs1, gweight = gwr.bisquare, coords = coords, hatmatrix = FALSE) 错误：(adapt <= 1) 不正确

我究竟做错了什么？

目标：我想在坐标对之间创建一个相异矩阵。我想使用这个矩阵作为输入来使用 Moran's I (LISA) 和后者在地理加权回归 (GWR) 中计算局部空间集群。

问题：我知道我可以dnearneigh{spdep}用来计算距离矩阵。但是，我想使用我已经估计的多边形之间的旅行时间。在实践中，我认为这就像输入一个相异矩阵，它根据另一个特征告诉多边形之间的距离/差异。我尝试将矩阵输入到dnearneigh{spdep}，但出现错误Error: ncol(x) == 2 is not TRUE

有什么建议么？下面有一个可重现的示例：

编辑：再深入一点，我想我可以使用mat2listw{spdep}，但我仍然不确定它是否保持矩阵和多边形之间的对应关系。如果我添加row.names = T它会返回一个错误row.names wrong length:(

可重现的例子

计算简单的距离矩阵

现在如何使用我的diss_matrix_invers这里？

我在R使用该spgwr库时运行了地理加权回归 (GWR)，现在我想返回准全局 R2（模型的拟合）。我已经深入研究了结果，summary(gwr_model)但我还没有找到提取这个值的方法。任何想法？

可重现的例子

我正在R使用该spgwr库进行地理加权回归 (GWR)。我知道可以使用 检索每个观察值的局部系数和标准误差gwr_fit$SDF。

现在我如何使用这些信息来确定哪些局部系数具有统计显着性，以便我可以将它们绘制在地图上？

可重现的例子

我在 R 中使用 spgwr 包运行了地理加权回归 (GWR)，并希望提取 AICc 的值并将其单独保存。

然而，这似乎是不可能的。那么，还有哪些其他方法可以提取或计算 GWR 的 AICc？

或者，如果这不可能，我们如何比较 GWR 模型和 OLS 模型？anova() 不是很有帮助，因为它不报告 p 值。

谢谢！

我正在使用 R 3.3.0，pachage 'spgwr'

我创建了一个空间点数据框，如图所示：

原始数据框 (points_df) 具有十进制度的 latlong 坐标。

我没有定义任何投影信息，现在 points_map 的 projargs 参数是 NA。我将此空间点数据框用于 GWR 带宽计算：

它已经运行了好几天，到目前为止的结果显示了两个带宽估计：

因为我没有直接输入long-lat坐标，而是没有投影的Spatial Points Data Frame坐标，所以我的问题是带宽的单位是多少？

函数信息表明，如果您使用longlat=TRUE，则​​单位为Km，但如果您有 SpatialPoints 对象，则 longlat=FALSE并且该值取自对象本身。我将衷心感谢您的帮助。谢谢

我找到了这个用于地理加权回归 (GWR) 的包，但我不知道如何安装它，以便在使用 Jupiter-notebook 时可以访问它？

我有一个大型栅格数据集，已将其转换为 SpatialPointsDataFrame 并从那里转换为常规数据框。我需要使用 spgwr 包运行地理加权回归 (GWR)。我已经使用 gwr.sel() 函数成功计算了 GWR 带宽，但是当我尝试使用 gwr() 运行 GWR 时，我收到一条错误消息，指出 R 无法分配大小为 25 Gb 的向量。我在想有可能只对整个数据框的随机子样本运行局部 GWR 回归，同时将所有其他点保留为回归的邻域值。

这是我到目前为止所拥有的，但我不确定如何在子样本上运行 GWR 而不会摆脱其他所有点。

我正在学习使用 R 通过使用 gwr.sel() 和 gwr() 函数来编写简单的地理加权回归。当我阅读 gwr() 的官方文档时，链接： gwr 定义

我发现有一个名为“权重”的属性，它的定义如下：

它包含在原始函数定义签名中：

我在这个文档或谷歌结果上都没有找到这个属性的直接使用示例，专家可以教我 GWR 过程中使用的参数是什么吗？提前致谢！

我正在尝试指定用于运行地理加权回归的集群数量。我有一堆已转换为数据框的 9 个栅格。我是否收到以下错误（最后）并且不知道如何解释错误，因为我能说的一切看起来都是正确的。示例 df.notmissing：

我想使用 spgwr 库运行 gwr。没有 clusterR 的 gwr 如下所示：

我创建了以下代码以与 clusterR 一起使用：

我在使用 clusterR 时遇到的错误：