Skip to contents

Computes risk scores for compensating muscles based on biomechanical overuse principles, integrating network topology, duration of compensation, and loading intensity.

Usage

mskCompensationRiskScore(
  compensation_result,
  hg = NULL,
  duration_weeks = 0,
  load_intensity = c("low", "moderate", "high")
)

Arguments

compensation_result

An "MSKCompensation" object from mskDetectCompensation().

hg

An MSKHypergraph object (NULL loads default).

duration_weeks

Numeric, how long compensation has been occurring.

load_intensity

Character, one of "low", "moderate", "high".

Value

An S3 object of class "MSKCompensationRisk" with:

per_muscle_risk

Data.frame of per-muscle risk scores

overall_risk

Numeric overall risk score

overall_category

Character risk category

highest_risk_muscle

Name of highest-risk muscle

recommendation

Clinical action recommendation

Risk Model

Per compensating muscle: risk = z_excess * degree_normalized * duration_factor * load_factor where z_excess = max(0, z_score - z_threshold), degree_normalized = hyperedgeDegree / mean_degree, duration_factor = 1 + log(1 + duration_weeks), load_factor = intensity multiplier (0.5, 1.0, 1.5).

Examples

if (FALSE) { # \dontrun{
risk <- mskCompensationRiskScore(compensation_result,
  duration_weeks = 4, load_intensity = "moderate")
} # }