cn.mops returning only "NA" as copy number
1
0
Entering edit mode
rc16955 • 0
@rc16955-12563
Last seen 7.1 years ago

Dear all,

I've been trying to use cn.mops to identify copy number variants among 10 individual of my (non-human) study species. Each individual was subjected to Illumina whole-genome sequencing, and the reads were mapped to the reference using Bowtie2. I have mostly been following the cn.mops manual:

>BAMFiles<-list.files(pattern=".bam$")

>bamDataRanges<-getReadCountsFromBAM(BAMFiles,sampleNames=paste("sample",1:10),mode="paired")

>resCNMOPS<-cn.mops(bamDataRanges)

All of this went through without error. However, when I try to view the data (with cnvs or cnvr), everything in the "logical" column, which I believe should contain the estimates for copy number simply says "<NA>" What could be the cause of this?

I am able to view plots of the data with plot and segplot just fine.

 

Thank you

cn.mops • 1.1k views
ADD COMMENT
0
Entering edit mode
@gunter-klambauer-5426
Last seen 3.3 years ago
Austria
This is the default behaviour. The function "calcIntegerCopyNumbers" replaces the NA values with integer copy numbers. Regards, Günter On 2017-03-11 00:15, rc16955 [bioc] wrote: > Activity on a post you are following on support.bioconductor.org > <https: support.bioconductor.org=""> > > User rc16955 <https: support.bioconductor.org="" u="" 12563=""/> wrote > Question: cn.mops returning only "NA" as copy number > <https: support.bioconductor.org="" p="" 93721=""/>: > > Dear all, > > I've been trying to use cn.mops to identify copy number variants among > 10 individual of my (non-human) study species. Each individual was > subjected to Illumina whole-genome sequencing, and the reads were > mapped to the reference using Bowtie2. I have mostly been following > the cn.mops manual: > > >BAMFiles<-list.files(pattern=".bam$") > > >bamDataRanges<-getReadCountsFromBAM(BAMFiles,sampleNames=paste("sample",1:10),mode="paired") > > >resCNMOPS<-cn.mops(bamDataRanges) > > All of this went through without error. However, when I try to view > the data (with cnvs or cnvr), everything in the "logical" column, > which I believe should contain the estimates for copy number simply > says "<na>" What could be the cause of this? > > I am able to view plots of the data with plot and segplot just fine. > > Thank you > > ------------------------------------------------------------------------ > > Post tags: cn.mops > > You may reply via email or visit > cn.mops returning only "NA" as copy number > -- Günter Klambauer, PhD Institute of Bioinformatics Johannes Kepler University, Linz, Austria http://www.bioinf.jku.at/people/klambauer/
ADD COMMENT

Login before adding your answer.

Traffic: 854 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