Errors in masaslin3 presence/absence results

Downstream analysis and statistics MaAsLin
1

Hello,

I am trying to use masaslin3 for the first time and when I run the example code from masaslin3 function, I have very strange things in the presence/absence results. I get the same error for all taxa: “non-integer number of successes in a binomial glm!” (I am not sure of the exact message since my output is in French). Thus, no qval_individual is computed and the presence/absence results are not being used. I have the same issue when running masaslin3 on other datasets with a simple grouping factor. Would you have an idea of where it could come from (maybe a specific package to update)? Below an example of the problematic output and my R session information.

image
image1920×609 177 KB

image
image1920×896 137 KB

Thank you for your help!
Aurélie

2

Hi Aurélie,

This is a funny problem! During the data augmentation step to improve logistic regression, the model uses small, non-integer weights for the augmented points. For historical reasons, the weights are assumed to be integers, but they work just fine even if they’re not integers. However, a warning is thrown, so we catch it (the English version) and get rid of it. However, because the error thrown is in French, the automatic catching doesn’t work. Two possible solutions:

  1. If it’s not that hard to switch R to English (or install the English version of MaAsLin 3 and lme4 if that’s a possibility), that should fix the issue.
  2. If that doesn’t work, you could just ignore those errors and manually calculate the q-values based on applying p.adjust to the p-values.

Will

3

Hi Will,

Indeed, this is a funny problem. Using Sys.setenv(LANGUAGE=‘en’) before running my code works. Thank you!

Aurélie