polyRAD icon indicating copy to clipboard operation
polyRAD copied to clipboard

PolyRAD VCF2RADdata Error - "Either possiblePloidies or taxaPloidy must consist entirely of even numbers" - Is this a bug?

Open YI-164 opened this issue 6 months ago • 2 comments

Dear PolyRAD developers/community, I am encountering an error when using the VCF2RADdata function in PolyRAD and would like to confirm whether this is expected behavior or a potential bug.

Error message:
Either possiblePloidies or taxaPloidy must consist entirely of even numbers.

My setup:

Sample ploidy: Taro (Colocasia esculenta) - mixture of diploid samples and triploid samples Variant calling: FreeBayes with triploid settings Sample information: CSV file containing ploidy levels for each sample

Code used:

rsample_info <- read.csv("ploidy.csv", stringsAsFactors = FALSE)

# Create taxaPloidy vector
ploidy_vector <- sample_info$Ploidy
names(ploidy_vector) <- sample_info$SampleName

# Extract possible ploidies
mydata <- VCF2RADdata(
  file = "satoimo_filtered.vcf",
  possiblePloidies = list(2,3),
  taxaPloidy = ploidy_vector,
  min.ind.with.reads = 50,
  min.ind.with.minor.allele = 2,
  expectedAlleles = 500000,
  yieldSize = 500000,
  refgenome = "Taro_Lachesis_assembly_Chr.fa"
)

Question: Is this error expected behavior when working with triploid organisms? The error suggests that PolyRAD requires even ploidy numbers, but taro is naturally triploid (3n). Should I:

Modify my approach to work with triploid data in PolyRAD? Use different parameters or preprocessing steps? Is this a limitation of the current PolyRAD version for odd-ploidy organisms?

Any guidance on handling triploid organisms in PolyRAD would be greatly appreciated. Thank you for your time and assistance.


I am not good at English, so I used Claude to make this message.

YI-164 avatar Aug 04 '25 04:08 YI-164