Sequential Monte Carlo with kernel embedded mappings: The mapping particle filter
- Autores
- Pulido, Manuel Arturo; Leeuwen, Peter Jan van
- Año de publicación
- 2019
- Idioma
- inglés
- Tipo de recurso
- artículo
- Estado
- versión publicada
- Descripción
- In this work, a novel sequential Monte Carlo filter is introduced which aims at an efficient sampling of the state space. Particles are pushed forward from the prediction to the posterior density using a sequence of mappings that minimizes the Kullback-Leibler divergence between the posterior and the sequence of intermediate densities. The sequence of mappings represents a gradient flow based on the principles of local optimal transport. A key ingredient of the mappings is that they are embedded in a reproducing kernel Hilbert space, which allows for a practical and efficient Monte Carlo algorithm. The kernel embedding provides a direct means to calculate the gradient of the Kullback-Leibler divergence leading to quick convergence using well-known gradient-based stochastic optimization algorithms. Evaluation of the method is conducted in the chaotic Lorenz-63 system, the Lorenz-96 system, which is a coarse prototype of atmospheric dynamics, and an epidemic model that describes cholera dynamics. No resampling is required in the mapping particle filter even for long recursive sequences. The number of effective particles remains close to the total number of particles in all the sequence. Hence, the mapping particle filter does not suffer from sample impoverishment.
Fil: Pulido, Manuel Arturo. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Nordeste. Instituto de Modelado e Innovación Tecnológica. Universidad Nacional del Nordeste. Facultad de Ciencias Exactas Naturales y Agrimensura. Instituto de Modelado e Innovación Tecnológica; Argentina. University of Reading; Reino Unido. Universidad Nacional del Nordeste. Facultad de Ciencias Exactas y Naturales y Agrimensura. Departamento de Física; Argentina
Fil: Leeuwen, Peter Jan van. University of Reading; Reino Unido - Materia
-
STEIN GRADIENT DESCENT
SEQUENTIAL BAYES
SWAM OPTIMIZATION
OPTIMAL TRANSPORT
KERNEL EMBEDDING
PARTICLE FLOWS - Nivel de accesibilidad
- acceso abierto
- Condiciones de uso
- https://creativecommons.org/licenses/by/2.5/ar/
- Repositorio
.jpg)
- Institución
- Consejo Nacional de Investigaciones Científicas y Técnicas
- OAI Identificador
- oai:ri.conicet.gov.ar:11336/105970
Ver los metadatos del registro completo
| id |
CONICETDig_9c1b0087c485192d1b6573ba9228f288 |
|---|---|
| oai_identifier_str |
oai:ri.conicet.gov.ar:11336/105970 |
| network_acronym_str |
CONICETDig |
| repository_id_str |
3498 |
| network_name_str |
CONICET Digital (CONICET) |
| spelling |
Sequential Monte Carlo with kernel embedded mappings: The mapping particle filterPulido, Manuel ArturoLeeuwen, Peter Jan vanSTEIN GRADIENT DESCENTSEQUENTIAL BAYESSWAM OPTIMIZATIONOPTIMAL TRANSPORTKERNEL EMBEDDINGPARTICLE FLOWShttps://purl.org/becyt/ford/1.2https://purl.org/becyt/ford/1In this work, a novel sequential Monte Carlo filter is introduced which aims at an efficient sampling of the state space. Particles are pushed forward from the prediction to the posterior density using a sequence of mappings that minimizes the Kullback-Leibler divergence between the posterior and the sequence of intermediate densities. The sequence of mappings represents a gradient flow based on the principles of local optimal transport. A key ingredient of the mappings is that they are embedded in a reproducing kernel Hilbert space, which allows for a practical and efficient Monte Carlo algorithm. The kernel embedding provides a direct means to calculate the gradient of the Kullback-Leibler divergence leading to quick convergence using well-known gradient-based stochastic optimization algorithms. Evaluation of the method is conducted in the chaotic Lorenz-63 system, the Lorenz-96 system, which is a coarse prototype of atmospheric dynamics, and an epidemic model that describes cholera dynamics. No resampling is required in the mapping particle filter even for long recursive sequences. The number of effective particles remains close to the total number of particles in all the sequence. Hence, the mapping particle filter does not suffer from sample impoverishment.Fil: Pulido, Manuel Arturo. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Nordeste. Instituto de Modelado e Innovación Tecnológica. Universidad Nacional del Nordeste. Facultad de Ciencias Exactas Naturales y Agrimensura. Instituto de Modelado e Innovación Tecnológica; Argentina. University of Reading; Reino Unido. Universidad Nacional del Nordeste. Facultad de Ciencias Exactas y Naturales y Agrimensura. Departamento de Física; ArgentinaFil: Leeuwen, Peter Jan van. University of Reading; Reino UnidoAcademic Press Inc Elsevier Science2019-05info:eu-repo/semantics/articleinfo:eu-repo/semantics/publishedVersionhttp://purl.org/coar/resource_type/c_6501info:ar-repo/semantics/articuloapplication/pdfapplication/pdfapplication/pdfhttp://hdl.handle.net/11336/105970Pulido, Manuel Arturo; Leeuwen, Peter Jan van; Sequential Monte Carlo with kernel embedded mappings: The mapping particle filter; Academic Press Inc Elsevier Science; Journal of Computational Physics; 396; 5-2019; 400-4150021-99911090-2716CONICET DigitalCONICETenginfo:eu-repo/semantics/altIdentifier/url/https://www.sciencedirect.com/science/article/pii/S0021999119304681info:eu-repo/semantics/altIdentifier/doi/10.1016/j.jcp.2019.06.060info:eu-repo/semantics/openAccesshttps://creativecommons.org/licenses/by/2.5/ar/reponame:CONICET Digital (CONICET)instname:Consejo Nacional de Investigaciones Científicas y Técnicas2026-04-08T11:38:24Zoai:ri.conicet.gov.ar:11336/105970instacron:CONICETInstitucionalhttp://ri.conicet.gov.ar/Organismo científico-tecnológicoNo correspondehttp://ri.conicet.gov.ar/oai/requestdasensio@conicet.gov.ar; lcarlino@conicet.gov.arArgentinaNo correspondeNo correspondeNo correspondeopendoar:34982026-04-08 11:38:24.892CONICET Digital (CONICET) - Consejo Nacional de Investigaciones Científicas y Técnicasfalse |
| dc.title.none.fl_str_mv |
Sequential Monte Carlo with kernel embedded mappings: The mapping particle filter |
| title |
Sequential Monte Carlo with kernel embedded mappings: The mapping particle filter |
| spellingShingle |
Sequential Monte Carlo with kernel embedded mappings: The mapping particle filter Pulido, Manuel Arturo STEIN GRADIENT DESCENT SEQUENTIAL BAYES SWAM OPTIMIZATION OPTIMAL TRANSPORT KERNEL EMBEDDING PARTICLE FLOWS |
| title_short |
Sequential Monte Carlo with kernel embedded mappings: The mapping particle filter |
| title_full |
Sequential Monte Carlo with kernel embedded mappings: The mapping particle filter |
| title_fullStr |
Sequential Monte Carlo with kernel embedded mappings: The mapping particle filter |
| title_full_unstemmed |
Sequential Monte Carlo with kernel embedded mappings: The mapping particle filter |
| title_sort |
Sequential Monte Carlo with kernel embedded mappings: The mapping particle filter |
| dc.creator.none.fl_str_mv |
Pulido, Manuel Arturo Leeuwen, Peter Jan van |
| author |
Pulido, Manuel Arturo |
| author_facet |
Pulido, Manuel Arturo Leeuwen, Peter Jan van |
| author_role |
author |
| author2 |
Leeuwen, Peter Jan van |
| author2_role |
author |
| dc.subject.none.fl_str_mv |
STEIN GRADIENT DESCENT SEQUENTIAL BAYES SWAM OPTIMIZATION OPTIMAL TRANSPORT KERNEL EMBEDDING PARTICLE FLOWS |
| topic |
STEIN GRADIENT DESCENT SEQUENTIAL BAYES SWAM OPTIMIZATION OPTIMAL TRANSPORT KERNEL EMBEDDING PARTICLE FLOWS |
| purl_subject.fl_str_mv |
https://purl.org/becyt/ford/1.2 https://purl.org/becyt/ford/1 |
| dc.description.none.fl_txt_mv |
In this work, a novel sequential Monte Carlo filter is introduced which aims at an efficient sampling of the state space. Particles are pushed forward from the prediction to the posterior density using a sequence of mappings that minimizes the Kullback-Leibler divergence between the posterior and the sequence of intermediate densities. The sequence of mappings represents a gradient flow based on the principles of local optimal transport. A key ingredient of the mappings is that they are embedded in a reproducing kernel Hilbert space, which allows for a practical and efficient Monte Carlo algorithm. The kernel embedding provides a direct means to calculate the gradient of the Kullback-Leibler divergence leading to quick convergence using well-known gradient-based stochastic optimization algorithms. Evaluation of the method is conducted in the chaotic Lorenz-63 system, the Lorenz-96 system, which is a coarse prototype of atmospheric dynamics, and an epidemic model that describes cholera dynamics. No resampling is required in the mapping particle filter even for long recursive sequences. The number of effective particles remains close to the total number of particles in all the sequence. Hence, the mapping particle filter does not suffer from sample impoverishment. Fil: Pulido, Manuel Arturo. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Nordeste. Instituto de Modelado e Innovación Tecnológica. Universidad Nacional del Nordeste. Facultad de Ciencias Exactas Naturales y Agrimensura. Instituto de Modelado e Innovación Tecnológica; Argentina. University of Reading; Reino Unido. Universidad Nacional del Nordeste. Facultad de Ciencias Exactas y Naturales y Agrimensura. Departamento de Física; Argentina Fil: Leeuwen, Peter Jan van. University of Reading; Reino Unido |
| description |
In this work, a novel sequential Monte Carlo filter is introduced which aims at an efficient sampling of the state space. Particles are pushed forward from the prediction to the posterior density using a sequence of mappings that minimizes the Kullback-Leibler divergence between the posterior and the sequence of intermediate densities. The sequence of mappings represents a gradient flow based on the principles of local optimal transport. A key ingredient of the mappings is that they are embedded in a reproducing kernel Hilbert space, which allows for a practical and efficient Monte Carlo algorithm. The kernel embedding provides a direct means to calculate the gradient of the Kullback-Leibler divergence leading to quick convergence using well-known gradient-based stochastic optimization algorithms. Evaluation of the method is conducted in the chaotic Lorenz-63 system, the Lorenz-96 system, which is a coarse prototype of atmospheric dynamics, and an epidemic model that describes cholera dynamics. No resampling is required in the mapping particle filter even for long recursive sequences. The number of effective particles remains close to the total number of particles in all the sequence. Hence, the mapping particle filter does not suffer from sample impoverishment. |
| publishDate |
2019 |
| dc.date.none.fl_str_mv |
2019-05 |
| dc.type.none.fl_str_mv |
info:eu-repo/semantics/article info:eu-repo/semantics/publishedVersion http://purl.org/coar/resource_type/c_6501 info:ar-repo/semantics/articulo |
| format |
article |
| status_str |
publishedVersion |
| dc.identifier.none.fl_str_mv |
http://hdl.handle.net/11336/105970 Pulido, Manuel Arturo; Leeuwen, Peter Jan van; Sequential Monte Carlo with kernel embedded mappings: The mapping particle filter; Academic Press Inc Elsevier Science; Journal of Computational Physics; 396; 5-2019; 400-415 0021-9991 1090-2716 CONICET Digital CONICET |
| url |
http://hdl.handle.net/11336/105970 |
| identifier_str_mv |
Pulido, Manuel Arturo; Leeuwen, Peter Jan van; Sequential Monte Carlo with kernel embedded mappings: The mapping particle filter; Academic Press Inc Elsevier Science; Journal of Computational Physics; 396; 5-2019; 400-415 0021-9991 1090-2716 CONICET Digital CONICET |
| dc.language.none.fl_str_mv |
eng |
| language |
eng |
| dc.relation.none.fl_str_mv |
info:eu-repo/semantics/altIdentifier/url/https://www.sciencedirect.com/science/article/pii/S0021999119304681 info:eu-repo/semantics/altIdentifier/doi/10.1016/j.jcp.2019.06.060 |
| dc.rights.none.fl_str_mv |
info:eu-repo/semantics/openAccess https://creativecommons.org/licenses/by/2.5/ar/ |
| eu_rights_str_mv |
openAccess |
| rights_invalid_str_mv |
https://creativecommons.org/licenses/by/2.5/ar/ |
| dc.format.none.fl_str_mv |
application/pdf application/pdf application/pdf |
| dc.publisher.none.fl_str_mv |
Academic Press Inc Elsevier Science |
| publisher.none.fl_str_mv |
Academic Press Inc Elsevier Science |
| dc.source.none.fl_str_mv |
reponame:CONICET Digital (CONICET) instname:Consejo Nacional de Investigaciones Científicas y Técnicas |
| reponame_str |
CONICET Digital (CONICET) |
| collection |
CONICET Digital (CONICET) |
| instname_str |
Consejo Nacional de Investigaciones Científicas y Técnicas |
| repository.name.fl_str_mv |
CONICET Digital (CONICET) - Consejo Nacional de Investigaciones Científicas y Técnicas |
| repository.mail.fl_str_mv |
dasensio@conicet.gov.ar; lcarlino@conicet.gov.ar |
| _version_ |
1861931235679404032 |
| score |
13.018236 |