Fast approximate likelihood evaluation for stable VARFIMA processes

For VARFIMA models with sub-Gaussian stable errors, we present fast approximate likelihood computation by using a multivariate preconditioned conjugate gradient (MPCG) algorithm, and Monte Carlo integration over unobserved variables. We illustrate our approach on daily average temperatures measured at several US cities.