[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