gap: muvar – R documentation

Become an expert in R — Interactive courses, Cheat Sheets, certificates and more!

muvar

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

Description

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.

Arguments

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

Details

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

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.

References

Sham P (1998). Statistics in Human Genetics. Arnold

Note

Adapted from an earlier C program written for the above book

Author(s)

Jing Hua Zhao

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

Genetic Analysis Package

v1.2.3-1

GPL (>= 2)

Authors

Jing Hua Zhao and colleagues with inputs from Kurt Hornik and Brian Ripley

Initial release

2021-4-21