MutSig2CV
MutSig2CV copied to clipboard
Error when I run MutSig2CV using hg38
Hi, when I run MutSig2CV using hg38 files like #10. There are some errors as below:
1.000000 [ 1946 377 85 13 3 6 ]
04:51:57 04:51:57 19260/19260 ZZZ3 0 0 CV 1.000000 CL 1.000000 FN 1.000000 [ 1947 377 85 13 3 6 ]
22 genes with q<=0.1
Saving results... Reference to non-existent field 'longname'.
Error in getfield (line 36)
f = s.(deblank(strField)); % deblank field name
Error in keep_fields (line 16)
f = getfield(S,flds{i});
Error in MutSig_2CV_v3_11_core (line 1178)
pgmc = keep_fields(M.gene, {'name', 'longname', 'chr', 'gene_start', 'gene_end',
'tot_exon_len', 'gc'});
Error in MutSig2CV (line 51)
G = MutSig_2CV_v3_11_core(args{:});
Despite these errors, some result files are generated:
final_analysis_set.maf MutSig_version.txt sample_sig_gene_table.txt
mutcateg_discovery.txt patient_counts_and_rates.txt sig_genes.txt
mutcategs.txt results.mat
Can this error be ignored, or are some result files not obtained? Thanks!
mutsig2cv only support hg19