R如何每月回滚12个月并进行总结

2 回答

• 0

  使用magrittr表示法定义函数 count_ids ，该函数采用逗号分隔ID的向量并输出唯一ID的数量 .

  然后将其转换为宽动态范围对象，每个组一列聚 Contract 时聚合，以便使用 toString 连接特定日期/组中的所有ID . 在其上运行 rollapply 以计算过去12个月中的唯一ID .

  最后将数据转换回长格式数据框 .

  请注意， read.zoo 将使用 NA 填写任何缺失的年/月计数，只要每个组中没有缺少年/月 . 显示的数据没有每个组中缺少的年/月，事实上，几乎没有任何缺失值，所以它应该工作 . （取消注释第一行以安装动物园的开发版本 - 这是必需的 . ）

  # install.packages("zoo", repos="http://R-Forge.R-project.org")
  library(magrittr)
  library(zoo)
  
  # . is a vector of comma-separated strings of IDs
  count_ids <- . %>%
      textConnection %>%
      scan(what = "", sep = ",", quiet = TRUE, strip.white = TRUE) %>%
      unique %>%
      length
  
  foo <- read.csv("foo_complete.csv")
  foo %>%
      read.zoo(index = "Date", split = "Group", aggregate = toString) %>%
      rollapplyr(12, count_ids) %>%
      fortify.zoo(melt = TRUE, names = c("Date", "Group", "Count"))
  

  如果一个宽动物园对象可以作为输出（参见下面的输出），那么我们可以省略最后一行 .

  仅

  动物园

  如果我们希望在没有magrittr的情况下表达上述内容，我们可以写：

  library(zoo)
  
  count_ids <- function(x) {
      s <- scan(textConnection(x), what = "", sep = ",", quiet = TRUE, strip.white = TRUE)
      length(unique(s))
  }
  
  foo <- read.csv("foo_complete.csv")
  z <- read.zoo(foo, index = "Date", split = "Group", aggregate = toString)
  z.roll <- rollapplyr(z, 12, count_ids)
  fortify.zoo(z.roll, melt = TRUE, names = c("Date", "Group", "Count"))
  

  rollapplyr输出

  rollapplyr 行的输出如下（最后一行的输出太长而无法显示）：

  > foo <- read.csv("foo_complete.csv")
  > foo %>%
  +     read.zoo(index = "Date", split = "Group", aggregate = toString) %>%
  +     rollapplyr(12, count_ids) 
             Min sport  Xox
  2016-07-01  17   137 2165
  2016-08-01  17   142 2273
  2016-09-01  19   149 2361
  2016-10-01  19   154 2429
  2016-11-01  21   156 2532
  2016-12-01  22   161 2578
  2017-01-01  22   163 2585
  2017-02-01  25   163 2594
  2017-03-01  25   167 2632
  2017-04-01  25   168 2648
  2017-05-01  26   166 2659
  2017-06-01  23   166 2686
  2017-07-01  24   162 2740
  2017-08-01  27   163 2758
  2017-09-01  28   162 2788
  2017-10-01  28   165 2831
  2017-11-01  28   158 2786
  2017-12-01  29   153 2688
  

  更新

  • 现在需要动物园的开发版本

  • 固定情况，其中ID出现在多个月内

  • 将 summarize 语句合并到 read.zoo 聚合参数中，以便它现在使用magrittr和动物园

  • 在 rollapply 中添加右对齐

  • 显示 rollapplyr 的输出

  回复于 2024-05-17T04:30:29+08:00
• 2

  我曾经用rubridate包更好地处理日期 . 我修改了你的filter子句并添加了两个for循环，猜测你的第一个日期是1.1.2015 . 其余的只是你的示例代码 .

  library(lubridate)
  result_list <- list()
  for(y in 2015:2017){
    for(m in 1:12){
      d <- as.Date(paste0(y,"-",m, '-01'), format = '%Y-%m-%d')
      d2 <-  d %m-% months(12)
      result_list[[paste0(y,'_',m)]] <- foo %>%
        filter(Date >= d2 & Date < d) %>%
        group_by(ID, Group) %>%
        count() %>%
        group_by(Group) %>%
        count
  
    }
  }
  

  我把结果放到一个列表中 . 也许您必须修改它以更容易地绘制结果 .

  回复于 2024-05-17T04:30:29+08:00