[adegenet-commits] r344 - pkg/R

noreply at r-forge.r-project.org noreply at r-forge.r-project.org
Wed Jun 3 00:51:36 CEST 2009


Author: jombart
Date: 2009-06-03 00:51:35 +0200 (Wed, 03 Jun 2009)
New Revision: 344

Modified:
   pkg/R/seqTrack.R
Log:
Slight tweaks on 'eaten from inside' ending.


Modified: pkg/R/seqTrack.R
===================================================================
--- pkg/R/seqTrack.R	2009-06-02 22:46:33 UTC (rev 343)
+++ pkg/R/seqTrack.R	2009-06-02 22:51:35 UTC (rev 344)
@@ -389,7 +389,7 @@
 ## 3) uncomment, adapt, and test code for missing data
 ##
 optimize.seqTrack <- function(nstep, seq.names, seq.dates, W, thres, optim=c("min","max"),
-                              prox.mat=NULL, mu0, seq.length, step.size=1000,
+                              prox.mat=NULL, mu0, seq.length, step.size=1e3,
                               rMissDate=.rUnifTimeSeq, ...){
 
     ## CHECKS ##
@@ -469,12 +469,12 @@
 
     ## DEFAULT CASE: NO MISSING DATES
     if(!any(isMissDate)){
-        ## >> one step of 'step.size simulations', all with same prior << ##
         ## dates initialisation, taken from initial prior
         newDates <- sapply(1:N, function(i)
                            .rTimeSeq(n=step.size, mu0=mu0[i], L=seq.length[i], maxNbDays=RANGE.DATES))
         newDates <- t(newDates)*24*3600 + seq.dates
 
+        ## >> one step of 'step.size simulations', all with same prior << ##
         for(i in 1:nstep){
             ## >> each step contains 'step.size' iterations << ##
             for(j in 1:step.size){



More information about the adegenet-commits mailing list