Search
Question: How to take sample peak from each chromosome in peak bed files ?
0
gravatar for Jurat Shahidin
11 months ago by
Italy
Jurat Shahidin60 wrote:

Hi,

I used data set from Encode consortium for my package development, due to size of actual peak files are rather big, I can't use these data set for my package use. Because actual size of Bioconductor package resulted from R CMD build must be less than 4Mb on disk, I have to use rather small peak file as an example data for my package . In Encode sample's data set, each peak files contains around 100,000 peaks each. How can I obtain small example data of Chip-seq peak files ? Is there any approach getting ideal example data from actual peak files for package use? Thanks in advance :)

Here is quick insight of my actual peak files in my packages :

myPkg
- inst
     - extdata
             - wgEncodeOpenChromChipK562CmycAlnRep1.bed
             - wgEncodeOpenChromChipK562CmycAlnRep2.bed
             - wgEncodeOpenChromChipK562CmycAlnRep3.bed
             - wgEncodeSydhTfbsK562CmycIfna6hStdAlnRep1.bed
             - wgEncodeSydhTfbsK562CmycIfna6hStdAlnRep2.bed
             - wgEncodeSydhTfbsK562CmycIfna30StdAlnRep1.bed
             - wgEncodeSydhTfbsK562CmycIfna30StdAlnRep2.bed
             - wgEncodeSydhTfbsK562CmycIfng6hStdAlnRep1.bed
             - wgEncodeSydhTfbsK562CmycIfng6hStdAlnRep2.bed
             - wgEncodeSydhTfbsK562CmycIggrabAlnRep1.bed
             - wgEncodeSydhTfbsK562CmycIggrabAlnRep2.bed
             - wgEncodeSydhTfbsK562CmycStdAlnRep1.bed
             - wgEncodeSydhTfbsK562CmycStdAlnRep2.bed
- R

 

Edit :

I am looking at Chip-seq for transcription factors.  how can I get small peak files such as each bed file has around 1000 peaks ? 

Should I take sample from each chromosome in each peak file ? How to get sample ? Any approach to get ideal, small example data from Chip-seq peak files? Thanks a lot

Best regards :

Jurat

ADD COMMENTlink modified 11 months ago • written 11 months ago by Jurat Shahidin60
0
gravatar for federico.comoglio
11 months ago by
Switzerland
federico.comoglio100 wrote:

Hi Jurat,

this discussion might be of help: https://www.biostars.org/p/85063/

Federico

ADD COMMENTlink written 11 months ago by federico.comoglio100

Hi Federico :

Thanks for this post from Biostar. Not familiar with using BEDOPS tools on windows. How can I facilitate the process of taking samples from each bed files ? I think  it will take some time to get familiar with BEDOPS tools and try the solution. I need something fast, general approach instead. Any idea ?

Best regards :

Jurat

ADD REPLYlink written 11 months ago by Jurat Shahidin60
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.2.0
Traffic: 141 users visited in the last hour