Norming Bias — T Scores Under Skewed Distributions

Many psychological tests (SCL-90-R, BDI, PHQ-9, BAI) assess symptom values on ordinal scales (e.g. 0–4). Sum scores across several items are strongly right-skewed in healthy norming samples: most people have low values, few have high ones. Even though this distribution is not normally distributed, linear T scores are often computed:

This transformation implicitly assumes that X is normally distributed. This leads to systematic biases.

Four right-skewed distribution types can be fitted to raw scores scaled to 0–40. For each raw score, the tool compares the T-norm-implied percentile rank (under the normal-distribution assumption) with the true empirical percentile rank from the chosen distribution.

Gamma (α, β): the rate parametrization is used here (as in R: dgamma(x, shape=α, rate=β)). Shape α determines the shape and skew: skewness = 2/√α. For α < 1, the distribution is extremely right-skewed (monotonically decreasing); as α → ∞ it approaches the normal distribution. Rate β is the reciprocal scale parameter (β = 1/scale): it compresses (large β) or stretches (small β) the distribution. Mean M = α/β, variance σ² = α/β², PDF: f(x) = β^α · x^(α−1) · e^(−βx) / Γ(α). Example: α=1.5, β=0.20 → M=7.5, SD=6.1, γ₁=1.63 — typical for a healthy SCL-90 norm sample (10 items, range 0–40). In practice: α controls the asymmetry, β the location (inversely proportional to the mean).

Log-Normal (μ_log, σ_log): results when the logarithm of the measured value is normally distributed. Skewness rises rapidly with σ_log. Common for measurements that cannot go negative.

Exponential (λ): a special case of the gamma distribution (α=1). Memoryless process, maximum skewness = 2. All raw scores equally likely to decrease — an extreme model.

Weibull (k, λ): a two-parameter distribution that takes very different shapes depending on the shape parameter k. k < 1: strictly monotonically decreasing, decreasing hazard rate (like exponential, but even more extreme). k = 1: identical to Exponential(1/λ). 1 < k < 3.6: right-skewed, bell-shaped — relevant for survival times and reaction times. k > 3.6: approximately normally distributed. Mean = λ · Γ(1 + 1/k), variance = λ² · [Γ(1+2/k) − Γ(1+1/k)²].

Ex-Gaussian (μ, σ, τ): a convolution of a normal distribution N(μ,σ²) and an exponential Exp(1/τ). Common for reaction-time data. Allows separate control over the normally distributed base component and the exponential outlier term (τ). Skewness ≈ 2τ³/(σ²+τ²)^(3/2).

Exact CDFs: for gamma, log-normal, exponential, and Weibull, the tool uses analytically exact cumulative distribution functions (e.g. the regularized incomplete gamma function), which show no discretization error near x=0. Only for the Ex-Gaussian distribution is numerical integration used (trapezoid rule, 1200 steps).

Continuous approximation: the tool models all distributions as continuous. In diagnostic practice (e.g. SCL-90-R), the raw score 0 is often a discrete block with a massive point probability (floor effect). This cliff probability at x=0 is not represented by the tool — it shows the continuous model as an idealization.

Here M̄ is the mean and SD the standard deviation of the empirical distribution — the same values used for the linear T transformation. The difference empir. PR − nominal PR shows the misclassification.

The area-true T norming (Lienert & Raatz, 1998; also: "area transformation" or "normalized standard scores") solves the skew problem by converting each raw score into a T score via its empirical percentile rank — not via a linear formula:

Here F_emp(x) is the empirical cumulative frequency of raw score x and Φ⁻¹ is the quantile function (inverse) of the standard normal distribution. The result: whoever sits at the 84th percentile of the norm sample always gets T_F = 60 — regardless of whether the raw-score distribution is skewed or not.

Comparing both methods:

The "Area Transf." toggle in the tool adds three additional visual elements:

① Green density curve in the first panel: a quantile-based normal approximation — a normal distribution whose parameters are robustly estimated from the empirical quantiles: center = empirical median, spread = (PR84 − PR16)/2. Important: this is not the density of the area transformation itself (which, in raw-score space, would be identical to the empirical density), but a heuristic visualization: it shows what a normal distribution would look like under robust parametrization — making the difference from the linearly-normed orange curve visible.

② Green dashed vertical lines (T_F = 40, 50, 60, 70) in the density plot: they mark exactly the raw-score positions that, under area-true norming, produce T scores 40, 50, 60, 70 — i.e. the empirical PR quantiles 15.9%, 50%, 84.1%, and 97.7%. Comparison with the red lines (linear) shows the shift.

③ Table columns "T Area" and "Δ T": for each linear T score (T = 30–80), the area-true T score at the same raw score is computed and the difference reported. Positive Δ T: the linear norming underestimates the clinical elevation — the patient is more pathological than T_linear suggests.

Lienert, G.A. & Raatz, U. (1998). Testaufbau und Testanalyse (6th ed.). Beltz PVU. [Ch. 7: Norming]
Franke, G.H. (2002). SCL-90-R. Symptom-Checkliste von L.R. Derogatis (2nd ed.). Beltz Test.
Blom, G. (1958). Statistical estimates and transformed beta-variables. Almqvist & Wiksell.