AE3: AE model using nuclear family trios

In gap: Genetic Analysis Package

AE model using nuclear family trios

Description

AE model using nuclear family trios

Usage

AE3(model, random, data, seed = 1234, n.sim = 50000, verbose = TRUE)

Arguments

model	a linear mixed model formula, see example below.
random	random effect, see exampe below.
data	data to be analyzed.
seed	random number seed.
n.sim	number of simulations.
verbose	a flag for printing out results.

Details

This function is adapted from example 7.1 of Rabe-Hesketh et al. (2008). It also procides heritability estimate and confidence intervals.

Value

The returned value is a list containing:

• lme.result the linear mixed model result.

• h2 the heritability estimate.

• CL confidence intervals.

Note

Adapted from f.mbf.R from the paper.

Author(s)

Jing Hua Zhao

References

\insertRef

rh08gap

Examples

## Not run: 
require(gap.datasets)
AE3(bwt ~ male + first + midage + highage + birthyr,
    list(familyid = pdIdent(~var1 + var2 + var3 -1)), mfblong)

## End(Not run)

gap documentation built on Sept. 11, 2024, 5:36 p.m.