Differential Metabolite analysis by limma
Examples
result <- mlimma(meta_dat, group)
result
#> # A tibble: 219 × 8
#> logFC AveExpr t P.Value adj.P.Val B logP name
#> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <chr>
#> 1 2.86 22.6 9.25 1.52e-10 0.0000000332 14.1 7.48 C02045
#> 2 2.44 26.2 7.80 6.83e- 9 0.000000748 10.4 6.13 C00267
#> 3 -1.82 27.1 -6.80 1.10e- 7 0.00000622 7.64 5.21 C00073
#> 4 -3.78 20.9 -6.79 1.14e- 7 0.00000622 7.61 5.21 C05674
#> 5 -2.20 21.4 -6.58 2.07e- 7 0.00000907 7.02 5.04 C00255
#> 6 -2.37 21.6 -6.45 2.98e- 7 0.0000109 6.66 4.96 C00242
#> 7 -3.00 22.1 -6.33 4.25e- 7 0.0000133 6.31 4.88 C02985
#> 8 -2.80 22.0 -6.28 4.91e- 7 0.0000134 6.17 4.87 C03413
#> 9 -3.11 26.7 -6.15 7.01e- 7 0.0000170 5.82 4.77 C00106
#> 10 -6.77 17.1 -6.06 9.17e- 7 0.0000191 5.55 4.72 C00010
#> # ℹ 209 more rows