|
a |
|
b/scripts/analysis_fly.py |
|
|
1 |
import os, sys |
|
|
2 |
|
|
|
3 |
import pandas as pd |
|
|
4 |
from loguru import logger |
|
|
5 |
|
|
|
6 |
sys.path.append('..') |
|
|
7 |
|
|
|
8 |
from pyMultiOmics.constants import * |
|
|
9 |
from pyMultiOmics.mapping import Mapper |
|
|
10 |
from pyMultiOmics.common import set_log_level_info |
|
|
11 |
from pyMultiOmics.analysis import * |
|
|
12 |
from pyMultiOmics.query import * |
|
|
13 |
from pyMultiOmics.pipelines import * |
|
|
14 |
|
|
|
15 |
|
|
|
16 |
DATA_FOLDER = os.path.abspath(os.path.join('..', 'notebooks', 'test_data', 'fly_data')) |
|
|
17 |
gene_data = pd.read_csv(os.path.join(DATA_FOLDER, 'flyatlas_data.csv'), index_col='Identifier') |
|
|
18 |
gene_design = pd.read_csv(os.path.join(DATA_FOLDER, 'flyatlas_design.csv'), index_col='sample') |
|
|
19 |
|
|
|
20 |
compound_data = pd.read_csv(os.path.join(DATA_FOLDER, 'fly_metabolomics_no_dupes.csv'), index_col='Identifier') |
|
|
21 |
compound_design = pd.read_csv(os.path.join(DATA_FOLDER, 'fly_df_design.csv'), index_col='sample') |
|
|
22 |
|
|
|
23 |
set_log_level_info() |
|
|
24 |
|
|
|
25 |
m = Mapper(DROSOPHILA_MELANOGASTER, metabolic_pathway_only=True, include_related_chebi=True) \ |
|
|
26 |
.set_gene(gene_data, gene_design) \ |
|
|
27 |
.set_compound(compound_data, compound_design) \ |
|
|
28 |
.build() |
|
|
29 |
|
|
|
30 |
ap = AnalysisPipeline(m) |