muvar                  package:gap                  R Documentation

_M_e_a_n_s _a_n_d _v_a_r_i_a_n_c_e_s _u_n_d_e_r _1- _a_n_d _2- _l_o_c_u_s (_b_i_a_l_l_e_l_i_c) _Q_T_L _m_o_d_e_l

_D_e_s_c_r_i_p_t_i_o_n:

     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.

_U_s_a_g_e:

     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)

_A_r_g_u_m_e_n_t_s:

  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

_V_a_l_u_e:

     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.

_R_e_f_e_r_e_n_c_e_s:

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

_N_o_t_e:

     Adapted from an earlier C program written for the above book

_A_u_t_h_o_r(_s):

     Jing hua Zhao

_E_x_a_m_p_l_e_s:

     ## 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)

