Mean distance [intra] and Mean distance [inter] are now output in F_S…

…T mode for TN93
veg · Apr 7, 2016 · 1e4f5ea · 1e4f5ea
1 parent 171dc22
commit 1e4f5ea
Showing 1 changed file with 3 additions and 1 deletion.
diff --git a/src/TN93.cpp b/src/TN93.cpp
@@ -503,7 +503,9 @@ int main (int argc, const char * argv[])
              intra_pop =(mean[0]+mean[1])/(weighted_counts[0] + weighted_counts[1]),
              pi_D = mean[2]/weighted_counts[2] - intra_pop;
 
-      (*outStream) << "\t\"Mean distance " << keys[3] << "\" : " << meta << ',' << endl;
+      (*outStream) << "\t\"Mean distance [intra] "  << "\" : " << intra_pop << ',' << endl;
+      (*outStream) << "\t\"Mean distance [inter] "  << "\" : " << mean[2]/weighted_counts[2] << ',' << endl;
+      (*outStream) << "\t\"Mean distance [combined]"  << "\" : " << meta << ',' << endl;
       (*outStream) << "\t\"F_ST\" : " << pi_D/(intra_pop+pi_D) << ',' << endl;
 
       for (unsigned long k = 0; k < 3; k++) {