Bresnahan-Reiss (1991) -- minor fix to data generation

* now I use np.sum(P>0, axis=1) instead of a cryptic searchsorted + apply_along_axis line
Merge remote-tracking branch 'origin/master'