/cobage

Constraint-Based Model Analysis of Gene Essenciality

Primary LanguagePython

cobage

Authors Agustin Maria Pardo
Github AgustinPardo
Email agustinmpardo@gmail.com

Function inputs

Usage

essen_test(model_tb, dic_return, dataset_name, dataset_excel, growth_thresh_mult)

FN, FP, TN, TP = essen_test(model_iEK_griffin, "Yes both", "griffin", griffin_excel, grow_thresh)

model_tb : COBRA model

dic_return : "Yes" or "Yes both" options select which dictionaries return.

dataset_name : Type of dataset. "griffin" or "loerger"

dataset_excel : dataset file

growth_thresh_mult : FBA growth threshold

Griffin

Inputs:

Models

model_griffin : iEK1011_griffinEssen_media.json

Essenciality set

griffin_file : ppat.1002251.s002.xlsx

#griffin essenciality threshold grif_thres = 0.1

Loerger

Inputs:

Models

model_loerger : iEK1011_deJesusEssen_media.json

Essenciality set

loerger_file : mbo002173137st3.xlsx

#loerger essenciality threshold

ES being near 0

NE being near the mean

GD approximately 1/10 the mean

GA 5 times the mean

Output

Count of FN, FP, TN, TP

Experimental data Griffin
Model predictions iEK1011_griffinEssen_media.json True (Growth) False (No Growth)
Positive (Growth) 579 161
Negative (No Growth) 29 229
Total 998

Accuracy (Exactitud):

Matthew Correlation Coefficient (MCC):

Dictionaries of FN, FP, TN, TP