[adegenet-commits] r347 - pkg/R

noreply at r-forge.r-project.org noreply at r-forge.r-project.org
Wed Jun 3 01:46:35 CEST 2009


Author: jombart
Date: 2009-06-03 01:46:32 +0200 (Wed, 03 Jun 2009)
New Revision: 347

Modified:
   pkg/R/seqTrack.R
Log:
Some bugs remain, have to check that tomorrow.


Modified: pkg/R/seqTrack.R
===================================================================
--- pkg/R/seqTrack.R	2009-06-02 23:32:34 UTC (rev 346)
+++ pkg/R/seqTrack.R	2009-06-02 23:46:32 UTC (rev 347)
@@ -388,8 +388,9 @@
 ## 2) VECTORIZE mu0 and seq.length, recycle if needed with a warning
 ## 3) uncomment, adapt, and test code for missing data
 ##
-optimize.seqTrack <- function(nstep=10, seq.names, seq.dates, W, thres=NULL, optim=c("min","max"),
-                              prox.mat=NULL, mu0, seq.length, step.size=1e3,
+optimize.seqTrack <- function(nstep=10, step.size=1e3,
+                              seq.names, seq.dates, W, thres=NULL, optim=c("min","max"),
+                              prox.mat=NULL, mu0, seq.length,
                               rMissDate=.rUnifTimeSeq, ...){
 
     ## CHECKS ##
@@ -513,7 +514,7 @@
 
                 ## new dates
                 newDates <- apply(date, 1, function(vec)
-                                  sample(vec, size=step.size, replace=TRUE, prob=w)
+                                  sample(vec, size=step.size, replace=TRUE, prob=w))
                 newDates <- t(newDates)
 
                 ## re-initialize posterior distributions



More information about the adegenet-commits mailing list