| muvar {gap} | R Documentation |
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.
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 )
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. |
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.
Adapted from an earlier C program written for the above book.
Jing Hua Zhao
Sham P (1998). Statistics in Human Genetics. Arnold
## 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)