Skip to main content

Table 2 Type I error of GAMLSS-BEZI and LMAS on real microbiome data

From: metamicrobiomeR: an R package for analysis of microbiome relative abundance data using zero-inflated beta GAMLSS and meta-analysis across studies using random effects models

Taxonomic level

GAMLSS-BEZI

LMAS

GAMLSS-BEZI

LMAS

GAMLSS-BEZI

LMAS

 

Alpha level = 0.01 (median (IQR))

Alpha level = 0.05 (median (IQR))

Alpha level = 0.1 (median (IQR))

Cross-sectional microbiome data

 Phylum (5 taxa)

0.010 (0.007, 0.017)

0.007 (0.003, 0.010)

0.043 (0.043, 0.050)

0.040 (0.033, 0.043)

0.100 (0.093, 0.113)

0.090 (0.073, 0.090)

 Family (33 taxa)

0.000 (0.000, 0.003)

0.000 (0.000, 0.007)

0.007 (0.000, 0.043)

0.033 (0.007, 0.050)

0.070 (0.003, 0.103)

0.083 (0.053, 0.107)

Longitudinal microbiome data

 Phylum (5 taxa)

0.007 (0.002, 0.012)

0.010 (0.008, 0.013)

0.047 (0.030, 0.060)

0.067 (0.063, 0.080)

0.110 (0.075, 0.123)

0.117 (0.113, 0.132)

 Family (33 taxa)

0.003 (0.000, 0.008)

0.010 (0.007, 0.013)

0.043 (0.036, 0.053)

0.050 (0.043, 0.064)

0.097 (0.082, 0.110)

0.107 (0.089, 0.117)

  1. GAMLSS-BEZI Generalized Additive Models for Location, Scale and Shape (GAMLSS) with a zero inflated beta (BEZI) family, LMAS linear model with arcsin square root transformation (implemented in the software MaAsLin); IQR interquartile range. For longitudinal data, subject random intercepts were added to the models