scholarly journals Spatiotemporal blocking of the bouncy particle sampler for efficient inference in state-space models

2021 ◽  
Vol 31 (5) ◽  
Author(s):  
Jacob Vorstrup Goldman ◽  
Sumeetpal S. Singh

AbstractWe propose a novel blocked version of the continuous-time bouncy particle sampler of Bouchard-Côté et al. (J Am Stat Assoc 113(522):855–867, 2018) which is applicable to any differentiable probability density. This alternative implementation is motivated by blocked Gibbs sampling for state-space models (Singh et al. in Biometrika 104(4):953–969, 2017) and leads to significant improvement in terms of effective sample size per second, and furthermore, allows for significant parallelization of the resulting algorithm. The new algorithms are particularly efficient for latent state inference in high-dimensional state-space models, where blocking in both space and time is necessary to avoid degeneracy of MCMC. The efficiency of our blocked bouncy particle sampler, in comparison with both the standard implementation of the bouncy particle sampler and the particle Gibbs algorithm of Andrieu et al. (J R Stat Soc Ser B Stat Methodol 72(3):269–342, 2010), is illustrated numerically for both simulated data and a challenging real-world financial dataset.

2018 ◽  
Vol 51 (1) ◽  
pp. 530-535 ◽  
Author(s):  
Santhosh Kumar Varanasi ◽  
Phanindra Jampana

2018 ◽  
Vol 75 (5) ◽  
pp. 691-703 ◽  
Author(s):  
Timothy J. Miller ◽  
Saang-Yoon Hyun

State-space models explicitly separate uncertainty associated with unobserved, time-varying parameters from that which arises from sampling the population. The statistical aspects of formal state-space models are appealing and these models are becoming more widely used for assessments. However, treating natural mortality as known and constant across ages continues to be common practice. We developed a state-space, age-structured assessment model that allowed different assumptions for natural mortality and the degree of temporal stochasticity in abundance. We fit a suite of models where natural mortality was either age-invariant or an allometric function of mass and interannual transitions of abundance were deterministic or stochastic to observations on Gulf of Maine – Georges Bank Acadian redfish (Sebastes fasciatus). We found that allowing stochasticity in the interannual transition in abundance was important and estimating age-invariant natural mortality was sufficient. A simulation study showed low bias in annual biomass estimation when the estimation and simulation model matched and the Akaike imformation criterion accurately measured relative model performance, but it was important to allow simulated data sets to include the stochasticity in interannual transitions of abundance-at-age.


2010 ◽  
Vol 30 (2) ◽  
pp. 192-215 ◽  
Author(s):  
Alexander Shkolnik ◽  
Michael Levashov ◽  
Ian R. Manchester ◽  
Russ Tedrake

A motion planning algorithm is described for bounding over rough terrain with the LittleDog robot. Unlike walking gaits, bounding is highly dynamic and cannot be planned with quasi-steady approximations. LittleDog is modeled as a planar five-link system, with a 16-dimensional state space; computing a plan over rough terrain in this high-dimensional state space that respects the kinodynamic constraints due to underactuation and motor limits is extremely challenging. Rapidly Exploring Random Trees (RRTs) are known for fast kinematic path planning in high-dimensional configuration spaces in the presence of obstacles, but search efficiency degrades rapidly with the addition of challenging dynamics. A computationally tractable planner for bounding was developed by modifying the RRT algorithm by using: (1) motion primitives to reduce the dimensionality of the problem; (2) Reachability Guidance, which dynamically changes the sampling distribution and distance metric to address differential constraints and discontinuous motion primitive dynamics; and (3) sampling with a Voronoi bias in a lower-dimensional “task space” for bounding. Short trajectories were demonstrated to work on the robot, however open-loop bounding is inherently unstable. A feedback controller based on transverse linearization was implemented, and shown in simulation to stabilize perturbations in the presence of noise and time delays.


Sign in / Sign up

Export Citation Format

Share Document