peaksNearBDP(myPeakList, mart,AnnotationData, MaxDistance=5000, PeakLocForDistance = c("start", "middle", "end"), FeatureLocForDistance = c("TSS", "middle", "start", "end","geneEnd"))
peaksWithBDP
percentPeaksWithBDP
n.peaks
n.peaksWithBDP
if (interactive())
{
data(myPeakList)
data(TSS.human.NCBI36)
annotatedBDP = peaksNearBDP(myPeakList[1:6,],
AnnotationData=TSS.human.NCBI36,
MaxDistance=5000,
PeakLocForDistance = "middle",
FeatureLocForDistance = "TSS")
c(annotatedBDP$percentPeaksWithBDP, annotatedBDP$n.peaks,
annotatedBDP$n.peaksWithBDP)
}
Run the code above in your browser using DataCamp Workspace