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
CONICET Digital (CONICET)
Institución
Consejo Nacional de Investigaciones Científicas y Técnicas
OAI Identificador
oai:ri.conicet.gov.ar:11336/105970

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