[Pomp-commits] r20 - pkg/tests

noreply at r-forge.r-project.org noreply at r-forge.r-project.org
Fri Jul 25 16:04:42 CEST 2008


Author: kingaa
Date: 2008-07-25 16:04:42 +0200 (Fri, 25 Jul 2008)
New Revision: 20

Modified:
   pkg/tests/ou2-mif.Rout.save
Log:
no need to set rw.sd for parameters that aren't estimated

Modified: pkg/tests/ou2-mif.Rout.save
===================================================================
--- pkg/tests/ou2-mif.Rout.save	2008-07-25 14:03:17 UTC (rev 19)
+++ pkg/tests/ou2-mif.Rout.save	2008-07-25 14:04:42 UTC (rev 20)
@@ -1,6 +1,6 @@
 
-R version 2.6.1 (2007-11-26)
-Copyright (C) 2007 The R Foundation for Statistical Computing
+R version 2.7.1 (2008-06-23)
+Copyright (C) 2008 The R Foundation for Statistical Computing
 ISBN 3-900051-07-0
 
 R is free software and comes with ABSOLUTELY NO WARRANTY.
@@ -54,17 +54,20 @@
 +                pars=c('alpha.1','alpha.4'),ivps=c('x1.0','x2.0'),
 +                rw.sd=c(
 +                  x1.0=5,x2.0=5,
-+                  alpha.1=0.1,alpha.2=0,alpha.3=0,alpha.4=0.1,
-+                  sigma.1=0,sigma.2=0,sigma.3=0,
-+                  tau=0
++                  alpha.1=0.1,alpha.4=0.1
 +                  ),
-+                alg.pars=list(Np=1000,var.factor=1,ic.lag=10,cooling.factor=0.95),
++                alg.pars=list(
++                  Np=1000,
++                  var.factor=1,
++                  ic.lag=10,
++                  cooling.factor=0.95
++                  ),
 +                max.fail=100
 +                )
 > mif.fit <- continue(mif.fit,Nmif=70,max.fail=100)
 > toc <- Sys.time()
 > print(toc-tic)
-Time difference of 1.953113 mins
+Time difference of 1.542085 mins
 > cat("PF estimated log likelihood at MIF MLE\n")
 PF estimated log likelihood at MIF MLE
 > print(pfilter(mif.fit)$loglik,digits=4)



More information about the pomp-commits mailing list