approx_prod_stats function

Approximate summary statistics for a product of phenotypes and a set of predictors