construct_pearson function

Compute Pearson Residuals for ARD matrix and fitted model