README.md 1.58 KB
Newer Older
Nicole Bussola's avatar
Nicole Bussola committed
1
### Integrative Network Fusion (INF)
Nicole Bussola's avatar
Nicole Bussola committed
2
![INF pipeline ](figs/INF_pipeline.jpeg)
3

Alessia Marcolini's avatar
Alessia Marcolini committed
4
5
**Setup**
```bash
Alessia Marcolini's avatar
Alessia Marcolini committed
6
7
git clone https://gitlab.fbk.eu/MPBA/INF
cd INF
Alessia Marcolini's avatar
Alessia Marcolini committed
8
9
10
11
12
13
14
15
16
conda env create -f env.yml -n inf
conda activate inf
```

To install the R dependencies (not in conda channels), run the following command via the R prompt:
```bash
install.packages("TunePareto")
```

Alessia Marcolini's avatar
Alessia Marcolini committed
17
To install `mlpy`, follow the instructions [here](https://gitlab.fbk.eu/MPBA/mlpy). 
18

Alessia Marcolini's avatar
Alessia Marcolini committed
19

20
21
22
23
24
To install `bootstrapped`:
```bash
pip install bootstrapped
```

25
26
27
28
29
30
31
32
33
**Input files**

* omics layer 1 data: samples x features, tab-separated, with row & column names
* omics layer 2 data: same as above (**samples must be in the same order as the first file**)
* omics layers 1+2 data: the juxtaposition of the above two files
* labels file: one column, just the labels, no header (**same order as the data files**)

**Example run**

Marco Chierici's avatar
Marco Chierici committed
34
35
36
37
38
39
40
41
The INF pipeline is implemented as a Snakefile.

The following directory tree is required:

* {datafolder}/{dataset}/{layer1}_{layer2}_{tr,ts}.txt
* {datafolder}/{dataset}/labels_{target}_{tr,ts}.txt
* {datafolder}/{dataset}/{layer1,layer2}_{tr,ts}.txt
* {outfolder}/{dataset}/{target}/{juxt,rSNF,rSNFi,single}/ _(these will be created if not present)_
42

Marco Chierici's avatar
Marco Chierici committed
43
All the {variables} can be specified either in a config.yaml file or on the command line; for example:
44

Marco Chierici's avatar
Marco Chierici committed
45
46
```{python}
snakemake --config datafolder="data" dataset="breast" target="ER" layer1="gene" layer2="cnv"
47
```
Marco Chierici's avatar
Marco Chierici committed
48
49
50
51
52

A maximum number of cores can also be set:

```{python}
snakemake [--config etc.] --cores 12
53
54
```

Marco Chierici's avatar
Marco Chierici committed
55
The pipeline can be "dry-run" using the `-n` flag:
56

Marco Chierici's avatar
Marco Chierici committed
57
58
59
```{python}
snakemake --cores 12 -n
```