Questions tagged [particle-filter]

Particle filters (or sequential Monte Carlo) is a form of genetic simulation algorithm used for filtering problems in signal analysis and time series analysis.

Filter by
Sorted by
Tagged with
1
vote
0answers
42 views

How does Particle Filters work?

I'm trying to figure out how particle filter works. Assume that I have selected propability function called $a \sim Gauss(\mu, \sigma)$. We call it proposial (Gaussian) Distribution. Then we have ...
0
votes
0answers
22 views

Do I need to know the distribution of the noise before I'm using Monte Carlo Sampling?

I'm going to use Particle Filter, which is a Monte Carlo Sampling. My simple question is: Do I need to know the distribution of the noise before I'm using Monte Carlo Sampling? Or can I just use a ...
1
vote
0answers
24 views

Computing Posterior via Particle Filter

I have a question regarding the computation of the posterior using particle filter. I start reporting details from https://cse.sc.edu/~terejanu/files/tutorialMC.pdf Consider random processes $(Q_t)_{t ...
2
votes
0answers
24 views

Do Particle Filters actually approximate the posterior distribution?

Im reading a tutorial paper about particle filters (Link) in which it is stated that as the number of samples tends to infinity the approximated posterior density given by $p(x_k|z_{1:k}) \approx \...
0
votes
0answers
17 views

How to handle additional measurement uncertainty from non-fixed sensor with Kalman filter?

Most examples I've seen of Kalman filters assume a measurement w.r.t. a global reference frame; i.e. direct measurement of location, and therefore direct measurement of a state. Consider the scenario ...
0
votes
0answers
8 views

Seeking more information on this Bayesian inference method

