Summarize and print results of linear mixed model fitted with gremlin.
## S3 method for class 'gremlin'summary(object,...)## S3 method for class 'summary.gremlin'print(x, digits = max(3, getOption("digits")-3),...)
Arguments
object, x: An object of class gremlin or summary.gremlin .
...: Additional arguments to be passed to control the output.
digits: An integer used for number formatting with signif() .
Returns
A list of class summary.gremlin or a printed value to the screen with no return values.
logLik: Model log-likelihood.
formulae: Function call and model fixed, random, and residual formulae.
runtime: A numeric of class difftime containing the length of time to run the model. See how this is handled in update.gremlin.
lambda: A logical indicating if the model was transformed to the variance ratio, or lambda scale.
residQuants: A named vector listing summary output for the model residuals.
varcompSummary: Table of variance components and approximate standard errors (calculated from the inverse of the average information matrix). If a (co)variance component is fixed or at the boundary of its parameter space then an NA is returned for the standard error and a column with constraint types is added to the table. Alternative methods (e.g., profile likelihood CIs) should be pursued for obtaining uncertainties associated with fixed or boundary parameters.
varcompSampCor: A matrix containing the sampling correlations of the (co)variance components. Note this is on the underlying nu
scale that the model is fitting.
coefficients: Table of fixed effects and standard errors (calculated from the corresponding diagonal elements of the inverse of the coefficient matrix, transformed where necessary).
Examples
grS <- gremlin(WWG11 ~ sex -1, random =~ sire, data = Mrode11)summary(grS)