sim {aster} | R Documentation |
Data on life history traits for four years and five fitness components
data(sim)
Loads nine objects.
The objects beta.true
, mu.true
, phi.true
, and
theta.true
are the simulation truth parameter values in
different parametrizations.
resp ~ varb + 0 + z1 + z2 + I(z1^2) + I(z1*z2) + I(z2^2)
.
The objects fam
, pred
, and vars
specify the aster model graphical and probabilistic structure.
The objects ladata
and redata
are the simulated data
in two forms "wide"
and "long"
in the terminology
of the reshape
function.
y
, z1
,
z2
used for Lande-Arnold type estimation of fitness landscape.
y
is the response, fitness, and z1
and z1
are
predictor variables, phenotypes.resp
, z1
,
z2
, varb
, id
, root
used for aster type estimation of fitness landscape.
resp
is the response, containing all components of fitness,
and z1
and z1
are predictor variables, phenotypes.
varb
is a factor whose levels are are elements of vars
indicating which elements of resp
go with which nodes of the
aster model graphical structure. The variables z1
and z2
have been set equal to zero except when grep("nseed", varb)
is
TRUE
. For the rationale see Section 3.2 of TR 669 referenced
below.
Geyer, C. J and Shaw, R. G. (2008) Supporting Data Analysis for a talk to be given at Evolution 2008. Technical Report No. 669. School of Statistics, University of Minnesota. http://www.stat.umn.edu/geyer/aster/.
Geyer, C. J and Shaw, R. G. (2009) Hypothesis Tests and Confidence Intervals Involving Fitness Landscapes fit by Aster Models. Technical Report No. 671. School of Statistics, University of Minnesota. http://www.stat.umn.edu/geyer/aster/.
data(sim) out6 <- aster(resp ~ varb + 0 + z1 + z2 + I(z1^2) + I(z1*z2) + I(z2^2), pred, fam, varb, id, root, data = redata) summary(out6) lout <- lm(y ~ z1 + z2 + I(z1^2) + I(z1*z2) + I(z2^2), data = ladata) summary(lout)