abline(fit1, col='green')
abline(fit.comp, col='blue')
mtext(paste(' Rho =', round(rho, 5)), at=-1)
mtext(paste(' Rho2 =', round(rho2, 5)), at=6)
Impute missing data, while forcing correlation coefficient to remain the same