AE3: AE model using nuclear family trios

Description Usage Arguments Value References Note Author(s) Examples

Description

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

Usage

1
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

References

Rabe-Hesketh S, Skrondal A, Gjessing HK. Biometrical modeling of twin and family data using standard mixed model software. Biometrics 2008, 64:280-288

Note

Adapted from f.mbf.R from the paper

Author(s)

Jing Hua Zhao

Examples

1
2
3
4
5
## Not run: 
AE3(bwt ~ male + first + midage + highage + birthyr,
    list(familyid = pdIdent(~var1 + var2 + var3 -1)), mfblong)

## End(Not run)


Search within the gap package
Search all R packages, documentation and source code

Questions? Problems? Suggestions? or email at ian@mutexlabs.com.

Please suggest features or report bugs with the GitHub issue tracker.

All documentation is copyright its authors; we didn't write any of that.