add covariances computation on reduced.r