2

我是分析的新手，所以我不知道我所做的是否正确。我有很多时空数据，我想绘制一个时空样本变异函数。

我的问题是我不太了解variogramST(package: gstat) 是如何工作的。是spacelag公里还是什么？这就是我现在所做的：

var_st <- variogramST (datast~1, data=DF, tunits="days", tlags=seq(0, 93, by= 31), 
                       boundaries=seq(0, 100, by=25), na.omit=T, assumeRegular=F)

我的第二个问题是：

它是否在精确距离 25、50、75、100 处绘制变异函数？例如，如果两个城市有距离例如 12 正在variogramST考虑吗？或者它只是考虑精确距离为 25、50、75、100 的观察对？此外，我有每月的数据。我读到，在这种情况下，最好的选择是选择tunits="days"，我选择tlags=seq(0, 93, by= 31)，对吗？我不知道如何处理它。

非常感谢所有愿意回答我的人！

4

1 回答 1

0

  1. 空间距离单位取决于输入数据的坐标参考系，详情请参阅 sp::spDistsN1。

  2. 经验变异函数依赖于分箱 - 距离类。空间和时间上的观测按它们的近似空间距离分组（例如，0 公里到 25 公里；25 公里到 50 公里，……根据cutoff和width）boundaries。随着时间的推移，它还取决于您提供给 variogramST 的对象的类别：

    • 对于 STFDF，它们tlags被评估为 STFDF 时隙的 id，在您的情况下，观察同一时间 id，有 31、62 和 93 个 id 差异 - 很可能不是您想要的，假设时隙中的每个 id ofDF指新月份，需要使用 0:3。这个论点tunit在这里没有作用。
    • 如果数据对象是 STIDF，则时间距离也被分箱到分隔范围，其中tlags用作具有时间单位的类之间的间隔tunit。

有关将时空数据表示为 STFDF、STSDF 和 STIDF 的详细信息，请参阅 package spacetime。

高温高压

于 2020-03-26T15:17:33.837 回答