estimate rho