CNVkit segment command RuntimeError: Subprocess command failed #901
Description
I am using CNVkit to call somatic CNAs from tumor-normal paired WES data. The BAM files were generated using GATK pipeline and the somatic VCF file was called using Mutect2 pipeline. However, I am getting the following error while running cnvkit segmentation command.
(cnvkit) iit@lab:/P03-TD$ cnvkit.py segment P03-TD.cnr -o P03-TD.cns -m cbs --smooth-cbs --vcf P03.somatic.vcf --sample-id P03-TD --normal-id P03-BD --min-variant-depth 10
Selected test sample P03-TD and control sample P03-BD
Loaded 3536 records; skipped: 0 somatic, 9 depth
VCF normal sample's genotypes are all 0/0 or missing; inferring genotypes from allele frequency instead
Skipping 99 additional somatic records based on T/N genotypes
Kept 31 heterozygous of 3437 VCF records
Segmenting with method 'cbs', significance threshold 0.0001, in 1 processes
Re-segmenting on variant allele frequency
.
.
Re-segmenting on variant allele frequency
Traceback (most recent call last):
File "/home/iit/miniconda3/envs/cnvkit/lib/python3.10/site-packages/cnvlib/core.py", line 23, in call_quiet
proc = subprocess.run(args, check=True, capture_output=True)
File "/home/iit/miniconda3/envs/cnvkit/lib/python3.10/subprocess.py", line 526, in run
raise CalledProcessError(retcode, process.args,
subprocess.CalledProcessError: Command '('Rscript', '--no-restore', '--no-environ', '/tmp/tmpgodv9oax')' returned non-zero exit status 1.
The above exception was the direct cause of the following exception:
Traceback (most recent call last):
File "/home/iit/miniconda3/envs/cnvkit/bin/cnvkit.py", line 10, in
sys.exit(main())
File "/home/iit/miniconda3/envs/cnvkit/lib/python3.10/site-packages/cnvlib/cnvkit.py", line 10, in main
args.func(args)
File "/home/iit/miniconda3/envs/cnvkit/lib/python3.10/site-packages/cnvlib/commands.py", line 1022, in _cmd_segment
results = segmentation.do_segmentation(
File "/home/iit/miniconda3/envs/cnvkit/lib/python3.10/site-packages/cnvlib/segmentation/init.py", line 81, in do_segmentation
rets = list(
File "/home/iit/miniconda3/envs/cnvkit/lib/python3.10/site-packages/cnvlib/segmentation/init.py", line 126, in _ds
return _do_segmentation(*args)
File "/home/iit/miniconda3/envs/cnvkit/lib/python3.10/site-packages/cnvlib/segmentation/init.py", line 209, in _do_segmentation
seg_out = core.call_quiet(
File "/home/iit/miniconda3/envs/cnvkit/lib/python3.10/site-packages/cnvlib/core.py", line 30, in call_quiet
raise RuntimeError(
RuntimeError: Subprocess command failed:
$ Rscript --no-restore --no-environ /tmp/tmpgodv9oax
Command '('Rscript', '--no-restore', '--no-environ', '/tmp/tmpgodv9oax')' returned non-zero exit status 1.
(cnvkit) iit@lab:P03-TD$
I have tried many ways but still getting this error. R package 'DNACopy' is successfully installed in my system.
Please let me know how to solve this!