compareComplex             package:ScISI             R Documentation

_A _f_u_n_c_t_i_o_n _t_o _c_o_m_p_a_r_e _t_w_o _b_i_p_a_r_t_i_t_e _g_r_a_p_h _m_a_t_r_i_c_e_s

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

     The 'compareComplex' function takes two bipartite graph matrix
     representations and calculates three statistics between all pairs
     of complexes, C-i and K-j: (1) the intersection between C-i and
     K-j, (2) the proteins in C-i and not in K-j, (3) the proteins in
     K-j and not in C-i. These stats are later used to calculate
     Jaccard and Dice-Sorenson Coefficients as well as probability
     distributions of a protein p in a complex C.

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

     compareComplex(BGMat1, BGMat2)

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

  BGMat1: The first bipartite graph matrix

  BGMat2: The second bipartite graph matrix

_D_e_t_a_i_l_s:

     The main point to remember is that we don't compare complexes
     within the same matrix. We only compare pairs of complexes from
     the BGMat1 and BGMat2. The 'runCompareComplex' function must be
     called first since all the labelled vertices present in BGMat1
     must be present in the BGMat2 and vice versa so a true comparison
     can be made.

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

     A list with the following entries:

intersect: A matrix of pairwise intersections of the complexes of
          BGMat1 and BGMat2

 cminusk: A matrix of the setdiff of the complexes of BGMat1 and those
          of BGMat2

 kminusc: A matric of the setdiff of the complexes of BGMat2 and those
          of BGMat1

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

     Tony Chiang

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

     go = getGOInfo(wantAllComplexes = FALSE)
     mips = getMipsInfo(wantALLComplexes = FALSE)
     goM = createGOMatrix(go)
     mipsM = createMipsMatrix(mips)
     runCompareComplex(goM, mipsM)

