Data on life history traits for the brown ambrosia aphid Uroleucon rudbeckiae
data(aphid)A data frame with records for 18 insects. Data are already in “long” format; no need to reshape.
Response vector.
Categorical. Gives node of graphical model corresponding
to each component of resp. See details below.
All ones. Root variables for graphical model.
Categorical. Indicates individual plants.
The levels of varb indicate nodes of the graphical model to which
the corresponding elements of the response vector resp belong.
This is the typical “long” format produced by the R reshape
function. For each individual, there are several response variables.
All response variables are combined in one vector resp.
The variable varb indicates which “original” variable
the number was for. The variable id indicates which individual
the number was for. The levels of varb, which are the names
of the “original” variables are the following.
S1 through S13 are Bernoulli: one if alive, zero if dead.
B2 through B9 are conditionally Poisson: the number of
offspring in the corresponding time period. Some variables in the
original data that were zero have been deleted.
These data were published in the following, where they were analyzed by non-aster methods.
Lenski, R. E. and Service, P. M. (1982). The statistical analysis of population growth rates calculated from schedules of survivorship and fecunidity. Ecology, 63, 655-662. tools:::Rd_expr_doi("10.2307/1936785").
These data are reanalyzed by aster methods in the following.
Shaw, R. G., Geyer, C. J., Wagenius, S., Hangelbroek, H. H., and Etterson, J. R. (2008) Unifying life history analyses for inference of fitness and population growth. American Naturalist, 172, E35-E47. tools:::Rd_expr_doi("10.1086/588063").
data(aphid)
### wide version
aphidw <- reshape(aphid, direction = "wide", timevar = "varb",
v.names = "resp", varying = list(levels(aphid$varb)))
Run the code above in your browser using DataLab