male_heights <- rnorm(100, mean = 175.4, sd = 7)
female_heights <- rnorm(100, mean = 161.74, sd = 6)
height_data <- data.frame(
height = c(male_heights, female_heights),
gender = factor(rep(c("Men", "Women"), each = 100))
)
ggplot(height_data, aes(x = gender, y = height, fill = gender)) +
geom_violin(trim = F, alpha = 0.6) +
stat_summary(fun = mean, geom = "point", shape = 23, size = 1,
color = "black", fill = "white") +
labs(title = "Findings of Average Height by Gender",
x = "Gender", y = "Height in Centimeters") +
scale_fill_manual(values = c("blue", "pink")) +
theme_minimal()