[Pomp-commits] r506 - pkg/inst/doc

noreply at r-forge.r-project.org noreply at r-forge.r-project.org
Wed Jun 1 15:13:07 CEST 2011


Author: kingaa
Date: 2011-06-01 15:13:07 +0200 (Wed, 01 Jun 2011)
New Revision: 506

Modified:
   pkg/inst/doc/intro_to_pomp.Rnw
   pkg/inst/doc/intro_to_pomp.pdf
   pkg/inst/doc/nlf-multi-short.rda
Log:
- abbreviate one results file


Modified: pkg/inst/doc/intro_to_pomp.Rnw
===================================================================
--- pkg/inst/doc/intro_to_pomp.Rnw	2011-06-01 12:44:00 UTC (rev 505)
+++ pkg/inst/doc/intro_to_pomp.Rnw	2011-06-01 13:13:07 UTC (rev 506)
@@ -1122,7 +1122,14 @@
   pars <- my.params
   pars["log.r"] <- log.r.vals[j]
   for(k in 1:4) {
-    fvals[j,k] <- nlf(my.pomp, start=pars, nasymp=5000, est=NULL, lags=lags[[k]], eval.only=TRUE)
+    fvals[j,k] <- nlf(
+                      my.pomp,
+                      start=pars,
+                      nasymp=5000,
+                      est=NULL,
+                      lags=lags[[k]],
+                      eval.only=TRUE
+                      )
   }
 }
 @ 
@@ -1135,7 +1142,14 @@
   pars["log.K"] <- log.K.vals[j]
   pars["X.0"] <- exp(log.K.vals[j])
   for(k in 1:4) {
-    fvals2[j,k] <- nlf(my.pomp, start=pars, nasymp=5000, est=NULL, lags=lags[[k]], eval.only=TRUE)
+    fvals2[j,k] <- nlf(
+                       my.pomp,
+                       start=pars,
+                       nasymp=5000,
+                       est=NULL,
+                       lags=lags[[k]],
+                       eval.only=TRUE
+                       )
   }
 }
 @ 
@@ -1217,7 +1231,7 @@
   load(binary.file)
 } else {
 <<nlf-multi-short>>
-  save(lags,nreps,ndata,fvals,new.pomp,file=binary.file)
+  save(lags,nreps,ndata,fvals,file=binary.file)
 }
 @ 
 The last line above expresses the objective function as the geometric mean (quasi)likelihood per data point. 
@@ -1298,6 +1312,7 @@
 <<>>=
 apply(pars,2,sd)
 @ 
+[CORRECT THIS COMMENT, OR EXPUNGE?]  
 The bootstrap standard errors are, as expected, slightly higher than the asymptotic estimates. 
 
 The code above implements a ``resampling cases'' approach to bootstrapping the data set to which the intermediate autoregressive model is fitted. 

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

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



More information about the pomp-commits mailing list