Question: CopywriteR can't finish calculating copy number
0
3.0 years ago by
thestaroceanster0 wrote:

Hello.

I'm trying to use copywriteR for calculating copy number across samples genome. My problem is that the program seem to be  running for  a short period of time until it suddenly stop and quit to cmd command line.

After showing the message about my job would be ran on multiple cpu, there is no other warning or error message. The last lines of log file shows like follow:

INFO [2016-05-31 00:35:50] CopywriteR has finished filtering low-quality and anomalous reads in the following (unique) samples:

Any ideas?

bam copywriter • 454 views
modified 3.0 years ago by t.kuilman140 • written 3.0 years ago by thestaroceanster0
Answer: CopywriteR can't finish calculating copy number
1
3.0 years ago by
t.kuilman140
Netherlands
t.kuilman140 wrote:

Hi,

To find out what is going on it would be very helpful if you could provide me with the exact commands / script you use to run CopywriteR, the exact common line output, the content of the CopywriteR.log file you find in your destination folder, and a sessionInfo(). One thing we have run into in the past is that if you use very large bam files, the countBam() function of the Rsamtools library does not work anymore (the output is an integer, and integer overflow can occur with large bams). You could test the bams that you want to process using CopywriteR using

library(Rsamtools)

for (bam.file in vector.with.file.paths.of.your.bam.files) {
countBam(bam.file)
}

Thanks for your answer. I also consider the possibility of too large bam files since this program ran fine with previous (and smaller) bam file. The problem seems to solve itself after I ran each sample per time but each run it required the system to be reboot (which I figure out after serveral failure). I think the bam file size could be the source of this problem. I will test the bam files to see if it is indeed the problem later.  Thanks for your help again.