CopywriteR can't finish calculating copy number
1
0
Entering edit mode
@thestaroceanster-10072
Last seen 21 months ago

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:

[1] "P1194_neoplasm_reclapse.final_chr_sorted_properreads.bam"
[2] "P1194_NC.final_chr_sorted_properreads.bam"              

Any ideas?

copywriter bam • 708 views
ADD COMMENT
1
Entering edit mode
@tkuilman-6868
Last seen 17 months ago
Netherlands

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)
}

Please find more information in this previous post Error using countBam (Rsamtools package).

ADD COMMENT
0
Entering edit mode

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.
 

ADD REPLY
0
Entering edit mode

Ok, I hope this solves the issue for you. Let me know how things go.

ADD REPLY

Login before adding your answer.

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