The cross section of options holds great promise for identifying return distributions and risk premia, but estimating dynamic option valuation models with latent state variables is challenging when using large option panels. We propose a particle Markov Chain Monte Carlo framework with a novel filtering approach and illustrate our method by estimating index option pricing models. Estimates of variance risk premiums, variance mean reversion, and higher moments differ from the literature. We show that these differences are due to the composition of the option sample. Restricting the option sample’s maturity dimension has the strongest impact on parameter inference and option fit in these models.

We propose no-arbitrage term structure models with volatility factors that follow GARCH processes. The models’ tractability is similar to canonical affine term structure models, but they fit yield volatility much better, especially for long-maturity yields. This improvement does not come at the expense of a deterioration in yield fit. Because of the improved volatility fit, the model performs substantially better in pricing Treasury futures options. We conclude that the specification of the volatility factors is critical. Modeling volatility as a function of (lagged) squared innovations to factors improves on models where volatility is a linear function of the factors.