add canonical covariate scores to output