功能“rollapply”的更快替代方案

2023-11-25

我需要对包含大约 7,000 行和 11,000 列的 xts 数据运行滚动窗口函数。我做了以下事情:

require(PerformanceAnalytics)
ssd60<-rollapply(wddxts,width=60,FUN=function(x) SemiDeviation(x),by.column=TRUE)

我等了12个小时,但计算还没有完成。但是,当我尝试使用小数据集时,如下所示:

sample<-wddxts[,1:5]
ssd60<-rollapply(sample,width=60,FUN=function(x) SemiDeviation(x),by.column=TRUE)

计算在 60 秒内完成。我在配备 Intel i5-2450M CPU、Windows 7 操作系统和 12 GB RAM 的计算机上运行它们。

如果有更快的方法在大型 xts 数据集上执行上述计算,有人可以建议我吗?


如果可以的话,将它们转换为动物园对象。rollapply.zoo比更有效rollapply.xts(在这种情况下。我不确定一般情况下哪个更有效):

R> require(PerformanceAnalytics)
R> set.seed(21)
R> x <- .xts(rnorm(7000,0,0.01), 1:7000)
R> system.time({
+   r <- rollapply(x, 60, SemiDeviation, by.column=TRUE, fill=NA)
+ })
   user  system elapsed 
  9.936   0.111  10.075 
R> system.time({
+   z <- rollapplyr(as.zoo(x), 60, SemiDeviation, by.column=TRUE, fill=NA)
+ })
   user  system elapsed 
  1.950   0.010   1.964 
本文内容由网友自发贡献,版权归原作者所有,本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容,请联系:hwhale#tublm.com(使用前将#替换为@)

