getpvalue              package:GraphAT              R Documentation

_F_u_n_c_t_i_o_n _t_o _o_b_t_a_i_n _P _v_a_l_u_e_s _f_r_o_m _t_h_e _E_d_g_e _p_e_r_m_u_t_a_t_i_o_n _a_n_d _N_o_d_e
_p_e_r_m_u_t_a_t_i_o_n _t_e_s_t_s _r_e_s_p_e_c_t_i_v_e_l_y

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

     The function takes as inputs two adjacency matrices. Let X denote
     the observed number of  edges in common between the two adjacency
     matrices.  To test the significance of the correlation between the
     two data sources, the function performs N random edge permutations
     and random node permutations respectively. For each permutation
     test, the function outputs the proportion of N realizations that
     resulted in X edges or more at the intersection of the two
     datasources

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

     getpvalue(act.mat, nonact.mat, num.iterations = 1000)

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

 act.mat: Adjacency matrix corresponding to first data source. That is,
          the i,j element of this matrix is 1 if data source one
          specifies a functional link between genes i and j 

nonact.mat: Adjacency matrix corresponding to first data source. That
          is, the i,j element of this matrix is 1 if data source two
          specifies a functional link between genes i and j

num.iterations: Number of realizations from random edge (node)
          permutation to be obtained

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

     We note that the first adjacency matrix, denoted act.mat is the
     data source that is permutated with respect to edges or notes

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

     A vector of length 2, where the first element is the P value from
     Random Edge Permutation and the second element is the P value from
     Random Node Permutation

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

     Raji Balasubramanian rbalasub@hsph.harvard.edu

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

     'permEdgesM2M', 'permNodesM2M', 'makeClustM'

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

     act.mat <- matrix(0,3,3)
     act.mat[2,1] <- 1 
     act.mat[3,1] <- 1
     nonact.mat  <- matrix(0,3,3)
     nonact.mat[2,1] <- 1
     nonact.mat[3,2] <- 1
     p.val <- getpvalue(act.mat, nonact.mat, num.iterations = 100)
     print(p.val)

