zoo相关内容

将每日气象周转换为标准气象周(R

我在使用xts、zoo或lubridate包将每日数据转换为每周数据时遇到了许多问题。没有一个答案适合我的问题。我已经尝试了以下代码 library(zoo) library(lubridate) library(xts) library(tidyverse) #Calculation for multistation set.seed(123) df ..
发布时间:2022-08-03 21:36:02 其他开发

在 data.table 中使用带有 .SD 的 na.locf 函数(zoo 包)

我正在尝试填写所有 NA,不包括第 1 列和第 4 列的前两个 NA,以及第 2 列和第 3 列的三个 NA,以及最近的非 NA 值.这是我的数据和代码: hh 但是,我正在寻找将 lapply 与 .SD 一起使用,因为每种类型都有两列以上.这可能吗? 解决方案 试试 setDT(hh)[, lapply(.SD, function(x) na.locf(x, na.rm=FA ..
发布时间:2022-01-13 19:38:08 其他开发

有没有一种快速的方法在 data.table 中运行滚动回归?

我在 R 中运行滚动回归,使用存储在 data.table 中的数据. 我有一个工作版本,但是它感觉就像一个黑客——我真的在使用我从 zoo 包中知道的东西,而 没有任何魔法code>data.table ...因此,感觉比应该的要慢. 结合 Joshua 的建议 - 下面 - 使用 lm.fit 而不是 lm 可以将速度提高约 12 倍. (修订)示例代码: 要求(动物园)需 ..
发布时间:2022-01-13 18:51:51 其他开发

按类别填写数据表中的缺失(空白) - 向后和向前

我正在处理我 11 年以上临床实践的大量账单记录数据集.相当多的行缺少转诊医生.但是,使用一些规则我可以​​很容易地填写它们,但不知道如何在 R 下的 data.table 中实现它.我知道 zoo 包中有 na.locf 之类的东西,data.table 包中的自滚动连接.我看到的例子太简单了,对我没有帮助. 这里有一些虚构的数据来指导你(作为 dput ASCII 文本表示) 结构( ..
发布时间:2022-01-13 18:51:33 其他开发

自适应移动平均线 - R 中的顶级性能

我正在寻找 R 中滚动/滑动窗口函数方面的一些性能提升.这是一项非常常见的任务,可用于任何有序的观察数据集.我想分享一些我的发现,也许有人可以提供反馈以使其更快. 重要的一点是,我关注的是 align="right" 和自适应滚动窗口的情况,所以 width 是一个向量(与我们的观察向量长度相同).如果我们将 width 作为标量,那么在 zoo 和 TTR 包中已经有了非常完善的功能,这些功能 ..
发布时间:2022-01-13 18:42:51 其他开发

R中的xts回归

是否有使用以下类型的 xts 对象运行回归的实用程序: lm(y ~ lab(x, 1) + lag(x, 2) + lag(x,3), data=as.data.frame(coredata(my_xts))) 其中 my_xts 是一个 xts 对象,其中包含一个 x 和一个 y.问题的关键是有没有办法避免做一堆滞后和合并以使 data.frame 具有所有滞后?我认为包 dyn 适用于 ..
发布时间:2022-01-11 10:10:52 其他开发

合并两个不同时间粒度的时间序列

我有两个包含时间序列数据的 CSV 文件.我想将两者合并到一个文件中.File1 具有 1 分钟间隔的周期性数据.File2 具有非周期性的事件触发数据.File2 中数据的时间戳可能与 File1 中的数据一致,也可能不一致.我想合并这两个数据集以创建一个数据集,其时间戳是 Data1 和 Data2 的并集.对于两者都不通用的时间戳,我希望将相应数据集的缺失条目指示为 NA. 这是 Fi ..
发布时间:2022-01-11 09:53:21 其他开发

将每小时的动物园时间序列输入函数 stl()

在你问之前,是的,我需要展示这么多数据.stl() 需要两个周期的数据.在这种情况下,一个周期是 24 个值,因此 stl() 至少需要 48 个值. 另外,来自 stl() 帮助: "....这应该是频率大于一的 "ts" 类的对象...." 我正在升级一些旧的计算,以便我的数据采用动物园格式.到目前为止,我已经升级了每月和每日数据,没有任何明显的问题,尽管存在 stl() 速 ..
发布时间:2022-01-11 09:52:38 其他开发

na.locf 但不做尾随 NA

我有以下时间序列 >y 直接 na.locf 给我这个: >na.locf(y)[,1]2011-09-04 不适用2011-09-05 不适用2011-09-06 32011-09-07 42011-09-08 42011-09-09 62011-09-10 72011-09-11 82011-09-12 82011-09-13 8 我如何做到这一点? [,1]2011-09-04 ..
发布时间:2022-01-11 09:50:11 其他开发

趋势长度 - 面板数据

我有一个平衡良好的面板数据集,其中包含 NA 观察结果.我将使用 LOCF,并想知道每个面板中有多少连续的 NA,然后再进行观察.LOCF 是一个过程,其中可以使用“最后一次观察结转"“填充"缺失值.这对一些时间序列应用程序来说是有意义的;也许我们有 5 分钟增量的天气数据:对缺失观测值的良好猜测可能是 5 分钟前进行的观测. 显然,在一个面板中将观察结果提前一小时比在同一面板中将相同的观察 ..
发布时间:2022-01-11 09:48:58 其他开发

有效地从 1 个数据帧中的多个时间序列的开头和结尾删除缺失值

使用 R,我正在尝试从包含多个时间序列的数据框的开头和结尾修剪 NA 值.我已经使用 for 循环和 zoo 包实现了我的目标,但正如预期的那样,它在大型数据帧上效率极低. 我的数据框看起来像这样,包含 3 列,每个时间序列由它的唯一 ID 标识.在这种情况下,AAA、B 和 CCC. id 日期值AAA 2010/01/01 不适用AAA 2010/02/01 34AAA 2010/03 ..
发布时间:2022-01-11 09:44:18 其他开发

R/zoo:‘order.by’中的索引条目不是唯一的

我有一个 .csv 文件,其中包含 4 列数据,对应一列日期/时间,间隔为一分钟.缺少一些时间戳,所以我试图生成缺少的日期/时间并在 Y 列中为它们分配 NA 值.我以前使用其他格式完全相同的 .csv 文件完成了此操作,没有任何问题.代码是: #读取csv文件har10 = read.csv(fpath, header=TRUE);# 设置日期har10$HAR.TS 警告信息:在动物园(ha ..
发布时间:2022-01-11 09:44:12 其他开发

R:在时间序列中填充缺失的日期?

我有一个缺少天数的动物园时间序列.为了填补它并有一个连续的系列,我这样做...... 我从开始到结束生成一个 chron 日期时间序列. 我将我的系列与这个合并. 我使用 na.locf 将 NA 替换为 las 观察. 我删除了 syntetic chron 序列. 我可以更轻松地做同样的事情吗?也许有一些与频率相关的索引函数? 解决方案 如果使用带有索引的 ..
发布时间:2022-01-11 09:28:03 其他开发