aiFun: Sampling (co)variances and correlations In nadiv: Functions to construct non-additive genetic relatedness matrices

Description

This function returns the sampling covariances, variances, and correlations of the random effects fitted in an mixed model solved using the Average Information algorithm

Usage

 `1` ```aiFun(model = NULL, AI.vec = NULL, inverse = TRUE, Dimnames = NULL) ```

Arguments

 `model ` A model object returned by a call to the `asreml` function. `AI.vec ` A numeric vector of the Average Information matrix. The order must be the row-wise lower triangle of the matrix (including the diagonal). `inverse ` A logical indicating whether the elements of the Average Information matrix are being provided. If FALSE, the Average Information matrix (and not its inverse) is being supplied. `Dimnames ` A vector of characters if names are desired for the output (co)variance/correlation matrix. If not specified, either the default labels from the `asreml` object will be used or the rows and columns will be un-labeled.

Details

The inverse of the Average Information matrix provides the sampling (co)variance of each random term in the mixed model. If a model from the ASReml-R function is supplied (`model` is not NULL), this function extracts the inverse of the AI matrix from an ASReml-R model and organizes it so that the sampling covariances between random terms are located below the diagonal, the sampling variances of random terms are located along the diagonal, and the sampling correlations between random terms are located above the diagonal. The order of the variances along the diagonal is the same as the order entered in the random section of the `asreml` function. This is also the same order as the rows of a call to the summary function, `summary(model)\$varcomp`.

If `model` is NULL then `AI.vec` should contain the vector of values from an Average Information matrix. The function will then reconstruct this matrix, invert it, and supply the sampling (co) variances and correlations for the random terms in the model as described above. Note, either `model` or `AI.vec` must be supplied, but not both.

Value

A matrix of k x k dimensions is returned, if k is the number of random effects estimated in the model. Sampling covariances are below the diagonal, variances along the diagonal, and correlations are above the diagonal. If `Dimnames` is specified, the row and column names are assigned according the vector of names in this argument.

Note

The vector of `Dimnames` should match the same order of variance components specified in the model.

Author(s)

[email protected]

References

Gilmour, A.R., Gogel, B.J., Cullis, B.R., & Thompson, R. 2009. ASReml User Guide Release 3.0. VSN International Ltd., Hemel Hempstead, UK.

Examples

 ``` 1 2 3 4 5 6 7 8 9 10 11 12 13``` ``` ## Not run: library(asreml) ginvA <- asreml.Ainverse(warcolak)\$ginv ginvD <- makeD(warcolak[,1:3])\$listDinv warcolak\$IDD <- warcolak\$ID warcolak.mod <- asreml(trait1 ~ 1, random = ~ped(ID) + giv(IDD), ginverse = list(ID = ginvA, IDD = ginvD), data = warcolak) summary(warcolak.mod)\$varcomp aiFun(model = warcolak.mod, Dimnames = c("Va", "Vd", "Ve"), inverse = TRUE) ## End(Not run) output <- c(7.3075921, 7.0635161, 12.3423380, 1.9539486, 2.7586340, 0.6626111) aiFun(AI.vec = output, inverse = FALSE, Dimnames = c("Va", "Vd", "Ve")) ```

nadiv documentation built on May 31, 2017, 2:38 a.m.