findSigLevelTrad           package:KCsmart           R Documentation

_F_i_n_d _s_i_g_n_i_f_i_c_a_n_c_e _l_e_v_e_l

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

     Method to find the cutoff at which gains and losses are considered
     significant using permutations

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

     findSigLevelTrad(data, observedSpm, n = 1, p = 0.05, maxmem = 1000)

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

    data: aCGH data in the same format as used for 'calcSpm' 

observedSpm: A sample point matrix as produced by 'calcSpm'

       n: Number of permutations 

       p: Alpha level for significance 

  maxmem: This parameter controls memory usage, set to lower value to
          lower memory consumption 

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

     The number of permutations needed for reliable results depends on
     the data and can not be determined beforehand. As a general
     rule-of-thumb around 100 permutations should be used for 'quick
     checks' and around 2000 permutations for more rigorous testing.

     p is the uncorrected alpha level, the method corrects for multiple
     testing internally using simple Bonferroni correction. See the
     referenced publication for more details.

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

     A list with the cutoffs corresponding to the given alpha level

    pos : The cutoff for the gains

    neg : The cutoff for the losses'

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

     Jorma de Ronde

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

     'plotScaleSpace'

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

     data(hsSampleData)
     data(hsMirrorLocs)

     spm1mb <- calcSpm(hsSampleData, hsMirrorLocs)

     sigLevel1mb <- findSigLevelTrad(hsSampleData, spm1mb, n=3)

     plot(spm1mb, sigLevels=sigLevel1mb)
     plotScaleSpace(list(spm1mb), list(sigLevel1mb), type='g')

