Skip to content

Impossible values in varQ matrix #84

@elcortegano

Description

@elcortegano

I am running fastStructure 1.0 isntalled from bioconda. The software is being run with the --full option to save the file with the variances of the Q matrix (*varQ).

Inspecting the file. however, I have noted that the results in there do not make sense. This is one example of row for a file with K = 11:

0.090920  0.090920  0.090920  0.090920  0.090920  0.090920  0.090920  0.090920  0.090920  0.090920  64250.090802

The 11th field is too high for the variance of a probability. I could assume here that 6425 as been wrongfully appended to 0.090802. But this does not solve all issues.

There is also the problem about how unlikely is that all variances are the same across clusters. The situation shown above repeats in several situations.

Other rows show an even more complicated picture, e.g.:

32748.272252  0.090920  16221.312401  0.090920  0.090920  0.090920  0.090920  0.090920  15280.687988  0.090920  0.090920

Here, the issue of concateanting numbers no longer applies. It is impossible for me to deduce what the real numbers could be.

Finally, I have also noted that all rows sum to approximately the same value of "variance", which is suspicious.

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions