init_mclust: Find initial cliques using mclust on the estimated...

In Rmomal/VEMtree: Network inference from Species counTs with missing actORs

Description

Find initial cliques using mclust on the estimated correlation matrix

Usage

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
#-- clique found by mclust
clique_mclust=init_mclust(Sigma_hat, r=1)
clique_mclust

Rmomal/VEMtree documentation built on Dec. 17, 2020, 5:17 p.m.