-
Notifications
You must be signed in to change notification settings - Fork 1
/
README.Rmd
100 lines (69 loc) · 4.51 KB
/
README.Rmd
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
<!-- README.md is generated from README.Rmd. Please edit that file -->
```{r setup, include = FALSE}
knitr::opts_chunk$set(
collapse = TRUE,
comment = "#>",
fig.path = "man/figures/README-",
out.width = "100%",
dpi=300,fig.width=7,
fig.keep="all"
)
```
# network ABC <img src="man/figures/logo.png" align="right" width="200"/>
----------------
# Network Reverse Engineering with Approximate Bayesian Computation
## Frédéric Bertrand and Myriam Maumy-Bertrand
<!-- badges: start -->
[](https://lifecycle.r-lib.org/articles/stages.html)
[](https://www.repostatus.org/#active)
[](https://github.com/fbertran/networkABC/actions)
[](https://app.codecov.io/gh/fbertran/networkABC?branch=master)
[](https://cran.r-project.org/package=networkABC)
[](https://cran.r-project.org/package=networkABC)
[](https://github.com/fbertran/networkABC)
[](https://zenodo.org/badge/latestdoi/169994404)
<!-- badges: end -->
The goal of `networkABC` is to provide **an inference tool based on approximate Bayesian computation** to decipher network data and assess the strength of their inferred links.
We developed an inference tool based on approximate Bayesian computation to decipher network data and assess the strength of the inferred links between network's actors.
It is a new multi-level approximate Bayesian computation (ABC) approach.
* At the first level, the method captures the global properties of the network, such as scale-freeness and clustering coefficients,
* whereas the second level is targeted to capture local properties, including the probability of each couple of genes being linked.
Up to now, Approximate Bayesian Computation (ABC) algorithms have been scarcely used in that setting and, due to the computational overhead, their application was limited to a small number of genes. On the contrary, our algorithm was made to cope with that issue and has low computational cost.
It can be used, for instance, for elucidating gene regulatory network, which is an important step towards understanding the normal cell physiology and complex pathological phenotype. Reverse-engineering consists in using gene expressions over time or over different experimental conditions to discover the structure of the gene network in a targeted cellular process. The fact that gene expression data are usually noisy, highly correlated, and have high dimensionality explains the need for specific statistical methods to reverse engineer the underlying network.
The package was accepted for presentation at the the useR! 2021 international conference. A technical note for the package was created and published on the website of the conference. It can be accessed here: [https://user2021.r-project.org/participation/technical_notes/t164/technote/](https://user2021.r-project.org/participation/technical_notes/t164/technote/).
This website and these examples were created by F. Bertrand and M. Maumy-Bertrand.
## Installation
You can install the released version of networkABC from [CRAN](https://CRAN.R-project.org) with:
```{r, eval = FALSE}
install.packages("networkABC")
```
You can install the development version of networkABC from [github](https://github.com) with:
```{r, eval = FALSE}
devtools::install_github("fbertran/networkABC")
```
## Examples
The `network_gen` allows to simulate networks with given clustering coefficient.
```{r networkgen}
set.seed(314)
library(networkABC)
resnet <- network_gen(10,1)
resnet
```
The `abc` function performs the network reverse engineering with Approximate Bayesian Computation.
```{r abc}
set.seed(314)
M10<-matrix(rnorm(30),10,3)
result<-abc(data=M10)
```
The `showHp` function plots the hub probabilities.
```{r showHp}
showHp(result)
```
The `showNp` function plots the neighbourhood probabilities.
```{r showNp}
showNp(result)
```
The `showNetwork` function plots the final network.
```{r showNetwork}
showNetwork(result,.2)
```