ld2                   package:gllm                   R Documentation

_E_s_t_i_m_a_t_e _l_i_n_k_a_g_e _d_i_s_e_q_u_i_l_i_b_r_i_u_m _b_e_t_w_e_e_n _t_w_o _c_o_d_o_m_i_n_a_n_t _a_u_t_o_s_o_m_a_l _l_o_c_i

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

     Fits a log-linear model for allelic association between two
     codominant autosomal loci.  Measures of LD are odds ratios.

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

     ld2(locus1, locus2) 

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

  locus1: is a character vector containing the genotypes at the first
          locus, or a RxC contingency table of genotype counts.

  locus2: is a character vector containing the genotypes at the second
          locus.

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

      m0: 

      m1: 

      m2: 

      m3: 

      m4: 

     {estimating LD and HWD coefficient(s)}

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

     MNS<-matrix(c(91,32,5,147,78,17,85,75,7), nr=3)
     colnames(MNS)<-c("S/S","S/s","s/s")
     rownames(MNS)<-c("M/M","M/N","N/N")
     class(MNS)<-"table"
     print(MNS)
     res<-ld2(MNS)
     print(res)

