added script for plotting frequentist estimatse