Hello everyone,
I use currently the package PLNmodels for my metabarcoding data set,
But , Impossible to obtain the Rsquare when I use the function PLN and I didn't find more information for add this parameters.
(See to the link the slide 23) or the script exemple just below.
Thank you for ur help 
library(PLNmodels)
packageVersion("PLNmodels")
data(mollusk)
str(mollusk, max.level = 1)
mollusk$Abundance %>%
dplyr::select(1:10) %>%
head() %>% DT::datatable(fillContainer = FALSE)
t(log(1 + mollusk$Abundance)) %>%
corrplot::corrplot(
is.corr = FALSE,
addgrid.col = NA,
tl.cex = .5,
cl.pos = "n"
)
mollusc <- prepare_data(
counts = mollusk$Abundance,
covariates = mollusk$Covariate,
{{ offset = "TSS" }} # look at other possible choices
)
M00_mollusc <- PLN(Abundance ~ 1, mollusc)
M00_mollusc
https://mia-paris.pages.mia.inra.fr/formation_abondance_reseau/tutoriels/PLN/mollusk.html#36