* TST: Added test for the failing sampler

* BUG: Make Beta.random to use the clipped random number generator

* DOC: Added to release notes

* Fix the PR link in the release notes

* FIX: use scipy.stats and change ref_rand to point to clipped_beta_rvs

* FIX: Use np.maximum and np.minimum to work with scalars and arrays