diff --git a/multipartite_fungustree.R b/multipartite_fungustree.R
index 77079bb..3dd8c5f 100644
--- a/multipartite_fungustree.R
+++ b/multipartite_fungustree.R
@@ -1,4 +1,5 @@
-library(GREMLINS) # Load the custom GREMLINS
+# library(GREMLINS) # Load the custom GREMLINS
+devtools::load_all("../GREMLINS/")
 library(sbm)
 
 data(fungusTreeNetwork)
@@ -20,7 +21,10 @@ v_distrib <- c("gaussian", "bernoulli")
 namesFG <- c("Tree", "Fungus")
 givenclassif <- list(tree_geo_sbm$memberships, fungus_tree_sbm$memberships[["row"]])
 
-geo_multi_noinit <- multipartiteBM(list_Net = list_Net, v_distrib = v_distrib, namesFG = namesFG, givenclassif = NULL)
+geo_multi_noinit <- multipartiteBM(list_Net = list_Net, v_distrib = v_distrib, namesFG = namesFG, givenclassif = NULL, initBM = TRUE, verbose = TRUE)
+geo_multi_noinit_nobm <- multipartiteBM(list_Net = list_Net, v_distrib = v_distrib, namesFG = namesFG, givenclassif = NULL, initBM = FALSE, verbose = TRUE)
+geo_multi_noinit_nobm_kmin <- multipartiteBM(list_Net = list_Net, v_distrib = v_distrib, namesFG = namesFG, givenclassif = NULL, initBM = FALSE, verbose = TRUE, v_Kmin = c(7, 2))
+geo_multi_noinit_kmin <- multipartiteBM(list_Net = list_Net, v_distrib = v_distrib, namesFG = namesFG, givenclassif = NULL, initBM = TRUE, verbose = TRUE, v_Kmin = c(7, 2))
 geo_multi_init <- multipartiteBM(list_Net = list_Net, v_distrib = v_distrib, namesFG = namesFG, givenclassif = givenclassif)
 
 dataR6 <- formattingData(list_Net, v_distrib)