Diff of /scripts/analysis_fly.py [000000] .. [7d5693]

Switch to unified view

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)