Sıralı bir medyanı bulmak için idiomatik yöntem

Bir ordinalin medyanını bulmam gerekiyor (yani, sipariş edilen faktör).

Bunu yapmak için standart kütüphanede bir yöntem bulamadım, bu yüzden aşağıdaki tıknaz çözümü buldum:

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])
}

Bu işte idiomatik çözüm ne olurdu?

2
quantile (x, 0.5, type = c (1,3)) sayısal olmayan işlem x
katma yazar James, kaynak

1 cevap

Bunu biraz basitleştirebilirsiniz (ve sıralı 'nın sıralı faktörlerin sınıfı olduğunu unutmayın, bu nedenle bunu yalnızca median (o) ile çağırabilirsiniz. code> değişkeninizdir):

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))
}

Kullanımı:

median(ordered(c("A", "B", "C"))) 
median(ordered(c("A", "B", "A", "B")))
5
katma
@csgillespie: Kabul edildi. Uyarıyı ekledim.
katma yazar Richie Cotton, kaynak
Orta noktayı "bölmek" gerektiğinde büyük olasılıkla bir uyarı eklemelisiniz. Örneğin, c ("A", "A", "B", "B") ortanca değeri nedir
katma yazar csgillespie, kaynak