search

cruk-mi / mesa

mesa package for Methylation Enrichment Sequencing Analysis
9 stars 3 forks source link

`plotRegionsHeatmap` fails if `regionsToOverlap` has multiple regions that overlap the same 300bp window. #18

Open paddyharker1005 opened 3 months ago

paddyharker1005 commented 3 months ago

If the regions given to regionsToOverlap has some regions that overlap the same 300bp window, plotRegionsHeatmap will give the following error:

Error in '.rowNamesDF<-'(x, value = value) : duplicate 'row.names' are not allowed

To reproduce this error run qseaSet %>% plotRegionsHeatmap(regionsToOverlap = mesa::hg38UltraStableProbes) with any qseaSet. There are some hyperStable probes which overlap the same windows.

I believe this is an issue with the plyranges::join_overlap_left() function in getWindowAnnotation(), leading to multiple matches with the same 300bp window.