[Pomp-commits] r92 - in pkg: . data inst inst/doc tests
noreply at r-forge.r-project.org
noreply at r-forge.r-project.org
Mon Apr 13 16:27:27 CEST 2009
Author: kingaa
Date: 2009-04-13 16:27:26 +0200 (Mon, 13 Apr 2009)
New Revision: 92
Modified:
pkg/DESCRIPTION
pkg/data/ou2.rda
pkg/inst/ChangeLog
pkg/inst/doc/compiled_code_in_pomp.pdf
pkg/inst/doc/intro_to_pomp.pdf
pkg/tests/ou2-mif.Rout.save
Log:
Modified: pkg/DESCRIPTION
===================================================================
--- pkg/DESCRIPTION 2009-04-13 12:21:14 UTC (rev 91)
+++ pkg/DESCRIPTION 2009-04-13 14:27:26 UTC (rev 92)
@@ -1,8 +1,8 @@
Package: pomp
Type: Package
Title: Statistical inference for partially observed Markov processes
-Version: 0.22-7
-Date: 2009-04-10
+Version: 0.22-8
+Date: 2009-04-12
Author: Aaron A. King, Edward L. Ionides, Carles Martinez Breto, Steve Ellner, Bruce Kendall
Maintainer: Aaron A. King <kingaa at umich.edu>
Description: Inference methods for partially-observed Markov processes
Modified: pkg/data/ou2.rda
===================================================================
(Binary files differ)
Modified: pkg/inst/ChangeLog
===================================================================
--- pkg/inst/ChangeLog 2009-04-13 12:21:14 UTC (rev 91)
+++ pkg/inst/ChangeLog 2009-04-13 14:27:26 UTC (rev 92)
@@ -1,5 +1,6 @@
2009-04-13 kingaa
+ * [r91] inst/ChangeLog:
* [r90] R/mif.R, man/mif.Rd, tests/ou2-mif.R: rewrite mif methods
to use 'mif.internal'.
this obviates the need for a user-visible '.ndone' argument.
Modified: pkg/inst/doc/compiled_code_in_pomp.pdf
===================================================================
(Binary files differ)
Modified: pkg/inst/doc/intro_to_pomp.pdf
===================================================================
(Binary files differ)
Modified: pkg/tests/ou2-mif.Rout.save
===================================================================
--- pkg/tests/ou2-mif.Rout.save 2009-04-13 12:21:14 UTC (rev 91)
+++ pkg/tests/ou2-mif.Rout.save 2009-04-13 14:27:26 UTC (rev 92)
@@ -69,7 +69,7 @@
> mif.fit <- continue(mif.fit,Nmif=70,max.fail=100)
> toc <- Sys.time()
> print(toc-tic)
-Time difference of 38.9845 secs
+Time difference of 1.054774 mins
> cat("PF estimated log likelihood at MIF MLE\n")
PF estimated log likelihood at MIF MLE
> print(pfilter(mif.fit)$loglik,digits=4)
@@ -123,4 +123,12 @@
> fit <- continue(fit,Nmif=2,ivps=c("x1.0"),pars=c("alpha.1"))
Warning message:
mif warning: the variable(s) alpha.4, x2.0 have positive random-walk SDs specified, but are included in neither ‘pars’ nor ‘ivps’. These random walk SDs are ignored.
+> s <- coef(fit)
+> s[2] <- 0.01
+> fit <- mif(fit,Nmif=10,start=s)
+> fit <- mif(ou2,Nmif=3,rw.sd=c(alpha.1=0.1,alpha.4=0.1),alg.pars=list(Np=1000,cooling.factor=0.98,var.factor=1,ic.lag=2))
+> fit <- continue(fit,Nmif=5,alg.pars=list(Np=2000,cooling.factor=0.98,var.factor=1,ic.lag=2))
+> fit <- continue(fit,ivps=c("x1.0"),rw.sd=c(alpha.1=0.1,alpha.4=0.1,x1.0=5,x2.0=5),Nmif=3)
+Warning message:
+mif warning: the variable(s) x2.0 have positive random-walk SDs specified, but are included in neither ‘pars’ nor ‘ivps’. These random walk SDs are ignored.
>
More information about the pomp-commits
mailing list