diff --git a/multipartite_taxtable.R b/multipartite_taxtable.R
index b95655e..d0fc7a0 100644
--- a/multipartite_taxtable.R
+++ b/multipartite_taxtable.R
@@ -40,3 +40,11 @@ multipartite_fit <- estimateMultipartiteSBM(listSBM = lbm_list, estimOptions = l
 multipartite_covariates_fit <- estimateMultipartiteSBM(listSBM = lbm_list_covariates, estimOptions = list(initBM = FALSE))
 
 plot(multipartite_covariates_fit)
+
+
+
+lbm_list <- lapply(seq_along(per_taxa_networks), function(idx) {
+    defineSBM(per_taxa_networks[[idx]], model = "poisson", dimLabels = c(row = paste0("Rank", idx + 1), col = "Sample"))
+})
+
+multipartite_fit <- estimateMultipartiteSBM(listSBM = lbm_list, estimOptions = list(initBM = FALSE))