I’m not well versed in R but am using the GLLVM package to assess the influence of season and water mass on plankton distribution – when I try to make the coefficient plots using coeffplot.gllvm() I get a plot per factor level eg; season_factorSPRING or water_factorNeritic but I want a coeff plot for just the factors (season + water mass). Also, the one level from each factor is missing (I assume this has become a reference level the others are based off but is this correct?) – is there a way to do this? I have been following [https://jenniniku.github.io/gllvm/articles/vignette1.html#model-with-environmental-variables] as a guide and the only difference I can think of that is causing this issue is my data is categorical while the AntTraits dataset is numerical?