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
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.
# NOT RUN {# the default 1-locus modelmuvar(n.loci=1,y1=c(0,1,1),p1=0.5)
# the default 2-locus modelmuvar(n.loci=2,y12=c(1,1,1,1,1,0,0,0,0),p1=0.99,p2=0.9)
# }