[R-sig-ME] Arctic Nematode Growth Data

Murray Jorgensen maj at stats.waikato.ac.nz
Thu Mar 26 01:47:11 CET 2009


I thought some list members might appreciate the following data set as 
an example suitable for nonlinear mixed effects modelling. It is data 
from Dennis Proctor's PhD thesis which I came by 30 years ago but have 
only just fully entered onto the computer. Unfortunately I have lost 
touch with Dennis but I have no doubt that he would not object to fair 
usage provided attribution was made.

At each age there are 3 replicate groups of nematodes measured. 
Apparently the length measurement was destructive and no nematodes are 
multiply measured.

The R code below reads the data into a data frame as well as creating 
some other structures of interest. I am well aware that constructing a 
vector by concatenation inside a loop is regarded as poor R style, but I 
could not remember what the approved alternative was!

# Arctic Nematode Growth Data from
# Proctor, D. L. C. (1979)  "Energy flow through free-
# living soil nematodes in High Arctic terrestrial
# communities. PhD thesis, Dept of Entomology, U. Alberta.
reps = scan()
1 4 21 7 2 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
5 20 6 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
2 15 10 7 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 1 8 10 13 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 4 9 10 4 6 2 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 1 7 8 13 4 4 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 1 2 3 14 5 5 3 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 1 1 15 13 2 3 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 1 7 12 17 1 2 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 1 0 0 7 13 7 2 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 2 1 12 9 9 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 3 1 2 6 13 7 2 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 1 0 2 2 1 1 4 5 0 2 5 5 4 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 1 12 5 5 5 2 3 4 0 0 0 0 0 0 0 0 0
0 0 0 0 1 0 1 0 0 1 0 4 7 6 4 7 2 1 0 0 0 0 0 0 0 0 0
0 0 0 0 2 1 1 1 1 0 1 2 5 1 1 4 3 2 4 3 0 0 0 0 0 0 0
0 1 0 0 1 0 0 0 1 2 4 2 6 0 1 4 4 4 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 1 0 1 1 0 1 0 0 0 0 0 6 9 9 4 2 0 0 0 0
0 0 0 0 0 1 0 0 1 0 0 1 2 2 1 1 2 0 3 6 11 3 2 0 0 0 0
1 0 0 0 0 1 0 0 0 0 1 1 3 3 1 4 1 3 6 5 1 3 1 0 0 0 0
0 0 0 0 0 0 0 1 0 0 0 0 0 1 1 0 0 0 0 4 8 8 7 4 1 0 0
0 0 0 0 0 0 0 1 0 0 0 2 0 0 1 1 2 6 7 6 5 1 0 0 0 0 0
0 0 0 0 1 0 2 0 0 0 0 0 0 0 0 1 1 2 7 12 7 2 0 1 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 2 2 8 9 8 5 1
0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 1 0 2 11 12 6 2 1 0
0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 3 3 0 1 5 3 7 2 5 2 0
0 0 0 0 0 0 0 0 0 0 0 1 0 1 0 0 0 0 6 7 12 6 0 0 1 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 6 12 6 9 1 0 0
0 0 0 0 0 0 0 0 0 0 0 1 1 0 0 1 1 0 1 3 1 3 6 8 6 3 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 2 3 7 16 5 5 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 2 14 11 2 2 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 1 5 7 11 9 1 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 3 3 5 8 6 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 2 2 6 3 3 7 6 2 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 2 1 2 7 9 10 4 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 1 2 1 10 14 4 1 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 2 3 12 6 4 1 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 2 4 6 8 6 2 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 2 1 6 5 8 5 4 1 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 3 3 2 6 6 6 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 4 0 3 5 11 6 5 1 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 2 4 7 10 4 2 1 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 5 4 8 11 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 3 3 8 10 2 6 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 2 1 1 5 9 5 3 0 0

len = seq(0.22, by=0.04,length=27) # class mid-points
repage = rep(seq(4,60,4),rep(3,15))
fr = matrix(reps,nrow=45,ncol=27,byrow = TRUE)
repsize = apply(fr,1,sum)
age = rep(repage,repsize)    # age of nemetode in days
repl = rep(1:45,repsize)     # replicate number
leng = NULL
for ( row in 1:45 ) {
    leng = c(leng, len[rep(1:27,fr[row,])])
    }
# leng is rounded length in mm
nematodes = data.frame(age,repl,leng)



-- 
Dr Murray Jorgensen      http://www.stats.waikato.ac.nz/Staff/maj.html
Department of Statistics, University of Waikato, Hamilton, New Zealand
Email: maj at waikato.ac.nz                                Fax 7 838 4155
Phone  +64 7 838 4773 wk    Home +64 7 825 0441   Mobile 021 0200 8350




More information about the R-sig-mixed-models mailing list