-
Notifications
You must be signed in to change notification settings - Fork 5
Description
I am using SNP data in regulatory regions to calculate enrichments. I am using the R package.
My Input data (hg38) looks like this:
GRanges object with 1800 ranges and 0 metadata columns:
seqnames ranges strand
[1] chr1 51851558-51851559 *
[2] chr1 51800569-51800570 *
[3] chr1 51825311-51825312 *
[4] chr1 51832478-51832479 *
[5] chr1 51845209-51845210 *
However, I am not getting any results.
category nb.overlaps random.average mapped.peaks.ratio effect.size p.significance p.value q.significance q.value
TAL1 TAL1 0 0 0 NaN NA NA NA NA
GABPA GABPA 0 0 0 NaN NA NA NA NA
FOXP1 FOXP1 0 0 0 NaN NA NA NA NA
NME2 NME2 0 0 0 NaN NA NA NA NA
NRF1 NRF1 0 0 0 NaN NA NA NA NA
CTCF CTCF 0 0 0 NaN NA NA NA NA
e.significance e.value
TAL1 NA NA
GABPA NA NA
FOXP1 NA NA
NME2 NA NA
NRF1 NA NA
CTCF NA NA
I am not sure what I am doing wrong and would appreciate any feedback.