kbyl                   package:gap                   R Documentation

_L_D _s_t_a_t_i_s_t_i_c_s _f_o_r _t_w_o _m_u_l_t_i_a_l_l_e_l_i_c _l_o_c_i

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

     LD statistics for two multiallelic loci

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

     kbyl(n1,n2,h,n,optrho=2)

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

      n1: number of alleles at marker 1

      n2: number of alleles at marker 2

       h: a vector of haplotype frequencies

       n: number of haplotypes

  optrho: type of contingency table association, 0=Pearson,
          1=Tschuprow, 2=Cramer (default)

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

     The returned value is a list containing:

      n1: the number of alleles at marker 1

      n2: the number of alleles at marker 2

       h: the haplotype frequency vector

       n: the number of haplotypes

   VarDp: variance of D'

Dijtable: table of Dij

Dmaxtable: table of Dmax

Dijptable: table of Dij'

VarDijtable: table of variances for Dij

VarDijptable: table of variances for Dij'

      x2: the Chi-squared statistic

    seX2: the variance of x2

     rho: the measure of association

     seR: the standard error of rho

  optrho: the method for calculating rho

  klinfo: the Kullback-Leibler information

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

     Bishop YMM, Fienberg SE, Holland PW (1975) Discrete Multivariate
     Analysis - Theory and Practice, The MIT press

     Cramer H (1946) Mathematical Methods of Statistics. Princeton
     Univ. Press

     Zapata C, Carollo C, Rodriquez S (2001) Sampleing variance and
     distribution of the D' measure of overall gametic disequilibrium
     between multiallelic loci. Ann. Hum. Genet. 65: 395-406

_N_o_t_e:

     adapted from 2ld.c

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

     Jing hua Zhao

_S_e_e _A_l_s_o:

     'tbyt'

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

     ## Not run: 
     # example of two SNPs
     h <- c(0.442356,0.291532,0.245794,0.020319)
     n <- 481*2
     kbyl(2,2,h,n)
     ## End(Not run)

