Optimizing relativistic energy density functionals: covariance analysis