3

我有一个如下所示的 data.frames 列表：

df=data.frame(
data_id=rep(LETTERS[1:10],each=1),
data_value=c(1,2,2,3,3,2,3,1,1,3))
df2=data.frame(
data_id=rep(LETTERS[1:10],each=1),
data_value=c(2,1,3,1,1,1,2,1,2,1))
df3=data.frame(
data_id=rep(LETTERS[1:10],each=1),
data_value=c(2,2,3,3,1,2,2,1,2,3))
df.list <- list(df, df2, df3)

单个 data.frame 如下所示：

         data_id    data_value
1        A          1
2        B          2
3        C          2
4        D          3
5        E          3
6        F          2
7        G          3
8        H          1
9        I          1
10       J          3

我想计算每个唯一值在 data_value 中出现的频率。我可以做这个：

for(i in 1:length(df.list)){
    daply(df.list[[i]], .(df.list[[i]]$data_value), nrow) -> freq
}

这给了我频率计数（在这种情况下只是最后一个，对于 df3）：

1 2 3 
2 5 3

我的实际数据集要大得多，所以我不能在这里发布。但是，它具有完全相同的结构。问题是，当我尝试获取实际数据集的频率计数时，我收到以下错误消息：

Error in dim(out_array) <- out_dim : dims [product 0] do not match the length of object [1]

谁能告诉我我需要从哪里开始解决这个问题？我不明白 'dim()' 的来源和作用。非常感谢。

4

1 回答 1

8

实际上，您可以做得比这更好，将 for 循环替换为 a laply，这意味着输入是一个列表，输出是一个矩阵/数组。

o <- laply(df.list, function(x) {
    table(x$data_value)
})
> o
#      1 2 3
# [1,] 3 3 4
# [2,] 6 3 1
# [3,] 2 5 3

为了检查你的错误原因，当你尝试这个时会发生什么？

o <- llply(df.list, function(x) {
    table(x$data_value)
})

编辑：为了使错误更容易理解，让我们创建这个data.frame：

d1 <- data.frame(a=1:4)
d2 <- data.frame(a=1:5)
d3 <- data.frame(a=1:6)
d4 <- data.frame(a=1:7)

dl <- list(d1,d2,d3,d4)

现在运行laply：

laply(dl, function(x) table(x$a))
# Error: Results must have the same dimensions.

为什么？要看到这一点，让我们打印它：

> laply(dl, function(x) print(table(x$a)))

# 1 2 3 4 
# 1 1 1 1 
# 
# 1 2 3 4 5 
# 1 1 1 1 1 
# 
# 1 2 3 4 5 6 
# 1 1 1 1 1 1 
# 
# 1 2 3 4 5 6 7 
# 1 1 1 1 1 1 1 

# Error: Results must have the same dimensions.

你看到问题了吗？每行中的元素数量不同。您不能有一个矩阵（除非您将具有较小元素的那些附加到等于行）。

相反，使用列表，以便它们成为列表的元素，以后可以使用[[number]]语法访问。

llply(dl, function(x) table(x$a))

# [[1]]
# 
# 1 2 3 4 
# 1 1 1 1 
# 
# [[2]]
# 
# 1 2 3 4 5 
# 1 1 1 1 1 
# 
# [[3]]
# 
# 1 2 3 4 5 6 
# 1 1 1 1 1 1 
# 
# [[4]]
# 
# 1 2 3 4 5 6 7 
# 1 1 1 1 1 1 1

希望这可以解决问题。

于 2013-02-13T16:18:50.070 回答