I came across an algorithm for performing Bayesian inference, see section 3.2 of this paper. Their approach is outlined as (where the "belief state" is a distribution on a discrete set of ...
0
votes
0answers
14 views

How do particle filters diversify the population?

I am trying to understand how to compute Bayesian filter with particle filter algorithms. My background is more from MCMC, importance sampling, nested sampling. I find the notations a bit confusing. ...
1
vote
0answers
24 views

Particle Gibbs Sampler For Regime-Switching Nonlinear Gaussian SSM

I'm reading this paper on using a non-linear Gaussian SSM for measuring regime-switching leverage effect using stock market data. I'm using it as jump-off point for an undergraduate paper. My advisor ...
1
vote
0answers
32 views

Particle filter/SMC - dynamic rotation in ICA (independent component analysis)

I struggle with the applicability of the bootstrap particle filter within dynamic rotations in independent component analysis. To be clear, suppose the following: $$Y_{t} = R(\delta_{t})\epsilon_{t}$$ ...
1
vote
1answer
47 views

Particle filter for likelihood evaluation

I struggle to implement a particle filter to evaluate the likelihood of a textbook example. I got the following process: $x_t = \alpha + \beta x_{t-1}/(1+x_{t-1}^2) + w_t$ where $w_t \sim \mathcal{N}(...
2
votes
0answers
90 views

Carré du champ operator is a quadratic variation

Let $X_t$ be a real valued Markov process (starting at $x$) with generator $L$. Let $\Gamma(f)$ denote Carré du champ operator i.e. $L(f^2) - 2f \cdot L (f)$. As far as I know under suitable ...
0
votes
0answers
53 views

How to update particles' weights and what are some ensemble evaluation metrics for a multivariate particle filter?

I am trying to figure out how to perform a multivariate (standard/bootstrap) particle filter to assimilate observations into a hydrological model. Most of the literature (in my field at least) is ...
0
votes
1answer
447 views

Sample from a distribution and plot in python

I am trying to understand Particle Filter and Importance Sampling Principle from a UniFreiburg Course and this USNA document on particle filters. Simultaneously, I am also trying to write a document ...
0
votes
1answer
18 views

Particle Filter for navigation through known map

I have some issues with understanding the Particle Filter for navigation through a known map. So, consider a situation where I want to write a Particle filter to navigate through a maze or a map that ...
0
votes
0answers
25 views

Deriving Hyperparameter updates in Online Interactive Collaborative Filtering

I've been going through "Online Interactive Collaborative Filtering Using Multi-Armed Bandit with Dependent Arms" by Wang et al. and am unable to understand how the update equations for the ...
2
votes
0answers
37 views

Particle Filtering: Derivation that mean of weights is the marginal likelihood

I see everywhere the following (for the Bootstrap Filter) $$ p(y_t \mid y_{1:t-1}) \approx \frac{1}{N} \sum_{i=1}^N W(x_{0:t}^i) $$ where $W(x_{0:t}^i)$ are the normalized weights defined as $$W(x_{...
0
votes
0answers
33 views

Monte Carlo probability approximation vs Histogram

I am trying to learn the sequential Monte Carlo method (particle filter) in data assimilation. In this method, the aim is to approximate the CDF of the target variable having a random sample of the ...
2
votes
0answers
50 views

Particle Filter for structural credit risk model

Kwon (2012)* proposes a structural credit risk model where the asset value process and the noise are estimated based on the observed equity prices: $S$ - equity prices $V$ - value of the assets $Z$ - ...
1
vote
0answers
34 views

Modelling Ball Movement using Delayed Measurements with Known Latency

I am a hobby programmer currently developing an algorithm to combine measurements of a dynamically moving ball position (and velocity) from multiple robots. Each robot measure and calculate the ...
0
votes
0answers
33 views

Simplying Bayes Theorem expression: SIS particle filter posteriori

In the book Beyond the Kalman Filter: Particle Filters for Tracking Applications on page 39 the weight update equation for the particle filter is derived. The derivations begins by introducing the ...
1
vote
1answer
51 views

Calculating the observation density

In the context of link. For a state space model $$ x_{k+1} = f(x_{k}, u_k, w_k) $$ $$ y_k = H x_k + v_k $$ where the measurement function is assumed linear and Gaussian and the state transition is ...
3
votes
1answer
142 views

Particle filter: Evaluating Optimal importance density

NOTE I posted this in the math stack exchange but I realized this may be the more appropriate place, old post here. I'm not sure if I should delete one of them so I just linked them in both? I am ...
0
votes
0answers
96 views

Evaluating Likelihood in Bootstrap Particle Filter

I am currently struggling with an attempt to apply a bootstrap particle filter to a linear, Gaussian state-space model $$s_t=A\,s_{t-1}+B\,\nu_t\qquad\text{( transition equation )}$$ $$\qquad z_t=C\,...
0
votes
1answer
216 views

Bootstrap Particle Filter (Gordon, Salmond, Smith, 2003) - Importance Weights

So, my endeavor to apply the is just for my own edificationI am currently struggling with an attempt to apply a bootstrap particle filter (Gordon, Salmond, Smith, 2003) to a linear, Gaussian state-...
1
vote
1answer
28 views

extension of sequential probability ratio tests to particle filters?

I've been wondering if there are extensions of the sequential probability ratio test to account for particle filters. I ask because, in my research, I'm working with distributions that cannot be ...
0
votes
1answer
113 views

Finite grid approximation to the Bayesian filtering problem

I need some hints for solving Ecercise 4.4 from Bayesian Filtering & Smoothing by Simo Särkkä: Select a finite interval in the state space, say, $x \in [-10, 10]$ and discretize it evenly to N ...
1
vote
1answer
54 views

Predictions after SMC

I have a statistical model given by $$ y_t\sim p(y_t|x_t, \theta)\\ x_t\sim p(x_t|x_{t-1},\theta)\\ \theta\sim p(\theta) $$ where $y$ is the only observed component. Using a sequential Monte Carlo ...
0
votes
1answer
58 views

What are the differences between Bayesian filters and adaptive filters?

I am learning about state estimation and I am having difficulty understanding the difference between Bayesian filters such as Kalman filter and particle filters compared to adaptive filters. According ...
0
votes
0answers
45 views

Particle filter for diagnosis

I have two annual measurements taken on medical images depicting a lung cancer tumor 's condition. I have likelihood function that taken in the measurement values and estimates malignancy of the tumor....
1
vote
0answers
110 views

Posterior as prior for correlated parameters [closed]

I want to use the posterior distribution of the model parameters $\theta$ given data in the time frame $[0,t]$ days, $P(\theta|y_{0:t})$; as a prior for the parameters in the time frame $[t+1, t+n]$ ...
2
votes
1answer
182 views

sequential Monte Carlo sampler, why the extended space and backward kernel?

Hello cross validated, I am currently studying sequential Monte Carlo samplers. My current understanding is as follows: We are interested in the marginal distribution of some sequence of joint ...
2
votes
1answer
104 views

Kernel for MCMC moves in sequential monte carlo

I'm trying to understand how to employ MCMC moves in a sequential Monte Carlo procedure for estimating static parameters as in the setting described by Chopin. He proposes, for example, the usage of a ...
2
votes
0answers
177 views

boostrap particle filter marginal likelihood

I want to calculate the marginal likelihood $p(y|\Theta)$ of the parameters of a Markov state space model with unknown parameters $\Theta$ that I am trying to estimate the marginal likelihood (...
3
votes
1answer
445 views

Why is it necessary to perform resampling step in particle filtering (or sequential monte carlo)?

I read the Wikipedia page on particle filter, it says that during 'prediction-updating', the samples from the distribution are weighted by a likelihood that represents the probability of that particle ...
1
vote
0answers
30 views

Inference for Maximum Likelihood Estimator Using Particle Filter

How does one compute standard errors for the MLE when using a particle filter approximation to the likelihood? I know that the estimator is asymptotically normal and that the variance-covariance ...
1
vote
1answer
132 views

Particle filter - expectations

I've recently been implementing some particle filter algorithms and I've realized there is a small detail I might have been doing incorrectly. Unfortunately the descriptions of the algorithms in ...
0
votes
1answer
194 views

Sequential monte carlo, resampling

In particle filters when one is doing sequential importance sampling, the quantity of interest that is being approximated is usually a weighted sum: $$\hat x_t = \sum_{i=1}^M \Bigl [f(v^{(i)}_{t}) \...
2
votes
1answer
132 views

Are the Sequential Monte Carlo algorithm invariant to the step at which we resample?

In a usual textual description (according to SMC in Practice book ) of a SMC algorithm for State-Space models, we usually expand the particles according to the distribution from the transition ...
2
votes
2answers
118 views

A doubt on the formula for updating the weights in Sequential Importance Sampling in a State-Space model

Let $x_{0:t}^{(i)}$ be the states from time $0$ to $t$ from sample $i$. Similarly for the observations $y_{1:t}$. The normalized weights are updated according to Where does the term $p(y_t|x_t^{(i)})...
1
vote
0answers
34 views

Bayes filter with delayed measurements

I have some straight and curve pieces with numbers, they are used to build tracks (of $5$ lanes) for my cars (figure $1$), I can send commands to the cars using an SDK on the Raspberry (set the speed ...
3
votes
1answer
109 views

Particle filter maximum likelihood with a discrete (Bernoulli) state variable, non-smooth loglikelihood

My model looks like this \begin{align} \begin{split} dY_{t} & = \sigma_{t} dW_{t} + Z_t dN_t \\ d\lambda_t & = \alpha(\lambda_\infty - \lambda_t)dt + \beta dN_t \end{split} \end{...
3
votes
1answer
204 views

Importance weight of conditioned particle in conditional SMC

In a generic particle filter, I understand the importance weights for each particle are calculated as $w_t^s \propto w_{t-1}^s \frac{p(y_t \mid z_t^s) p(z_t^s \mid z_{t-1}^s)}{q(z_t^s \mid z_{t-1}^s, ...
4
votes
2answers
1k views

Understanding Sequential Importance Sampling and Particle Filtering

I am struggling with SIS for particle filtering in the following aspect: In particle filtering (as per this book), the objective is to estimate the full posterior $p( x_{0:k} \mid y_{1:k} )$ rather ...
0
votes
1answer
285 views

notation conditional normal distribution [duplicate]

I'm describing parameter search using a particle filter, for which I use West M. (1993) Approximating Posterior Distributions by Mixture. On page 8 of the document, he states "and $p(\theta)$ is (...
3
votes
0answers
60 views

Soft Question: What background do I need to understand Feynmann Kac Formulae by Pierre Del Moral?

I am attempting to understand Sequential Monte Carlo(SMC) deeply, but with little theoretical background on probability theory and stochastic processes. Usually, the 'statistics' perspective of markov ...
1
vote
1answer
944 views

How to calculate importance weights for update step of an SIR (Sequential Importance Resampling) Particle filter?

I understand that one may use a particle filter to solve the filtering problem (estimating the hidden state of a system which can be described as a Hidden Markov Model). If I have a system where I ...
4
votes
2answers
547 views

Resample-move algorithm for Sequential Monte Carlo

I'm reading about the resample-move strategy, originally by Gilks and Berzuini, but my question will use the slightly more verbose description from the review of Doucet and Johansen, section 4.4, PDF ...
1
vote
0answers
18 views

Stochastic models: freedom in choosing error type?

I have a general question concerning probabilistic models. Imagine that I have a system with a high-dimensional state $x$ with a pdf $p(x)$. The state $x$ is time-dependent, and a propagating ...
5
votes
1answer
572 views

Stochastic volatility: particle filter vs Metropolis-Hastings

In many of the papers on particle filter I've read (e.g. Douc, Moulines and Olsson, 2007), stochastic volatility is a common example to show that a newly-proposed filter is working. At the same time, ...
0
votes
1answer
110 views

Multi-Target Tracking Filters

I am trying to solve a multi-target tracking problem, which is in some parts different to some filters I have already researched such as the PHD filter. I am asking for advise which filters to start ...