[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