Bayesian parameter inference for partially observed stochastic volterra equations

In this article we consider Bayesian parameter inference for a type of partially observed stochastic Volterra equation (SVE). SVEs are found in many areas such as physics and mathematical finance. In the latter field they can be used to represent long memory in unobserved volatility processes. In many cases of practical interest, SVEs must be time-discretized and then parameter inference is based upon the posterior associated to this time-discretized process. Based upon recent studies on time-discretization of SVEs (e.g. Richard et al. in Stoch Proc Appl 141:109–138, 2021) we use Euler–Maruyama methods for the afore-mentioned discretization. We then show how multilevel Markov chain Monte Carlo (MCMC) methods (Jasra et al. in SIAM J Sci Comp 40:A887–A902, 2018) can be applied in this context. In the examples we study, we give a proof that shows that the cost to achieve a mean square error (MSE) of \(\mathcal {O}(\epsilon ^2)\) $\mathcal{O}({ϵ}^2)$ , \(\epsilon >0\) $ϵ>0$ , is \(\mathcal {O}(\epsilon ^{-\tfrac{4}{2H+1}})\) $\mathcal{O}({ϵ}^{-\frac{4}{2H+1}})$ , where H is the Hurst parameter. If one uses a single level MCMC method then the cost is \(\mathcal {O}(\epsilon ^{-\tfrac{2(2H+3)}{2H+1}})\) $\mathcal{O}({ϵ}^{-\frac{2(2H+3)}{2H+1}})$ to achieve the same MSE. We illustrate these results in the context of state-space and stochastic volatility models, with the latter applied to real data.