[adegenet-commits] r405 - pkg/R

[noreply at r-forge.r-project.org]

Author: jombart
Date: 2009-06-16 17:06:36 +0200 (Tue, 16 Jun 2009)
New Revision: 405

Modified:
   pkg/R/seqTrack.R
Log:
best consensus... my ass... 


Modified: pkg/R/seqTrack.R
===================================================================
--- pkg/R/seqTrack.R	2009-06-16 12:43:10 UTC (rev 404)
+++ pkg/R/seqTrack.R	2009-06-16 15:06:36 UTC (rev 405)
@@ -771,8 +771,8 @@
 
 
     if(mode=="best"){
-        toKeep <- which.max(listres$valsim)
-        nbDraws <- sum(listres$valsim > (max(listres$valsim) - 1e-10 )) -1
+        toKeep <- which.min(listres$valsim)
+        nbDraws <- sum(listres$valsim < (min(listres$valsim) + 1e-10 )) -1
         cat("\nThere were\n",nbDraws, "draws.\n")
 
         res$ances <- listres$ances[,toKeep]