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...