library(ggplot2)

age_data <- c(
  1,3,3,3,5,6,6,6,7,8,8,9,9,9,9,9,
  10,10,10,10,10,10,11,11,11,11,11,11,
  12,12,12,12,12,12,13,13,13,13,13,
  14,14,14,14,14,14,15,15,15,15,15,15,15,15,
  16,16,16,16,17,17,17,17,17,17,18,18,18,18,
  19,19,19,20,21,22,22,24,24
)

# Create data
data <- data.frame(
  Age = age_data
)

# Histogram
# Change binwidth number to see how number of bins changes the graph
ggplot(data, aes(x = Age)) +
  geom_histogram(binwidth = 2, fill = "#D44803", color = "black") +
  labs(
    title = "Histogram: Age Distribution with PikBioStat",
    x = "Age",
    y = "Count"
  ) +
  theme_minimal() +
  theme(panel.grid = element_blank())