fixed bug that counted missing data symbol when counting total number…

… of alleles at each locus.
brannala · Apr 5, 2024 · a320c4c · a320c4c
1 parent 245f8c2
commit a320c4c
Showing 1 changed file with 4 additions and 1 deletion.
diff --git a/src/main.cpp b/src/main.cpp
@@ -1332,7 +1332,10 @@ void checkDataSize(void)
 	{
 	  std::string regex_str = "\\s+"; 
 	  std::vector<std::string> list1 = split(aline, regex_str);
-	  alleleLabels[list1[2]].insert(list1[3]);
+	  /* bug fix Apr 4 2024 don't count missing data symbol "0" as an allele */
+ 	  if(list1[3]!="0")
+	    alleleLabels[list1[2]].insert(list1[3]);
+	  if(list1[4]!="0")
 	  alleleLabels[list1[2]].insert(list1[4]);
 	}
   for (auto all = alleleLabels.begin(); all != alleleLabels.end(); all++) {