muvar: Means and variances under 1- and 2- locus (biallelic) QTL...

In gap: Genetic Analysis Package

Means and variances under 1- and 2- locus (biallelic) QTL model

Description

Means and variances under 1- and 2- locus (biallelic) QTL model

Usage

muvar(
  n.loci = 1,
  y1 = c(0, 1, 1),
  y12 = c(1, 1, 1, 1, 1, 0, 0, 0, 0),
  p1 = 0.99,
  p2 = 0.9
)

Arguments

n.loci	number of loci, 1=single locus, 2=two loci.
y1	the genotypic means of aa, Aa and AA.
y12	the genotypic means of aa, Aa and AA at the first locus and bb, Bb and BB at the second locus.
p1	the frequency of the lower allele, or the that for the first locus under a 2-locus model.
p2	the frequency of the lower allele at the second locus.

Details

Function muvar() gives means and variances under 1-locus and 2-locus QTL model (simple); in the latter case it gives results from different avenues. This function is included for experimental purpose and yet to be generalized.

Value

Currently it does not return any value except screen output; the results can be kept via R's sink() command or via modifying the C/R codes.

Note

Adapted from an earlier C program written for the above book.

Author(s)

Jing Hua Zhao

References

\insertRef

sham97gap

Examples

## Not run: 
# the default 1-locus model
muvar(n.loci=1,y1=c(0,1,1),p1=0.5)

# the default 2-locus model
muvar(n.loci=2,y12=c(1,1,1,1,1,0,0,0,0),p1=0.99,p2=0.9)

## End(Not run)

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