plot_heterogeneity_p function

Heterogeneity p-plot