Calibrating the smile with multivariate time-changed brownian motion and the esscher transform

In this study, we investigate two multivariate time-changed Brownian motion option pricing models in which the connection between the historical measure P and the risk-neutral measure Q is given by the Esscher transform. The models incorporate skewness, kurtosis and more complex dependence structures among stocks log-returns than the simple correlation matrix. The two models can be seen as a multivariate extension of the normal inverse Gaussian (NIG) model and the variance gamma (VG) model, respectively. We discuss two possible approaches to estimate historical asset returns and calibrate univariate option implied volatilities. While the first approach considers only time series of log-returns, the second approach makes use of both time series of log-returns and univariate observed volatility surfaces. To calibrate the models, there is no need of liquid multivariate derivative quotes.