Description Usage Arguments Details Value Value Warning References See Also Examples
View source: R/SelectionGradients.R
glamx
is used to fit generalized linear models, specified by error distributions and link functions as denoted by family
, using approaches based on the quantitative framework established by Lande and Arnold (1983). Statistical methods are based on glm
for linear models and glmer
for linear mixed effects models. An option to employ the Janzen and Stern (1998) correction factor for logistic regression models is available via JS = TRUE
. Model formulae are constructed such that regression coefficients and standard errors for quadratic terms do NOT need to be doubled (Stinchcombe et al. 2008). *Note* Only the OLS is currently supported in this version of the function.
1 |
|
Fitness measure. Gaussian fitness types should be use relative fitness, which is calculated as the absolute fitness for each individual |
|
Phenotypic traits. |
|
Choice of whether analyses will be output for linear selection, nonlinear selection, or both. |
The Lande-Arnold Method is based on the 1983 paper by Russell Lande and Stevan Arnold, entitled "The measurement of selection on correlated characters". Their method involves applying ordinary least-squares (OLS) regression to estimate selection gradients.
The function returns an object of classes "glam
", "lm
", and "glm
."
GL
These analyses are currently only available for longitudinal data. Selection gradients for cross-sectional data must be evalated using matrix algebra rather than OLS regressions (Lande and Arnold 1983).
Lande R, Arnold SJ. 1983. The measurement of selection on correlated characters. Evolution 37(6): 1210-1226. http://www.jstor.org/stable/2408842
1 2 3 4 5 6 7 8 9 | # use the BumpusMales data set
data(BumpusMales)
# Define the input
fitness <- BumpusMales$w
z <- BumpusMales[,3:11]
# Calculate the selection gradients using glam
mod1 <- glam(fitness, z, method = "linear")
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.