Iterated filtering
From Wikipedia, the free encyclopedia
Iterated filtering algorithms are a tool for maximum likelihood inference on partially observed dynamical systems. Stochastic perturbations to the unknown parameters are used to explore the parameter space. Applying sequential Monte Carlo (the particle filter) to this extended model results in the selection of the parameter values that are more consistent with the data. Appropriately constructed procedures, iterating with successively diminished perturbations, converge to the maximum likelihood estimate.[1][2][3] Iterated filtering methods have so far been used most extensively to study infectious disease transmission dynamics. Case studies include cholera,[4][5] Ebola virus,[6] influenza,[7][8][9][10] malaria,[11][12][13] HIV,[14] pertussis,[15][16] poliovirus[17] and measles.[5][18] Other areas which have been proposed to be suitable for these methods include ecological dynamics[19][20] and finance.[21][22]
The perturbations to the parameter space play several different roles. Firstly, they smooth out the likelihood surface, enabling the algorithm to overcome small-scale features of the likelihood during early stages of the global search. Secondly, Monte Carlo variation allows the search to escape from local minima. Thirdly, the iterated filtering update uses the perturbed parameter values to construct an approximation to the derivative of the log likelihood even though this quantity is not typically available in closed form. Fourthly, the parameter perturbations help to overcome numerical difficulties that can arise during sequential Monte Carlo.
Overview
The data are a time series collected at times . The dynamic system is modeled by a Markov process which is generated by a function in the sense that
where is a vector of unknown parameters and is some random quantity that is drawn independently each time is evaluated. An initial condition at some time is specified by an initialization function, . A measurement density completes the specification of a partially observed Markov process. We present a basic iterated filtering algorithm (IF1)[1][2] followed by an iterated filtering algorithm implementing an iterated, perturbed Bayes map (IF2).[3][23]
Procedure: Iterated filtering (IF1)
- Input: A partially observed Markov model specified as above; Monte Carlo sample size ; number of iterations ; cooling parameters and ; covariance matrix ; initial parameter vector
- for to
- draw for
- set for
- set
- for to
- draw for
- set for
- set for
- draw such that
- set and for
- set to the sample mean of , where the vector has components
- set to the sample variance of
- set
- Output: Maximum likelihood estimate
Variations
- For IF1, parameters which enter the model only in the specification of the initial condition, , warrant some special algorithmic attention since information about them in the data may be concentrated in a small part of the time series.[1]
- Theoretically, any distribution with the requisite mean and variance could be used in place of the normal distribution. It is standard to use the normal distribution and to reparameterise to remove constraints on the possible values of the parameters.
- Modifications to the IF1 algorithm have been proposed to give superior asymptotic performance.[24][25]
Procedure: Iterated filtering (IF2)
- Input: A partially observed Markov model specified as above; Monte Carlo sample size ; number of iterations ; cooling parameter ; covariance matrix ; initial parameter vectors
- for to
- set for
- set for
- for to
- draw for
- set for
- set for
- draw such that
- set and for
- set for
- Output: Parameter vectors approximating the maximum likelihood estimate,
Software
"pomp: statistical inference for partially-observed Markov processes" : R package.
References
Wikiwand - on
Seamless Wikipedia browsing. On steroids.