diff --git a/multipartite_fungustree.R b/multipartite_fungustree.R index 3dd8c5f..45bcc37 100644 --- a/multipartite_fungustree.R +++ b/multipartite_fungustree.R @@ -23,8 +23,13 @@ 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, 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) +# Les kmins donnent le même résultat ICL -449.21 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_kmin <- multipartiteBM(list_Net = list_Net, v_distrib = v_distrib, namesFG = namesFG, givenclassif = givenclassif, v_Kmin = c(7, 2)) # Ne s'autorise pas à descendre +## + +# Le meilleur ICL est ici : -408.05 avec une classif de base issue de blockmodels indépendants geo_multi_init <- multipartiteBM(list_Net = list_Net, v_distrib = v_distrib, namesFG = namesFG, givenclassif = givenclassif) dataR6 <- formattingData(list_Net, v_distrib)