Find initial cliques using mclust on the estimated correlation matrix

init_mclust(Sigma, r, n.noise = NULL)

Arguments

Sigma	A covariance matrix estimate (p x p).
r	Number of missing actors.
n.noise	Quantity of noise for mclust, set to 3 x p by default.

Value

A list of size r, with initial cliques of neighbors for each missing actor.

Examples

data=generate_missing_data(n=100,p=10,r=1,type="scale-free", plot=FALSE)
EMtree::draw_network(data$G,layout="nicely",curv=0,btw_rank=1,groupes=c(rep(1,10),2),
 nodes_size=c(5,5),nodes_label=1:11, pal_nodes= c("#adc9e0","#e7bd42"),
 pal_edges = "#31374f")$G
PLNfit<-norm_PLN(data$Y)
Sigma_hat=PLNfit$sigma_O
#-- original true clique
data$TC
#> [[1]]
#> [1] 1 2 3 4 8 9
#> 
#-- clique found by mclust
clique_mclust=init_mclust(Sigma_hat, r=1)
#> Warning: `data_frame()` is deprecated as of tibble 1.1.0.
#> Please use `tibble()` instead.
#> This warning is displayed once every 8 hours.
#> Call `lifecycle::last_warnings()` to see where this warning was generated.
clique_mclust
#> [[1]]
#> [1] 2 4 5 8
#>