Barplots avec plusieurs groupes de facteurs et moyenne de variable entre ces facteurs

Question

J'essaie de créer un graphique à barres qui montre les salaires horaires moyens des travailleurs syndiqués et non syndiqués regroupés par célibataires ou mariés regroupés par diplômé ou non. Alors que j'ai réussi à construire un barplot passable avec deux groupes de facteurs, je ne peux pas comprendre comment le faire avec trois groupes de facteurs. Les exemples que j'ai vus qui ont trois facteurs ne tiennent compte que du nombre de fréquences, donc je ne suis pas sûr de savoir comment incorporer la moyenne d'une autre variable dans tous les facteurs de la parcelle. Ce que je cherche à créer quelque chose qui ressemble à ceci (créé en Stata): Average Hourly Wage by Union Status, Marital Status, and College Graduation Mon code ressemble à ceci:

levelbar = tapply(wage, list(as.factor(union), as.factor(married), 
as.factor(collgrad)), mean) 
par(mfrow = c(1, 2)) 
barplot(levelbar, beside = TRUE) 
barplot(t(levelbar), beside = TRUE) 

Quand je lance cela, cependant, je reçois l'erreur:

Error in barplot.default(levelbar, beside = TRUE) : 
'height' must be a vector or a matrix 

Toute aide à ce sujet serait appréciée. Je suis sûr que ggplot pourrait être utile ici, mais je n'ai pas beaucoup d'expérience dans l'utilisation de ce paquet.

Answer 1

Voici un exemple reproductible utilisant ggplot et le jeu de données intégré Titanic. Notez que nous calculons les moyennes d'abord et utilisons stat = identity pour nous assurer que nous obtenons ceux dans l'intrigue.

# Format the Titanic dataframe 
Titanic_df <- Titanic %>% as_tibble() 

# Make Class, Sex, Age, and Survived factors 
for (col in c("Class", "Sex", "Age", "Survived")) { 
    Titanic_df[[col]] <- factor(Titanic_df[[col]]) 
} 

# Get by group means 
means <- Titanic_df %>% 
    group_by(Class, Sex, Survived) %>% 
    summarise(
    mean_n = mean(n) 
) 

# Plot: facets are the Classes, bar colors are the two Sexes, and the groupings in each facet are Survived vs. Not Survived 
ggplot(data = means) + 
    geom_bar(aes(x = Survived, y = mean_n, fill = Sex), stat = "identity", position = "dodge") + 
    facet_wrap(~ Class)