You signed in with another tab or window. Reload to refresh your session.You signed out in another tab or window. Reload to refresh your session.You switched accounts on another tab or window. Reload to refresh your session.Dismiss alert
Came across an issue where using seqGDS2SNP() on a "SeqVarGDSClass" object housing haploid variant calls that had been filtered. Effectively the process alters the genotypes, with the resulting "SNPGDSFileClass" object's genotype matrix not matching at all to the source data (even missing values are replaced). Extracting the genotype array from the source "SeqVarGDSClass" object works fine, and so does extraction of the dosage matrix, both of these match each other (with the dosage value accounting for REF allele count as it should). However it seems that when seqGDS2SNP() parses through the "SeqVarGDSClass" genotype array it produces something else altogether.
I also tried to produce the "SNPGDSFileClass" object with the dosage = T option, and this resulted in:
Error in seqApply(gdsfile, dosage, as.is = gGeno, .progress = verbose, : The GDS node "annotation/format/DS/data" does not exist.
This may be an unrelated issue, but I thought it would be worth noting.
The text was updated successfully, but these errors were encountered:
Came across an issue where using seqGDS2SNP() on a "SeqVarGDSClass" object housing haploid variant calls that had been filtered. Effectively the process alters the genotypes, with the resulting "SNPGDSFileClass" object's genotype matrix not matching at all to the source data (even missing values are replaced). Extracting the genotype array from the source "SeqVarGDSClass" object works fine, and so does extraction of the dosage matrix, both of these match each other (with the dosage value accounting for REF allele count as it should). However it seems that when seqGDS2SNP() parses through the "SeqVarGDSClass" genotype array it produces something else altogether.
I also tried to produce the "SNPGDSFileClass" object with the dosage = T option, and this resulted in:
Error in seqApply(gdsfile, dosage, as.is = gGeno, .progress = verbose, : The GDS node "annotation/format/DS/data" does not exist.
This may be an unrelated issue, but I thought it would be worth noting.
The text was updated successfully, but these errors were encountered: