parmigene (version 1.1.1)

aracne.a: Algorithm for the Reconstruction of Accurate Cellular Networks

Description

A function that implements the ARACNE algorithm for the reconstruction of gene interaction networks (additive model).

Usage

aracne.a(mi, eps=0.05)

Value

A square weighted adjacency matrix of the inferred network.

Arguments

mi

matrix of the mutual information.

eps

a positive numeric value used to remove the weakest edge of each triple of nodes.

Details

This algorithm considers each triple of edges independently and removes the weakest one if

$$MI(i;j) < MI(j;k) - \varepsilon$$

and

$$MI(i;j) < MI(i;k) - \varepsilon$$

By default, the function uses all the available cores. You can set the actual number of threads used to N by exporting the environment variable OMP_NUM_THREADS=N.

References

Adam A. Margolin, Ilya Nemenman, Katia Basso, Chris Wiggins, Gustavo Stolovitzky, Riccardo Dalla Favera, and Andrea Califano. Aracne : An algorithm for the reconstruction of gene regulatory networks in a mammalian cellular context. BMC Bioinformatics, 2006.

See Also

aracne.m

clr

mrnet

Examples

mat <- matrix(rnorm(1000), nrow=10)
mi  <- knnmi.all(mat)
grn <- aracne.a(mi, 0.05)