Assignment 7

Due Date

Friday, October 26, 2012

Data Source

The data for this assignment is Pyrausta.txt, a tab-delimited text file.

Background

The data are from Beall (1940) and concern the European corn borer, Pyrausta nubilalis. Table 1 tabulates the number of egg masses found on individual corn plants for three independent samples of corn plants from three consecutive years.

Table 1   Frequency of egg masses on individual corn plants by year
Number of egg masses found Number of corn plants
1936 1937 1938
0 26 37 14
1 19 14 21
2 7 2 11
3 1 3 7
4 3 0 2

Questions

  1. Fit one or more Poisson regression models to test whether the mean number of egg masses is the same or varied by year.
  2. Produce an appropriate graph of the raw data on which you superimpose the predictions of your final model. Your graph should clearly show the distribution of egg mass counts separately by year.
  3. Carry out a goodness of fit test separately by year (three tests). Because we're fitting a fairly minimal model and we don't have many categories to work with I suggest using a parametric bootstrap, i.e., a simulation-based approach.
    1. What do you conclude from the tests?
    2. If you find a significant lack of fit what observation(s) is(are) driving it?
  4. Carry out a goodness of fit test using all three years simultaneously. What do you conclude?

Hints

  1. The number of corn plants examined is not the same in each of the three years.
  2. To answer Question 4 you will need to construct your own parametric bootstrap in which you simulate data from each year separately but you combine the results and use them to calculate a single Pearson statistic. See lecture 17.

Cited references

Course Home Page


Jack Weiss
Phone: (919) 962-5930
E-Mail: jack_weiss@unc.edu
Address: Curriculum in Ecology, Box 3275, University of North Carolina, Chapel Hill, 27599
Copyright © 2012
Last Revised--October 18, 2012
URL: https://sakai.unc.edu/access/content/group/3d1eb92e-7848-4f55-90c3-7c72a54e7e43/public/docs/assignments/assign7.htm