Blockmodels, использующие r

Question

Я пытаюсь классифицировать некоторые биологические данные, используя библиотечные блочные модели.Я получил лучшую модель для моих данных (с 6 класса).Вопрос в том, как я могу получить метку каждого из данных (метки векторов).Это мой код (мат 132 * 132 массив):

my_model <- BM_poisson("SBM_sym",mat)
my_model$estimate()
which.max(my_model$ICL)
my_model$model_parameters[6]
my_model$plot_obs_pred(6)

Спасибо

Answer 1

Вот как я решил проблему:

library("network")
library("ramify")
mynet <- as.network(mat)
my_model <- BM_bernoulli("SBM",mat)
my_model$estimate()
#K ist most probable number of classes
K<-which.max(my_model$ICL)
#argmax of Z is the most probable classmemebership
classmembership <-argmax(my_model$memberships[[K]]$Z)
#set classmembership as attribute for plotting
network::set.vertex.attribute(mynet, "Class", classmembership)