Calculation Methods & Formulas
Complete documentation of all threshold calculations used, with exact formulas, scientific sources, and validation data.
Threshold Calculations (LT1/LT2)
Modified Dmax (Scientific)
Standard βGold standard for threshold determination. Finds the point with the maximum perpendicular distance to the connecting line, starting before the first significant lactate rise.
1. Find start point S:
If La[i] - La[i-1] β₯ 0.4 mmol/L β S = (P[i-1], La[i-1])
2. Calculate distance for each point P(x,y):
d(x,y) = |AΒ·x + BΒ·y + C| / β(AΒ² + BΒ²)
where: A = -(La_end - La_start) / (P_end - P_start)
B = 1
C = AΒ·P_start - La_start
3. LT2 = Point with max(d)Dmax (Original)
LegacyOriginal Dmax method by Cheng (1992). Uses the first data point as the start point, which can lead to underestimation.
Identical to Modified Dmax, but:
Start point S = (P[0], La[0]) β First data point
β οΈ With elevated stress lactate at the start, this can
underestimate the threshold by 10-20W.Can underestimate MLSS (203 Β± 53 W vs. MLSS 225 Β± 60 W)
Dickhuth (IAT)
Individual Anaerobic Threshold based on the lactate equivalent. The minimum of the equivalent indicates maximum metabolic efficiency.
1. Calculate lactate equivalent:
LE(i) = La(i) / P(i)
2. Find minimum:
i_min = argmin[LE(i)]
La_min = La(i_min)
3. Thresholds:
LT1 = Power at (La_min + 0.5 mmol/L)
LT2 = Power at (La_min + 1.5 mmol/L)Delta (+1.5 mmol/L)
Simple and practical method. Adds a fixed offset to the baseline lactate value.
1. Determine baseline:
La_baseline = min(La[0], La[1], La[2])
2. Thresholds:
LT1 = Power at (La_baseline + 0.4 mmol/L)
LT2 = Power at (La_baseline + 1.5 mmol/L)
3. Interpolate power:
LT2_P = P[j] + (LT2_La - La[j])/(La[j+1] - La[j]) Γ (P[j+1] - P[j])Mader (4.0 mmol/L)
LegacyFixed threshold at 4.0 mmol/L (OBLA). Historically widely used, but not individual.
LT2_Lactate = 4.0 mmol/L (fixed value)
LT2_P = P[j] + (4.0 - La[j])/(La[j+1] - La[j]) Γ (P[j+1] - P[j])
β οΈ The individual MLSS varies between 2.0-7.0 mmol/L!This method is NOT individual and systematically overestimates the anaerobic threshold in many athletes.
Log-Log
ExperimentalDmax calculation in logarithmic space. Originally developed for ventilatory thresholds.
1. Log transformation:
x' = logββ(P)
y' = logββ(La)
2. Dmax in Log space:
d'(x',y') = |A'Β·x' + B'Β·y' + C'| / β(A'Β² + B'Β²)
3. Back transformation:
LT2_P = 10^(x'_max)
LT2_La = 10^(y'_max)Interactive method comparison
LAB TIERLoad a saved test and see all 6 threshold methods side by side on the same lactate curve. This visualizes how strongly methods can diverge on your data.
β Unlock LabUncertainty & Validation
β V4The core innovation: honest measurement uncertainty instead of false precision. Every threshold receives a CI, every curve gets coverage, and every progression comparison is checked against measurement tolerance.
Bootstrap confidence intervals
β New in V4Bootstrap procedure for estimating LT1/LT2 measurement uncertainty. Two variants are used: non-parametric resampling with replacement for n β₯ 6, and parametric heteroscedastic Gaussian noise Ο(y) = 0.10 + 0.03Β·y for n = 4-5. Runs 500 deterministic seeded iterations for stable UI output.
Ο(y) = Ο_base + Ο_cv Β· y
Ο_base = 0.10 mmol/L (device floor)
Ο_cv = 0.03 (3% CV)
CIββ = [Pββ, Pββ] β Β±Half-widthGoodness of fit (RΒ²)
β New in V4Polynomial fit (cubic for n β₯ 5, quadratic for n = 4) used to quantify lactate-curve noise. Low RΒ² can indicate noisy measurements, stress lactate or a non-monotone curve.
RΒ² = 1 β (Ξ£(y_i β Ε·_i)Β²) / (Ξ£(y_i β Θ³)Β²)
Categories:
β₯ 95% β Excellent
β₯ 90% β Good
β₯ 80% β Acceptable
< 80% β NoisyCoverage score
β New in V4Checks whether the test actually reached LT2. RΒ² only measures noise: a 4-stage test at 40-60% of LT2 can have RΒ² β 1.0 and still miss LT2 completely. Coverage answers whether the curve truly brackets LT2.
Complete βΊ all 3 criteria are met:
1. Peak lactate β₯ 4.0 mmol/L
2. β₯ 2 stages > 2.5 mmol/L
3. Ξ(last 2 stages) β₯ 0.8 mmol/LProgression significance (error propagation)
β New in V4Compares two tests with error propagation: is the improvement significant relative to measurement uncertainty? This prevents athletes from interpreting noise as progress. Three levels are used: 95%, 68% and not significant.
Ο_combined = β(CI_aΒ² + CI_bΒ²)
|Ξ| > 1.96Β·Ο β sig. 95%
|Ξ| > 1.00Β·Ο β sig. 68%
otherwise β n.s.LT1-anchored training zones
β New in V4The upper boundary of Z2 is the measured LT1, not a fixed Coggan percentage. According to San MillΓ‘n & Brooks 2018, LT1 is the transition from fat-dominant to mixed substrate use, making it the correct Z2 boundary. The middle of Z3 is set between LT1 and LT2. No heuristic is needed.
lt1Frac = lt1 / lt2 (bike)
lt1Frac = lt2 / lt1 (run, speed ratio)
Z1_upper = 0.85 Β· lt1Frac
Z2_upper = lt1Frac β LT1
Z3_upper = (lt1Frac + 1) / 2
Z4_upper = 1.05 β +5% LT2MLSS calibration & method ranking
β New in V4Gold-standard validation: the user uploads MLSS constant-test values, the app compares all 7 methods against this anchor, and ranks them by absolute error. Across multiple comparison pairs this creates a per-athlete method recommendation. Lab tier.
For each method M:
error_M = |LT2_M β MLSS|
Ranking: ascending by mean(error_M)
across all (stepTest, mlssTest) pairs.VLamax Calculation
VLamax Sprint Test (Gold Standard)
Standard βPrecise determination of the maximal glycolytic rate through an all-out sprint.
VLamax = ΞLa / t_eff
where:
ΞLa = La_peak - La_rest [mmol/L]
t_eff = t_sprint - 5.5s [s]
5.5s = Alactic time constant (PCr depletion)
Example:
Sprint: 15s, Rest-La: 1.2, Peak-La: 12.8
t_eff = 15 - 5.5 = 9.5s
VLamax = (12.8 - 1.2) / 9.5 = 1.22 mmol/L/sVLamax Step Test Estimation
ExperimentalApproximation from step test data. Accuracy: Β±15-25%.
VLamax β 0.3 + (ΞLa/ΞP Γ 50) Γ 1.2
where:
ΞLa = La[end] - La[mid] [mmol/L]
ΞP = P[end] - P[mid] [W]
50 = Normalization factor (per 50W)
Range: 0.2 - 1.5 mmol/L/sThis is an estimate! For precise values use the Sprint Test.
Mathematical Foundations
Linear Interpolation
Simple interpolation between two known points.
y = yβ + (x - xβ)/(xβ - xβ) Γ (yβ - yβ)
For threshold calculation (Y β X):
x = xβ + (y - yβ)/(yβ - yβ) Γ (xβ - xβ)Fritsch-Carlson Monotone Spline
β New in V4Monotone cubic Hermite interpolation. It guarantees that no new extrema appear between data points, unlike Catmull-Rom splines, which can oscillate on steep or non-monotone curves. This matters for Dmax geometry because overshoot can create false inflection points.
1. Ξ΄_k = (y_{k+1} β y_k) / (x_{k+1} β x_k)
2. m_k = (Ξ΄_{k-1} + Ξ΄_k) / 2
3. if Ξ΄_{k-1}Β·Ξ΄_k β€ 0 β m_k = 0
4. if Ξ±Β² + Ξ²Β² > 9 β scale m_k, m_{k+1}
5. Hermite basis: hββ, hββ, hββ, hββPoint-to-Line Distance (Dmax)
Core formula for all Dmax-based methods.
Perpendicular distance of a point (x,y) to the line:
|(yβ-yβ)Β·x - (xβ-xβ)Β·y + xβΒ·yβ - yβΒ·xβ|
d(x,y) = βββββββββββββββββββββββββββββββββββββββββ
β[(yβ-yβ)Β² + (xβ-xβ)Β²]
where (xβ,yβ) = Start point, (xβ,yβ) = End pointReferences
1. Jamnick NA et al. (2018). Manipulating graded exercise test variables affects the validity of the lactate threshold.PLoS ONE
2. San MillΓ‘n I & Brooks GA (2018). Assessment of Metabolic Flexibility by Means of Measuring Blood Lactate, Fat, and Carbohydrate Oxidation Responses to Exercise in Professional Endurance Athletes and Less-Fit Individuals.Sports Med 48:467β479
3. Allen H & Coggan A (2019). Training and Racing with a Power Meter (3rd ed.). VeloPress.
4. Jeukendrup A (2014). A step towards personalized sports nutrition: carbohydrate intake during exercise.Sports Med 44(Suppl 1):S25β33
5. Fritsch FN & Carlson RE (1980). Monotone Piecewise Cubic Interpolation.SIAM J Numer Anal 17(2):238β246
6. Pyne DB et al. (2001). Relationships between repeated sprint testing, speed, and endurance.J Sci Med Sport 4(1):17β27
7. Bishop D et al. (1998). The validity of the modified Dmax method for determining the lactate threshold.Med Sci Sports Exerc 30(2):272β278
8. Dickhuth HH et al. (1999). Ventilatory, lactate-derived and catecholamine thresholds.Int J Sports Med
9. Mader A, Heck H (1986). A theory of the metabolic origin of anaerobic threshold.Int J Sports Med
10. Beneke R (2003). Methodological aspects of maximal lactate steady state β implications for performance testing.Eur J Appl Physiol 89(1):95β99