search

liusihan / seGMM

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

Error at GMM model #14

Open Lucianabarros opened 7 months ago

Lucianabarros commented 7 months ago

Hello,

I am trying to run seGMM for the first time and I got an error at the end. I was running only four samples to test. It generates several outputs but fails to write the last file.

Could you help solve this issue?

The analysis was ok and can I use the files?

Beginning to generate features at Fri Jan 19 15:58:39 2024

Collected feature of X chromosome heterozygosity Finish generate features of X chromosome heterozygosity at Fri Jan 19 15:58:41 2024

Collected feature of X mapping rate Finish generate features of X mapping rate at Fri Jan 19 16:01:01 2024

Collected feature of Y mapping rate Finish generate features of Y mapping rate at Fri Jan 19 16:01:01 2024

Collected feature of mean depth of SRY gene Finish generate features of mean depth of SRY gene at Fri Jan 19 16:01:03 2024

Combine features into a single file

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

Analysis complete for seGMM at Fri Jan 19 16:01:04 2024 Total time elapsed: 2.0m:25.00999999999999s

liusihan commented 7 months ago

Hi, it seems that the problem is caused by only four specific samples. To resolve this, I would recommend merging your current dataset with relevant reference data.