Question: create.set on MEDIPS
0
gravatar for oussama.badad
15 months ago by
oussama.badad10 wrote:

Hi all

i am trying to create set for my medipseq data , here is my command

AKM-ANF_MeDIP = MEDIPS.createSet(file = "/scratch/obadad/MEDIP-Seq_BAMs/AKM-ANF-mapped.sorted.bam", BSgenome = "BSgenome.Oeuropaea.IOGC.v1", extend = 300, shift = 0, uniq = 1e-3, window_size = 100, paired=T)

everything is good till it gets to Calculating short read coverage at genome wide windows...

All chromosomes in the reference BSgenome will be processed: [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" "13" "14" "15" [16] "16" "17" "18" "19" "20" "21" "22" "23" Reading bam alignment AKM-ANF-mapped.sorted.bam considering 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 using bam index Total number of imported short reads: 94574662 Extending reads... Creating GRange Object... Keep at most 2 read(s) mapping to the same genomic location Number of remaining reads: 70758499 Calculating genomic coordinates... Creating Granges object for genome wide windows... Calculating short read coverage at genome wide windows... Error in AKM - ANF_MeDIP = MEDIPS.createSet(file = "/scratch/obadad/MEDIP-Seq_BAMs/AKM-ANF-mapped.sorted.bam", : *object 'AKM' not found*

createset medips • 407 views
ADD COMMENTlink written 15 months ago by oussama.badad10
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1265 users visited in the last hour