We combine two models: S&P total returns vs volatility only or vs volatility and earnings yield. See the data and code in https://github.com/asarantsev/earnings-yield-annual-simulator
This research continues blog posts on modeling returns vs earnings yield and volatility and returns vs only volatility. But previous Python code combines fitting of this model and simulation. Here, we removed plots (quantile-quantile and autocorrelation function for residuals) and analysis of skewness, kurtosis, etc. We only keep fitting regression code necessary for actual simulations. Also, we combine these two options (with or without earnings yield) so a person can make a choice. We also present a choice between simulating innovations using Gaussian multivariate distribution or kernel density estimation.
We also change the models slightly: Previously, we normalized returns and earnings growth
by dividing them by annual volatility
Here, we instead model these as regression
Future research will include Shiller CAPE (with various averaging windows), the new valuation measure (also adapted for various averaging windows), and bond spreads.

Leave a comment