phylogenie 2.0.1__py3-none-any.whl → 2.0.3__py3-none-any.whl

phylogenie/__init__.py

@@ -0,0 +1,72 @@
+from phylogenie.generators import (
+    AliSimDatasetGenerator,
+    BDEITreeDatasetGenerator,
+    BDSSTreeDatasetGenerator,
+    BDTreeDatasetGenerator,
+    CanonicalTreeDatasetGenerator,
+    DatasetGenerator,
+    DatasetGeneratorConfig,
+    EpidemiologicalTreeDatasetGenerator,
+    FBDTreeDatasetGenerator,
+    TreeDatasetGeneratorConfig,
+)
+from phylogenie.io import load_fasta, load_newick
+from phylogenie.msa import MSA
+from phylogenie.skyline import (
+    SkylineMatrix,
+    SkylineMatrixCoercible,
+    SkylineParameter,
+    SkylineParameterLike,
+    SkylineVector,
+    SkylineVectorCoercible,
+    SkylineVectorLike,
+    skyline_matrix,
+    skyline_parameter,
+    skyline_vector,
+)
+from phylogenie.tree import Tree
+from phylogenie.treesimulator import (
+    Event,
+    get_BD_events,
+    get_BDEI_events,
+    get_BDSS_events,
+    get_canonical_events,
+    get_epidemiological_events,
+    get_FBD_events,
+    simulate_tree,
+)
+
+__all__ = [
+    "AliSimDatasetGenerator",
+    "BDEITreeDatasetGenerator",
+    "BDSSTreeDatasetGenerator",
+    "BDTreeDatasetGenerator",
+    "CanonicalTreeDatasetGenerator",
+    "DatasetGenerator",
+    "DatasetGeneratorConfig",
+    "EpidemiologicalTreeDatasetGenerator",
+    "FBDTreeDatasetGenerator",
+    "SkylineMatrix",
+    "SkylineMatrixCoercible",
+    "skyline_matrix",
+    "SkylineParameter",
+    "SkylineParameterLike",
+    "skyline_parameter",
+    "SkylineVector",
+    "SkylineVectorCoercible",
+    "SkylineVectorLike",
+    "skyline_vector",
+    "Tree",
+    "TreeDatasetGeneratorConfig",
+    "Event",
+    "get_BD_events",
+    "get_BDEI_events",
+    "get_BDSS_events",
+    "get_canonical_events",
+    "get_epidemiological_events",
+    "get_FBD_events",
+    "simulate_tree",
+    "load_fasta",
+    "load_newick",
+    "MSA",
+]

phylogenie/generators/__init__.py

@@ -4,11 +4,29 @@ from pydantic import Field
 
 from phylogenie.generators.alisim import AliSimDatasetGenerator
 from phylogenie.generators.dataset import DatasetGenerator
-from phylogenie.generators.trees import TreeDatasetGeneratorConfig
+from phylogenie.generators.trees import (
+    BDEITreeDatasetGenerator,
+    BDSSTreeDatasetGenerator,
+    BDTreeDatasetGenerator,
+    CanonicalTreeDatasetGenerator,
+    EpidemiologicalTreeDatasetGenerator,
+    FBDTreeDatasetGenerator,
+    TreeDatasetGeneratorConfig,
+)
 
 DatasetGeneratorConfig = Annotated[
     TreeDatasetGeneratorConfig | AliSimDatasetGenerator,
     Field(discriminator="data_type"),
 ]
 
-__all__ = ["DatasetGeneratorConfig", "DatasetGenerator"]
+__all__ = [
+    "DatasetGeneratorConfig",
+    "DatasetGenerator",
+    "AliSimDatasetGenerator",
+    "CanonicalTreeDatasetGenerator",
+    "EpidemiologicalTreeDatasetGenerator",
+    "FBDTreeDatasetGenerator",
+    "BDTreeDatasetGenerator",
+    "BDEITreeDatasetGenerator",
+    "BDSSTreeDatasetGenerator",
+]

phylogenie/generators/configs.py

@@ -38,4 +38,6 @@ SkylineParameterConfig = ScalarConfig | SkylineParameterModel
 SkylineVectorConfig = (
     str | pgt.Scalar | list[SkylineParameterConfig] | SkylineVectorModel
 )
