Ordinamento di un boxplot basato sul valore mediano

Question

Vorrei utilizzare R per creare una serie di grafici a scatole che sono ordinati per valore mediano. Supponiamo poi eseguo:

boxplot(cost ~ type) 

Questo mi darebbe alcuni grafici a scatole sono stati costo è mostrato sull'asse y e la categoria tipo è visibile sul asse x:

-----  ----- 
    |   | 
[ ]  | 
    |  [ ] 
    |   | 
-----  ----- 
    A   B 

Tuttavia, quello che ho' d come sono le figure del boxplot ordinate dal valore mediano più alto al più basso. Il mio sospetto è che quello che devo fare è cambiare le etichette del tipo (A o B) per indicare numericamente quale è il valore medio più basso e più alto, ma mi chiedo se c'è un modo più intelligente per risolvere il problema.

Answer 1

Check out ?reorder. L'esempio sembra essere quello che vuoi, ma ordinato nell'ordine opposto. Ho modificato -count nella prima riga in basso per ordinare nell'ordine desiderato.

bymedian <- with(InsectSprays, reorder(spray, -count, median)) 
    boxplot(count ~ bymedian, data = InsectSprays, 
      xlab = "Type of spray", ylab = "Insect count", 
      main = "InsectSprays data", varwidth = TRUE, 
      col = "lightgray") 

Answer 2

Sì, questa è l'idea:

> set.seed(42)      # fix seed  
> DF <- data.frame(type=sample(LETTERS[1:5], 100, replace=TRUE), 
+     cost=rnorm(100)) 
> 
> boxplot(cost ~ type, data=DF) # not ordered by median 
> 
> # compute index of ordered 'cost factor' and reassign   
> oind <- order(as.numeric(by(DF$cost, DF$type, median)))  
> DF$type <- ordered(DF$type, levels=levels(DF$type)[oind]) 
> 
> boxplot(cost ~ type, data=DF) # now it is ordered by median 

Answer 3

Beware of valori mancanti, si deve aggiungere na.rm = TRUE per farlo funzionare. In caso contrario, il codice semplicemente non funziona. Mi ci sono volute ore per scoprirlo.

bymedian <- with(InsectSprays, reorder(spray, -count, median, **na.rm = TRUE**) 
    boxplot(count ~ bymedian, data = InsectSprays, 
      xlab = "Type of spray", ylab = "Insect count", 
      main = "InsectSprays data", varwidth = TRUE, 
      col = "lightgray")