This helper function is called when doing Metabolites wide association analysis. It reports the results of linear regression models to study the association of a test variable to each metabolites individually and corrected for the covariates indicated.
do_metabowas(
phen,
dat,
test_variable = "age",
covariates = c("sex"),
adj_method = "BH",
quiet = TRUE
)results= the results of the MetaboWAS (estimate, tstatistics, pvalue, BH corrected pvalue)
phenotypes data.frame
metabolites data.frame
the variable to be investigated
the covariates that you want to add
correction method.
if FALSE it will plot the amount of people avaialble