r - 使用plyr按类别计算最频繁的级别

Question

我想使用下面的代码使用plyr按类别计算最常见的因子级别。数据框b显示请求的结果。为什么c$mlevels只有值“数字”？

require(plyr)
set.seed(0)
a <- data.frame(cat=round(runif(100, 1, 3)),
                levels=factor(round(runif(100, 1, 10))))
mode <- function(x) names(table(x))[which.max(table(x))]
b <- data.frame(cat=1:3,
                mlevels=c(mode(a$levels[a$cat==1]),
                       mode(a$levels[a$cat==2]),
                       mode(a$levels[a$cat==3])))
c <- ddply(a, .(cat), summarise,
           mlevels=mode(levels))

Answer 1

当您使用 summarise , plyr在检查 base 中的函数之前，似乎“看不到”全局环境中声明的函数:

我们可以使用 Hadley 的方便 pryr 来检查这一点。包裹。您可以通过以下命令安装它:

library(devtools)
install_github("pryr")


require(pryr)
require(plyr)
c <- ddply(a, .(cat), summarise, print(where("mode")))
# <environment: namespace:base>
# <environment: namespace:base>
# <environment: namespace:base>

基本上，它不读/不知道/见 您的 mode功能。有两种选择。第一个是@AnandaMahto 的建议，我也会这样做，并建议您坚持下去。另一种选择是不使用 summarise并使用 function(.) 调用它使 mode您的全局环境中的功能是“可见的”。

c <- ddply(a, .(cat), function(x) mode(x$levels))
#   cat V1
# 1   1  6
# 2   2  5
# 3   3  9

为什么这样做？

c <- ddply(a, .(cat), function(x) print(where("mode")))
# <environment: R_GlobalEnv>
# <environment: R_GlobalEnv>
# <environment: R_GlobalEnv>

因为正如您在上面看到的，它读取位于 global environment 中的函数。 .

> mode # your function
# function(x)
#     names(table(x))[which.max(table(x))]
> environment(mode) # where it sits
# <environment: R_GlobalEnv>

与:

> base::mode # base's mode function
# function (x) 
# {
#     some lines of code to compute mode
# }
# <bytecode: 0x7fa2f2bff878>
# <environment: namespace:base>

Here's an awesome wiki在 environments来自哈德利，如果您有兴趣进一步阅读/探索。