月度数据的季节性汇总
我有数据框每个x,y点. 我正在尝试获得季节性的总体. 我需要计算季节性手段,即 对于冬季平均值(12月,1月,2月);对于春季平均值(3月,4月,5月),夏季平均值(6月,7月,8月)和秋季平均值(9月,10月,11月,11月). 数据看起来类似于: set.seed(1) df
8 2024-01-19
编程技术问答社区
按周汇总(R)
在r中经常按月(以zoo为单位)汇总每日数据,使用类似的内容: result apply.weekly(zoo(1:10, as.Date("2010-01-01") + 1:10), mean) 2010-01-03 2010-01-10 2010-01-11 3 42 10
2 2024-01-19
编程技术问答社区
从时间中减去月份,格式为年-月
我想从给定的日期和每月格式的给定日期减去几个月. global_date = "2017-01" 我正在用zoo库将其转换为: as.yearmon(global_date) - 0.1 ,但它给了我2016年11月,我希望它为'201612' 我该如何在r? 中做到这一点 解决方案 由于我们想减去一个月,我们应该减去1/12 0.083而不是0.1 library(zoo) as.yearmon(global_date) - (1/12) #[1] "Dec 2016" 如果我们需要以上述格式输出 format(as.yearmon(global_date) - (1/12), "%Y%m") #[1] "201612" 其他解决方案 仅使用base R,仅对OP代码进行最小更改 format(as.Date(paste0(global_date, "-01")) - 0.1*10, "%Y%m") #[1] "201612"
6 2023-12-26
编程技术问答社区
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
编程技术问答社区
带置信区间的滚动回归(tidyverse)。
这与组在组中按组滚动回归整理? 再次考虑这个简单的示例 library(dplyr) library(purrr) library(broom) library(zoo) library(lubridate) mydata = data_frame('group' = c('a','a', 'a','a','b', 'b', 'b', 'b'), 'y' = c(1,2,3,4,2,3,4,5), 'x' = c(2,4,6,8,6,9,12,15), 'date' = c(ymd('2016-06-01', '2016-06-02', '2016-06-03', '2016-06-04', '2016-06-03', '2016-06-04', '2016-06-05','20
6 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
编程技术问答社区
矩阵中的R-环路
i有两个变量,第一个是包含230个数据的1D流量,第二个是2D温度矩阵(230*44219). 我试图找到每个流动值与相应44219温度之间的相关矩阵.这是我下面的代码. Houlgrave_flow_1981_2000 = window(Houlgrave_flow_average, start = as.Date("1981-11-15"),end = as.Date("2000-12-15")) > str(Houlgrave_flow_1981_2000) ‘zoo’ series from 1981-11-15 to 2000-12-15 Data: num [1:230] 0.085689 0.021437 0.000705 0 0.006969 ... Index: Date[1:230], format: "1981-11-15" "1981-12-15" "1982-01-15" "1982-02-15" ... Hulgrave_SST_1981_
4 2023-12-26
编程技术问答社区
对多列的时间序列进行线性回归
我有以下问题.我想以矩阵形式计算年度时间序列的回归.我总共有56个时间序列从我想检查的区域的网格点中提取,因此我每年有56个值.我将所有值绘制为数字中的点.现在,我想在此图中添加回归线,其中包含所有数据. 我的目标是计算整个矩阵的回归. library(zoo) pdf(file="/home/user/name.pdf", pointsize=20, onefile = FALSE, width=18, height=11, paper = "special") plot(mat.zoo[,1], pch=20, type="p", ylim=c(8,max(mat.zoo)), yaxt = "n", xaxt = "n", lwd = 1.5, main = "Some title", ylab = "ylabtext", xlab ="", col = "black") tt
8 2023-12-26
编程技术问答社区
R-没有 .fill_short_gaps 的文档
我正在查看动物园软件包的R源代码(许多功能非常有用).我注意到一个函数 .fill_short_gaps 使用了很多,但是我在动物园源代码或基本源代码中都找不到任何文档. . 这是内部功能吗?这个功能应该做什么? 解决方案 这是一个内部功能.源文件的第661版的检查注释说:"使用基本R编码样式约定作为内部非出口函数:.fill_short_gaps()而不是fillShortGaps()." 我在R-Forge上找到了来源: .fill_short_gaps()在该文件的底部. 由于该函数最近被重命名,因此您应该确保使用该参考的所有库都使用Zoo的兼容版本.
6 2023-12-26
编程技术问答社区
按日期分组的滚动平均数
我有一天的数据.每天总是具有1个以上的价值,并且日子不一定是连续的. 我想在特定的一天和滚动14天的平均值中为平均值创建一个计算.我尝试在R中这样做,但我没有太多运气.我想我很亲密. 基本上是按日期分组的均值和常规平均值. 非常感谢,因为我觉得自己很近,所以我会发疯! library(plyr) library(zoo) help=ddply(data, .(DATE), roll_avg14 =rollmean(score, 14, align="right", na.pad=TRUE ),mean_on_that_day = mean(score)) #help 如果您想与一个示例一起进行 - 将其复制为数据框架 structure(list(DATE = structure(c(16700, 16702, 16703, 16704, 16705, 16706, 16707, 16708, 16709, 16710, 16712,
8 2023-12-26
编程技术问答社区
Rollapply:by.column 有什么作用?
我在手册中读过by.column的描述,但我不明白如何使用它.请参阅下面: x=matrix(1:60,nrow=10) library('zoo') rollapply(x,3,mean,fill=NA,align="right",by.column=FALSE) [1] NA NA 27 28 29 30 31 32 33 34 当我使用by.column= FALSE:它适用于宽度(3)线滚动数mean(x[1:3,]) 现在,如果我使用by.column=TRUE,我得到: x=matrix(1:60,nrow=10) rollapply(x,3,mean,fill=NA,align="right",by.column=TRUE) [,1] [,2] [,3] [,4] [,5] [,6] [1,] NA NA NA NA NA NA [2,] NA NA NA NA
2 2023-12-26
编程技术问答社区
R: 为什么我的矢量会紊乱?
说我有这个动物园向量 mine mine[14] 01/15/70 1 > mine[5] 01/06/70 1 > mine[41] 02/11/70 1 它有效! 但是现在我尝试以不同的方式进行 zz
2 2023-12-26
编程技术问答社区
读取 CSV 时出现 "索引在数据行处有坏条目 "错误
我有一个.csv文件,其中包含以下格式的股票价格 Date Open High Low Close 3/7/2014 10:15 228.28 236.58 228.23 236.18 3/7/2014 11:15 236.23 241.27 236.09 241.14 3/7/2014 12:15 241.1 246.12 241.1 242.96 3/7/2014 13:15 242.84 243.92 242 242.32 当我运行以下脚本时: test
2 2023-12-26
编程技术问答社区
将季度年度转换为 R 中季度的最后日期
我使用as.Date(as.yearqtr(test[,1],format ="%qQ%Y"),frac =1)时会有一个问题,但是它返回错误,而Quater年迄今没有更改.错误是: error in as.yearqtr(as.numeric(x)) (list) object cannot be coerced to type 'double' 这是我在r. 中的数据框 在 TIME VALUE 1Q2019 1 2Q2019 2 3Q2019 3 4Q2019 4 理想的输出是 在 TIME VALUE 2019-03-31 1 2019-06-30 2 2019-09-30 3 2019-12-31 4 解决方案 我们可以使用动物园转换为日期,并使用frac获得本季度的最后日期.我们在zoo的合适格式中使用一些正则重新排列: df$TIME=as.Date(as.yearqtr(gsub("(\\d)(Q)
8 2023-12-26
编程技术问答社区
求每七天后的总和,但 R 中缺少天数
我的R DF如下: | date_entered | returning | new | | ------------ | --------- | --- | | 2021-06-02 | 0 | 14 | | 2021-06-03 | 12 | 8 | | 2021-06-04 | 8 | 0 | | 2021-06-05 | 7 | 15 | | 2021-06-07 | 10 | 4 | | 2021-06-08 | 12 | 22 | | 2021-06-09 | 4 | 15 | | 2021-06-10 | 0 | 28 | | 2021-06-12 | 22 | 16 | | 2021-06-15 | 20 | 5 | | 2021-06-16 |
2 2023-12-26
编程技术问答社区