Added the function of cvrcov() for covariate-by-covariate constrained randomization.
changed the function name cvcrand() to be cvrall() for constrained randomization.
Added the arguments check_validity to check validity by giving outputs of cluster pair coincidence matrix and descriptive statistics.