AE model using nuclear family trios

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.

Value

The returned value is a list containing:

• lme.result the linear mixed model result.

• h2 the heritability estimate.

• CL confidence intervals.

Details

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

Note

Adapted from f.mbf.R from the paper.

References

Rabe-Hesketh S, Skrondal A, Gjessing HK (2008). “Biometrical modeling of twin and family data using standard mixed model software.” Biometrics, 64(1), 280-8. doi:10.1111/j.1541-0420.2007.00803.x .

Author

Jing Hua Zhao

Examples

if (FALSE) { # \dontrun{
require(gap.datasets)
AE3(bwt ~ male + first + midage + highage + birthyr,
    list(familyid = pdIdent(~var1 + var2 + var3 -1)), mfblong)
} # }