I have followed the songbird tutorial (GitHub - biocore/songbird: Vanilla regression methods for microbiome differential abundance analysis) and generated the attached summarized results. Please, see below.
I ran 3 models for this analysis: 1st) A null model, 2nd) a 2nd model containing few important variables, and a 3rd) full model containing all the important variables. It looks like my full model was the best, judging by the lowest cv_error, lowest loss, and highest Pseudo-Q-square compared to the 2nd model (incomplete model) vs Null.
However, I have few questions regarding the interpretation of the qbeta results.
What is the interpretation of the first line of qbeta plots? Can you please explain in simple terms the y and x axis? I think the x-axis are the “interactions” but I didn’t quite understand the meaning of that in the tutorial.
How about the interpretation of the second line of qbeta plots containing the qbeta distributions? Are the numbers on the y-axis the interactions? What about the numbers in front of the distribution? The shape of the Full-model distribution seems normally distributed while the Null model or the incomplete model have several peaks. What does that tell you?
Thank you so much for your time and support!