1.II.13I

Statistical Modelling | Part II, 2007

This problem deals with data collected as the number of each of two different strains of Ceriodaphnia organisms are counted in a controlled environment in which reproduction is occurring among the organisms. The experimenter places into the containers a varying concentration of a particular component of jet fuel that impairs reproduction. Hence it is anticipated that as the concentration of jet fuel grows, the mean number of organisms should decrease.

The table below gives a subset of the data. The full dataset has n=70n=70 rows. The first column provides the number of organisms, the second the concentration of jet fuel (in grams per litre) and the third specifies the strain of the organism.

 number fuel  strain 82015800450.51270.50290.751151.25161.25181.5041.750...\begin{array}{lll} \text { number fuel } & \text { strain } \\ 82 & 0 & 1 \\ 58 & 0 & 0 \\ 45 & 0.5 & 1 \\ 27 & 0.5 & 0 \\ 29 & 0.75 & 1 \\ 15 & 1.25 & 1 \\ 6 & 1.25 & 1 \\ 8 & 1.5 & 0 \\ 4 & 1.75 & 0 \\ . & . & . \end{array}

Explain and interpret the RR commands and (slightly abbreviated) output below. In particular, you should describe the model being fitted, explain how the standard errors are calculated, and comment on the hypothesis tests being described in the summary.

The following R\mathrm{R} code fits two very similar models. Briefly explain the difference between these models and the one above. Motivate the fitting of these models in light of

Part II 2007 the summary from the fit of the one above.

Denote by H1,H2,H3H_{1}, H_{2}, H_{3} the three hypotheses being fitted in sequence above.

Explain the hypothesis tests, including an approximate test of the fit of H1H_{1}, that can be performed using the output from the following R\mathrm{R} code. Use these numbers to comment on the most appropriate model for the data.

>c(fit1$dev>c(f i t 1 \$ d e v, fit2$dev, fit3$dev)

[1] 84.5955786.37646118.9950384.59557 \quad 86.37646118 .99503

>qchisq(0.95,df=1)>\operatorname{qchisq}(0.95, \mathrm{df}=1)

[1] 3.8414593.841459

Typos? Please submit corrections to this page on GitHub.