gravatar for chualec

2 hours ago by

Hi,
I'm comparing two datasets using GRanges to see how many peaks are retained or leave after a certain treatment. However, when I do the code, it doesn't quite add up.
I have 2 datasets called "untreated" and "treated"

peaks_stay <- intersect(untreated, treated)

peaks_leave <- setdiff(untreated, treated)

length(peaks_stay)

16850

length(peaks_leave)

66873

length(untreated)

50883

I am not quite sure why there are more peaks identified using setdiff compared to the original. Can anyone help me with this?

Thanks in advance!

link

modified 2 hours ago

written
2 hours ago
by

chualec0



Source link