我有一个数据框。我想创建一个频率表,用“组”显示bin频率。如果有一个包含0个实体的bin,我希望它显示该bin中有0个实体。
如果我使用table()
函数,我会得到数据框中所有二进制数的频率计数,但不是“组”。它也没有告诉我,例如,我在第1组Bin 3中没有任何行。我也查看了tabulate()
,但这似乎不是我需要的。不知怎的,我需要告诉它实际上可能的箱子是什么。
以下是一些示例代码。
df = as.data.frame(rbind(c(1,1.2), c(1,1.4), c(1,2.1), c(1,2.5), c(1,2.7), c(1,4.1), c(2,1.6), c(2,4.5), c(2,4.3), c(2,4.8), c(2,4.9)))
colnames(df) = c("Group", "Value")
df.in = split(df, df$Group)
FindBin = function(df){
maxbin = max(ceiling(df$Value),na.rm=TRUE)+1 #what is the maximum bin value.
bin = seq(from=0, to=maxbin, by=1) #Specify your bins: 0 to the maximum value by increments of 1
df$bin_index = findInterval(df$Value, bin, all.inside = TRUE) #Determine which bin the value is in
return(df)
}
df.out = lapply(names(df.in), function(x) FindBin(df.in[[x]]))
df.out2 = do.call(rbind.data.frame, df.out) #Row bind the list of dataframes to one dataframe
df.out2的输出如下所示:
Group Value bin_index
1 1 1.2 2
2 1 1.4 2
3 1 2.1 3
4 1 2.5 3
5 1 2.7 3
6 1 4.1 5
7 2 1.6 2
8 2 4.5 5
9 2 4.3 5
10 2 4.8 5
11 2 4.9 5
除了上面的输出之外,我还想看一下我的结果的摘要输出,如下所示:
Group Bin Freq
1 1 0
1 2 2
1 3 3
1 4 0
1 5 1
2 1 0
2 2 1
2 3 0
2 4 0
2 5 4
有什么想法吗?
答案 0 :(得分:2)
table
不会为您的第一个问题执行您想要的操作:
df$bin_index <- factor(df$bin_index, levels=1:5)
table(df[, c("Group", "bin_index")])
# bin_index
# Group 1 2 3 4 5
# 1 0 2 3 0 1
# 2 0 1 0 0 4
它显示了bin 3,group 2的0
条目(我认为这就是你的意思,组1中的bin 3有行)。此外,通过设置因子水平,我也能够显示bin_index 1。关于第二个问题,请使用melt
:
library(reshape2)
melt(table(df[, c("Group", "bin_index")]))
# Group bin_index value
# 1 1 1 0
# 2 2 1 0
# 3 1 2 2
# 4 2 2 1
# 5 1 3 3
# 6 2 3 0
# 7 1 4 0
# 8 2 4 0
# 9 1 5 1
# 10 2 5 4