How to “reconstruct” my observed data from factor scores and loadings using 'fa' function in psych (R)?

Question

I\'m running factor analysis with a \'varimax\' or \'promax\' rotation and covar=TRUE. I would like to "reconstruct" back my data by matrix multiplication of the l