Description Usage Arguments Details Value Filtering failures Author(s) References See Also
Modified version of the Liu and West (2001) algorithm.
1 2 3 4 5 6 7 8 | ## S4 method for signature 'data.frame'
bsmc2(data, Np, smooth = 0.1, tol = 1e-17,
max.fail = 0, params, rprior, rinit, rprocess, dmeasure, partrans, ...,
verbose = getOption("verbose", FALSE))
## S4 method for signature 'pomp'
bsmc2(data, Np, smooth = 0.1, tol = 1e-17,
max.fail = 0, ..., verbose = getOption("verbose", FALSE))
|
data |
either a data frame holding the time series data, or an object of class ‘pomp’, i.e., the output of another pomp calculation. |
Np |
number of particles |
smooth |
Kernel density smoothing parameter.
The compensating shrinkage factor will be |
tol |
positive numeric scalar;
particles with likelihood less than |
max.fail |
integer; the maximum number of filtering failures allowed (see below).
If the number of filtering failures exceeds this number, execution will terminate with an error.
By default, |
params |
optional; named numeric vector of parameters.
This will be coerced internally to storage mode |
rprior |
optional; prior distribution sampler, specified either as a C snippet, an R function, or the name of a pre-compiled native routine available in a dynamically loaded library.
For more information, see here.
Setting |
rinit |
simulator of the initial-state distribution.
This can be furnished either as a C snippet, an R function, or the name of a pre-compiled native routine available in a dynamically loaded library.
Setting |
rprocess |
simulator of the latent state process, specified using one of the rprocess plugins.
Setting |
dmeasure |
evaluator of the measurement model density, specified either as a C snippet, an R function, or the name of a pre-compiled native routine available in a dynamically loaded library.
Setting |
partrans |
optional parameter transformations, constructed using Many algorithms for parameter estimation search an unconstrained space of parameters.
When working with such an algorithm and a model for which the parameters are constrained, it can be useful to transform parameters.
One should supply the |
... |
additional arguments supply new or modify existing model characteristics or components.
See When named arguments not recognized by |
verbose |
logical; if |
bsmc2
uses a version of the original algorithm (Liu \& West 2001), but discards the auxiliary particle filter.
The modification appears to give superior performance for the same amount of effort.
Samples from the prior distribution are drawn using the rprior
component.
This is allowed to depend on elements of params
, i.e., some of the elements of params
can be treated as “hyperparameters”.
Np
draws are made from the prior distribution.
An object of class ‘bsmcd_pomp’. The following methods are avaiable:
plot
produces diagnostic plots
as.data.frame
puts the prior and posterior samples into a data frame
If the degree of disagreement between model and data becomes sufficiently large, a “filtering failure” results.
A filtering failure occurs when, at some time point, none of the Np
particles is compatible with the data.
In particular, if the conditional likelihood of a particle at any time is below the tolerance value tol
, then that particle is considered to be uninformative and its likelihood is taken to be zero.
A filtering failure occurs when this is the case for all particles.
A warning is generated when this occurs unless the cumulative number of failures exceeds max.fail
, in which case an error is generated.
Michael Lavine, Matthew Ferrari, Aaron A. King, Edward L. Ionides
Liu, J. and M. West. Combining Parameter and State Estimation in Simulation-Based Filtering. In A. Doucet, N. de Freitas, and N. J. Gordon, editors, Sequential Monte Carlo Methods in Practice, pages 197-224. Springer, New York, 2001.
Other particle filter methods: cond.logLik
,
eff.sample.size
, filter.mean
,
filter.traj
, mif2
,
pfilter
, pmcmc
,
pred.mean
, pred.var
Other pomp parameter estimation methods: abc
,
kalman
, mif2
,
nlf
, pmcmc
,
pomp2-package
, probe.match
,
spect.match
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.