nadiv (version 2.9)

aiCI: Confidence Intervals for Variance Components

Description

Produces the 1-alpha Upper and Lower Confidence Limits for the variance components in an ASReml-R model.

Usage

aiCI(asr.model, Dimnames = NULL, alpha = 0.05)

Arguments

asr.model

Object from a call to asreml

Dimnames

A vector of characters if names are desired for the output. If not specified, the default labels from the asreml object will be used.

alpha

A numeric value indicating the level of Type I error for constructing the Confidence Intervals.

Value

• A matrix is returned with a row for each variance component. The three columns correspond to the Upper Confidence Limit, estimate from the asreml model, and Lower Confidence Limit for each variance component.

Details

Variances from the inverse of the Average Information matrix of an ASReml model are translated according to the varTrans function and used in constructing the 1-alpha Confidence Interval.

See Also

aiFun

Examples

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
    aiCI(warcolak.mod)