R code used to produce the graph shown below.

library(vcd) 
families <- c(3, 24, 104, 286, 670, 1033, 1343, 
              1112, 829, 478, 181, 45, 7) 
boys <- 0:12*families           # Count number of boys 
# Estimate probability of a boy 
phat <- sum(boys) / sum(families*12) 
rootogram(families, 6115*dbinom(0:12, size=12, p=phat), 
          points_gp = gpar(col = fgblue), 
          lines_gp = gpar(col = fgblue, lwd = 3)) 
# grid.ls(fullNames=TRUE) 
grid.edit("GRID.yaxis.20", gp=gpar(fontsize=20)) 
#grid.edit("GRID.text.21", gp=gpar(fontsize=22)) 
grid.edit("GRID.text.22", gp=gpar(fontsize=20)) 
grid.edit("GRID.text.23", gp=gpar(fontsize=20))  # Y akse 
#grid.edit("GRID.text.24", gp=gpar(fontsize=16))  # X akse