getzlab / MutSig2CV

MutSig2CV from Lawrence et al. 2014
Other
30 stars 8 forks source link

Error when I run MutSig2CV using hg38 #30

Open lh12565 opened 8 months ago

lh12565 commented 8 months ago

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!

xxx58 commented 3 months ago

mutsig2cv only support hg19