-
Notifications
You must be signed in to change notification settings - Fork 8
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
implement expression-limited fluxes (E-Flux like algorihm) #777
Conversation
PS let this be the final feature of the 1.x branch. I'll release 1.5, the rest of the 1.5 milestone will be applied as bugfixes to 1.5.x or better to 2.x directly. |
Codecov ReportPatch coverage:
Additional details and impacted files@@ Coverage Diff @@
## master #777 +/- ##
==========================================
+ Coverage 91.34% 91.39% +0.05%
==========================================
Files 82 86 +4
Lines 2196 2209 +13
==========================================
+ Hits 2006 2019 +13
Misses 190 190
☔ View full report in Codecov by Sentry. |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Nice, 🚀
E-Flux algorithm is closer described by: *Colijn, Caroline, Aaron Brandes, | ||
Jeremy Zucker, Desmond S. Lun, Brian Weiner, Maha R. Farhat, Tan-Yun Cheng, D. | ||
Branch Moody, Megan Murray, and James E. Galagan. "Interpreting expression data | ||
with metabolic flux models: predicting Mycobacterium tuberculosis mycolic acid | ||
production." PLoS computational biology 5, no. 8 (2009): e1000489*. |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
People love making questionable cobra methods :/
Create E-Flux-like bounds for a reaction that requires gene products given by | ||
`grr`, with expression "choke" data given by relative probabilities (between 0 | ||
and 1) in `relative_expression`. |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
interesting
isnothing(grr) && return 1.0 | ||
lup(g) = get(relative_expression, g, 1.0) | ||
1 - prod(1 - prod(lup.(gr)) for gr in grr) |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Ah, I thought this would have some stochasticity here
Closes #552
I'll make a separate port to the 2.0 branch later.