raphael-group / THetA

Tumor Heterogeneity Analysis (THetA) and THetA2 are algorithms that estimate the tumor purity and clonal/subclonal copy number aberrations directly from high-throughput DNA sequencing data. This repository includes the updated algorithm, called THetA2.
http://compbio.cs.brown.edu/projects/theta/
70 stars 33 forks source link

how to output all possible CNV #21

Open WRui opened 5 years ago

WRui commented 5 years ago

I install THeTA-0.7, when run RunBAFModel.py , I meet error code: python /WPSnew/wangrui/Software/THetA-0.7/python/RunBAFModel.py /WPSnew/wangrui/Project/Colon_Cancer/Exon/OtherPackageTest/CNVkit/P7_WXY_LM2R2_Tumor.tumor.snp_formatted.txt /WPSnew/wangrui/Project/Colon_Cancer/Exon/OtherPackageTest/CNVkit/P7_WXY_LM2R2_Tumor.normal.snp_formatted.txt /WPSnew/wangrui/Project/Colon_Cancer/Exon/OtherPackageTest/CNVkit/Result//P7/P7_WXY_LM2R2.theta2.interval_count /WPSnew/wangrui/Project/Colon_Cancer/Exon/OtherPackageTest/THetA2_V2/P7_WXY_LM2R2/P7_WXY_LM2R2.n2.results -P /WPSnew/wangrui/Project/Colon_Cancer/Exon/OtherPackageTest/THetA2_V2/P7_WXY_LM2R2/BAF/P7_WXY_LM2R2.BAF -O /WPSnew/wangrui/Project/Colon_Cancer/Exon/OtherPackageTest/THetA2_V2/P7_WXY_LM2R2/BAF --PLOT_OPTION ALL

Error: Traceback (most recent call last): File "/WPSnew/wangrui/Software/THetA-0.7/python/RunBAFModel.py", line 544, in main() File "/WPSnew/wangrui/Software/THetA-0.7/python/RunBAFModel.py", line 541, in main run_BAF_model(*kwargs) File "/WPSnew/wangrui/Software/THetA-0.7/python/RunBAFModel.py", line 121, in run_BAF_model plot_results(BAFVec, meansVec, posVec, chrmVec, NLLVec, chrmsToUse, plotOption, directory, prefix, plotDim) File "/WPSnew/wangrui/Software/THetA-0.7/python/RunBAFModel.py", line 242, in plot_results numberResults, fig, colors, plotNum=i) File "/WPSnew/wangrui/Software/THetA-0.7/python/RunBAFModel.py", line 148, in plot_single_result ax = fig.add_subplot(numberResults, 1, plotNum) File "/WPSnew/wangrui/Software/anaconda/lib/python2.7/site-packages/matplotlib/figure.py", line 1257, in add_subplot a = subplot_class_factory(projection_class)(self, args, *kwargs) File "/WPSnew/wangrui/Software/anaconda/lib/python2.7/site-packages/matplotlib/axes/_subplots.py", line 67, in init ).format(maxn=rowscols, num=num)) ValueError: num must be 1 <= num <= 1, not 0