使用小时应用功能将10分钟数据聚合为小时平均值失败
我有一个文件/时间数据及其在上述日期和时间的测量值.在一个月的过程中,每十分钟进行一次测量值,我最终试图进行时间序列分析.但是,在此之前,我想通过计算每60分钟的平均测量值将10分钟间隔汇总到每小时间隔.这是我的数据示例(总共有4319个观测值): Date/Time Value 2013-01-01 00:00:00 31,439999 2013-01-01 00:10:00 33,439999 2013-01-01 00:20:00 39 2013-01-01 00:30:00 35,279999 2013-01-01 00:40:00 37,279999 2013-01-01 00:50:00 32,199997 2013-01-01 01:00:00 35,279999 2013-01-01 01:10:00 38,1999
2 2024-01-20
编程技术问答社区
使用lapply和split执行以前的蜱虫聚合
我试图在过去三个月中解决这个问题.请帮助. 我有许多属于单个交易所的股票的刻度数据(价格和数量).每个股票在硬盘上都有自己的.rds文件.我有兴趣清理它: 通过中位数合并多个相同的时间戳 仅用于交换小时的子集 以前 我知道 高频软件包中的功能聚集 可以执行以前的刻度聚合操作.但是,该功能仅占一天的数据数据. 为了证明我正在使用原始tick数据(命名为贸易)的问题. dput(head(trade,50)) structure(c(54.7, 54.7, 54.5, 54.5, 54.5, 54.6, 54.6, 54.65, 54.65, 54.6, 54.65, 54.65, 54.65, 54.65, 54.7, 54.7, 54.8, 54.8, 54.85, 54.85, 54.85, 54.85, 54.8, 54.8, 54.8, 54.8, 54.65, 54.65, 54.8, 54.8, 54.8, 5
4 2024-01-19
编程技术问答社区
R-时间序列中的下一个最高值
一个相对简单的问题,但是我似乎找不到任何示例. 我有一个简单的外汇价格数据,该数据位于2列XTS对象中,称为subx1: Datetime, Price 2016-09-01 00:00:01, 1.11563 2016-09-01 00:00:01, 1.11564 2016-09-01 00:00:02, 1.11564 2016-09-01 00:00:03, 1.11565 ...等等. 我试图在下午2点之后第一次,价格高于下午2点之前,该价格在另一个对象的列中保存在名为daypeakxts $ the the.high and 中. daypeakxts样本是: Date, before2.High 2016-09-01, 1.11567 2016-09-02, 1.11987 这是我要做的事情的一个糟糕例子: subxresult dayp
8 2023-12-26
编程技术问答社区
在 R 中每周分割时间序列
我想每周在R中拆分xts/zoo时间序列.时区设置为"亚洲/加尔各答" Sys.setenv(TZ="Asia/Kolkata") library(xts) seqs head(weekdata[[2]],2) value 2016-01-04 05:30:00 156 2016-01-04 06:00:00 157 > head(weekdata[[3]],2) value 2016-01-11 05
4 2023-12-26
编程技术问答社区
在xts-object的一行中只对负值进行求和
我有一个XTS-Object,我只想在第一行中概括. 代码 sum(test[test dim(test)是216 39. 解决方案 有可能与此处的开放XTS问题有关----> 看来XTS对象不能以这种方式子集.相反,您应该使用coredata(),即将"核心数据"视为矩阵的XTS函数. 有点混乱,但是: sum(coredata(test)[1,][coredata(test)[1,]
6 2023-12-26
编程技术问答社区
用zoo在R中读取CSV
我有以下格式的CSV: TICKER,PER,DATE,TIME,CLOSE SYMBOL,1,20160104,1002,14180.0000000 SYMBOL,1,20160104,1003,14241.0000000 我想将其读成一个时间序列: f
8 2023-12-26
编程技术问答社区
滚动回归返回多个对象
我正在尝试基于示例来构建滚动回归函数>,但是除了返回预测值外,我还想返回一些滚动模型诊断(即系数,T值和Mabye R^2).我希望根据结果类型在离散对象中返回结果.上面链接中提供的示例可成功地创建滚动预测,但是我需要一些帮助包装并写出滚动模型诊断: 最后,我希望函数返回三(3)个对象: 预测 系数 t值 r^2 以下是代码: require(zoo) require(dynlm) ## Create Some Dummy Data set.seed(12345) x
4 2023-12-26
编程技术问答社区
绘制xts对象时,对点有效,对线无效
我已经合并了两个XTS对象,并希望将它们绘制在一个显示中.当我使用点(type =" p")时,这可以正常工作.但是,当我使用行(类型=" l")时发生问题时:仅在第二个系列不涵盖的索引区域中显示第一个系列.我希望这些线与"点"一样长.可重现的示例如下. 由于默认值和ggplot绘图命令都会发生这种情况,因此我怀疑这与某些时间序列数据有关. 这种行为的原因是什么?是否有适当的绘制此类数据的方法? ## Minimal example for Reproduction library(xts) library(ggplot) # create two artificial xts objects xts1
2 2023-12-26
编程技术问答社区
R 如何使用 stackApply 计算每年的季节平均值?
我想计算每年的季节平均值,而不是整个期间.我将自己的季节定义为:DJF(12月至1月),妈妈(3月 - 玛),JJA(6月至8月)和儿子(9月至9月). 灵感来自于解决 fredrick 的启发. ,我创建了一个代表季节的索引"组",然后我应用了命令" stackapply",但是该命令计算了整个期间的平均季节性价值.我解释说,获得的最后一层仅包含4个栅格,但是对于我的情况,我想计算"每年的季节平均值,因此,每年必须有4个射手,总共栅格扇形应有136个射手./p> 我的代码下方 感谢您的帮助 library(raster) set.seed(123) r
2 2023-12-26
编程技术问答社区
如何在ggplot2中绘制xts?
我得到了这个情节 使用此代码 library(xts) library(zoo) plot.new() par(mai=c(2, 2, 1, 1)) plot(nats[,2], ylim=c(0, 2.5), xlab = "", ylab="", main="", major.ticks="months", major.format="%d-%b-%Y", minor.ticks=FALSE, las=1, cex.axis=0.8) points(nats[rans,2], col="darkgrey", pch=16, cex=0.5) points(M1[rans], col="red", pch=16, cex=0.5) points(M2[rans], col="blue", pch=16, cex=0.5) points(M3[rans], col="green", pch=16, cex=0.5) legend(x="
6 2023-12-26
编程技术问答社区
在两列上使用 Rollapply
我正在尝试做类似的事情,我要求在这里不幸的是我无法解决. 这是我的数据框架(数据),一个时间序列: Date Price Vol 1998-01-01 200 0.3 1998-01-02 400 0.4 1998-01-03 600 -0.2 1998-01-04 100 0.1 ... 1998-01-20 100 0.1 1998-01-21 200 -0.4 1998-01-21 500 0.06 .... 1998-02-01 100 0.2 1998-02-02 200 0.4 1998-02-03 500 0.3 1998-02-04 100 0.1 etc. 我想告诉r, 以" Vol"的第一值,将其除以"价格"的20个值,然后 以"
2 2023-12-26
编程技术问答社区
将时间序列转换为带有 NA 的季度序列
我有一个多元zoo时间序列,其中包含NAS,我想将其转换为季度系列. df1 zoo.object stock1 stock2 2013-01-01 1 7 2013-02-05 2 NA 2013-03-12 3 9 2013-04-16 NA 10 2013-05-21 5 11 2013-06-25 6 NA 理想情况下,我想保留每个股票的较早数据. 我已经尝试了XTS软件包,但问题在于它与NAS删除了线条. > to.quarterly(zoo.
2 2023-12-26
编程技术问答社区
不规则的时间序列数据--能否使其规则化?
我的数据已编程为每5小时获取每5小时的信息,这意味着每天多个数据点.问题有时是数据记录器故障或电池死亡或其他任何东西,并且缺少数据或失败的尝试. 我所做的所有分析都是基于这些日期.我需要根据此日期/时间戳进行分类和过滤并选择所有内容.截至目前,我正在浏览"长期挂"的脚本,只是为了选择和子集数据.此外,我还有几个数据记录仪(个人),所以我也在批处理处理.例如,我想每天子集并仅使用一个数据点,而我 #Subset data by minimum time interval SubsetData = T MinInterval = 47 #minimum time interval (in hours) between consecutive locations ########### ##Setting the date and time to be recognized as date time loc$DateTime = as
2 2023-12-26
编程技术问答社区
用XTS/动物园r的一天中的一个小时绘图价值
我确实有一个像这样的时间序列(分钟值): "timestamp", value "2012-04-09 05:03:00",2 "2012-04-09 05:04:00",4 "2012-04-09 05:05:00",5 "2012-04-09 05:06:00",0 "2012-04-09 05:07:00",0 "2012-04-09 05:08:00",3 "2012-04-09 05:09:00",0 "2012-04-09 05:10:00",1 有一种简单的方法可以在一天中的一个小时内绘制这些值:x轴从1到24小时(OR 0和23).因此 - 在第5小时等上5:00到5:59之间的所有值等.这并不取决于日期,我只是对当天的混蛋感兴趣. 谢谢! 其他问题:我可以将其绘制为箱形图吗?现在是plot(df$hh,df$coredata.._data..)以动物园格式. 解决方案 使用x
2 2023-12-26
编程技术问答社区
XTS至周期返回不同的每周终点
我有一个问题,endpoints() xts中的功能(以及使用端点的to.weekly函数)有时在星期五返回,作为一周结束,有时会返回星期一. 我的数据集称为sp2. > head(sp2) [,1] 2012-01-09 1.78 2012-01-10 1.78 2012-01-11 1.77 2012-01-12 1.80 2012-01-13 1.77 2012-01-16 1.77 > tail(sp2) [,1] 2012-06-28 1.94 2012-06-29 1.92 2012-07-02 1.92 2012-07-03 1.90 2012-07-04 1.89 2012-07-05 1.89 > class(sp2) [1] "xts" "zoo" > weekdays(index(sp2)) [1] "Monday" "Tuesday" "Wednesday" "Thursday" "Fr
6 2023-12-26
编程技术问答社区