-SkylineMatrixConfig = str | pgt.Scalar | list[SkylineVectorConfig] | SkylineMatrixModel
+SkylineMatrixConfig = (
+    str | pgt.Scalar | list[SkylineVectorConfig] | SkylineMatrixModel | None
+)

phylogenie/generators/factories.py

@@ -148,7 +148,10 @@ def one_or_many_2D_scalars(
 
 def skyline_matrix(
     x: cfg.SkylineMatrixConfig, data: dict[str, Any]
-) -> SkylineMatrixCoercible:
+) -> SkylineMatrixCoercible | None:
+    if x is None:
+        return None
+
     if isinstance(x, str):
         e = _eval_expression(x, data)
         if tg.is_one_or_many_2D_scalars(e):

phylogenie/generators/trees.py

@@ -89,8 +89,8 @@ class CanonicalTreeDatasetGenerator(TreeDatasetGenerator):
     birth_rates: cfg.SkylineVectorConfig = 0
     death_rates: cfg.SkylineVectorConfig = 0
     removal_probabilities: cfg.SkylineVectorConfig = 0
-    migration_rates: cfg.SkylineMatrixConfig = 0
-    birth_rates_among_states: cfg.SkylineMatrixConfig = 0
+    migration_rates: cfg.SkylineMatrixConfig = None
+    birth_rates_among_states: cfg.SkylineMatrixConfig = None
     states: list[str] | None = None
 
     def _get_events(self, rng: Generator, data: dict[str, Any]) -> list[Event]:
@@ -116,8 +116,8 @@ class EpidemiologicalTreeDatasetGenerator(TreeDatasetGenerator):
     become_uninfectious_rates: cfg.SkylineVectorConfig = 0
     sampling_proportions: cfg.SkylineVectorConfig = 1
     removal_probabilities: cfg.SkylineVectorConfig = 1
