5 votes
Accepted

Hypergeometic test for the overlap of two set of genomic intervals (e.g. from Chip-seq data)

There is no such thing as a hypergeometric test, at least in statistical textbooks. It's a fisher test based on hypergeometric distribution. If it is chip-seq for the same target, i.e biological ...
  • 1,648
4 votes
Accepted

Occupancy of TFs with the target genes

I would ignore peak calling for this and instead compute enrichment of ChIP/input for the genome (e.g., with deepTools or presumably homer) and then plot it for the genes of interest individually (e.g....
  • 19.3k
3 votes

determine if ChIP-seq peaks are broad or narrow

I'd say mostly this is a question of understanding the underlying biology and the relevant literature. If it is not known in the literature whether a mark is peaky or broad, evidence might come from ...
  • 3,221
3 votes

Detecting broad peaks in sRNA-seq data

The sure why any software should be confined to only broad peaks, but there is two packages that should do what I thing you have in mind. Just to clarify, for smallRNAs peak detection is sometimes ...
2 votes

Score predefined ChIP-seq peaks with MACS2 or equivalent

You could try setting a high p-value threshold when you're calling peaks to retain the "non-significant" values. Something like: ...
  • 1,222
2 votes

Finding peaks and estimating cell population sizes in multi-dimensional flow cytometry data

Automated gating methods are gaining popularity over manual and time consuming analysis in say FlowJo, you should take a look at the openCyto bioconductor package this is a framework which builds on ...
  • 1,029
2 votes

Finding peaks and estimating cell population sizes in multi-dimensional flow cytometry data

There are quite a few algorithms developed for the automatic classification of multidimensional flow cytometry data, you can see a (not so recent?) review here: http://www.nature.com/nmeth/journal/v10/...
1 vote

annotation using ChIPseeker package error

that was solved by changing seq level: ...
  • 95
1 vote

Separating peaks of chip-seq with specific length

You can use rtracklayer from R bioconductor to import it, filter and do much more stuff, i use the code from this blog for ...
  • 1,648
1 vote
Accepted

Separating peaks of chip-seq with specific length

You can do it by awk. awk -v OFS="\t" '{ if ( ($3-$2) >= 500 && ($3-$2) <= 1000) { print } }' in.bed > out.bed Assuming there is no ...
  • 321
1 vote
Accepted

How to calculate the number of peaks that are upstream/downstream of some other peaks

You can use windowBed from bedtools with -c if you just want counts. Or use -wa -wb and then ...
  • 321
1 vote

differential analysis of chip-seq data

For sharing, you can test the proportion of overlapping peaks using bedtools intersect or the find.overlap function in ...
1 vote

Hypergeometic test for the overlap of two set of genomic intervals (e.g. from Chip-seq data)

Just to make sure, I took a quick look at chippeakanno, and it uses a hypergeometric test for gene enrichment, not necessarily overlap of two chip-seq results. As you mentioned, you can imagine in ...
1 vote

How to make Venn diagram for the Macs peak calling output of two data sets?

Venn diagrams can be helpful, but for peak calling they can be somewhat misleading. This topic is addressed in the DiffBind documentation, in the "Comparison of occupancy and affinity based analyses" ...
  • 1,222
1 vote

How to make Venn diagram for the Macs peak calling output of two data sets?

Check out the tool Intervene. It's exactly designed for this. Github link.
  • 21
1 vote
Accepted

How to get the intersection of peaks after peak calling using MACS2?

bedtools intersect with multiple files after -b performs pairwise intersections between -a ...
  • 1,222
1 vote
Accepted

Peak-calling using homer

You just need to strip the / from $d, for which there are a number of options: ...
  • 19.3k

Only top scored, non community-wiki answers of a minimum length are eligible