[adegenet-commits] r424 - pkg/R

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

Author: jombart
Date: 2009-06-28 15:17:47 +0200 (Sun, 28 Jun 2009)
New Revision: 424

Modified:
   pkg/R/seqTrack.R
Log:
Small fix to likelihood compute. Does not change values, but more elegant


Modified: pkg/R/seqTrack.R
===================================================================
--- pkg/R/seqTrack.R	2009-06-23 13:07:51 UTC (rev 423)
+++ pkg/R/seqTrack.R	2009-06-28 13:17:47 UTC (rev 424)
@@ -904,9 +904,9 @@
         nb.days <- abs(as.integer(anc.dates-dates))
         nb.mut <- x$weight
         mu <- mu0/365
-        mu <- mu*nb.days
+        ##mu <- mu*nb.days
 
-        res <- dbinom(nb.mut, size=seq.length, prob=mu)
+        res <- dbinom(nb.mut, size=seq.length*nb.days, prob=mu)
     } else{
         cat("Method not implemented.")
     }