[Pomp-commits] r773 - branches/mif2/tests

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

Author: nxdao2000
Date: 2012-08-13 18:44:22 +0200 (Mon, 13 Aug 2012)
New Revision: 773

Modified:
   branches/mif2/tests/ou2-mif2-mif.R
Log:
update with continue example

Modified: branches/mif2/tests/ou2-mif2-mif.R
===================================================================
--- branches/mif2/tests/ou2-mif2-mif.R	2012-08-13 16:34:16 UTC (rev 772)
+++ branches/mif2/tests/ou2-mif2-mif.R	2012-08-13 16:44:22 UTC (rev 773)
@@ -11,7 +11,7 @@
 guess2[c('x1.0','x2.0','alpha.2','alpha.3')] <- 1.2*guess1[c('x1.0','x2.0','alpha.2','alpha.3')]
 
 set.seed(64857673L)
-mif1 <- mif(ou2,Nmif=100,start=guess1,
+mif0 <- mif(ou2,Nmif=50,start=guess1,
 		pars=c('alpha.2','alpha.3'),ivps=c('x1.0','x2.0'),
 		rw.sd=c(
 				x1.0=.5,x2.0=.5,
@@ -19,14 +19,15 @@
 		transform=F,
 		Np=1000,
 		var.factor=1,
-		cooling.scalar=400,
 		ic.lag=10,
 		cooling.factor=0.95,
+		cooling.fraction=0.05,
 		option="mif2",
 		.ndone=0
-		
+
 )
-set.seed(64857673L)
+mif1<-continue(mif0,Nmif=50)
+
 mif2 <- mif(ou2,Nmif=100,start=guess1,
 		pars=c('alpha.2','alpha.3'),ivps=c('x1.0','x2.0'),
 		rw.sd=c(
@@ -40,8 +41,10 @@
 		max.fail=100,
 		option="mif",
 		.ndone=0
-		
+
 )  
+
 compare.mif(list(mif1,mif2))
 
 dev.off()
+