findRegions               package:ACME               R Documentation

_F_i_n_d _a_l_l _r_e_g_i_o_n_s _i_n _d_a_t_a _a_b_o_v_e _p-_v_a_l_u_e _t_h_r_e_s_h_o_l_d

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

     After the ACME calculation, each probe is associated with a
     p-value of enrichment.  However, one often wants the contiguous
     regions associated with runs of p-values above a given p-value
     threshold.

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

     findRegions(calcobj, thresh = 1e-04)

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

 calcobj: An 'aGFFCalc' object

  thresh: The p-value threshold

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

     Runs of p-values above the p-value threshold will be reported as
     one "region".  These can be used for downstream analyses, export
     to browsers, submitted for transcription factor binding enrichment
     analyses, etc.

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

     A data frame with these columns: 

  Length: The length of the region in probes

      TF: Either TRUE or FALSE; TRUE regions represent regions of
          enrichment while FALSE regions are the regions between the
          TRUE regions

StartInd: The starting Index of the region

  EndInd: The ending Index of the region

  Sample: The sample containing the region

Chromosome: The Chromosome of the region

   Start: The starting basepairof the region

     End: The ending basepair of the region

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

     Sean Davis <sdavis2@mail.nih.gov>

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

     'do.aGFF.calc', 'findClosestGene'

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

     data(example.agff)
     example.agffcalc <- do.aGFF.calc(example.agff,window=1000,thresh=0.9)
     foundregions <- findRegions(example.agffcalc,thresh=0.001)
     foundregions[1:6,]

