我需要找到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/