我需要找到R中的序数(即有序因子)的中位数。

我在标准库中找不到方法来执行此操作,因此我想出了以下笨拙的解决方案:

ordinal.median <- function(x){
        lbls <- levels(x)
        new.vars <- c(NA, 1:length(lbls))
        new.vars[1] <- median(as.numeric(x))
        return(factor(new.vars, labels=lbls, ordered=T)[1])
}


R中对此的惯用解决方案是什么?

最佳答案

您可以简化一下(请注意,ordered是序数因子的类,因此您可以仅使用median(o)来调用它,其中o是您的变量):

median.ordered <- function(x)
{
    levs <- levels(x)
    m <- median(as.integer(x))
    if(floor(m) != m)
    {
      warning("Median is between two values; using the first one")
      m <- floor(m)
    }
    ordered(m, labels = levs, levels = seq_along(levs))
}


用法:

median(ordered(c("A", "B", "C")))
median(ordered(c("A", "B", "A", "B")))

关于r - 求序数中位数的惯用方法,我们在Stack Overflow上找到一个类似的问题:https://stackoverflow.com/questions/7925102/

10-12 19:48