The goal of GREMLINS is to perform statistical analysis of multipartite networks through a block model approach.
Multipartite networks consist in the joint observation of several networks implying some common individuals. The individuals (or entities represented by nodes) at stake are partitioned into groups defined by their nature. In what follows, these groups will be referred to as .
Installation
You can install the released version of GREMLINS GitHub with:
#devtools::install_github("Demiperimetre/GREMLINS") library(GREMLINS)
Mathematical Background
A collection of networks
Assume that Q functional groups of individuals are at stake; Let nq be the number of individuals in the q-th functional group.
A multipartite network is a collection of networks: each network may be simple (relations inside a functional group) or bipartite (relations between individuals of two functional groups). We index the collection of networks by pairs of functional groups (q, q′).
The set E denotes the list of pairs of functional groups for which we observe an interaction network.
For any pair (q, q′) ∈ E, the interaction network is encoded in a matrix Xqq′ : $X^{qq’}_{ii’} 0 $ if there is an edge from unit i of functional group q to unit i′ of functional group q′, 0 otherwise. - If q ≠ q′, Xqq′ is said to be an . - Xqq is an : it is symmetric if the relation inside the functional group q is non-oriented, non-symmetric otherwise.
A block model
Assume that, each functional group q is divided into Kq blocks (or equivalently clusters). ∀q ∈ {1, …, Q} and $ i {1,,n_q}$, let Ziq be the latent random variable such that Ziq = k if individual i of functional group q belongs to cluster k. The random variables Ziq’s are assumed to be independent and such that: ∀k ∈ {1, …, Kq}, ∀q ∈ {1, …, Q}, ∀i ∈ {1, …, nq}:
with $\sum_{k=1}^{K_q}\pi^{q}_k=1$, ∀q ∈ {1, …, Q}. Let $\bZ = \left(Z^{q}_i\right)_{i\in \{1,\ldots,n_q\}, q \in \{1,\ldots,Q\}}$ denote the set of latent variables.
