def corr_equi(k_vars, rho):
 corr = np.empty((k_vars, k_vars))
 corr.fill(rho)
 corr[np.diag_indices_from(corr)] = 1
 return corr

def corr_equi ( k_vars , rho ) : corr = np . empty ( ( k_vars , k_vars ) ) corr . fill ( rho ) corr [ np . diag_indices_from ( corr ) ] = 1 return corr
create equicorrelated correlation matrix with rho on off diagonal parameters k_vars : int number of variables .