We need to use a function like pooled_cov from matchit to get the right covariance values. Build on sj514-fast-matches-greedy