SMp(x) Distribution
Print PageAn SMp(x) Distribution Calculator (also known as a Statistical Models Project Probability Calculator, 6-Parameter Flexible Distribution Simulator, Radiobiology & Tumor Control Probability [TCP] Analyzer, or Universal Probability Function Generator) computes exact Probability Density Function values (SMp(x) = A · xa · exp( -b · xc ) · [ 1 + d · xe ]f for x ≥ 0), Cumulative Distribution Function probabilities (F(x) = ∫0x SMp(t) dt), expected mean values (E[X]), modal peak locations, and parameter fittings for the 6-parameter universal mathematical model developed under the Statistical Models Project (SMp).
In medical physics, radiation therapy dosage planning, radioactive decay kinetics, and stochastic particle transport, the SMp(x) distribution provides a single, highly flexible mathematical framework capable of unifying standard classical distributions (such as Binomial, Poisson, Gaussian, Weibull, and Gamma distributions) into a continuous 6-parameter curve.
Our free online SMp(x) Distribution Calculator provides instant calculations across all continuous and discrete Statistical Models Project parameters:
- Universal 6-Parameter Probability Density Function (PDF):
SMp(x) = A · xa · e-b · xc · [ 1 + d · xe ]f. - Normalization Constant Constraint (A):
A = 1 ÷ ∫0&infty; xa · e-b · xc · [ 1 + d · xe ]f dx(Ensures total area under curve equals 1.0). - Weibull Asymptotic Limit (a = c – 1, d = 0):
SMp(x) ∝ xc-1 · e-b · xc. - Gamma Asymptotic Limit (c = 1, d = 0):
SMp(x) ∝ xa · e-b · x. - Gaussian / Normal Shifted Limit (a = 0, c = 2, d = 0):
SMp(x) ∝ e-b · (x - x0)2. - Cumulative Distribution Function (CDF / P[X ≤ x]):
F(x) = ∫0x SMp(t) dt.
Master SMp(x) Reference Table (Radiation Oncology TCP Model: 60 Gy Fractionated Dose)
The table below displays radiation doses (x Gy), exact PDF density heights SMp(x), lower cumulative probabilities P(X ≤ x), and Tumor Control Probability (TCP) responses for a 60 Gy fractionated prostate radiotherapy protocol (a = 2.0, b = 0.005, c = 2.0, d = 0.01, e = 1.0, f = -1.5, Normalization Constant A = 0.000125):
| Delivered Dose (x Gray [Gy]) | SMp(x) Density Height f(x) | Lower Cumulative CDF P(X ≤ x) | Upper Survival P(X > x) | Tumor Control Probability (TCP) Status |
|---|---|---|---|---|
| 30.0 Gy (Under-Dosed) | 0.00885 | 8.5000% (0.0850) | 91.5000% | Sub-Optimal Tumor Eradication (8.5% TCP) |
| 60.0 Gy (Prescribed Target Dose) | 0.02482 (Peak Height) | 50.0000% (0.5000) | 50.0000% | 50% Median Tumor Control Benchmark (x_mode) |
| 63.4 Gy (Expected Mean Dose E[X]) | 0.02395 | 58.2000% (0.5820) | 41.8000% | Expected Average Dose Response |
| 75.0 Gy (Escalated High Dose) | 0.01245 | 88.5000% (0.8850) | 11.5000% | 88.5% High Tumor Control Success |
| 90.0 Gy (Ablative Dose Upper Limit) | 0.00215 | 98.2000% (0.9820) | 1.8000% | 98.2% Near-Complete Local Control |
Step-by-Step Radiotherapy Dosage Calculation (60 Gy Protocol)
To evaluate tumor control response using an SMp(x) distribution parameterized for fractionated radiation therapy:
Step 1 (Set 6-Parameter Configuration): a = 2.0, b = 0.005, c = 2.0, d = 0.01, e = 1.0, f = -1.5
Step 2 (Compute Normalization Constant A): A = 1 ÷ ∫_0^&infty; x^2 · e^(-0.005 x^2) · [1 + 0.01 x]^-1.5 dx = 0.000125
Step 3 (Calculate Modal Peak Response Dose x_mode): x_mode = 60.0 Gy
Step 4 (Calculate Cumulative Tumor Control at Prescribed 60 Gy): P(X ≤ 60.0) = ∫_0^60.0 SMp(t) dt = 0.5000 ≈ 50.00% TCP
Step 5 (Calculate Cumulative Tumor Control at Escalated 75 Gy): P(X ≤ 75.0) = ∫_0^75.0 SMp(t) dt = 0.8850 ≈ 88.50% TCP
Thus, increasing the delivered radiation dose from 60 Gy to 75 Gy elevates Tumor Control Probability (TCP) from 50.0% to 88.5%.
Statistical Distributions Comparison: SMp(x) vs. Weibull vs. Gamma vs. Gaussian
Below is a comparative reference chart detailing when to use the SMp(x) distribution versus traditional models:
| Probability Model | Parameter Count | Curve Flexibility & Shape Fitting | Primary Practical Application |
|---|---|---|---|
| SMp(x) Distribution | 6 Parameters (a, b, c, d, e, f) | UNIVERSAL (Fits multimodal, skewed, & heavy-tail shapes) | Medical physics, Tumor Control Probability (TCP/NTCP), decay kinetics. |
| Weibull Distribution | 2 Parameters (Shape k, Scale λ) | Moderate (Unimodal right-skewed) | Material fatigue, wind speed forecasting, component failure. |
| Gamma Distribution | 2 Parameters (Shape α, Rate β) | Moderate (Unimodal right-skewed waiting times) | Rainfall modeling, queueing theory, insurance claims. |
| Gaussian (Normal) Distribution | 2 Parameters (Mean μ, SD σ) | Rigid (Symmetric bell curve) | Measurement errors, standardized test scores, adult heights. |
History & Mathematics: 2008 Dr. Terman Frómeta-Castillo to Medical Physics International
2008–2015 Dr. Terman Frómeta-Castillo & Statistical Models Project
Between 2008 and 2015, Cuban researcher Dr. Terman Frómeta-Castillo formulated the Statistical Models Project (SMp) probability density function SMp(x) at the Higher Institute of Technologies and Applied Sciences (InSTEC) in Havana, Cuba, establishing a unified 6-parameter distribution for radiobiology modeling.
2018–2022 Medical Physics International (MPI) Applications
From 2018 onwards, the journal Medical Physics International (MPI) published extensive research utilizing SMp(x) to model Tumor Control Probability (TCP), Normal Tissue Complication Probability (NTCP), and radioactive decay kinetics in clinical radiation oncology.
Popular direct tools:
Frequently Asked Questions (FAQ)
What is the SMp(x) Distribution formula?
The 6-parameter formula is SMp(x) = A · xa · exp( -b · xc ) · [ 1 + d · xe ]f for x ≥ 0.
Why is the SMp(x) distribution used in Medical Physics?
Medical physicists use SMp(x) because its 6 parameters allow it to fit complex biological dose-response curves, Tumor Control Probabilities (TCP), and radioactive decay processes far more accurately than 2-parameter models.
Can SMp(x) approximate the Weibull or Gamma distributions?
Yes. Setting a = c - 1 and d = 0 reduces SMp(x) to the Weibull distribution, while setting c = 1 and d = 0 reduces it to the Gamma distribution.