Learn R Programming

QTLRel (version 0.2-15)

genMatrix: Derive genetic matrices

Description

Derive genetic matrices from Jacquard condensed identity coefficients or genotypic data.

Usage

genMatrix(x)

Arguments

x

an object of cic or ibs, or genotypic data in a matrix or a data frame with each row representing an observation and each column a marker locus and entry being "AA", "AB", "BB" (or 1, 2, 3) without missing genotypes.

Value

AA

additive genetic matrix.

DD

dominance genetic matrix.

AD,HH,MH

other three genetic matrices (see Abney et. al. 2000).

ib

inbreeding coefficients.

References

Abney, M., M. S. McPeek, and C. Ober (2000). Estimation of variance components of quantitative traits in inbred populations. Am. J. Hum. Genet. 141, 629-650.

See Also

cic

Examples

Run this code
# NOT RUN {
data(miscEx)

ids<- sample(pedF8$id[300:500],20)

# }
# NOT RUN {
# get condensed identity coefficients
oo<- cic(pedF8, ids=ids, df=0)
ksp<- kinship(pedF8, ids=ids) # kinship coefficients only
# extract genetic matrices
gm<- genMatrix(oo)
sum((gm$AA-2*ksp)>1e-7) # same results
# }

Run the code above in your browser using DataLab