search

liusihan / seGMM

A new tool to infer sex from massively parallel sequencing data.
MIT License
13 stars 2 forks source link

Error #10

Open Abusamra opened 9 months ago

Abusamra commented 9 months ago

I run seGMM using command:seGMM -vcf /ddn/data/BDC/Analysis/NGS_Analysis/habusamra/230920_A00840_0041_AHGMT3DRX3/Alignment_Germline_Aziz/PBG-862-23_S3_HaplotypeCaller.vcf -i test.txt -a BAM -t WES -n 20 -o /ddn/data/BDC/Analysis/NGS_Analysis/habusamra/seGMM/

I got error in >> Running sample classification based on GMM model stage: Error in if (G[1] == 1) { : missing value where TRUE/FALSE needed Calls: Mclust -> eval -> eval -> mclustBIC Execution halted

liusihan commented 9 months ago

The information you provided suggests that the error may be due to having only one sample for classification.

To resolve this issue, you can 1) combine your data with the provided reference data in the reference forder or 2) generate a dataset comprising a minimum of 10 samples and execute seGMM without requiring reference data.