1

我有一个缺失值的不平衡季度面板数据集。我想A2从A1后续季度中减去变量。请注意，我不想得到 A2 的差异，而是从彼此中减去不同的变量。应为每个 uid 分别计算差异。除了像 1999 年第 4 季度和 2000 年第 1 季度这样的变化年份之外，还意味着随后会发生变化。

我真的不确定我是否应该在这里连接我的时间索引，因为像 zoo 这样的包只需要一个索引。但这不是这里的问题。这是一些示例数据：

structure(list(uid = c(1, 1, 1, 2, 2, 3, 3, 3), tndx = c(1999.4, 
2000.1, 2000.2, 1999.4, 2000.1, 2000.1, 2000.2, 2000.3), A1 = c(2, 
2, 2, 10, 11, 1, 1, 1), A2 = c(3, 3, 3, 14, 14, 2, 100, 2)), .Names = c("uid", 
"tndx", "A1", "A2"), row.names = c(NA, -8L), class = "data.frame")

# which results in
  uid   tndx A1  A2
1   1 1999.4  2   3
2   1 2000.1  2   3
3   1 2000.2  2   3
4   2 1999.4 10  14
5   2 2000.1 11  14
6   3 2000.1  1   2
7   3 2000.2  1 100
8   3 2000.3  1   2

如果您更喜欢分隔索引，请使用以下示例：

# Thx Andrie!
x2 <- data.frame(x, colsplit(x$tndx, "\\.", names=c("year", "qtr")))

有没有一个很好的方法来解决这个问题reshape2， plyr甚至是 base 还是你更愿意编写一个自定义函数？

请注意，也有可能uid只出现一次。显然，您无法计算滞后差异。我仍然需要检查并创建一个 NA 。

4

2 回答 2

2

我们在uidusing上拆分它，by并在对单个行的每组行进行操作的函数中uid，我们创建一个动物园对象z，，使用yearqtr类作为索引。然后我们将时间序列与一个空序列合并，该序列具有所有所需的季度，包括任何缺失的中间季度，zm并执行计算zz。最后我们转换为data.frame出路的形式：

library(zoo)
to.yearqtr <- function(x) as.yearqtr(trunc(x) + (10*(x-trunc(x))-1)/4)

DF <- do.call("rbind", by(x, x$uid, function(x) {
    # columns of x are: uid tndx A1 A2
    z <- zoo(x[c("A1", "A2")], to.yearqtr(x$tndx))
    zm <- merge(z, zoo(, seq(start(z), end(z), 1/4)))
    zz <- with(zm, cbind(zm, `A1 - A2 lag` = A1 - lag(A2, -1)))
    if (ncol(zz) <= ncol(z)) zz$`A1 - A2 lag` <- NA # append NA if col not added
    data.frame(uid = x[1, 1], tndx = time(zz), coredata(zz), check.names = FALSE)
}))

这给出了这个：

> DF
    uid    tndx A1 A2 result A1 - A2 lagged
1.1   1 1999 Q4  2  3     NA             NA
1.2   1 2000 Q1  2  2     NA             -1
1.3   1 2000 Q2  2  3     NA              0
2.1   2 1999 Q4  2  4     NA             NA
2.2   2 2000 Q1 NA NA     NA             NA
2.3   2 2000 Q2 NA NA     NA             NA
2.4   2 2000 Q3 NA NA     NA             NA
2.5   2 2000 Q4 NA NA     NA             NA
2.6   2 2001 Q1  3  4     NA             NA
3.1   3 2000 Q1  1  2     NA             NA
3.2   3 2000 Q2  1 NA     NA             -1
3.3   3 2000 Q3  1  2     NA             NA

编辑：根据进一步的讨论完全重新做了解决方案。请注意，这不仅会添加额外的列，还会将索引转换为"yearqtr"类并添加额外的缺失行。

编辑：函数中的一些小的简化by。

于 2011-11-10T13:58:43.320 回答
1

我并不完全清楚你想要什么，因为你没有包含“正确答案”。如果你想从另一个非滞后变量中减去一个滞后变量，你可以使用偏移的索引来做到这一点。（如果您希望将结果放回数据框中，则确实需要填充结果。

 x$A1lagA2 <- ave(x[, c("A1", "A2")], x$uid, FUN=function(z) {
            with(z, c(NA, A1[2:NROW(z)] -A2[1:(NROW(z)-1)]) ) } )[[1]]
 x
  uid   tndx A1  A2 A1lagA2
1   1 1999.4  2   3      NA
2   1 2000.1  2   3      -1
3   1 2000.2  2   3      -1
4   2 1999.4 10  14      NA
5   2 2000.1 11  14      -3
6   3 2000.1  1   2      NA
7   3 2000.2  1 100      -1
8   3 2000.3  1   2     -99

当参数是多列时，您确实会使用 ave() 得到烦人的重复额外列，但我只取了第一个。

于 2011-11-10T14:58:39.863 回答