[Pomp-commits] r784 - branches/mif2/inst/doc

noreply at r-forge.r-project.org noreply at r-forge.r-project.org
Tue Sep 25 22:23:42 CEST 2012


Author: nxdao2000
Date: 2012-09-25 22:23:42 +0200 (Tue, 25 Sep 2012)
New Revision: 784

Modified:
   branches/mif2/inst/doc/intro_to_pomp.Rnw
Log:
Need to change back the comment line

Modified: branches/mif2/inst/doc/intro_to_pomp.Rnw
===================================================================
--- branches/mif2/inst/doc/intro_to_pomp.Rnw	2012-09-25 20:22:49 UTC (rev 783)
+++ branches/mif2/inst/doc/intro_to_pomp.Rnw	2012-09-25 20:23:42 UTC (rev 784)
@@ -997,22 +997,22 @@
 <<>>=
 pf.truth <- pfilter(ricker,Np=1000,max.fail=50,seed=1066L)
 pf.guess <- pfilter(ricker,params=guess,Np=1000,max.fail=50,seed=1066L)
-pf.mf <- pfilter(mf,Np=1000,seed=1066L)
-pf.pm <- pfilter(pm,Np=1000,max.fail=10,seed=1066L)
+#pf.mf <- pfilter(mf,Np=1000,seed=1066L)
+#pf.pm <- pfilter(pm,Np=1000,max.fail=10,seed=1066L)
 pb.mf <- probe(mf,nsim=1000,probes=plist,seed=1066L)
 res <- rbind(
              cbind(guess=guess,truth=coef(ricker),MLE=coef(mf),PM=coef(pm)),
              loglik=c(
                pf.guess$loglik,
-               pf.truth$loglik,
-               pf.mf$loglik,
-               pf.pm$loglik
+               pf.truth$loglik
+               #pf.mf$loglik
+               #pf.pm$loglik
                ),
              synth.loglik=c(
                summary(pb.guess)$synth.loglik,
                summary(pb.truth)$synth.loglik,
-               summary(pb.mf)$synth.loglik,
-               summary(pm)$synth.loglik
+               summary(pb.mf)$synth.loglik
+             #  summary(pm)$synth.loglik
                )
              )
 



More information about the pomp-commits mailing list