(function from pyomo.contrib.parmest.graphics)
pyomo.contrib.parmest.graphics
Fit a Gaussian kernel-density distribution to theta values
theta_values (DataFrame) – Theta values, columns = variable names
scipy.stats.gaussian_kde distribution