[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