[Pomp-commits] r329 - pkg/tests

[noreply at r-forge.r-project.org]

Author: kingaa
Date: 2010-09-26 19:06:12 +0200 (Sun, 26 Sep 2010)
New Revision: 329

Modified:
   pkg/tests/ricker-probe.Rout.save
Log:
- improve tests for probes


Modified: pkg/tests/ricker-probe.Rout.save
===================================================================
--- pkg/tests/ricker-probe.Rout.save	2010-09-26 17:01:29 UTC (rev 328)
+++ pkg/tests/ricker-probe.Rout.save	2010-09-26 17:06:12 UTC (rev 329)
@@ -27,6 +27,28 @@
 > po <- ricker
 > pb <- probe(
 +             po,
++             probes=probe.median("y"),
++             nsim=1000,
++             seed=838775L
++             )
+> plot(pb)
+> summary(pb)
+$coef
+    log.r log.sigma   log.phi       N.0       e.0 
+ 3.800000 -1.203973  2.302585  7.000000  0.000000 
+
+$nsim
+[1] 1000
+
+$quantiles
+[1] 0.784
+
+$pvals
+[1] 0.3376623
+
+> 
+> pb <- probe(
++             po,
 +             probes=probe.marginal(
 +               var="y",
 +               transform=sqrt,
@@ -150,7 +172,7 @@
 +                               )
 +             )
    user  system elapsed 
- 10.072   0.011  10.084 
+ 10.083   0.030  10.115 
 > plot(pm)
 > 
 > cbind(truth=coef(ricker),est=coef(pm),guess=coef(po))