-    migration_rates: cfg.SkylineMatrixConfig = 0
-    reproduction_numbers_among_states: cfg.SkylineMatrixConfig = 0
+    migration_rates: cfg.SkylineMatrixConfig = None
+    reproduction_numbers_among_states: cfg.SkylineMatrixConfig = None
 
     def _get_events(self, rng: Generator, data: dict[str, Any]) -> list[Event]:
         return get_epidemiological_events(
@@ -142,8 +142,8 @@ class FBDTreeDatasetGenerator(TreeDatasetGenerator):
     turnover: cfg.SkylineVectorConfig = 0
     sampling_proportions: cfg.SkylineVectorConfig = 1
     removal_probabilities: cfg.SkylineVectorConfig = 0
-    migration_rates: cfg.SkylineMatrixConfig = 0
-    diversification_between_types: cfg.SkylineMatrixConfig = 0
+    migration_rates: cfg.SkylineMatrixConfig = None
+    diversification_between_types: cfg.SkylineMatrixConfig = None
 
     def _get_events(self, rng: Generator, data: dict[str, Any]) -> list[Event]:
         return get_FBD_events(

phylogenie/io.py

@@ -85,7 +85,7 @@ def load_fasta(
             else:
                 try:
                     time = float(id.split("|")[-1])
-                except ValueError:
+                except:
                     time = None
             chars = next(f).strip()
             sequences.append(Sequence(id, chars, time))

phylogenie/skyline/matrix.py

@@ -159,9 +159,13 @@ class SkylineMatrix:
     @overload
     def __getitem__(self, item: slice) -> "SkylineMatrix": ...
     @overload
-    def __getitem__(
-        self, item: tuple[int | slice, int | slice]
-    ) -> Union[SkylineParameter, SkylineVector, "SkylineMatrix"]: ...
+    def __getitem__(self, item: tuple[int, int]) -> SkylineParameter: ...
+    @overload
+    def __getitem__(self, item: tuple[int, slice]) -> SkylineVector: ...
+    @overload
+    def __getitem__(self, item: tuple[slice, int]) -> SkylineVector: ...
+    @overload
+    def __getitem__(self, item: tuple[slice, slice]) -> "SkylineMatrix": ...
     def __getitem__(
         self, item: int | slice | tuple[int | slice, int | slice]
     ) -> Union[SkylineParameter | SkylineVector, "SkylineMatrix"]:

phylogenie/treesimulator/events.py

@@ -55,7 +55,7 @@ class DeathEvent(Event):
 
 
 class MigrationEvent(Event):
-    def __init__(self, state: str, target_state: str, rate: SkylineParameterLike):
+    def __init__(self, rate: SkylineParameterLike, state: str, target_state: str):
         super().__init__(rate, state)
         self.target_state = target_state
 
@@ -68,8 +68,8 @@ class SamplingEvent(Event):
     def __init__(
         self,
         rate: SkylineParameterLike,
-        removal_probability: SkylineParameterLike,
         state: str | None = None,
+        removal_probability: SkylineParameterLike = 0,
     ):
         super().__init__(rate, state)
         self.removal_probability = skyline_parameter(removal_probability)
@@ -85,8 +85,8 @@ def get_canonical_events(
     birth_rates: SkylineVectorCoercible = 0,
     death_rates: SkylineVectorCoercible = 0,
     removal_probabilities: SkylineVectorCoercible = 0,
-    migration_rates: SkylineMatrixCoercible = 0,
-    birth_rates_among_states: SkylineMatrixCoercible = 0,
+    migration_rates: SkylineMatrixCoercible | None = None,
+    birth_rates_among_states: SkylineMatrixCoercible | None = None,
     states: list[str] | None = None,
 ) -> list[Event]:
     N = 1 if states is None else len(states)
@@ -95,31 +95,32 @@ def get_canonical_events(
     death_rates = skyline_vector(death_rates, N)
     sampling_rates = skyline_vector(sampling_rates, N)
     removal_probabilities = skyline_vector(removal_probabilities, N)
-    if migration_rates and N == 1:
-        raise ValueError(f"Migration rates require multiple states (got {states}).")
-    if birth_rates_among_states and N == 1:
-        raise ValueError(
-            f"Birth rates among states require multiple states (got {states})."
-        )
 
     events: list[Event] = []
     for i in range(N):
         state = None if states is None else states[i]
         events.append(BirthEvent(birth_rates[i], state, state))
         events.append(DeathEvent(death_rates[i], state))
-        events.append(SamplingEvent(sampling_rates[i], removal_probabilities[i], state))
-        if N > 1:
-            migration_rates = skyline_matrix(migration_rates, N, N - 1)
-            birth_rates_among_states = skyline_matrix(
-                birth_rates_among_states, N, N - 1
-            )
-            assert states is not None
-            assert state is not None
+        events.append(SamplingEvent(sampling_rates[i], state, removal_probabilities[i]))
+
+    if states is not None and migration_rates is not None:
+        migration_rates = skyline_matrix(migration_rates, N, N - 1)
+        for i, state in enumerate(states):
+            for j, other_state in enumerate([s for s in states if s != state]):
+                events.append(MigrationEvent(migration_rates[i, j], state, other_state))
+    elif migration_rates is not None:
+        raise ValueError(f"Migration rates require states to be provided.")
+
+    if states is not None and birth_rates_among_states is not None:
+        birth_rates_among_states = skyline_matrix(birth_rates_among_states, N, N - 1)
+        for i, state in enumerate(states):
             for j, other_state in enumerate([s for s in states if s != state]):
-                events.append(MigrationEvent(state, other_state, migration_rates[i][j]))
                 events.append(
-                    BirthEvent(birth_rates_among_states[i][j], state, other_state)
+                    BirthEvent(birth_rates_among_states[i, j], state, other_state)
                 )
+    elif birth_rates_among_states is not None:
+        raise ValueError(f"Birth rates among states require states to be provided.")
+
     return [event for event in events if event.rate]
 
 
@@ -128,8 +129,8 @@ def get_epidemiological_events(
     reproduction_numbers: SkylineVectorCoercible = 0,
     become_uninfectious_rates: SkylineVectorCoercible = 0,
     removal_probabilities: SkylineVectorCoercible = 1,
-    migration_rates: SkylineMatrixCoercible = 0,
-    reproduction_numbers_among_states: SkylineMatrixCoercible = 0,
+    migration_rates: SkylineMatrixCoercible | None = None,
+    reproduction_numbers_among_states: SkylineMatrixCoercible | None = None,
     states: list[str] | None = None,
 ) -> list[Event]:
     N = 1 if states is None else len(states)
@@ -138,20 +139,20 @@ def get_epidemiological_events(
     become_uninfectious_rates = skyline_vector(become_uninfectious_rates, N)
     sampling_proportions = skyline_vector(sampling_proportions, N)
     removal_probabilities = skyline_vector(removal_probabilities, N)
-    if N == 1 and reproduction_numbers_among_states:
-        raise ValueError(
-            f"Reproduction numbers among states require multiple states (got {states})."
-        )
-    reproduction_numbers_among_states = (
-        skyline_matrix(reproduction_numbers_among_states, N, N - 1) if N > 1 else 0
-    )
 
     birth_rates = reproduction_numbers * become_uninfectious_rates
     sampling_rates = become_uninfectious_rates * sampling_proportions
-    birth_rates_among_states = (
-        reproduction_numbers_among_states * become_uninfectious_rates
-    )
     death_rates = become_uninfectious_rates - removal_probabilities * sampling_rates
+    birth_rates_among_states = None
+    if states is None and reproduction_numbers_among_states is not None:
+        raise ValueError(
+            f"Reproduction numbers among states require states to be provided."
+        )
+    elif reproduction_numbers_among_states is not None:
+        birth_rates_among_states = (
+            skyline_matrix(reproduction_numbers_among_states, N, N - 1)
+            * become_uninfectious_rates
+        )
 
     return get_canonical_events(
         states=states,
@@ -169,32 +170,30 @@ def get_FBD_events(
     turnover: SkylineVectorCoercible = 0,
     sampling_proportions: SkylineVectorCoercible = 1,
     removal_probabilities: SkylineVectorCoercible = 0,
-    migration_rates: SkylineMatrixCoercible = 0,
-    diversification_between_types: SkylineMatrixCoercible = 0,
+    migration_rates: SkylineMatrixCoercible | None = None,
+    diversification_between_types: SkylineMatrixCoercible | None = None,
     states: list[str] | None = None,
-):
+) -> list[Event]:
     N = 1 if states is None else len(states)
 
     diversification = skyline_vector(diversification, N)
     turnover = skyline_vector(turnover, N)
     sampling_proportions = skyline_vector(sampling_proportions, N)
     removal_probabilities = skyline_vector(removal_probabilities, N)
-    if N == 1 and diversification_between_types:
-        raise ValueError(
-            f"Diversification rates among states require multiple states (got {states})."
-        )
-    diversification_between_types = (
-        skyline_matrix(diversification_between_types, N, N - 1) if N > 1 else 0
-    )
 
     birth_rates = diversification / (1 - turnover)
     death_rates = turnover * birth_rates
-    sampling_rates = (
-        sampling_proportions
-        * death_rates
-        / (1 - removal_probabilities * sampling_proportions)
-    )
-    birth_rates_among_states = diversification_between_types + death_rates
+    sampling_rates_dividend = 1 - removal_probabilities * sampling_proportions
+    sampling_rates = sampling_proportions * death_rates / sampling_rates_dividend
+    birth_rates_among_states = None
+    if states is None and diversification_between_types is not None:
+        raise ValueError(
+            f"Diversification rates among states require states to be provided."
+        )
+    elif diversification_between_types is not None:
+        birth_rates_among_states = (
+            skyline_matrix(diversification_between_types, N, N - 1) + death_rates
+        )
 
     return get_canonical_events(
         states=states,

{phylogenie-2.0.1.dist-info → phylogenie-2.0.3.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: phylogenie
-Version: 2.0.1
+Version: 2.0.3
 Summary: Generate phylogenetic datasets with minimal setup effort
 Author: Gabriele Marino
 Author-email: gabmarino.8601@gmail.com
@@ -23,7 +23,7 @@ Description-Content-Type: text/markdown
 ---
 
 [![AliSim](https://img.shields.io/badge/Powered%20by-AliSim-orange?style=flat-square)](https://iqtree.github.io/doc/AliSim)
-[![PyPI package](https://badge.fury.io/py/phylogenie.svg)](https://pypi.org/project/phylogenie/)
+[![PyPI version](https://img.shields.io/pypi/v/phylogenie)](https://pypi.org/project/phylogenie/)
 [![PyPI downloads](https://shields.io/pypi/dm/phylogenie)](https://pypi.org/project/phylogenie/)
 
 Phylogenie is a [Python](https://www.python.org/) package designed to easily simulate phylogenetic datasets—such as trees and multiple sequence alignments (MSAs)—with minimal setup effort. Simply specify the distributions from which your parameters should be sampled, and Phylogenie will handle the rest!
@@ -53,9 +53,6 @@ Phylogenie comes packed with useful features, including:
 - **Arithmetic operations on parameters** 🧮  
   Perform flexible arithmetic operations between parameters directly within the config file.
 
-- **Modular and extendible architecture** 🧩  
-  Easily add new simulation backends as needed.
-
 ## 📦 Installation
 Phylogenie requires [Python](https://www.python.org/) 3.10 to be installed on your system. There are several ways to install Python and managing different Python versions. One popular option is to use [pyenv](https://github.com/pyenv/pyenv).
 
@@ -83,7 +80,7 @@ It includes a collection of thoroughly commented configuration files, organized
 
 For quick start, pick your favorite config file and run Phylogenie with:
 ```bash
-phylogenie examples/<config_file>.yaml
+phylogenie examples/config_file.yaml
 ```
 This command will create the output dataset in the folder specified inside the configuration file, including data directories and metadata files for each dataset split defined in the config.
 

{phylogenie-2.0.1.dist-info → phylogenie-2.0.3.dist-info}/RECORD

@@ -1,28 +1,28 @@
-phylogenie/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
-phylogenie/generators/__init__.py,sha256=kHDCQUSCl2O-wtPxcYEc1d4NTh5wjt-Oci9bYVNbhIY,440
+phylogenie/__init__.py,sha256=1w_0H9lg7hI3b-NLjKuzc34GbJJGyjLq9LrlogecTzI,1759
+phylogenie/generators/__init__.py,sha256=zsOxy28-9j9alOQLIgrOAFfmM58NNHO_NEtW-KXQXAY,888
 phylogenie/generators/alisim.py,sha256=Psh_fAHHrJBCd1kqVtz_BdoVVmNvieHJRMmapQncJhM,2288
-phylogenie/generators/configs.py,sha256=2NctdpKhOK7uuTrau18QCG1OqAAmmTl_c6QYwFnRSU0,1108
+phylogenie/generators/configs.py,sha256=uiqPvdGSwmItQmhRMT-s1xG7DVnfpwaDTZ_wztL1Pwo,1123
 phylogenie/generators/dataset.py,sha256=Q8GsO6b4792VACvwle26QAsbeiyOGUqID2Q3-_urVS4,2543
-phylogenie/generators/factories.py,sha256=AxLNuvZKnrSSUltr--kJ1zr-JjObD-KIip31ETSN0wk,6952
-phylogenie/generators/trees.py,sha256=CXFXeYLYjgJe7tDwCaYrtHBVLl8Wn83in_pDM_2oH6M,8717
+phylogenie/generators/factories.py,sha256=lJ517JC0T8fVu_2NBNC1GUbMPQfRNeNpO5rGdxubQOA,6998
+phylogenie/generators/trees.py,sha256=pSqnYslDMPl4OVsIUFgQu_Rtc2htlBVwhIWnYgSAr3s,8735
 phylogenie/generators/typeguards.py,sha256=6WIQIy6huHok4vYqA9ym6g2kvyyXPggBONWwbiHtDiY,925
-phylogenie/io.py,sha256=lDLWz-7fsw85HCpVjtpGpw2oWxf5Aurd3_jTx2Gt_pU,2865
+phylogenie/io.py,sha256=ZXlofnSh7FX5UJiP0svRHrTraMSNgKa1GiAv0bMz7jU,2854
 phylogenie/main.py,sha256=vtvSpQxBNlYABoFQ25czl-l3fIr4QRo3svWVd-jcArw,1170
 phylogenie/msa.py,sha256=JDGyZUsAq6-m-SQjoCDjAkAZIxfgyl_PDIhdYn5HOow,2064
 phylogenie/py.typed,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
 phylogenie/skyline/__init__.py,sha256=7pF4CUb4ZCLzNYJNhOjpuTOLTRhlK7L6ugfccNqjIGo,620
-phylogenie/skyline/matrix.py,sha256=S9CTsiK8hh4hXEnbsnPFnIXVNUFL9Oe3MMBKWRNEsSE,8951
+phylogenie/skyline/matrix.py,sha256=BnIreEyCymWzo5jiFtLWXE33ck0pbOFNAfcJW7S2m2s,9147
 phylogenie/skyline/parameter.py,sha256=CJ5OEyRQG2Tg1WJWQ1IpfX-6hjJv80Zj8lMoRke5nnQ,4648
 phylogenie/skyline/vector.py,sha256=becZlWLeT12mvpDC5MLGXZPKlgmbkY-DEcuhdQpnsSc,7135
 phylogenie/tree.py,sha256=nNUQrTtUqgc-Z6CHflqROZ8W2ZL1YWGG7xeWmcKwZUU,1684
 phylogenie/treesimulator/__init__.py,sha256=koyVZxKUeMtY1S2WUgqbeIthoGbdUmx9vP4svFi4a6U,459
-phylogenie/treesimulator/events.py,sha256=4TJpf9v9Cx-8g886Ilw2fbjSCuSkIoQ_xfOkm_a4ex8,9356
+phylogenie/treesimulator/events.py,sha256=RfFodTxMAPUtI3ZkDmw5ekXCaRGnQ8Sy3tjR4oJ8CW4,9636
 phylogenie/treesimulator/gillespie.py,sha256=p2bjwc6iyaf4mMI-xregJwOnXkjT2tIhxXbcGncQ5b0,2293
 phylogenie/treesimulator/model.py,sha256=pPl4JYnBIHjWeSgpIypaDOes11tA4bFAGuMRArXbsVo,3577
 phylogenie/typeguards.py,sha256=WBOSJSaOC8VDtrYoA2w_AYEXTpyKdCfmsM29KaKXl3A,1350
 phylogenie/typings.py,sha256=O1X6lGKTjJ2YJz3ApQ-rYb_tEJNUIcHdUIeYlSM4s5o,500
-phylogenie-2.0.1.dist-info/LICENSE.txt,sha256=NUrDqElK-eD3I0WqC004CJsy6cs0JgsAoebDv_42-pw,1071
-phylogenie-2.0.1.dist-info/METADATA,sha256=d019w9ZXrEKhJlf2EaXDpySllJXPky4O9Dr7Fdd5jYA,5555
-phylogenie-2.0.1.dist-info/WHEEL,sha256=sP946D7jFCHeNz5Iq4fL4Lu-PrWrFsgfLXbbkciIZwg,88
-phylogenie-2.0.1.dist-info/entry_points.txt,sha256=Rt6_usN0FkBX1ZfiqCirjMN9FKOgFLG8rydcQ8kugeE,51
-phylogenie-2.0.1.dist-info/RECORD,,
+phylogenie-2.0.3.dist-info/LICENSE.txt,sha256=NUrDqElK-eD3I0WqC004CJsy6cs0JgsAoebDv_42-pw,1071
+phylogenie-2.0.3.dist-info/METADATA,sha256=VoJ3oxfg0YcD-uXcYE4iE-v2vXh3TQRc_Sp8aEN-pmY,5456
+phylogenie-2.0.3.dist-info/WHEEL,sha256=sP946D7jFCHeNz5Iq4fL4Lu-PrWrFsgfLXbbkciIZwg,88
+phylogenie-2.0.3.dist-info/entry_points.txt,sha256=Rt6_usN0FkBX1ZfiqCirjMN9FKOgFLG8rydcQ8kugeE,51
+phylogenie-2.0.3.dist-info/RECORD,,