bootstrapMI            package:maigesPack            R Documentation

_C_a_l_c_u_l_a_t_e _b_o_o_t_s_t_r_a_p _p-_v_a_l_u_e_s _f_o_r _m_u_t_u_a_l _i_n_f_o_r_m_a_t_i_o_n (_M_I) _m_e_a_s_u_r_e_s

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

     This function takes a numerical matrix (or two vectors) and
     calculates bootstrapped (by permutation) p-values to test if the
     mutual information value is equal to zero. If the first argument
     is a matrix, the p-values are calculated between all pairs of rows
     of the matrix.

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

     bootstrapMI(x, y=NULL, bRep, ret="p-value")

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

       x: numerical matrix or vector to be analysed. If a vector, the
          argument 'y' must be informed.

       y: numerical vector. Must be informed if 'x' is a vector. If 'x'
          is a matrix, this argument is ignored. Defaults to NULL.

    bRep: number of permutation to be done in the test.

     ret: character string with the value to return. Must be 'p-value'
          (default) for the usual p-value or 'max', to return the
          maximum absolute correlation value obtained by the
          permutation.

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

     The method implemented in this function is proposed by Butte and
     Kohane (2000). The MI value is calculated using the function 'MI'.

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

     The result of this function is a square matrix (length equal to
     the number of rows of 'x') if 'x' is a matrix or a numerical value
     if 'x' and 'y' are vectors. The result is the p-values or maximum
     MI values calculated by permutation tests.

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

     Gustavo H. Esteves <gesteves@vision.ime.usp.br>

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

     Butte, A.J. and Kohane, I.S. Mutual information relevance
     networks: functional genomic clustering using pairwise entropy
     measurements. In Pacific Symposium on Biocomputing, 5, 415-426,
     2000 (<URL:
     http://psb.stanford.edu/psb-online/proceedings/psb00/>)

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

     'MI'

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

     x <- runif(50, 0, 1)
     y <- rbeta(50, 1, 2)
     bootstrapMI(x, y, bRep=100)

     z <- matrix(rnorm(100, 0, 1), 4, 25)
     bootstrapMI(z, bRep=100)

