Using R

The data frame Geissler in the vcdExtra package contains the complete data from Geissler's (1889) tabulation of family sex composition in Saxony. The table below gives the number of boys in families of size 11.

(a) Read these data into R

(b) Following Example 3.13, use goodfit() to fit the binomial model and plot the results. Is there an indication that the binomial does not fit these data?

(c) Diagnose the form of the distribution using the methods described in Section 3.4.

(d) Try fitting the negative binomial distribution, and use distplot() to diagnose whether the negative binomial is a reasonable fit.