所有可能的分档的频率计数

时间:2014-02-03 21:31:00

标签: r frequency

我有一个数据框。我想创建一个频率表,用“组”显示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

有什么想法吗?

1 个答案:

答案 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