gap (version 1.1-16)

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

Value

command or via modifying the C/R codes.

References

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

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)

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)