MutSig2CV icon indicating copy to clipboard operation
MutSig2CV copied to clipboard

Error when I run MutSig2CV using hg38

Open lh12565 opened this issue 1 year ago • 1 comments

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!

lh12565 avatar Feb 01 '24 06:02 lh12565

mutsig2cv only support hg19

xxx58 avatar Jun 18 '24 07:06 xxx58