问题标签 [time-series]

嘿，我正在尝试在 statistica 中计算余弦分析，但不知道如何计算。我需要计算昼夜节律数据的 MESOR、AMPLITUDE 和 ACROPHASE。

http://www.wepapers.com/Papers/73565/Cosinor_analysis_of_accident_risk_using__SPSS%27s_regression_procedures.ppt

有一个链接显示了如何做到这一点，公式等，但它并没有给我太多帮助。有谁知道它的代码，无论是在 statistica 还是 SPSS 中？

我真的需要完成这项工作，因为这是一篇重要的论文

对于金融时间序列，如每日股票价格或日内数据，首选哪些时间序列包？xts，普通动物园，或 timeSeries 或其他什么？我同时使用 xts 和 zoo，但有时不确定只使用 xts，或者有时 zoo 具有更轻开销的优势；另外，我记得 Rmetrics 对所有这些软件包的评论论文，其中声称 xts 甚至无法完成他们所做的一些测试。但是我现在找不到那张纸了。

我刚刚开始使用 R - 并且在一些时间序列数据上遇到了一些僵局。

我有一个“动物园”格式的时间序列（日期和值），我想在事件发生时用叉号进行注释。这些事件是不规则的，采用 csv 格式（只是日期，有时会重复）。

我已经设法将日期等读入 R 接受的格式 - 但我似乎无法获得一种方法来绘制主要时间序列，并在顶部注释次要事件？

更新： 对不起，我之前错过了 - 以下是我正在使用的数据类型：

我现在想比较推文的数量（对于搜索词），但我只有以下形式的参差不齐的数据：

有些是重复的。到目前为止，我已经设法编写了一个 python 脚本来进行一些清理（即日期元组，出现次数），但我希望我可以使用 R 处理原始数据？

非常感谢

我的数据如下所示：

每个组名都是一个数据系列。由于我想分别绘制每个系列，因此我将它们分开如下：

预计到达时间：

生成的动物园对象如下所示：

因此，如果我想根据时间绘制 dist.mean 并为每个系列包含等于 +/- dist.sd 的误差线：

• 如何结合 A、B、C dist.mean 和 dist.sd？
• 如何制作条形图，或者更好的是结果对象的折线图？

我有一个矩阵，每列代表一个随时间变化的特征。我需要找到给定窗口大小的这些值的移动平均值。

有没有类似 MATLAB 中的函数？

我想编写一个代码，涉及使用 ETS 和自回归模型自动选择最佳复合模型。我应该根据什么标准进行选择？

此外，如果我使用 auto.arima 函数从 R 中的预测包中推断 AR 项的数量和相应的系数，我的输入序列是否必须是固定的？还是会自动选择 d 的值，从而返回非平稳模型？

谢谢， 帕尼

我最近一直在处理大型数据集（超过 40 万行）。到目前为止，我一直在使用 XTS 格式，它适用于包含几万个元素的“小型”数据集。

现在项目增长了，在检索数据库数据并将其放入 XTS 时，R 只会崩溃。

据我了解，R 应该能够拥有大小不超过 2^32-1 个元素（或根据版本为 2^64-1）的向量。因此，我得出的结论是 XTS 可能有一些限制，但我在文档中找不到答案。（也许我对我对理论上可能的向量大小的理解有点过分自信）。

总而言之，我想知道是否：

1. XTS 确实有大小限制
2. 您认为处理大型时间序列的最聪明方法是什么？（我正在考虑将分析分成几个较小的数据集）。
3. 我没有收到错误消息，R 只是自动关闭。这是一种已知的行为吗？

解决方案

1. 与 R 相同，它取决于所使用的内存类型（64 位、32 位）。无论如何，它非常大。
2. 分块数据确实是一个好主意，但它不是必需的。
3. 此问题来自R 2.11.0 中的错误，该错误已在 R 2.11.1 中解决。长日期向量（这里是 XTS 的索引）存在问题。

我是 stackoverflow 的新手，对 R 也很陌生，但我搜索了很长时间，但找不到以下问题的答案。

我有许多数据文件是温度与时间序列的关系。我将 CSV 作为 ZOO 对象导入，然后转换为 XTS。一个正确的文件看起来像这样，一个小时和一个半小时的读数：

但有些时间值稍有偏差 - 即 23:59:00 而不是 00:00:00，或 00:29:00 而不是 00:30:00。

我想纠正这些时间序列，因为微小的差异对我的分析并不重要，我最终想要合并文件，所以每个时间序列都需要有相同的时间。

我想要一个命令，可以只说“将时间序列向前移动 1 分钟，但不要更改数据列（例如 S_21）。我gsub()在更容易的更改方面有一些运气，并考虑使用复杂的正则表达式来更改数据在转换为 ZOO 或 XTS 之前。我已经阅读过lag()，diff()但它们似乎相对于时间序列移动了数据值；如果我错了，请纠正我。

任何解决此问题的帮助将不胜感激。

我有以下问题：我有一个包含超过 10000 个条目的时间序列，我想对每个条目执行一些计算。仅此一项不会有问题，但我需要获得最后一个计算值才能获得下一个值。我需要的一个非常简单的形式如下所示：

Val(n) = Val(n-1) + (time-series-entry / 2) （或类似的东西！）

我不知道如何管理这个。简单地做这样的事情：

不会工作，因为不能（至少我不知道怎么做！）得到最后的计算值。然后我想：好吧，让我们使用 Loop-Recur。这将为我提供与时间序列条目相对应的值，但对于下一个我将不得不再次进行所有计算。迭代是正确的，但它不起作用，因为该函数有副作用。

所以我被困在这里。如果有人能给我一个提示，那就太好了。

我正在寻找能够将大量时间序列数据映射到强度图的工具或 API 的建议/示例。

数据包括国家、系列和年份的维度。这是一个示例http://spreadsheets.google.com/ccc?key=t9ZwziZAgy768ZTXDEg8Maw&authkey=CPn0pdoH&hl=en_GB&ui=1