model1.py now outputs data to plot predictive density

plot_predicted_alphas.R also handles a single distribution correctly