功能“rollapply”的更快替代方案 的相关文章

  • picker输入字体或背景颜色

    我在闪亮的仪表板中使用 pickerInput 这很好 除了一个问题 背景颜色和字体颜色太相似 使得过滤器选择难以阅读 有什么办法可以改变背景或字体颜色吗 如果可能的话 我想继续使用 pickerInput 但如果有一个带有 selectI
  • 绘制 Cox 回归的 Kaplan-Meier 图

    我使用 R 中的以下代码设置了一个 Cox 比例风险模型来预测死亡率 添加协变量 A B 和 C 只是为了避免混淆 即年龄 性别 种族 但我们真正对预测变量 X 感兴趣 X 是一个连续变量 cox model lt coxph Surv t
  • 如何在 R 中的 for 循环内将值存储在向量中

    我正在开始使用 R 但我对以下问题感到非常沮丧 我试图将 for 循环内完成的某些计算的值存储到我之前定义的向量中 问题是如何进行索引 因为for循环迭代代码的次数取决于用户的输入 所以变量i不一定要从1开始 它可以从80开始 for举个例
  • 更新 R6 对象实例中的方法定义

    如何更新 R6 类实例的方法定义 正如我所期望的 S3 使用当前的方法定义 对于 R5 参考类 我可以使用 myInstance myInstance copy 在 R6 中 我尝试了 myInstance myInstance clone
  • sapply - 保留列名称

    我试图总结数据集中许多不同列 变量 的平均值 标准差等 我已经编写了自己的汇总函数 以准确返回我需要和正在使用的内容sapply立即将此函数应用于所有变量 它工作正常 但是返回的数据帧没有列名 我似乎甚至无法使用列号引用重命名它们 也就是说
  • 在 R 的 for 循环中创建动态命名对象并分配动态值

    我正在尝试创建一套动态命名的新对象 例如 temp2015 使用 for 循环 并存储动态值 具体来说 其他对象的名称 例如 Y2015 和 for 循环中使用的值 例如 2015 在动态命名的新对象中 我不确定为什么下面的代码不起作用 Y
  • R - 重塑 - 熔化错误

    我正在尝试融化数据框 但出现了这个奇怪的错误 有什么想法吗 str zx7 data frame 519 obs of 5 variables calday new Date format 2011 01 03 2011 01 04 201
  • 在r中的某个阈值处破坏 cumsum() 函数

    例如我有以下代码 cumsum 1 100 我想打破它 如果一个元素 i 1 大于3000 我怎样才能做到这一点 因此 而不是这个结果 1 1 3 6 10 15 21 28 36 45 55 66 78 91 105 120 136 15
  • 在 RMarkdown 输出到 PDF 时缩进而不添加项目符号点或编号

    之前有人问过如何在没有项目符号的情况下缩进文本 RMarkdown 中的点 但这是针对 HTML 输出的 在 RMarkdown 中缩进而不添加项目符号点或数字 https stackoverflow com questions 47087
  • 使用选定因子水平的值向 ggplot-barchart 添加水平线

    在这个情节中 df lt data frame factor as factor c rep A 3 rep B 3 Treatment c rep c A B C 2 values runif 6 0 1 ggplot df aes Tr
  • dplyr:连接中的 NSE (by)

    我很难弄清楚如何使用 dplyr left join 和 NSE 连接两个表 问题是我无法为 by 提供正确的值 我想我现在已经找到了解决方案 但感觉我正在以一种额外复杂的方式来做 因此 如果您知道更简单 更优雅的解决方案 请告诉我 这就是
  • R:按组,测试一个变量的每个值是否存在于另一个变量中

    我有一个数据框架 结构如下 a lt c 1 1 1 2 2 2 3 3 3 3 4 4 b lt c 1 2 3 1 2 3 1 2 3 4 1 2 c lt c NA NA 2 NA 1 1 NA NA 1 1 NA NA df lt
  • 使用点阵个性化 R 上显示的 X 轴值

    我收集了大量包含日期 客户端及其 NFS 使用情况的数据 我正在使用lattice R包进行绘图 正如对超级用户的建议 https superuser com questions 523195 plot custom log data on
  • 如何将旋转的 NetCDF 转换回正常的纬度/经度网格?

    我有一个带有旋转坐标的 NetCDF 文件 我需要将其转换为正常的纬度 经度坐标 经度为 180到180 纬度为 90到90 library ncdf4 nc open dat nf 对于尺寸 它显示 1 5 variables exclu
  • StatET调试工具

    我想我只是很密集 但我似乎无法弄清楚如何在 Eclipse 中的 R 中使用调试工具 StatET 插件 有人有关于这个主题的任何提示或教程吗 StatET 2 00 现在对高级 可视化调试提供实验性支持 需要 Eclipse 3 6 或
  • 获取包含矩阵行内最大值的列名称,该矩阵在数组内包含单独的最大值

    例如给出 dim1 lt c P PO C T dim2 lt c LL RR R Y dim3 lt c Jerry1 Jerry2 Jerry3 Q lt array 1 48 c 4 4 3 dimnames list dim1 di
  • 条件和分组 mutate dplyr

    假设我有以下每个抽屉库存增加的数据 gt socks year drawer nbr sock total 1990 1 2 1991 1 2 1990 2 3 1991 2 4 1990 3 2 1991 3 1 我想要一个二进制变量来标
  • 如何在R中分离两个图?

    每当我运行这段代码时 第一个图就会简单地覆盖前一个图 R中有没有办法分开得到两个图 plot pc title main abc xlab xx ylab yy plot pcs title main sdf xlab sdf ylab x
  • 斯皮尔曼相关性和联系

    我正在一小组配对排名上计算斯皮尔曼的 rho 斯皮尔曼因处理领带不当而闻名 例如 取2组8个排名 即使两组中有6个是平局 相关性仍然很高 gt cor test c 1 2 3 4 5 6 7 8 c 0 0 0 0 0 0 7 8 met
  • 闪亮井板宽度

    library shiny library shinydashboard ui lt dashboardPage dashboardHeader dashboardSidebar dashboardBody wellPanel tags d

