Merge branch 'develop' into fix-zfpkm

Author: JoshLoecker

main/como/rnaseq_preprocess.py

@@ -459,6 +459,7 @@ async def _create_config_df(  # noqa: C901
 
         fragment_label = f"{context_name}_{label}_fragment_size.txt"
         frag_paths = [p for p in aux_lookup["fragment"].values() if p.name == fragment_label]
+        print(f"HERE: {prep}")
         if not frag_paths and prep.lower() != RNAType.TRNA.value.lower():
             logger.warning(f"No fragment file for '{label}'; defaulting to 100 bp (needed for zFPKM).")
             mean_frag = 100.0

stubs/cobra/__init__.pyi

@@ -0,0 +1,5 @@
+from cobra import flux_analysis as flux_analysis, io as io, medium as medium, sampling as sampling, summary as summary
+from cobra.core import Configuration as Configuration, DictList as DictList, Gene as Gene, Metabolite as Metabolite, Model as Model, Object as Object, Reaction as Reaction, Solution as Solution, Species as Species
+from cobra.util import show_versions as show_versions
+
+__version__: str

stubs/cobra/core/__init__.pyi

@@ -0,0 +1,10 @@
+from cobra.core.configuration import Configuration as Configuration
+from cobra.core.dictlist import DictList as DictList
+from cobra.core.gene import GPR as GPR, Gene as Gene
+from cobra.core.group import Group as Group
+from cobra.core.metabolite import Metabolite as Metabolite
+from cobra.core.model import Model as Model
+from cobra.core.object import Object as Object
+from cobra.core.reaction import Reaction as Reaction
+from cobra.core.solution import Solution as Solution, get_solution as get_solution
+from cobra.core.species import Species as Species

stubs/cobra/core/configuration.pyi

@@ -0,0 +1,32 @@
+import pathlib
+import types
+from numbers import Number
+
+from _typeshed import Incomplete
+
+from .singleton import Singleton
+
+__all__ = ["Configuration"]
+
+
+class Configuration(metaclass=Singleton):
+    tolerance: float
+    lower_bound: Incomplete
+    upper_bound: Incomplete
+    processes: Incomplete
+    max_cache_size: Incomplete
+    cache_expiration: Incomplete
+
+    def __init__(self, **kwargs) -> None: ...
+    @property
+    def solver(self) -> types.ModuleType: ...
+    @solver.setter
+    def solver(self, value) -> None: ...
+    @property
+    def bounds(self) -> tuple[Number | None, Number | None]: ...
+    @bounds.setter
+    def bounds(self, bounds: tuple[Number | None, Number | None]) -> None: ...
+    @property
+    def cache_directory(self) -> pathlib.Path: ...
+    @cache_directory.setter
+    def cache_directory(self, path: pathlib.Path | str) -> None: ...

stubs/cobra/core/dictlist.pyi

@@ -0,0 +1,35 @@
+from .object import Object as Object
+from typing import Any, Callable, Iterable, Iterator, Pattern
+
+class DictList(list):
+    def __init__(self, *args) -> None: ...
+    def has_id(self, id: Object | str) -> bool: ...
+    def get_by_id(self, id: Object | str) -> Object: ...
+    def list_attr(self, attribute: str) -> list: ...
+    def get_by_any(self, iterable: list[str | Object | int]) -> list: ...
+    def query(self, search_function: str | Pattern | Callable, attribute: str | None = None) -> DictList: ...
+    def append(self, entity: Object) -> None: ...
+    def union(self, iterable: Iterable[Object]) -> None: ...
+    def extend(self, iterable: Iterable[Object]) -> None: ...
+    def __sub__(self, other: Iterable[Object]) -> DictList: ...
+    def __isub__(self, other: Iterable[Object]) -> DictList: ...
+    def __add__(self, other: Iterable[Object]) -> DictList: ...
+    def __iadd__(self, other: Iterable[Object]) -> DictList: ...
+    def __reduce__(self) -> tuple[type['DictList'], tuple, dict, Iterator]: ...
+    def index(self, id: str | Object, *args) -> int: ...
+    def __contains__(self, entity: str | Object) -> bool: ...
+    def __copy__(self) -> DictList: ...
+    def insert(self, index: int, entity: Object) -> None: ...
+    def pop(self, *args) -> Object: ...
+    def add(self, x: Object) -> None: ...
+    def remove(self, x: str | Object) -> None: ...
+    def reverse(self) -> None: ...
+    def sort(self, cmp: Callable = None, key: Callable = None, reverse: bool = False) -> None: ...
+    def __getitem__(self, i: int | slice | Iterable | Object | DictList) -> DictList | Object: ...
+    def __setitem__(self, i: slice | int, y: list | Object) -> None: ...
+    def __delitem__(self, index: int) -> None: ...
+    def __getslice__(self, i: int, j: int) -> DictList: ...
+    def __setslice__(self, i: int, j: int, y: list | Object) -> None: ...
+    def __delslice__(self, i: int, j: int) -> None: ...
+    def __getattr__(self, attr: Any) -> Any: ...
+    def __dir__(self) -> list: ...

stubs/cobra/core/formula.pyi

@@ -0,0 +1,15 @@
+from .object import Object as Object
+from _typeshed import Incomplete
+
+element_re: Incomplete
+
+class Formula(Object):
+    formula: Incomplete
+    elements: Incomplete
+    def __init__(self, formula: str | None = None, **kwargs) -> None: ...
+    def __add__(self, other_formula: Formula | str) -> Formula: ...
+    def parse_composition(self) -> None: ...
+    @property
+    def weight(self) -> float: ...
+
+elements_and_molecular_weights: Incomplete

stubs/cobra/core/gene.pyi

@@ -0,0 +1,58 @@
+import sympy.logic.boolalg as spl
+from ..util import resettable as resettable
+from ..util.util import format_long_string as format_long_string
+from .dictlist import DictList as DictList
+from .species import Species as Species
+from _typeshed import Incomplete
+from ast import AST, BoolOp, Expression, Module, Name, NodeTransformer, NodeVisitor
+from sympy import Symbol
+from typing import Iterable
+
+logger: Incomplete
+keywords: Incomplete
+keyword_re: Incomplete
+number_start_re: Incomplete
+replacements: Incomplete
+
+class GPRWalker(NodeVisitor):
+    gene_set: Incomplete
+    def __init__(self, **kwargs) -> None: ...
+    def visit_Name(self, node: Name) -> None: ...
+    def visit_BoolOp(self, node: BoolOp) -> None: ...
+
+class GPRCleaner(NodeTransformer):
+    gene_set: Incomplete
+    def __init__(self, **kwargs) -> None: ...
+    def visit_Name(self, node: Name) -> Name: ...
+    def visit_BinOp(self, node: BoolOp) -> None: ...
+
+def parse_gpr(str_expr: str) -> tuple: ...
+
+class Gene(Species):
+    def __init__(self, id: str = None, name: str = '', functional: bool = True) -> None: ...
+    @property
+    def functional(self) -> bool: ...
+    @functional.setter
+    @resettable
+    def functional(self, value: bool) -> None: ...
+    def knock_out(self) -> None: ...
+
+class GPR(Module):
+    body: list | None
+    def __init__(self, gpr_from: Expression | Module | AST = None, **kwargs) -> None: ...
+    @classmethod
+    def from_string(cls, string_gpr: str) -> GPR: ...
+    @property
+    def genes(self) -> frozenset: ...
+    def update_genes(self) -> None: ...
+    def eval(self, knockouts: DictList | set | str | Iterable = None) -> bool: ...
+    def to_string(self, names: dict = None) -> str: ...
+    def copy(self) -> GPR: ...
+    def __copy__(self) -> GPR: ...
+    def as_symbolic(self, names: dict = None) -> spl.Or | spl.And | Symbol: ...
+    @classmethod
+    def from_symbolic(cls, sympy_gpr: spl.BooleanFunction | Symbol) -> GPR: ...
+    def __eq__(self, other) -> bool: ...
+
+def eval_gpr(expr: Expression | GPR, knockouts: DictList | set) -> bool: ...
+def ast2str(expr: Expression | GPR, level: int = 0, names: dict = None) -> str: ...

stubs/cobra/core/group.pyi

@@ -0,0 +1,16 @@
+from .object import Object as Object
+from _typeshed import Incomplete
+from typing import Iterable
+
+class Group(Object):
+    KIND_TYPES: Incomplete
+    def __init__(self, id: str, name: str = '', members: Iterable | None = None, kind: str | None = None) -> None: ...
+    def __len__(self) -> int: ...
+    @property
+    def members(self) -> set: ...
+    @property
+    def kind(self) -> str: ...
+    @kind.setter
+    def kind(self, kind: str) -> None: ...
+    def add_members(self, new_members: list) -> None: ...
+    def remove_members(self, to_remove: list) -> None: ...

stubs/cobra/core/metabolite.pyi

@@ -0,0 +1,32 @@
+from ..exceptions import OptimizationError as OptimizationError
+from ..util.solver import check_solver_status as check_solver_status
+from ..util.util import format_long_string as format_long_string
+from .formula import elements_and_molecular_weights as elements_and_molecular_weights
+from .species import Species as Species
+from _typeshed import Incomplete
+from cobra.core import Solution as Solution
+from cobra.summary.metabolite_summary import MetaboliteSummary as MetaboliteSummary
+from optlang.interface import Container as Container
+from pandas import DataFrame as DataFrame
+
+element_re: Incomplete
+
+class Metabolite(Species):
+    formula: Incomplete
+    compartment: Incomplete
+    charge: Incomplete
+    def __init__(self, id: str | None = None, formula: str | None = None, name: str | None = '', charge: float | None = None, compartment: str | None = None) -> None: ...
+    @property
+    def constraint(self) -> Container: ...
+    @property
+    def elements(self) -> dict[str, int | float] | None: ...
+    @elements.setter
+    def elements(self, elements_dict: dict[str, int | float]) -> None: ...
+    @property
+    def formula_weight(self) -> int | float: ...
+    @property
+    def y(self) -> float: ...
+    @property
+    def shadow_price(self) -> float: ...
+    def remove_from_model(self, destructive: bool = False) -> None: ...
+    def summary(self, solution: Solution | None = None, fva: float | DataFrame | None = None) -> MetaboliteSummary: ...

stubs/cobra/core/model.pyi

@@ -0,0 +1,91 @@
+import optlang
+import pandas as pd
+import types
+from ..medium import find_boundary_types as find_boundary_types, find_external_compartment as find_external_compartment, sbo_terms as sbo_terms
+from ..util.context import HistoryManager as HistoryManager, get_context as get_context, resettable as resettable
+from ..util.solver import add_cons_vars_to_problem as add_cons_vars_to_problem, assert_optimal as assert_optimal, check_solver as check_solver, interface_to_str as interface_to_str, remove_cons_vars_from_problem as remove_cons_vars_from_problem, set_objective as set_objective
+from ..util.util import AutoVivification as AutoVivification, format_long_string as format_long_string
+from .configuration import Configuration as Configuration
+from .dictlist import DictList as DictList
+from .gene import Gene as Gene
+from .group import Group as Group
+from .metabolite import Metabolite as Metabolite
+from .object import Object as Object
+from .reaction import Reaction as Reaction
+from .solution import get_solution as get_solution
+from _typeshed import Incomplete
+from cobra import Solution as Solution
+from cobra.summary import ModelSummary as ModelSummary
+from cobra.util.solver import CONS_VARS as CONS_VARS
+from optlang.container import Container as Container
+from optlang.symbolics import Basic
+from types import ModuleType
+from typing import Iterable
+
+logger: Incomplete
+configuration: Incomplete
+
+class Model(Object):
+    genes: Incomplete
+    reactions: Incomplete
+    metabolites: Incomplete
+    groups: Incomplete
+    def __init__(self, id_or_model: str | Model | None = None, name: str | None = None) -> None: ...
+    @property
+    def solver(self) -> optlang.interface.Model: ...
+    @solver.setter
+    @resettable
+    def solver(self, value: str | ModuleType) -> None: ...
+    @property
+    def tolerance(self) -> float: ...
+    @tolerance.setter
+    def tolerance(self, value: float) -> None: ...
+    @property
+    def compartments(self) -> dict: ...
+    @compartments.setter
+    def compartments(self, value: dict) -> None: ...
+    @property
+    def medium(self) -> dict[str, float]: ...
+    @medium.setter
+    def medium(self, medium: dict[str, float]) -> None: ...
+    def copy(self) -> Model: ...
+    def add_metabolites(self, metabolite_list: list | Metabolite) -> None: ...
+    def remove_metabolites(self, metabolite_list: list | Metabolite, destructive: bool = False) -> None: ...
+    def add_boundary(self, metabolite: Metabolite, type: str = 'exchange', reaction_id: str | None = None, lb: float | None = None, ub: float | None = None, sbo_term: str | None = None) -> Reaction: ...
+    def add_reactions(self, reaction_list: Iterable[Reaction]) -> None: ...
+    def remove_reactions(self, reactions: str | Reaction | list[str | Reaction], remove_orphans: bool = False) -> None: ...
+    def add_groups(self, group_list: str | Group | list[Group]) -> None: ...
+    def remove_groups(self, group_list: str | Group | list[Group]) -> None: ...
+    def get_associated_groups(self, element: Reaction | Gene | Metabolite) -> list[Group]: ...
+    def add_cons_vars(self, what: list['CONS_VARS'] | tuple['CONS_VARS'], **kwargs) -> None: ...
+    def remove_cons_vars(self, what: list['CONS_VARS'] | tuple['CONS_VARS']) -> None: ...
+    @property
+    def problem(self) -> optlang.interface: ...
+    @property
+    def variables(self) -> Container: ...
+    @property
+    def constraints(self) -> Container: ...
+    @property
+    def boundary(self) -> list[Reaction]: ...
+    @property
+    def exchanges(self) -> list[Reaction]: ...
+    @property
+    def demands(self) -> list[Reaction]: ...
+    @property
+    def sinks(self) -> list[Reaction]: ...
+    def slim_optimize(self, error_value: float | None = ..., message: str | None = None) -> float: ...
+    def optimize(self, objective_sense: str | None = None, raise_error: bool = False) -> Solution: ...
+    def repair(self, rebuild_index: bool = True, rebuild_relationships: bool = True) -> None: ...
+    @property
+    def objective(self) -> optlang.Objective: ...
+    @objective.setter
+    def objective(self, value: dict[Reaction, 'Container'] | str | int | Reaction | optlang.interface.Objective | Basic) -> None: ...
+    @property
+    def objective_direction(self) -> str: ...
+    @objective_direction.setter
+    @resettable
+    def objective_direction(self, value: str) -> None: ...
+    def summary(self, solution: Solution | None = None, fva: pd.DataFrame | float | None = None) -> ModelSummary: ...
+    def __enter__(self) -> Model: ...
+    def __exit__(self, type: type[BaseException] | None, value: BaseException | None, traceback: types.TracebackType | None) -> None: ...
+    def merge(self, right: Model, prefix_existing: str | None = None, inplace: bool = True, objective: str = 'left') -> Model: ...