kinemotion 0.28.0__py3-none-any.whl → 0.29.1__py3-none-any.whl

kinemotion/core/cmj_metrics_validator.py

@@ -0,0 +1,717 @@
+"""CMJ metrics validation using physiological bounds.
+
+Comprehensive validation of Counter Movement Jump metrics against
+biomechanical bounds, cross-validation checks, and consistency tests.
+
+Provides severity levels (ERROR, WARNING, INFO) for different categories
+of metric issues.
+"""
+
+from dataclasses import dataclass, field
+from enum import Enum
+
+from kinemotion.core.cmj_validation_bounds import (
+    AthleteProfile,
+    CMJBounds,
+    MetricBounds,
+    MetricConsistency,
+    RSIBounds,
+    TripleExtensionBounds,
+    estimate_athlete_profile,
+)
+
+
+class ValidationSeverity(Enum):
+    """Severity level for validation issues."""
+
+    ERROR = "ERROR"  # Metrics invalid, likely data corruption
+    WARNING = "WARNING"  # Metrics valid but unusual, needs review
+    INFO = "INFO"  # Normal variation, informational only
+
+
+@dataclass
+class ValidationIssue:
+    """Single validation issue."""
+
+    severity: ValidationSeverity
+    metric: str
+    message: str
+    value: float | None = None
+    bounds: tuple[float, float] | None = None
+
+
+@dataclass
+class ValidationResult:
+    """Complete validation result for CMJ metrics."""
+
+    issues: list[ValidationIssue] = field(default_factory=list)
+    status: str = "PASS"  # "PASS", "PASS_WITH_WARNINGS", "FAIL"
+    athlete_profile: AthleteProfile | None = None
+    rsi: float | None = None
+    height_flight_time_consistency: float | None = None  # % error
+    velocity_height_consistency: float | None = None  # % error
+    depth_height_ratio: float | None = None
+    contact_depth_ratio: float | None = None
+
+    def add_error(
+        self,
+        metric: str,
+        message: str,
+        value: float | None = None,
+        bounds: tuple[float, float] | None = None,
+    ) -> None:
+        """Add error-level issue."""
+        self.issues.append(
+            ValidationIssue(
+                severity=ValidationSeverity.ERROR,
+                metric=metric,
+                message=message,
+                value=value,
+                bounds=bounds,
+            )
+        )
+
+    def add_warning(
+        self,
+        metric: str,
+        message: str,
+        value: float | None = None,
+        bounds: tuple[float, float] | None = None,
+    ) -> None:
+        """Add warning-level issue."""
+        self.issues.append(
+            ValidationIssue(
+                severity=ValidationSeverity.WARNING,
+                metric=metric,
+                message=message,
+                value=value,
+                bounds=bounds,
+            )
+        )
+
+    def add_info(
+        self,
+        metric: str,
+        message: str,
+        value: float | None = None,
+    ) -> None:
+        """Add info-level issue."""
+        self.issues.append(
+            ValidationIssue(
+                severity=ValidationSeverity.INFO,
+                metric=metric,
+                message=message,
+                value=value,
+            )
+        )
+
+    def finalize_status(self) -> None:
+        """Determine final pass/fail status based on issues."""
+        has_errors = any(
+            issue.severity == ValidationSeverity.ERROR for issue in self.issues
+        )
+        has_warnings = any(
+            issue.severity == ValidationSeverity.WARNING for issue in self.issues
+        )
+
+        if has_errors:
+            self.status = "FAIL"
+        elif has_warnings:
+            self.status = "PASS_WITH_WARNINGS"
+        else:
+            self.status = "PASS"
+
+
+class CMJMetricsValidator:
+    """Comprehensive CMJ metrics validator."""
+
+    def __init__(self, assumed_profile: AthleteProfile | None = None):
+        """Initialize validator.
+
+        Args:
+            assumed_profile: If provided, validate against this specific profile.
+                            Otherwise, estimate from metrics.
+        """
+        self.assumed_profile = assumed_profile
+
+    def validate(self, metrics: dict) -> ValidationResult:
+        """Validate CMJ metrics comprehensively.
+
+        Args:
+            metrics: Dictionary with CMJ metric values
+
+        Returns:
+            ValidationResult with all issues and status
+        """
+        result = ValidationResult()
+
+        # Estimate athlete profile if not provided
+        if self.assumed_profile:
+            result.athlete_profile = self.assumed_profile
+        else:
+            result.athlete_profile = estimate_athlete_profile(metrics)
+
+        profile = result.athlete_profile
+
+        # PRIMARY BOUNDS CHECKS
+        self._check_flight_time(metrics, result, profile)
+        self._check_jump_height(metrics, result, profile)
+        self._check_countermovement_depth(metrics, result, profile)
+        self._check_concentric_duration(metrics, result, profile)
+        self._check_eccentric_duration(metrics, result, profile)
+        self._check_peak_velocities(metrics, result, profile)
+
+        # CROSS-VALIDATION CHECKS
+        self._check_flight_time_height_consistency(metrics, result)
+        self._check_velocity_height_consistency(metrics, result)
+        self._check_rsi_validity(metrics, result, profile)
+
+        # CONSISTENCY CHECKS
+        self._check_depth_height_ratio(metrics, result)
+        self._check_contact_depth_ratio(metrics, result)
+
+        # TRIPLE EXTENSION ANGLES
+        self._check_triple_extension(metrics, result, profile)
+
+        # Finalize status
+        result.finalize_status()
+
+        return result
+
+    def _check_flight_time(
+        self, metrics: dict, result: ValidationResult, profile: AthleteProfile
+    ) -> None:
+        """Validate flight time."""
+        flight_time = metrics.get("flight_time")
+        if flight_time is None:
+            return
+
+        bounds = CMJBounds.FLIGHT_TIME
+
+        if not bounds.is_physically_possible(flight_time):
+            if flight_time < bounds.absolute_min:
+                result.add_error(
+                    "flight_time",
+                    f"Flight time {flight_time:.3f}s below frame rate resolution limit",
+                    value=flight_time,
+                    bounds=(bounds.absolute_min, bounds.absolute_max),
+                )
+            else:
+                result.add_error(
+                    "flight_time",
+                    f"Flight time {flight_time:.3f}s exceeds elite human capability",
+                    value=flight_time,
+                    bounds=(bounds.absolute_min, bounds.absolute_max),
+                )
+        elif bounds.contains(flight_time, profile):
+            result.add_info(
+                "flight_time",
+                f"Flight time {flight_time:.3f}s within expected range for {profile.value}",
+                value=flight_time,
+            )
+        else:
+            # Outside expected range but physically possible
+            expected_min, expected_max = self._get_profile_range(profile, bounds)
+            result.add_warning(
+                "flight_time",
+                f"Flight time {flight_time:.3f}s outside typical range "
+                f"[{expected_min:.3f}-{expected_max:.3f}]s for {profile.value}",
+                value=flight_time,
+                bounds=(expected_min, expected_max),
+            )
+
+    def _check_jump_height(
+        self, metrics: dict, result: ValidationResult, profile: AthleteProfile
+    ) -> None:
+        """Validate jump height."""
+        jump_height = metrics.get("jump_height")
+        if jump_height is None:
+            return
+
+        bounds = CMJBounds.JUMP_HEIGHT
+
+        if not bounds.is_physically_possible(jump_height):
+            if jump_height < bounds.absolute_min:
+                result.add_error(
+                    "jump_height",
+                    f"Jump height {jump_height:.3f}m essentially no jump (noise)",
+                    value=jump_height,
+                    bounds=(bounds.absolute_min, bounds.absolute_max),
+                )
+            else:
+                result.add_error(
+                    "jump_height",
+                    f"Jump height {jump_height:.3f}m exceeds human capability",
+                    value=jump_height,
+                    bounds=(bounds.absolute_min, bounds.absolute_max),
+                )
+        elif bounds.contains(jump_height, profile):
+            result.add_info(
+                "jump_height",
+                f"Jump height {jump_height:.3f}m within expected range for {profile.value}",
+                value=jump_height,
+            )
+        else:
+            expected_min, expected_max = self._get_profile_range(profile, bounds)
+            result.add_warning(
+                "jump_height",
+                f"Jump height {jump_height:.3f}m outside typical range "
+                f"[{expected_min:.3f}-{expected_max:.3f}]m for {profile.value}",
+                value=jump_height,
+                bounds=(expected_min, expected_max),
+            )
+
+    def _check_countermovement_depth(
+        self, metrics: dict, result: ValidationResult, profile: AthleteProfile
+    ) -> None:
+        """Validate countermovement depth."""
+        depth = metrics.get("countermovement_depth")
+        if depth is None:
+            return
+
+        bounds = CMJBounds.COUNTERMOVEMENT_DEPTH
+
+        if not bounds.is_physically_possible(depth):
+            if depth < bounds.absolute_min:
+                result.add_error(
+                    "countermovement_depth",
+                    f"Countermovement depth {depth:.3f}m essentially no squat",
+                    value=depth,
+                    bounds=(bounds.absolute_min, bounds.absolute_max),
+                )
+            else:
+                result.add_error(
+                    "countermovement_depth",
+                    f"Countermovement depth {depth:.3f}m exceeds physical limit",
+                    value=depth,
+                    bounds=(bounds.absolute_min, bounds.absolute_max),
+                )
+        elif bounds.contains(depth, profile):
+            result.add_info(
+                "countermovement_depth",
+                f"Countermovement depth {depth:.3f}m within expected range for {profile.value}",
+                value=depth,
+            )
+        else:
+            expected_min, expected_max = self._get_profile_range(profile, bounds)
+            result.add_warning(
+                "countermovement_depth",
+                f"Countermovement depth {depth:.3f}m outside typical range "
+                f"[{expected_min:.3f}-{expected_max:.3f}]m for {profile.value}",
+                value=depth,
+                bounds=(expected_min, expected_max),
+            )
+
+    def _check_concentric_duration(
+        self, metrics: dict, result: ValidationResult, profile: AthleteProfile
+    ) -> None:
+        """Validate concentric duration (contact time)."""
+        duration = metrics.get("concentric_duration")
+        if duration is None:
+            return
+
+        bounds = CMJBounds.CONCENTRIC_DURATION
+
+        if not bounds.is_physically_possible(duration):
+            if duration < bounds.absolute_min:
+                result.add_error(
+                    "concentric_duration",
+                    f"Concentric duration {duration:.3f}s likely phase detection error",
+                    value=duration,
+                    bounds=(bounds.absolute_min, bounds.absolute_max),
+                )
+            else:
+                result.add_error(
+                    "concentric_duration",
+                    f"Concentric duration {duration:.3f}s likely includes standing phase",
+                    value=duration,
+                    bounds=(bounds.absolute_min, bounds.absolute_max),
+                )
+        elif bounds.contains(duration, profile):
+            result.add_info(
+                "concentric_duration",
+                f"Concentric duration {duration:.3f}s within expected range for {profile.value}",
+                value=duration,
+            )
+        else:
+            expected_min, expected_max = self._get_profile_range(profile, bounds)
+            result.add_warning(
+                "concentric_duration",
+                f"Concentric duration {duration:.3f}s outside typical range "
+                f"[{expected_min:.3f}-{expected_max:.3f}]s for {profile.value}",
+                value=duration,
+                bounds=(expected_min, expected_max),
+            )
+
+    def _check_eccentric_duration(
+        self, metrics: dict, result: ValidationResult, profile: AthleteProfile
+    ) -> None:
+        """Validate eccentric duration."""
+        duration = metrics.get("eccentric_duration")
+        if duration is None:
+            return
+
+        bounds = CMJBounds.ECCENTRIC_DURATION
+
+        if not bounds.is_physically_possible(duration):
+            result.add_error(
+                "eccentric_duration",
+                f"Eccentric duration {duration:.3f}s outside physical limits",
+                value=duration,
+                bounds=(bounds.absolute_min, bounds.absolute_max),
+            )
+        elif bounds.contains(duration, profile):
+            result.add_info(
+                "eccentric_duration",
+                f"Eccentric duration {duration:.3f}s within expected range for {profile.value}",
+                value=duration,
+            )
+        else:
+            expected_min, expected_max = self._get_profile_range(profile, bounds)
+            result.add_warning(
+                "eccentric_duration",
+                f"Eccentric duration {duration:.3f}s outside typical range "
+                f"[{expected_min:.3f}-{expected_max:.3f}]s for {profile.value}",
+                value=duration,
+                bounds=(expected_min, expected_max),
+            )
+
+    def _check_peak_velocities(
+        self, metrics: dict, result: ValidationResult, profile: AthleteProfile
+    ) -> None:
+        """Validate peak eccentric and concentric velocities."""
+        # Eccentric
+        ecc_vel = metrics.get("peak_eccentric_velocity")
+        if ecc_vel is not None:
+            bounds = CMJBounds.PEAK_ECCENTRIC_VELOCITY
+            if not bounds.is_physically_possible(ecc_vel):
+                result.add_error(
+                    "peak_eccentric_velocity",
+                    f"Peak eccentric velocity {ecc_vel:.2f} m/s outside limits",
+                    value=ecc_vel,
+                    bounds=(bounds.absolute_min, bounds.absolute_max),
+                )
+            elif bounds.contains(ecc_vel, profile):
+                result.add_info(
+                    "peak_eccentric_velocity",
+                    f"Peak eccentric velocity {ecc_vel:.2f} m/s within range for {profile.value}",
+                    value=ecc_vel,
+                )
+            else:
+                expected_min, expected_max = self._get_profile_range(profile, bounds)
+                result.add_warning(
+                    "peak_eccentric_velocity",
+                    f"Peak eccentric velocity {ecc_vel:.2f} m/s outside typical range "
+                    f"[{expected_min:.2f}-{expected_max:.2f}] for {profile.value}",
+                    value=ecc_vel,
+                    bounds=(expected_min, expected_max),
+                )
+
+        # Concentric
+        con_vel = metrics.get("peak_concentric_velocity")
+        if con_vel is not None:
+            bounds = CMJBounds.PEAK_CONCENTRIC_VELOCITY
+            if not bounds.is_physically_possible(con_vel):
+                if con_vel < bounds.absolute_min:
+                    result.add_error(
+                        "peak_concentric_velocity",
+                        f"Peak concentric velocity {con_vel:.2f} m/s insufficient to leave ground",
+                        value=con_vel,
+                        bounds=(bounds.absolute_min, bounds.absolute_max),
+                    )
+                else:
+                    result.add_error(
+                        "peak_concentric_velocity",
+                        f"Peak concentric velocity {con_vel:.2f} m/s exceeds elite capability",
+                        value=con_vel,
+                        bounds=(bounds.absolute_min, bounds.absolute_max),
+                    )
+            elif bounds.contains(con_vel, profile):
+                result.add_info(
+                    "peak_concentric_velocity",
+                    f"Peak concentric velocity {con_vel:.2f} m/s within range for {profile.value}",
+                    value=con_vel,
+                )
+            else:
+                expected_min, expected_max = self._get_profile_range(profile, bounds)
+                result.add_warning(
+                    "peak_concentric_velocity",
+                    f"Peak concentric velocity {con_vel:.2f} m/s outside typical range "
+                    f"[{expected_min:.2f}-{expected_max:.2f}] for {profile.value}",
+                    value=con_vel,
+                    bounds=(expected_min, expected_max),
+                )
+
+    def _check_flight_time_height_consistency(
+        self, metrics: dict, result: ValidationResult
+    ) -> None:
+        """Verify jump height is consistent with flight time."""
+        flight_time = metrics.get("flight_time")
+        jump_height = metrics.get("jump_height")
+
+        if flight_time is None or jump_height is None:
+            return
+
+        # h = g * t^2 / 8
+        g = 9.81
+        expected_height = (g * flight_time**2) / 8
+        error_pct = abs(jump_height - expected_height) / expected_height
+
+        result.height_flight_time_consistency = error_pct
+
+        if error_pct > MetricConsistency.HEIGHT_FLIGHT_TIME_TOLERANCE:
+            result.add_error(
+                "height_flight_time_consistency",
+                f"Jump height {jump_height:.3f}m inconsistent with flight time {flight_time:.3f}s "
+                f"(expected {expected_height:.3f}m, error {error_pct*100:.1f}%)",
+                value=error_pct,
+                bounds=(0, MetricConsistency.HEIGHT_FLIGHT_TIME_TOLERANCE),
+            )
+        else:
+            result.add_info(
+                "height_flight_time_consistency",
+                f"Jump height and flight time consistent (error {error_pct*100:.1f}%)",
+                value=error_pct,
+            )
+
+    def _check_velocity_height_consistency(
+        self, metrics: dict, result: ValidationResult
+    ) -> None:
+        """Verify peak velocity is consistent with jump height."""
+        velocity = metrics.get("peak_concentric_velocity")
+        jump_height = metrics.get("jump_height")
+
+        if velocity is None or jump_height is None:
+            return
+
+        # h = v^2 / (2*g)
+        g = 9.81
+        expected_velocity = (2 * g * jump_height) ** 0.5
+        error_pct = abs(velocity - expected_velocity) / expected_velocity
+
+        result.velocity_height_consistency = error_pct
+
+        if error_pct > MetricConsistency.VELOCITY_HEIGHT_TOLERANCE:
+            error_msg = (
+                f"Peak velocity {velocity:.2f} m/s inconsistent with "
+                f"jump height {jump_height:.3f}m (expected {expected_velocity:.2f} "
+                f"m/s, error {error_pct*100:.1f}%)"
+            )
+            result.add_warning(
+                "velocity_height_consistency",
+                error_msg,
+                value=error_pct,
+                bounds=(0, MetricConsistency.VELOCITY_HEIGHT_TOLERANCE),
+            )
+        else:
+            result.add_info(
+                "velocity_height_consistency",
+                f"Peak velocity and jump height consistent (error {error_pct*100:.1f}%)",
+                value=error_pct,
+            )
+
+    def _check_rsi_validity(
+        self, metrics: dict, result: ValidationResult, profile: AthleteProfile
+    ) -> None:
+        """Validate Reactive Strength Index."""
+        flight_time = metrics.get("flight_time")
+        concentric_duration = metrics.get("concentric_duration")
+
+        if (
+            flight_time is None
+            or concentric_duration is None
+            or concentric_duration == 0
+        ):
+            return
+
+        rsi = flight_time / concentric_duration
+        result.rsi = rsi
+
+        if not RSIBounds.is_valid(rsi):
+            if rsi < RSIBounds.MIN_VALID:
+                result.add_error(
+                    "rsi",
+                    f"RSI {rsi:.2f} below physiological minimum (likely error)",
+                    value=rsi,
+                    bounds=(RSIBounds.MIN_VALID, RSIBounds.MAX_VALID),
+                )
+            else:
+                result.add_error(
+                    "rsi",
+                    f"RSI {rsi:.2f} exceeds physiological maximum (likely error)",
+                    value=rsi,
+                    bounds=(RSIBounds.MIN_VALID, RSIBounds.MAX_VALID),
+                )
+        else:
+            expected_min, expected_max = RSIBounds.get_rsi_range(profile)
+            if expected_min <= rsi <= expected_max:
+                result.add_info(
+                    "rsi",
+                    f"RSI {rsi:.2f} within expected range [{expected_min:.2f}-{expected_max:.2f}] "
+                    f"for {profile.value}",
+                    value=rsi,
+                )
+            else:
+                result.add_warning(
+                    "rsi",
+                    f"RSI {rsi:.2f} outside typical range [{expected_min:.2f}-{expected_max:.2f}] "
+                    f"for {profile.value}",
+                    value=rsi,
+                    bounds=(expected_min, expected_max),
+                )
+
+    def _check_depth_height_ratio(
+        self, metrics: dict, result: ValidationResult
+    ) -> None:
+        """Check countermovement depth to jump height ratio."""
+        depth = metrics.get("countermovement_depth")
+        jump_height = metrics.get("jump_height")
+
+        if (
+            depth is None or jump_height is None or depth < 0.05
+        ):  # Skip if depth minimal
+            return
+
+        ratio = jump_height / depth
+        result.depth_height_ratio = ratio
+
+        if ratio < MetricConsistency.DEPTH_HEIGHT_RATIO_MIN:
+            result.add_warning(
+                "depth_height_ratio",
+                f"Jump height {ratio:.2f}x countermovement depth: "
+                f"May indicate incomplete squat or standing position detection error",
+                value=ratio,
+                bounds=(
+                    MetricConsistency.DEPTH_HEIGHT_RATIO_MIN,
+                    MetricConsistency.DEPTH_HEIGHT_RATIO_MAX,
+                ),
+            )
+        elif ratio > MetricConsistency.DEPTH_HEIGHT_RATIO_MAX:
+            result.add_warning(
+                "depth_height_ratio",
+                f"Jump height only {ratio:.2f}x countermovement depth: "
+                f"Unusually inefficient (verify lowest point detection)",
+                value=ratio,
+                bounds=(
+                    MetricConsistency.DEPTH_HEIGHT_RATIO_MIN,
+                    MetricConsistency.DEPTH_HEIGHT_RATIO_MAX,
+                ),
+            )
+        else:
+            result.add_info(
+                "depth_height_ratio",
+                f"Depth-to-height ratio {ratio:.2f} within expected range",
+                value=ratio,
+            )
+
+    def _check_contact_depth_ratio(
+        self, metrics: dict, result: ValidationResult
+    ) -> None:
+        """Check contact time to countermovement depth ratio."""
+        contact = metrics.get("concentric_duration")
+        depth = metrics.get("countermovement_depth")
+
+        if contact is None or depth is None or depth < 0.05:
+            return
+
+        ratio = contact / depth
+        result.contact_depth_ratio = ratio
+
+        if ratio < MetricConsistency.CONTACT_DEPTH_RATIO_MIN:
+            result.add_warning(
+                "contact_depth_ratio",
+                f"Contact time {ratio:.2f}s/m to depth ratio: Very fast for depth traversed",
+                value=ratio,
+                bounds=(
+                    MetricConsistency.CONTACT_DEPTH_RATIO_MIN,
+                    MetricConsistency.CONTACT_DEPTH_RATIO_MAX,
+                ),
+            )
+        elif ratio > MetricConsistency.CONTACT_DEPTH_RATIO_MAX:
+            result.add_warning(
+                "contact_depth_ratio",
+                f"Contact time {ratio:.2f}s/m to depth ratio: Slow for depth traversed",
+                value=ratio,
+                bounds=(
+                    MetricConsistency.CONTACT_DEPTH_RATIO_MIN,
+                    MetricConsistency.CONTACT_DEPTH_RATIO_MAX,
+                ),
+            )
+        else:
+            result.add_info(
+                "contact_depth_ratio",
+                f"Contact-depth ratio {ratio:.2f} s/m within expected range",
+                value=ratio,
+            )
+
+    def _check_triple_extension(
+        self, metrics: dict, result: ValidationResult, profile: AthleteProfile
+    ) -> None:
+        """Validate triple extension angles."""
+        angles = metrics.get("triple_extension")
+        if angles is None:
+            return
+
+        hip = angles.get("hip_angle")
+        if hip is not None:
+            if not TripleExtensionBounds.hip_angle_valid(hip, profile):
+                result.add_warning(
+                    "hip_angle",
+                    f"Hip angle {hip:.1f}° outside expected range for {profile.value}",
+                    value=hip,
+                )
+            else:
+                result.add_info(
+                    "hip_angle",
+                    f"Hip angle {hip:.1f}° within expected range for {profile.value}",
+                    value=hip,
+                )
+
+        knee = angles.get("knee_angle")
+        if knee is not None:
+            if not TripleExtensionBounds.knee_angle_valid(knee, profile):
+                result.add_warning(
+                    "knee_angle",
+                    f"Knee angle {knee:.1f}° outside expected range for {profile.value}",
+                    value=knee,
+                )
+            else:
+                result.add_info(
+                    "knee_angle",
+                    f"Knee angle {knee:.1f}° within expected range for {profile.value}",
+                    value=knee,
+                )
+
+        ankle = angles.get("ankle_angle")
+        if ankle is not None:
+            if not TripleExtensionBounds.ankle_angle_valid(ankle, profile):
+                result.add_warning(
+                    "ankle_angle",
+                    f"Ankle angle {ankle:.1f}° outside expected range for {profile.value}",
+                    value=ankle,
+                )
+            else:
+                result.add_info(
+                    "ankle_angle",
+                    f"Ankle angle {ankle:.1f}° within expected range for {profile.value}",
+                    value=ankle,
+                )
+
+    @staticmethod
+    def _get_profile_range(
+        profile: AthleteProfile, bounds: MetricBounds
+    ) -> tuple[float, float]:
+        """Get min/max bounds for specific profile."""
+        if profile == AthleteProfile.ELDERLY:
+            return (bounds.practical_min, bounds.recreational_max)
+        elif profile == AthleteProfile.UNTRAINED:
+            return (bounds.practical_min, bounds.recreational_max)
+        elif profile == AthleteProfile.RECREATIONAL:
+            return (bounds.recreational_min, bounds.recreational_max)
+        elif profile == AthleteProfile.TRAINED:
+            trained_min = (bounds.recreational_min + bounds.elite_min) / 2
+            trained_max = (bounds.recreational_max + bounds.elite_max) / 2
+            return (trained_min, trained_max)
+        elif profile == AthleteProfile.ELITE:
+            return (bounds.elite_min, bounds.elite_max)
+        return (bounds.absolute_min, bounds.absolute_max)