rrr.fit: Fitting reduced-rank regression with a specific rank
Description
Given a response matrix and a covariate matrix, this function fits reduced
rank regression for a specified rank. It reduces to singular value
decomposition if the covariate matrix is the identity matrix.