[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