compute_models_with_limma function

Compute lipid-specific regression models