问题标签 [biwavelet]

这个问题相当模糊，但有没有人在 R 中使用过 biwavelet 包并成功了？我有以下代码：

它返回错误：

image.default(x$t, yvals, t(zvals), zlim = zlims, ylim = rev(range(yvals)), 中的错误：'x' 和 'y' 值必须是有限且非缺失的

我不理解这个错误，因为我的数据不包含任何缺失值，而且它们都是有限的。

尝试重现示例时：

它很烦人，因此我不知道为什么我的数据在这个包中失败。上面示例中的数据与我的数据属于同一类，即“矩阵”。任何帮助或建议将不胜感激。

更改：

谁能建议为什么以下示例代码不起作用：

它会生成一个错误说明：

我将如何解决这个问题？

额外的：

作为对 Gavin Simpsons 的回答：如果我将数据保留为仅包含一个频率但更改时间向量，则代码可以正常工作。

在以下示例中：

如果我正在查看给定时间序列的昼夜周期，可以使用的最低分辨率是多少？此示例当前使用以 1 小时为间隔测量的数据，因此我一天有 24 个数据点。如果我使用 2,3 甚至 4 小时的时间步长，这是否同样重要？

还有一种方法可以计算在降低数据分辨率时结果的重要性是否会发生变化

我想生成一个结合了基本图形和 ggplot 图形的图形。以下代码使用 R 的基本绘图功能显示了我的图形：

产生

这些面板中的大多数看起来足以让我包含在我的报告中。但是，显示自相关的图需要改进。使用 ggplot 看起来好多了：

但是，由于 ggplot 不是基本图形，我们不能将 ggplot 与 layout 或 par(mfrow) 结合使用。如何用 ggplot 生成的自相关图替换从基础图形生成的自相关图？我知道如果我的所有数字都是用 ggplot 制作的，我可以使用 grid.arrange 但是如果在 ggplot 中只生成一个图，我该怎么做？

我有一个包含多个站点 30 年数据的每日时间序列。我正在使用biwaveletR 中的包来使用以下代码测试数据的周期性

运行代码时出现此错误

我使用的一些数据集每年都有一些缺失的每日价值。我是相当新的 R 和编程。使用biwavelet包时，似乎缺少值会导致这样的错误？我的问题有没有办法在 R 中解决这个问题？

我似乎找不到此页面上使用的 R 中的 wt() 函数：在 R 图形窗口中组合基础和 ggplot 图形

我在 ??wt 和谷歌搜索中找不到它。它必须在某个包中。上面的页面上没有图书馆声明。谢谢你的帮助。

我试图理解 R 中的交叉小波函数，但无法弄清楚如何使用双小波包将相位滞后箭头转换为时间滞后。例如：

这是我的两个时间序列。两者相同，但一个落后于另一个。箭头表示 45 度的相位角——显然向下或向上表示 90 度（同相或异相），所以我的解释是我看到的是 45 度的滞后。

我现在如何将其转换为时间滞后，即如何计算这些信号之间的时间滞后？

我在网上读到这只能针对特定波长进行（我认为这意味着一段时间？）。那么，假设我们对 365 周期感兴趣，并且信号之间的时间步长是一天，那么如何计算时间滞后呢？

image我对包中的 R 函数有疑问biwavelet。此函数用于绘制函数中的相关z值image。我正在绘制完全相关的数据集（R-squared=1 作为我的 z 值）。绘制图像时，彩色框会有一些变化。我对此感到困惑，因为我的z矩阵中的所有值都是 1。任何想法为什么会发生这种情况，以及如何解决它？

zvals是一个全为 1 的矩阵

更新：根据以下要求

以上是实际数据。下面是一个示例，但按预期绘制，而我的数据没有：

我正在使用 R 中的 biwavelet 包进行小波相干性分析。我有兴趣在每个时间步获得两个时间序列之间的相关性（rsq）。我正在使用

命令，它返回一个矩阵，其列数等于我的时间序列中的时间步数，但有 35 行。我不确定每一行表示什么，或者如何解释这个矩阵。任何帮助，将不胜感激。

编辑：来自 dput 时间序列 1（前 10 行）的数据 - 我在运行 wtc 函数时记录转换

时间序列 2（前 10 行）

该子集中的 rsq 矩阵

我正在使用双小波包的绘图功能，并且不断收到此消息

轮廓（）：圆形/长seglist - 设置选项（“max.contour.segments”）> 25000？

我在 R 中尝试了以下代码，但它不起作用

如何增加 R 中的 max.contour.segments。