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

noreply at r-forge.r-project.org noreply at r-forge.r-project.org
Tue Jan 8 12:06:53 CET 2013


Author: kingaa
Date: 2013-01-08 12:06:53 +0100 (Tue, 08 Jan 2013)
New Revision: 795

Modified:
   branches/mif2/tests/aou2-mif-fp.R
   branches/mif2/tests/ou2-mif2-mif.R
Log:
- reformat


Modified: branches/mif2/tests/aou2-mif-fp.R
===================================================================
--- branches/mif2/tests/aou2-mif-fp.R	2013-01-08 10:57:25 UTC (rev 794)
+++ branches/mif2/tests/aou2-mif-fp.R	2013-01-08 11:06:53 UTC (rev 795)
@@ -12,38 +12,38 @@
 
 set.seed(64857673L)
 mif1 <- mif(ou2,Nmif=100,start=guess1,
-		pars=c('alpha.2','alpha.3'),ivps=c('x1.0','x2.0'),
-		rw.sd=c(
-				x1.0=.5,x2.0=.5,
-				alpha.2=0.1,alpha.3=0.1),
-		transform=F,
-		Np=1000,
-		var.factor=1,
-		ic.lag=10,
-		cooling.factor=0.95,
-		cooling.fraction=0.05,
-		option="mif2",
-		.ndone=0,tol=10^-8,
-		verbose=F
-)
-#set.seed(64857673L)
+            pars=c('alpha.2','alpha.3'),ivps=c('x1.0','x2.0'),
+            rw.sd=c(
+              x1.0=.5,x2.0=.5,
+              alpha.2=0.1,alpha.3=0.1),
+            transform=F,
+            Np=1000,
+            var.factor=1,
+            ic.lag=10,
+            cooling.factor=0.95,
+            cooling.fraction=0.05,
+            option="mif2",
+            .ndone=0,tol=10^-8,
+            verbose=F
+            )
+##set.seed(64857673L)
 mif2 <- mif(ou2,Nmif=100,start=guess1,
-		pars=c('alpha.2','alpha.3'),ivps=c('x1.0','x2.0'),
-		rw.sd=c(
-				x1.0=0.5,x2.0=.5,
-				alpha.2=0.1,alpha.3=0.1),
-		transform=F,
-		cooling.scalar=430,
-		Np=1000,
-		var.factor=1,
-		ic.lag=10,
-		cooling.factor=0.95,
-		cooling.fraction=0.5,
-		max.fail=100,
-		option="mif",
-		.ndone=0,tol=10^-8,
-		verbose=F
-)  
+            pars=c('alpha.2','alpha.3'),ivps=c('x1.0','x2.0'),
+            rw.sd=c(
+              x1.0=0.5,x2.0=.5,
+              alpha.2=0.1,alpha.3=0.1),
+            transform=F,
+            cooling.scalar=430,
+            Np=1000,
+            var.factor=1,
+            ic.lag=10,
+            cooling.factor=0.95,
+            cooling.fraction=0.5,
+            max.fail=100,
+            option="mif",
+            .ndone=0,tol=10^-8,
+            verbose=F
+            )  
 compare.mif(list(mif1,mif2))
 
 dev.off()

Modified: branches/mif2/tests/ou2-mif2-mif.R
===================================================================
--- branches/mif2/tests/ou2-mif2-mif.R	2013-01-08 10:57:25 UTC (rev 794)
+++ branches/mif2/tests/ou2-mif2-mif.R	2013-01-08 11:06:53 UTC (rev 795)
@@ -12,39 +12,49 @@
 
 set.seed(64857673L)
 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,
-				alpha.2=0.1,alpha.3=0.1),
-		transform=F,
-		Np=1000,
-		var.factor=1,
-		ic.lag=10,
-		cooling.factor=0.95,
-		cooling.fraction=0.05,
-		option="mif2",
-		.ndone=0
+            pars=c('alpha.2','alpha.3'),ivps=c('x1.0','x2.0'),
+            rw.sd=c(
+              x1.0=.5,x2.0=.5,
+              alpha.2=0.1,alpha.3=0.1),
+            transform=F,
+            Np=1000,
+            var.factor=1,
+            ic.lag=10,
+            cooling.factor=0.95,
+            cooling.fraction=0.05,
+            option="mif2",
+            .ndone=0
 
-)
+            )
 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(
-				x1.0=0.5,x2.0=.5,
-				alpha.2=0.1,alpha.3=0.1),
-		transform=F,
-		Np=1000,
-		var.factor=1,
-		ic.lag=10,
-		cooling.factor=0.95,
-		max.fail=100,
-		option="mif",
-		.ndone=0
+            pars=c('alpha.2','alpha.3'),ivps=c('x1.0','x2.0'),
+            rw.sd=c(
+              x1.0=0.5,x2.0=.5,
+              alpha.2=0.1,alpha.3=0.1),
+            transform=F,
+            Np=1000,
+            var.factor=1,
+            ic.lag=10,
+            cooling.factor=0.95,
+            max.fail=100,
+            option="mif",
+            .ndone=0
+            )  
 
-)  
-
 compare.mif(list(mif1,mif2))
 
 dev.off()
 
+
+
+
+
+
+
+
+
+
+
+



More information about the pomp-commits mailing list