基于 R 中选定月份的子集数据

Question

Year         tot_precip
1/1/1989     0
1/2/1989    .23
1/3/1989    0
1/4/1989    .43
1/5/1989    0.254
1/6/1989    0
1/7/1989    0
1/8/1989    0
1/9/1989    0
1/10/1989   .21

我正在尝试根据某些月份对我的降水数据进行子集化。我做了一些研究。

tt=as.POSIXct(paste(prec$Year,prec$tot_precip), format="%m/%d/%Y")
datZoo <- zoo(prec[,-c(1,2)], tt)
month <- function (x) as.numeric(format(x, "%m"))
veranoIdx <- which(month(tt) %in% 6:8)
veranoZoo <- datZoo[veranoIdx]
veranoZoo

上面的代码可以根据某些月份（6 月、7 月 8 月）提取年份列，但我不知道如何提取这些日期的降水值（即 veranoZoo 系列）

Answer 1

假设Lines 如下注所示，将数据读入动物园对象。由于这显然是月度数据，因此最好使用"yearmon" 类作为索引。现在使用subset 获取六月、七月、八月的子集。

z <- read.zoo(text = Lines, header = TRUE, FUN = as.yearmon, format = "%d/%m/%Y")
subset(z, cycle(z) %in% 6:8)

给这个动物园系列：

Jun 1989 Jul 1989 Aug 1989 
       0        0        0 

这也可以：

z[cycle(z) %in% 6:8]

（如果您的数据在文件 myfile.dat 中，则在 read.zoo 语句中将 text=Lines 替换为 "myfile.dat"。）

注意：

Lines <- "Year tot_precip
1/1/1989 0
1/2/1989 .23
1/3/1989 0
1/4/1989 .43
1/5/1989 0.254
1/6/1989 0
1/7/1989 0
1/8/1989 0
1/9/1989 0
1/10/1989 .21"

Answer 2

考虑输入为：

dft <- read.table(header = TRUE, text = "Year         tot_precip

                  1/1/1989     0

                  1/2/1989    .23

                  1/3/1989    0

                  1/4/1989    .43

                  1/5/1989    0.254

                  1/6/1989    0

                  1/7/1989    0

                  1/8/1989    0

                  1/9/1989    0

                  1/10/1989   .21",stringsAsFactors=FALSE)

您可以使用tidyverse 函数并尝试：

dft %>% filter(month(Year) %in% c(6,7,8))

给出：

      Year tot_precip
1 1/6/1989          0
2 1/7/1989          0
3 1/8/1989          0