随机推荐

  • 有没有办法从 Google Fit API 读取每日静息心率?

    iOS Google Fit 应用程序将静息心率作为其指标之一 但是 我在 Android Google Fit 应用程序中没有看到它 我想使用 Google Fit 数据在 Android 应用程序中显示静息心率 然而 Android G
  • Android 的 Context.openFileOutput() 使用什么文件系统路径?

    我不明白为什么 Android 开发者文档中没有这个问题的答案 我发现他们总是令人沮丧 重新Context类上的openFileOutput 方法打开一个文件进行写入 它写入到什么内部存储文件路径 http developer androi
  • 浏览器中的 JavaScript 是否可能发生并发读/写读/写冲突?

    我遇到的情况是 我正在进行多次 比如四个 ajax 调用 使用 AngularJS http get 如果这很重要的话 并且我希望每次调用都回调并增加一个计数器 这样我就可以知道所有 四个 线程何时完成完全的 我担心的是 由于 JavaSc
  • GStreamer:如何连接动态打击垫

    我正在尝试使用 GStreamer 播放文件中的 MP4 视频 我已成功使用 playbin2 并在命令提示符下使用以下命令播放该文件 gst launch filesrc location bbb mp4 decodebin2 autov
  • 如何使我的 Rails 应用程序的每个独角兽工作者日志到不同的文件?

    如何让 Rails 应用程序的每个独角兽工作者写入不同的日志文件 原因 混合日志文件的问题 在默认配置中 Rails 会将其日志消息写入单个日志文件 log
  • 为什么我的返回类型毫无意义? [复制]

    这个问题在这里已经有答案了 我正在尝试使用返回类型const MyClass const 但是 我收到警告 警告 815 D 返回类型上的类型限定符毫无意义 这不是一个有效的类型吗 我想要一个不能更改的指针 并且我希望它指向的东西也不能被更
  • 如何使 gcc ftrapv 工作?

    我将下面的代码编译为 gcc Wall ftrapv test c 然而 运行生成的可执行文件总是打印 2147483648 这不是我所期望的 我正在运行 gcc 版本 4 6 3 Ubuntu Linaro 4 6 3 1ubuntu5
  • 如何在 symfony 1.4 中创建/使用自定义类和助手?

    将自定义库或辅助方法放入 symfony 中的最佳方法是什么 我在我的项目中使用了学说 我考虑放置的一个地方是 project root lib vendor MyClasses 但是 如果我想创建一个类或辅助函数 它将使用一些核心 sym
  • 使用 jQuery 解析 JSON 响应

    我正在我的一个应用程序中处理 JSON 响应 我已经使用 jsonp 成功建立了连接 但我无法解析我的回复 Code
  • 在 Windows 中将 python 与 Anaconda 结合使用

    我刚刚在计算机上安装了 Anaconda 因为我需要使用 Numpy 好吧 当我使用 python 时 由于某种原因 我必须位于与python exe当然 现在我想使用 Anaconda 我必须在Anaconda3 Scripts文件夹在哪
  • 无效的对象名称 - 存储过程

    我正在通过 SSMS 在 SQL Server 中创建一个存储过程 我已经编写了下面的存储过程 但是当我单击execute出现错误 消息 208 第 16 级 状态 6 程序 NewQuestion 第 11 行 对象名称 hgomez N
  • python 跟踪分段错误

    我正在从 python 开发 C 扩展 并且遇到了一些段错误 在开发过程中不可避免 我正在寻找一种方法来显示发生段错误的代码行 一个想法就像跟踪每一行代码 我该怎么做 如果你在linux上 在gdb下运行python gdb python
  • C++11 lambda 表达式末尾的括号

    我对 C 11 lambda 遇到的一些示例感到困惑 例如 include
  • Django:没有名为“app”的模块

    姜戈呕吐了 ImportError at store No module named store 但就在那里有调试消息和设置 INSTALLED APPS django contrib auth django contrib content
  • 是否无法使用 JSON.stringify 对错误进行字符串化?

    重现问题 我在尝试使用网络套接字传递错误消息时遇到了问题 我可以使用复制我面临的问题JSON stringify为了迎合更广泛的受众 node v0 10 15 gt var error new Error simple error mes
  • 最准确的编码检测器是什么? [关闭]

    Closed 此问题正在寻求书籍 工具 软件库等的推荐 不满足堆栈溢出指南 目前不接受答案 经过一番调查 我发现java界有一些编码检测项目 如果getEncoding in InputStreamReader不起作用 儒尼沃萨尔沙代 jc
  • 字体大小相对于用户的屏幕分辨率?

    我有一个流畅的网站 菜单占其宽度的 20 我希望正确测量菜单的字体大小 以便它始终适合框的宽度 并且永远不会换行到下一行 我正在考虑使用 em 作为单位 但它与浏览器的字体大小相关 因此当我更改分辨率时 字体大小保持不变 还尝试了分数和百分
  • Cookie 不与 Fetch 一起存储

    我已经阅读了我能找到的所有其他主题 但没有一个解决方案有效 我正在使用 React Redux Express 并尝试将 JWT 存储在 cookie 中 如下所示 https auth0 com blog 2015 09 28 5 ste
  • GSON 将特定字段的整数值动态转换为布尔值

    如何处理获取同名但不同类型的字段 我在同一请求中有时从 API 获取整数值 有时获取布尔值 我想知道当我得到这样的 Json 时如何处理 我创建了类型适配器 但它不起作用 我考虑过创建不同的 POJO 类 但这个问题不仅仅针对一个请求 由于
  • 功能“rollapply”的更快替代方案

    我需要对包含大约 7 000 行和 11 000 列的 xts 数据运行滚动窗口函数 我做了以下事情 require PerformanceAnalytics ssd60 lt rollapply wddxts width 60 FUN f