问题标签 [gwr]

For questions regarding programming in ECMAScript (JavaScript/JS) and its various dialects/implementations (excluding ActionScript). Note JavaScript is NOT the same as Java! Please include all relevant tags on your question; e.g., [node.js], [jquery], [json], [reactjs], [angular], [ember.js], [vue.js], [typescript], [svelte], etc.

0 投票
0 回答
290 浏览

r - 在 r 中的地理加权回归中获得新点的预测

我有一个观察点列表,我希望在这些点上估计地理加权回归。然后,我有一组单独且不同的点,我想使用这个估计的模型在这些点上进行预测。这是使用以下代码模拟的:

我的问题是 S.gwr 结构中唯一的拟合值是长度 100,并且是用于估计模型的 100 个点,而不是 pred.df 中我想要预测 y_hat 的 50 个点。谁能帮我指出没有为我提供 pred.df 中点的 50 个预测的问题?谢谢。

0 投票
1 回答
673 浏览

python - 面对 LinAlgError:矩阵是奇异的

社区,

语言障碍对我打击很大。不知何故,我设法用一个模型学习,但我无法预测,因为它现在说我的矩阵是奇异的——我不明白——它是唯一的另一个词吗?

无论如何,这里是重要的部分:

(457, 123)

(457, 2)

文件“...\venv3\lib\site-packages\spglm\iwls.py”,第 37 行,在 _compute_betas_gwr xtx_inv_xt = linalg.solve(xtx, xT)

文件“...\venv3\lib\site-packages\scipy\linalg\basic.py”,第 216 行,在求解 _solve_check(n, info)

文件“...\venv3\lib\site-packages\scipy\linalg\basic.py”,第 31 行,在 _solve_check raise LinAlgError('Matrix is single.')

numpy.linalg.LinAlgError:矩阵是奇异的。

.

在 X 和 coord 是数字(正数和少数负数,coord 是坐标经度和纬度)。

我试图使用的模型来自这个库:

在这里找到的文档

预测在这里

任何帮助将不胜感激!

最好的祝福

0 投票
1 回答
149 浏览

r - 地理加权逻辑回归不同维度

我正在尝试进行地理加权逻辑回归来量化我的数据集中的空间变化。但是,在运行 gwr 模型时,我收到输入数据和坐标具有不同维度的错误。

这是我用于荷兰边境的代码:

我的数据有关于病原体流行的二元结果 (0/1)。

接下来我进行了逻辑回归。这里没有问题。

我用平面 CRS 转换了空间 * 对象中的数据

但是一旦我运行这条线,我就会得到一个错误

gwr(pathogen ~ Age_Category, data = spt, bandwidth = bw, 中的错误:输入数据和坐标具有不同的维度

有人知道如何解决这个问题吗?先感谢您!

0 投票
0 回答
35 浏览

r - 地理加权套索 (pkg = gwrr) - Wi %*% X[-i, ] 中的错误:需要数字/复数矩阵/向量参数

我正在尝试在 packgae 中进行地理加权套索(GWL)gwrr。但是得到这个错误:

代码是:> col.gwl <- gwl.est(NPb~NAl2O3, locs, df, "exp")

当我检查错误时,在线答案是关于数据框中的因素/字符。还有矩阵as.matrix。但是我的数据框中的变量df都是numeric.

有人可以帮我解决这个问题吗?太感谢了!

0 投票
0 回答
410 浏览

r - R中的地理加权负二项式回归

我已经阅读了一些与地理加权回归 (GWR) 相关的参考资料,我发现一些 GWR 技术在 R 中可用,例如在 GWModel 的包中。那里有一些 GWR 类型,包括泊松族 (GWPR)。但是,我想知道 R 中是否有另一个包为 GWR 提供负二项式。一些论文将其称为 GWNBR。不幸的是,这种技术在 GWModel 包中不可用。因此,如果有人知道这一点,那么使用 GWNBR 完成一些任务可能会对我有很大帮助。

0 投票
1 回答
249 浏览

r - 如何测试 R 中的空间非平稳性以确定是否需要局部回归模型?

我有一个数据集,我为它实现了一个回归模型,并从中假设系数局部变化。如果给出了空间非平稳性,则运行局部回归模型是有意义的,在我的例子中是地理加权回归 (GWR)。

为了找出是否存在空间非平稳性,我知道 Koenker 检验可以在任何 GIS 软件中使用普通最小二乘 (OLS) 回归模型进行计算。但是对于这个项目,我正在使用 R 并且我需要找出是否给出了空间非平稳性。R中是否有任何统计测试可以给我这个结果?

我知道 R 中的 spgwr 包,但它只是实现 GWR,而不是测试是否需要 GWR。

0 投票
1 回答
36 浏览

r - 我正在使用 R 中的 ... 作为我正在制作的函数中的控制变量,如何将控制变量转换为列表?

我有这个函数可以进行地理加权回归,我有我的形状文件,x用于 x 变量和yy 变量,...用于控制变量。例如,我想先调用回归但有多个控制变量 - 我该怎么做?

我认为它需要使用列表,但我不确定如何执行此操作。

0 投票
1 回答
41 浏览

r - 如何允许 R 忽略我的函数的参数,而不跳到下一个参数?

我有这个功能

有 3 个 shapefile,我希望 R 允许 shape2 和 shape3 丢失。

虽然例如如果我使用 if(missing(shape2)) {} 然后输入:

R 将无法识别只有 2 个 shapefile 并且缺少第二个。

0 投票
1 回答
67 浏览

r - 如果我希望它是来自新对象的变量,如何指定 tm_fill()?

我正在尝试创建一个 R 函数,该函数将对用户从空间多边形数据框中指定的变量运行 GWR。运行该函数的最终结果是两个映射——一个自变量的值和一个来自 GWR 模型的系数值。第二张地图有问题。

我已经成功地为要可视化的系数创建了 GWR 模型和“结果”对象。

要创建 GWR 系数的可视化,我必须将 my 指定tm_fill()为“结果”对象中的一列,但我不知道如何做到这一点,以便该函数可以用于任何空间多边形数据框。到目前为止,我已经尝试使用该paste0()函数,如下所示:

但是我收到一个错误,说填充参数既不是颜色也不是有效的变量名。

我将不胜感激任何可以帮助我解决问题的提示。

0 投票
1 回答
66 浏览

r - 如何在 spgwr::ggwr() 中获取局部回归系数的标准误差?

我正在使用spgwr::ggwr()泊松模型和对数链接函数来拟合广义地理加权回归。结果提供了局部系数估计,但我错过了如何获得它们的标准误差(或 t 统计量)来计算伪 p 值。

下面是一个使用SpatialEpi::NYleukemia数据集的玩具示例:

有什么办法可以得到系数的标准误差?

谢谢!