Error with cn.MOPS run on exomeseq data
0
3
Entering edit mode
7.2 years ago
vakul.mohanty ▴ 260

Hello,

I'm having trouble with cn.MOPS. I run the following piece of code and receive an error

library(cn.mops)
bam = list.files(path = "/home/trantor/mounts/komurov_lab/SKCM_WES",pattern=".bam$",full.names=TRUE)
seg = read.table(file = "/home/trantor/mounts/komurov_lab/Ensemble_exon", sep="\t",as.is =TRUE)
gr = GRanges(seg[,1],IRanges(seg[,2],seg[,3]))
seqlevels(gr) = gsub("chr","",seqlevels(gr))
gr = disjoin(gr)
#specifing the chromosome to run the analysis on
gr = gr[seqnames(gr) == "1"]
count = getSegmentReadCountsFromBAM(bam, GR=gr, mode="paired")
res = exomecn.mops(count)

Error:

Normalizing...
Error in normalizeChromosomes(X, chr = chr, normType = normType, qu = normQu) :
  Some normalization factors are zero! Remove samples or chromosomes for which the average read count is zero, e.g. chromosome Y

My GRanges object is limited to chromosome 1 and the count object shows counts in every sample so I'm not quiet sure how the normalization factor is zero.

I will be grateful for help with this.

Thanking You,

Vakul

exomeseq cn.MOPS cnmops CNV • 2.0k views
ADD COMMENT

Login before adding your answer.

Traffic: 2479 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6