thetaHat = 15/50 via BernBetaExampleJAW.R