nadiv (version 2.18.0)

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)

Value

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

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.

Author

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

  if (FALSE) {
    library(asreml)
    ginvA <- ainverse(warcolak)
    ginvD <- makeD(warcolak[, 1:3])$listDinv
      attr(ginvD, "rowNames") <- as.character(warcolak[, 1])
      attr(ginvD, "INVERSE") <- TRUE
    warcolak$IDD <- warcolak$ID
    warcolak.mod <- asreml(trait1 ~ sex,
     random = ~ vm(ID, ginvA) + vm(IDD, ginvD), 
	data = warcolak) 
    summary(warcolak.mod)$varcomp
    aiCI(warcolak.mod)    
   }