[Pomp-commits] r596 - in pkg: . inst/doc

noreply at r-forge.r-project.org noreply at r-forge.r-project.org
Mon Jan 16 14:27:23 CET 2012


Author: kingaa
Date: 2012-01-16 14:27:23 +0100 (Mon, 16 Jan 2012)
New Revision: 596

Modified:
   pkg/DESCRIPTION
   pkg/inst/doc/gompertz-multi-mif.rda
   pkg/inst/doc/intro_to_pomp.Rnw
   pkg/inst/doc/intro_to_pomp.pdf
Log:
- fix SD/SE of pfilter log-likelihoods in 'intro_to_pomp'


Modified: pkg/DESCRIPTION
===================================================================
--- pkg/DESCRIPTION	2012-01-15 19:45:14 UTC (rev 595)
+++ pkg/DESCRIPTION	2012-01-16 13:27:23 UTC (rev 596)
@@ -2,7 +2,7 @@
 Type: Package
 Title: Statistical inference for partially observed Markov processes
 Version: 0.40-3
-Date: 2012-01-15
+Date: 2012-01-16
 Author: Aaron A. King, Edward L. Ionides, Carles Breto, Steve Ellner, Bruce Kendall, Helen Wearing, Matthew J. Ferrari, Michael Lavine, Daniel C. Reuman
 Maintainer: Aaron A. King <kingaa at umich.edu>
 URL: http://pomp.r-forge.r-project.org

Modified: pkg/inst/doc/gompertz-multi-mif.rda
===================================================================
(Binary files differ)

Modified: pkg/inst/doc/intro_to_pomp.Rnw
===================================================================
--- pkg/inst/doc/intro_to_pomp.Rnw	2012-01-15 19:45:14 UTC (rev 595)
+++ pkg/inst/doc/intro_to_pomp.Rnw	2012-01-16 13:27:23 UTC (rev 596)
@@ -693,10 +693,10 @@
 loglik.mif <- replicate(n=10,logLik(pfilter(mf[[1]],params=theta.mif,Np=10000)))
 bl <- mean(loglik.mif)
 loglik.mif.est <- bl+log(mean(exp(loglik.mif-bl)))
-loglik.mif.se <- sd(exp(loglik.mif-bl))/sqrt(length(loglik.mif))/exp(loglik.mif.est-bl)
+loglik.mif.se <- sd(exp(loglik.mif-bl))/exp(loglik.mif.est-bl)
 loglik.true <- replicate(n=10,logLik(pfilter(gompertz,params=theta.true,Np=10000)))
 loglik.true.est <- bl+log(mean(exp(loglik.true-bl)))
-loglik.true.se <- sd(exp(loglik.true-bl))/sqrt(length(loglik.true))/exp(loglik.true.est-bl)
+loglik.true.se <- sd(exp(loglik.true-bl))/exp(loglik.true.est-bl)
 @ 
 
 <<gompertz-multi-mif-eval,echo=F,results=hide>>=

Modified: pkg/inst/doc/intro_to_pomp.pdf
===================================================================
(Binary files differ)



More information about the pomp-commits mailing list