R - résumé dplyr sur des combinaisons de facteurs

Question

Si j'ai une trame de données simple avec 2 facteurs (a et b) avec 2 niveaux (1 et 2) et 1 variable (x), comment puis-je obtenir les valeurs médianes de x : médiane x sur chaque niveau de facteur a, chaque niveau de facteur b, et chaque combinaison de a * b?

library(dplyr)  
df <- data.frame(a = as.factor(c(1,1,1,1,1,1,1,1,2,2,2,2,2,2,2,2)), 
    b = as.factor(c(1,1,1,1,2,2,2,2,1,1,1,1,2,2,2,2)), 
    x = c(runif(16))) 

J'ai essayé différentes (beaucoup) de versions:

df %>% 
    group_by_(c("a", "b")) %>% 
    summarize(med_rate = median(df$x)) 

Les résultats devraient ressembler à ceci pour la médiane x de chaque niveau de facteur a:

une médiane
1 0,58811
2 0,53167

Et comme celui-ci pour la médiane x de chaque niveau de facteur b:

b médiane
1 0,60622
2 0,46096

Et comme celui-ci pour la médiane x pour chaque combinaison de a et b:

ab médiane
1 1 0,66745
1 2 0,34656
2 1 0,50903
2 2 0,55990

Merci à l'avance pour toute aide.

Answer 1

set.seed(123) ##make your example reproducible 
require(data.table) 
df <- data.table(a = as.factor(c(1,1,1,1,1,1,1,1,2,2,2,2,2,2,2,2)), 
      b = as.factor(c(1,1,1,1,2,2,2,2,1,1,1,1,2,2,2,2)), 
      x = c(runif(16))) 

df[, median(x), by = a] 
df[, median(x), by = b] 
df[, median(x), by = .(a,b)] 

Answer 2

Ce qui suit est pas très élégant, mais crée un data.frame unique qui répond à vos résultats escomptés.

Nous créons trois données data.frames (pour a, b et a * b) et de les combiner en un seul.

bind_rows(
    df %>% 
    group_by(a) %>% 
    rename(factor_g = a) %>% 
    summarize(med_rate = median(x)), 
    df %>% 
    group_by(b) %>% 
    rename(factor = b) %>% 
    summarize(med_rate = median(x)), 
    df %>% 
    # We create a column for grouping a*b 
    mutate(factor = paste(a, b)) %>% 
    group_by(factor) %>% 
    summarize(med_rate = median(x)) 
)