cov_nest function

Compute covariance or correlation matrix with treatments for clusters and missing values