动物园对象聚合

Question

尊敬的社区，

我收到的数据将在一个数据框中：

Var_1      Var_2         Date        VaR_3  VaR_4   VaR_5   Var_6
1           4       2010-01-18         7    apple    10    sweet
2           5       2010-07-19         8    orange   11    sour
3           6       2010-01-18         9    kiwi     12    juicy
...        ...      ...               ...   ...     ...    ... 

我想使用 zoo，因为它似乎是一个灵活的对象类。我刚从 R 开始，我试图阅读包的描述（小插图）。

问题：

1. 给定以上数据作为数据框，建议用哪种方法将完整的df转换成zoo对象，告诉zoo应该使用第三列作为日期列（日期可以在数据中出现多次）？
2. 如何使用 zoo 内置函数每月聚合所有其他列，除了第 4 列和第 6 列？ zoo 是否能够自动丢弃分类变量而只使用那些适合聚合的列？
3. 如何每月汇总所有数字列，针对第 4 列中的每个类别（第 6 列不应包括在内，因为它不是数字）。

感谢您的支持。

Answer 1

zoo 对象是时间序列，通常是数字向量或矩阵。似乎您真正拥有的是一堆不同的时间序列，其中第 5 列标识了它是哪个序列。即有苹果系列、橙系列、猕猴桃系列等，每个系列都有几列。

删除最后一列，因为它不是数字，使用第三列作为索引并在第 5 列上拆分：

# create test data
Lines <- "Var_1      Var_2         Date        VaR_3  VaR_4   VaR_5   Var_6
1           4       2010-01-18         7    apple    10    sweet
2           5       2010-07-19         8    orange   11    sour
3           6       2010-01-18         9    kiwi     12    juicy"
cat(Lines, "\n", file = "data.txt")

library(zoo)
z <- read.zoo("data.txt", header = TRUE, index = 3, split = "VaR_5",
  colClasses = c(Var_6 = "NULL"))

结果是：

> z
           Var_1.apple Var_2.apple VaR_3.apple VaR_5.apple Var_1.kiwi
2010-01-18           1           4           7          10          3
2010-07-19          NA          NA          NA          NA         NA
           Var_2.kiwi VaR_3.kiwi VaR_5.kiwi Var_1.orange Var_2.orange
2010-01-18          6          9         12           NA           NA
2010-07-19         NA         NA         NA            2            5
           VaR_3.orange VaR_5.orange
2010-01-18           NA           NA
2010-07-19            8           11

以上假设对于第 5 列的给定值，日期是唯一的。如果不是这种情况，则包含aggregate = mean 参数或aggregate 的其他值。

现在将其汇总为我们拥有的每月动物园系列：

aggregate(z, as.yearmon, mean)

也可以使用FUN = as.yearmon 参数直接将其转换为每月：

zm <- read.zoo("data.txt", header = TRUE, index = "Date", split = "VaR_4", 
  FUN = as.yearmon, colClasses = c(Var_6 = "NULL"), aggregate = mean)

请参阅?read.zoo、vignette("zoo-read")、?aggregate.zoo 以及其他小插曲和帮助文件。