Error in if (all(x <= minReadCount) & lambda <= minReadCount) { : missing value where TRUE/FALSE needed
0
0
Entering edit mode
s1812756 • 0
@s1812756-21082
Last seen 4.8 years ago

Dear Dr Klambauer,

I met an error when using singlecn.mops. At the modelling part, it runs well at the beginning, but when operating on NW_020228985.1, an error occurred as below, but it runs well when I try it using the same reference genome. Will it matter if the reference genome is not in the same directory?

......
Reference sequence:  NW_020228924.1
Reference sequence:  NW_020228957.1
Reference sequence:  NW_020228985.1
Error in if (all(x <= minReadCount) & lambda <= minReadCount) { : 
  missing value where TRUE/FALSE needed
Calls: singlecn.mops -> lapply -> FUN -> .singlecn.mops
Execution halted

Thank you Mengyuan Li

cnmops cn.mops error • 936 views
ADD COMMENT
0
Entering edit mode

You have NAs in x probably.

ADD REPLY
0
Entering edit mode

Thank you. And, I finally found that problem. The solution is to use "samtools view name.bam | awk '{print $3}' | uniq -c" to get the read count from the BAM file. And use the chromosome names in the second line as refSeqNames parameter

ADD REPLY
0
Entering edit mode

For what it's worth sequence names are available from within R using Rsamtools quickly and easily, for example, after running example(countBam) to get a path to a bam file fl

> seqnames(seqinfo(BamFile(fl)))
[1] "seq1" "seq2"
> idxstatsBam(fl)$seqnames
[1] seq1 seq2
Levels: seq1 seq2
> seqnames(seqinfo(BamFile(fl)))
[1] "seq1" "seq2"
ADD REPLY

Login before adding your answer.

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