findSigLevelFdr           package:KCsmart           R Documentation

_T_h_i_s _f_u_n_c_t_i_o_n _h_a_s _n_o_t _b_e_e_n _p_r_o_p_e_r_l_y _i_m_p_l_e_m_e_n_t_e_d _y_e_t

_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:

     findSigLevelFdr(data, observedSpm, n = 1, fdrTarget=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 

fdrTarget: Target False Discovery Rate (FDR) 

  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. 
     The FDR method is less conservatie than the p-value based approach
     since instead of controlling the family wise error rate (FWER,
     P(false positive > 1)) it controls the false discovery rate (FDR)
     (false positives / total number of called data points).

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

     A list with the cutoffs corresponding to the given FDR

    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')

