pyaging 0.6.0__tar.gz → 0.6.2__tar.gz

{pyaging-0.6.0 → pyaging-0.6.2}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.5
 Name: pyaging
-Version: 0.6.0
+Version: 0.6.2
 Summary: A Python-based compendium of GPU-optimized aging clocks.
 Project-URL: homepage, https://github.com/lucascamillomd/pyaging
 Project-URL: repository, https://github.com/lucascamillomd/pyaging

{pyaging-0.6.0 → pyaging-0.6.2}/src/pyaging/__init__.py

@@ -4,4 +4,4 @@ from . import data, logger, models, utils
 from . import predict as pred
 from . import preprocess as pp
 
-__version__ = "0.6.0"
+__version__ = "0.6.2"

{pyaging-0.6.0 → pyaging-0.6.2}/src/pyaging/models/_models.py

@@ -1725,14 +1725,11 @@ class PhenoAge(pyagingModel):
 
     def postprocess(self, x):
         """Convert the fitted log hazard to clinical Phenotypic Age in years."""
-        # log(1 - mortality) equals the negative Gompertz hazard exactly.
-        # Cancel the nested exp/log terms before finite precision rounds the
-        # mortality to zero or one, preserving the published constants.
-        # Levine 2018 supplementary methods, p. 2: gamma = 0.0076927.
-        # The 0.0192 on p. 1 is the Cox variable-selection penalty, not gamma.
-        gamma = 0.0076927
-        log_scale = math.log(0.00553 * math.expm1(120 * gamma) / gamma)
-        return 141.50225 + (x + log_scale) / 0.090165
+        # Kuo et al. (2021), section 4.3, doi:10.1111/acel.13376.
+        # Substitute the negative Gompertz hazard for log(1 - mortality).
+        # Cancel log(exp(xb)) before mortality can round to zero or one.
+        log_scaled_hazard = x + math.log(0.00553 * 1.51714 / 0.0076927)
+        return 141.50 + log_scaled_hazard / 0.09165
 
 
 class PhenoAgeSaoPaulo(pyagingModel):