Statistics & Probability Statistics & Probability major 100% Free Open Access
Chapter 3 • Theory & Derivations

Measurement Theory, Psychometrics, Reliability & Validity

Levels of measurement: nominal, ordinal, interval, and ratio scales (Stevens' taxonomy) and admissible statistical operations; Classical Test Theory (CTT): true score model X = T + E, assumptions regarding measurement error, and error propagation; Reliability estimation methods: test-retest reliability, parallel-forms reliability, and internal consistency; Cronbach's alpha (alpha): algebraic derivation from variance-covariance matrix, split-half Spearman-Brown prophecy formula, and limitations; Modern internal consistency metrics: McDonald's omega_total and omega_hierarchical via confirmatory factor analysis; Validity frameworks: content validity, criterion-related validity (concurrent and predictive), and construct validity; The Campbell-Fiske Multitrait-Multimethod (MTMM) Matrix: establishing convergent and discriminant validity.

§3.1 Stevens' Taxonomy of Measurement Scales: Nominal, Ordinal, Interval, and Ratio

In a landmark 1946 paper in Science, S. Smith Stevens established the foundational taxonomy of measurement scales, classifying scales according to the mathematical group of transformations under which the empirical information is invariant:

  1. Nominal Scale:
  • Permissible Transformations: Any one-to-one mapping (injective and surjective function) $f(x)$.
  • Mathematical Structure: Equivalence relations ($=$ or $\ne$). Numbers serve exclusively as qualitative labels or category identifiers (e.g., Biological Sex: $1 = \text{Male}, 2 = \text{Female}$; Blood Types: $1 = \text{A}, 2 = \text{B}, 3 = \text{AB}, 4 = \text{O}$).
  • Admissible Statistics: Mode, frequency counts, contingency tables, chi-square tests, entropy, and mutual information. Arithmetic operations (mean, variance) are completely meaningless.
  1. Ordinal Scale:
  • Permissible Transformations: Any strictly monotonic increasing transformation $f(x)$, where $x_1 > x_2 \implies f(x_1) > f(x_2)$.
  • Mathematical Structure: Strict order relations (@@MATH_TOKEN_0@@gt;$, @@MATH_TOKEN_1@@lt;$, $=$). The scale indicates relative rank order, but distances between adjacent points are non-equal and undefined (e.g., Cancer Staging: Stage I, II, III, IV; Education Level; Likert agreement categories).
  • Admissible Statistics: Median, percentiles, interquartile range, Spearman's rank correlation ($\rho$), Kendall's tau ($\tau$), and non-parametric rank tests (Mann-Whitney $U$, Kruskal-Wallis).
  1. Interval Scale:
  • Permissible Transformations: Positive linear (affine) transformations:
$$f(x) = a x + b, \quad a > 0$$
  • Mathematical Structure: Differences between numerical values represent equal empirical quantities (equal intervals), but the scale lacks an absolute, non-arbitrary physical zero point (e.g., Celsius and Fahrenheit temperature, standardized IQ scores, Gregorian calendar years). Because zero is arbitrary, ratios are meaningless: $20^\circ\mathrm{C}$ is not "twice as hot" as $10^\circ\mathrm{C}$ (in Fahrenheit, $20^\circ\mathrm{C} = 68^\circ\mathrm{F}$ and $10^\circ\mathrm{C} = 50^\circ\mathrm{F}$, where $68 / 50 = 1.36 \ne 2.0$).
  • Admissible Statistics: Arithmetic mean, variance, standard deviation, Pearson product-moment correlation ($r$), regression coefficients, and Analysis of Variance (ANOVA).
  1. Ratio Scale:
  • Permissible Transformations: Positive scalar similarity transformations:
$$f(x) = a x, \quad a > 0$$
  • Mathematical Structure: Possesses an absolute, physically meaningful, non-arbitrary zero point representing the total absence of the measured attribute (e.g., Kelvin temperature, mass in kilograms, length in meters, monetary income, reaction time in milliseconds).
  • Admissible Statistics: Geometric mean, harmonic mean, coefficient of variation ($s/\bar{x}$), percentage ratios, and all logarithmic and power transformations.

A pervasive controversy in behavioral research is the statistical treatment of Likert response items (e.g., 5-point scales: Strongly Disagree to Strongly Agree). While individual Likert items are strictly ordinal, psychometricians frequently sum or average multiple items into multi-item composite scale scores, treating the resulting composite as interval-level data. This practice is mathematically justified only when the composite scale satisfies the axioms of modern Item Response Theory or conjoint measurement.

§3.2 Classical Test Theory (CTT): The True Score Model $X = T + E$

Classical Test Theory, formulated by Charles Spearman (1904) and synthesized by Frederic Lord and Melvin Novick (1968), serves as the foundational mathematical model of psychometric measurement.

Let $X$ denote the observed score obtained by an individual on a measurement instrument. CTT posits that $X$ is the linear sum of an unobservable true score $T$ and a random measurement error $E$:

$$X = T + E$$

The true score $T$ is defined mathematically as the expected value of observed scores across hypothetical infinitely repeated, independent administrations of the instrument to the identical individual:

$$T \equiv \mathbb{E}[X]$$

From this definition, the fundamental axioms of Classical Test Theory follow:

  1. Axiom of Zero Expected Error: The expected value of measurement error across the population is zero:
$$\mathbb{E}[E] = \mathbb{E}[X - T] = \mathbb{E}[X] - \mathbb{E}[T] = 0$$
  1. Axiom of Error Independence: Measurement error is uncorrelated with the true score:
$$\mathrm{Cov}(T, E) = 0 \implies \rho_{TE} = 0$$
  1. Axiom of Parallel Error Uncorrelatedness: Measurement errors on two distinct, independent measurement occasions or distinct parallel test forms are uncorrelated:
$$\mathrm{Cov}(E_1, E_2) = 0 \implies \rho_{E_1 E_2} = 0$$

From these foundational axioms, the decomposition of observed score variance follows directly:

$$\mathrm{Var}(X) = \mathrm{Var}(T + E) = \mathrm{Var}(T) + \mathrm{Var}(E) + 2\mathrm{Cov}(T, E)$$

Since $\mathrm{Cov}(T, E) = 0$:

$$\sigma_X^2 = \sigma_T^2 + \sigma_E^2$$

The total observed phenotypic variance across a population decomposes into two mutually orthogonal components: true score variance ($\sigma_T^2$) reflecting genuine individual differences in the latent attribute, and error variance ($\sigma_E^2$) reflecting random noise.

§3.3 Derivation and Estimation of Reliability

In Classical Test Theory, reliability ($\rho_{XX'}$) is defined formally as the ratio of true score variance to total observed score variance:

$$\rho_{XX'} \equiv \frac{\sigma_T^2}{\sigma_X^2} = \frac{\sigma_T^2}{\sigma_T^2 + \sigma_E^2} = 1 - \frac{\sigma_E^2}{\sigma_X^2}$$

Reliability is a normalized index bounded strictly between $0.0$ and $1.0$:

  • $\rho_{XX'} = 1.0$: The instrument exhibits perfect precision; error variance $\sigma_E^2 = 0$, and observed scores are identical to true scores ($X = T$).
  • $\rho_{XX'} = 0.0$: The instrument measures pure random noise; true score variance $\sigma_T^2 = 0$, and observed differences reflect chance fluctuations.

Because true scores $T$ and errors $E$ are latent and cannot be directly observed, $\sigma_T^2$ cannot be calculated directly. Reliability must be estimated empirically through parallel tests. Two test forms $X_1$ and $X_2$ are defined as strictly parallel if:

  1. They share identical true scores for every individual: $T_1 = T_2 = T$.
  2. They exhibit identical error variances: $\sigma_{E_1}^2 = \sigma_{E_2}^2 = \sigma_E^2$.

Theorem: The Pearson correlation between two strictly parallel tests equals the reliability of either test:

$$\rho_{X_1 X_2} = \frac{\mathrm{Cov}(X_1, X_2)}{\sigma_{X_1} \sigma_{X_2}} = \frac{\mathrm{Cov}(T + E_1, T + E_2)}{\sigma_X^2} = \frac{\mathrm{Var}(T) + \mathrm{Cov}(T, E_2) + \mathrm{Cov}(E_1, T) + \mathrm{Cov}(E_1, E_2)}{\sigma_X^2}$$

By the CTT axioms, all covariance terms involving error vanish:

$$\rho_{X_1 X_2} = \frac{\sigma_T^2}{\sigma_X^2} = \rho_{XX'}$$

Practical Reliability Estimation Paradigms:

  1. Test-Retest Reliability: Administering the identical test instrument to the same cohort on two occasions separated by a time interval $\Delta t$. Captures temporal stability, but is vulnerable to carry-over testing effects, memory recall, and true developmental maturation.
  2. Parallel (Alternate) Forms Reliability: Administering two distinct, psychometrically matched forms of the test to the same cohort. Eliminates memory recall, but creating truly parallel forms is intellectually challenging and expensive.
  3. Internal Consistency Reliability: Assessing the degree of inter-correlation and homogeneity among the items comprising a single administration of an instrument (Split-Half, Cronbach's Alpha, McDonald's Omega).

The Spearman-Brown Prophecy Formula

When an internal consistency estimate is computed by splitting a test into two halves ($X_1, X_2$) and calculating the split-half correlation $r_{12}$, the resulting correlation reflects the reliability of a half-length test. Charles Spearman and William Brown independently derived the prophecy formula to calculate the reliability of a test lengthened or shortened by a factor of $k$:

$$\rho_{kk} = \frac{k \rho_{11}}{1 + (k - 1)\rho_{11}}$$

For doubling a half-test to full length ($k = 2$):

$$r_{\text{full}} = \frac{2 r_{12}}{1 + r_{12}}$$

§3.4 Cronbach's Alpha ($\alpha$): Algebraic Derivation from Covariance Matrix and Limitations

The most ubiquitous metric of internal consistency in empirical science is Cronbach's Alpha, formulated by Lee J. Cronbach in 1951.

Consider a composite test scale $X$ composed of $k$ distinct items:

$$X = \sum_{j=1}^k Y_j$$

The total observed scale variance is the sum of the elements in the $k \times k$ item variance-covariance matrix $\boldsymbol{\Sigma}$:

$$\sigma_X^2 = \mathrm{Var}\left(\sum_{j=1}^k Y_j\right) = \sum_{j=1}^k \sigma_j^2 + \sum_{j \ne l} \sigma_{jl} = \sum_{j=1}^k \sigma_j^2 + 2 \sum_{j < l} \mathrm{Cov}(Y_j, Y_l)$$

Cronbach's Alpha is defined as:

$$\alpha \equiv \frac{k}{k - 1} \left( 1 - \frac{\sum_{j=1}^k \sigma_j^2}{\sigma_X^2} \right) = \frac{k}{k - 1} \left( \frac{\sum_{j \ne l} \sigma_{jl}}{\sigma_X^2} \right)$$

When all items are standardized to unit variance ($\sigma_j^2 = 1$), the formula simplifies to the Standardized Alpha:

$$\alpha_{\text{standardized}} = \frac{k \bar{r}}{1 + (k - 1)\bar{r}}$$

where $\bar{r}$ is the mean off-diagonal inter-item correlation.

Essential $\tau$-Equivalence and Limitations of $\alpha$:

Cronbach's alpha equals the true population reliability $\rho_{XX'}$ if and only if the items satisfy the strict assumption of essential $\tau$-equivalence:

  1. Every item measures the identical single underlying latent factor with equal factor loadings:
$$Y_j = \mu_j + \lambda \xi + E_j \quad (\lambda_1 = \lambda_2 = \dots = \lambda_k = \lambda)$$
  1. Error terms are mutually uncorrelated: $\mathrm{Cov}(E_j, E_l) = 0$ for all $j \ne l$.

In real-world psychometric data, the assumption of equal loadings ($\tau$-equivalence) is almost universally violated: items vary in discrimination power (congeneric models: $\lambda_j \ne \lambda_l$). Under congeneric measurement, Cronbach's alpha is a mathematical lower bound to true reliability:

$$\alpha \le \rho_{XX'}$$

Furthermore, $\alpha$ suffers from two critical vulnerabilities:

  1. Inflation via Test Length: Notice from the formula that as the number of items $k \to \infty$, $\alpha \to 1.0$ even if the inter-item correlation $\bar{r}$ is modest. A bloated 50-item scale with redundant questions will yield $\alpha > 0.95$, while a pristine 4-item scale may yield $\alpha \approx 0.78$ despite higher psychometric quality.
  2. Fallacy of Unidimensionality: High alpha does not imply that a scale is unidimensional. A multidimensional scale measuring two completely distinct constructs (e.g., Depression and Anxiety) will easily produce $\alpha > 0.85$ if the total number of items is sufficiently large.

§3.5 Modern Composite Reliability: McDonald's $\omega_{\text{total}}$ and Factor Analytic Models

To overcome the severe limitations of Cronbach's alpha, modern psychometrics relies on structural equation modeling and confirmatory factor analysis (CFA) to calculate McDonald's Omega ($\omega_{\text{total}}$), introduced by Roderick P. McDonald (1999).

Consider a congeneric single-factor measurement model where each manifest item $Y_j$ loads onto a general latent construct $\xi$:

$$Y_j = \nu_j + \lambda_j \xi + \delta_j, \quad j = 1, \dots, k$$

where:

  • $\nu_j$ is the item intercept.
  • $\lambda_j$ is the unconstrained item factor loading.
  • $\xi$ is the latent factor with standardized variance $\mathrm{Var}(\xi) = 1$.
  • $\delta_j$ is the unique measurement error with variance $\theta_j = \mathrm{Var}(\delta_j)$.

The total observed composite score is $X = \sum_{j=1}^k Y_j$. The true score variance is the variance attributable to the shared latent factor:

$$\sigma_T^2 = \mathrm{Var}\left( \sum_{j=1}^k \lambda_j \xi \right) = \left( \sum_{j=1}^k \lambda_j \right)^2 \mathrm{Var}(\xi) = \left( \sum_{j=1}^k \lambda_j \right)^2$$

The total observed scale variance is:

$$\sigma_X^2 = \left( \sum_{j=1}^k \lambda_j \right)^2 + \sum_{j=1}^k \theta_j$$

McDonald's Omega Total is defined as:

$$\omega_{\text{total}} \equiv \frac{\left( \sum_{j=1}^k \lambda_j \right)^2}{\left( \sum_{j=1}^k \lambda_j \right)^2 + \sum_{j=1}^k \theta_j}$$

For standardized factor loadings (where $\lambda_j^2 + \theta_j = 1 \implies \theta_j = 1 - \lambda_j^2$):

$$\omega_{\text{total}} = \frac{\left( \sum_{j=1}^k \lambda_j \right)^2}{\left( \sum_{j=1}^k \lambda_j \right)^2 + \sum_{j=1}^k (1 - \lambda_j^2)}$$

McDonald's $\omega_{\text{total}}$ provides an unbiased estimate of composite reliability under congeneric models. Methodologists strongly advise reporting $\omega_{\text{total}}$ alongside or in place of Cronbach's $\alpha$.

§3.7 Validity Frameworks: Content, Criterion, and Construct Validity

While reliability addresses whether a measurement instrument is free from random error (precision and consistency), validity addresses whether the instrument actually measures the intended target construct (truth and meaning).

In modern psychometric standards (AERA, APA, NCME), validity is conceptualized as an integrated, unitary construct supported by distinct lines of empirical and theoretical evidence:

  1. Content-Related Validity:

The extent to which the items comprising an instrument represent a comprehensive, unbiased sample of the entire theoretical domain of the construct. Evaluated qualitatively and quantitatively via Expert Panels using the Content Validity Ratio (CVR, Lawshe 1975):

$$\mathrm{CVR} = \frac{n_e - N/2}{N/2}$$

where $N$ is total panel experts and $n_e$ is the number rating an item as "essential."

  1. Criterion-Related Validity:

The extent to which scores on the instrument correlate with an external, objective gold-standard behavioral outcome or physiological criterion:

  • Concurrent Validity: The criterion is measured simultaneously with the instrument (e.g., correlating a rapid 5-minute cognitive screen with a comprehensive 4-hour neuropsychological battery).
  • Predictive Validity: The instrument predicts a future real-world outcome (e.g., SAT scores predicting first-year university GPA, or cardiovascular risk scores predicting myocardial infarction within 10 years).
  1. Construct Validity:

The overarching validation framework establishing that the instrument functions as a faithful operationalization of the latent theoretical construct within a nomological network (evaluated via factor analysis, MTMM matrices, and experimental manipulations).

§3.8 The Campbell-Fiske Multitrait-Multimethod (MTMM) Matrix

In 1959, Donald T. Campbell and Donald W. Fiske revolutionized construct validation by introducing the Multitrait-Multimethod (MTMM) Matrix.

Campbell and Fiske recognized that any observed measurement score reflects a compound of trait variance (the target construct) and method variance (the specific measurement modality, such as self-report survey, peer evaluation, or physiological assay):

$$\text{Observed Score} = \text{Trait} + \text{Method} + \text{Error}$$

If two self-report scales correlate at $r = 0.60$, the correlation may reflect shared method bias (e.g., common rater response style, social desirability) rather than genuine construct commonality.

An MTMM study measures at least two distinct traits ($T_1, T_2, T_3$) using at least two distinct measurement methods ($M_1, M_2, M_3$), generating four distinct correlation quadrants:

  1. Monotrait-Monomethod (Reliability Diagonal): Correlations between identical traits measured by identical methods (test-retest reliability). These must be the highest values in the matrix ($r \ge 0.80$).
  2. Monotrait-Heteromethod (Validity Diagonal): Correlations between identical traits measured by different methods (e.g., Self-Report Anxiety vs. Clinical Interview Anxiety). This directly quantifies Convergent Validity: values must be statistically significant and substantial ($r \ge 0.50$).
  3. Heterotrait-Monomethod: Correlations between different traits measured by the same method (e.g., Self-Report Anxiety vs. Self-Report Depression). High values reveal Common Method Bias.
  4. Heterotrait-Heteromethod: Correlations between different traits measured by different methods. These quantify baseline independence and must be the lowest correlations in the matrix.

Discriminant validity is confirmed if and only if the convergent validity coefficients (Monotrait-Heteromethod) are substantially higher than both the Heterotrait-Monomethod and Heterotrait-Heteromethod coefficients.

§3.9 Item Response Theory (IRT)

Classical Test Theory operates at the total test score level, suffering from two fatal circularities:

  1. Item difficulty is sample-dependent: items appear "easy" if administered to an advanced sample and "hard" if administered to a remedial sample.
  2. Examinee ability is test-dependent: examinees score "high" on an easy test and "low" on a hard test.

Modern psychometrics resolves these limitations through Item Response Theory (IRT), which models the probability of an individual's response to an individual item as a mathematical function of the examinee's underlying latent trait level ($\theta \in (-\infty, +\infty)$) and intrinsic item parameters.

Major Dichotomous IRT Models:

  1. The One-Parameter Logistic (1PL / Rasch) Model:
$$P(Y_{ij} = 1 \mid \theta_i) = \frac{1}{1 + \exp(-(\theta_i - b_j))}$$

where $\theta_i$ is person ability and $b_j$ is item difficulty (the location on the $\theta$ scale where the probability of a correct response is exactly $0.50$). Under the Rasch model, raw sum scores are minimal sufficient statistics for $\theta$.

  1. The Two-Parameter Logistic (2PL) Model:
$$P(Y_{ij} = 1 \mid \theta_i) = \frac{1}{1 + \exp(-a_j(\theta_i - b_j))}$$

where $a_j$ is item discrimination (the slope of the Item Characteristic Curve at inflection point $b_j$). Items with higher $a_j$ differentiate sharply between examinees just below versus just above difficulty $b_j$.

  1. The Three-Parameter Logistic (3PL) Model:
$$P(Y_{ij} = 1 \mid \theta_i) = c_j + \frac{1 - c_j}{1 + \exp(-a_j(\theta_i - b_j))}$$

where $c_j \in [0, 1)$ is the pseudo-guessing parameter, representing the asymptotic lower bound probability that an examinee with $\theta \to -\infty$ answers a multiple-choice item correctly by blind guessing.

Information Functions:

In IRT, measurement precision is not a single static number (like CTT reliability $\rho_{XX'}$); precision varies continuously across the trait spectrum $\theta$. The Item Information Function (IIF) for a 2PL item is:

$$I_j(\theta) = a_j^2 P_j(\theta) [1 - P_j(\theta)]$$

By the conditional independence of items, the Test Information Function (TIF) is additive:

$$I(\theta) = \sum_{j=1}^k I_j(\theta)$$

The standard error of measurement is directly inversely related to information:

$$\mathrm{SE}(\theta) = \frac{1}{\sqrt{I(\theta)}}$$

This property enables Computerized Adaptive Testing (CAT): tailoring item administration dynamically to maximize information at the examinee's estimated trait level $\hat{\theta}$.

§3.10 Differential Item Functioning (DIF) and Measurement Invariance

A critical ethical and psychometric requirement in educational, psychological, and medical assessment is test fairness and measurement invariance. An instrument must measure the identical construct with identical calibration across diverse demographic subgroups (e.g., gender, race, language background).

Differential Item Functioning (DIF) occurs when examinees from different subgroups who possess the identical underlying trait level $\theta$ exhibit different probabilities of answering an item correctly:

$$P(Y_j = 1 \mid \theta, G = \text{Focal}) \ne P(Y_j = 1 \mid \theta, G = \text{Reference})$$

DIF is categorized into:

  • Uniform DIF: The item favors one group uniformly across all trait levels $\theta$ (a shift in difficulty parameter $b_j$).
  • Non-Uniform DIF: The item favors different groups at different trait levels (an interaction reflected in differing discrimination slopes $a_j$).

In multigroup structural equation modeling, measurement invariance is tested hierarchically:

  1. Configural Invariance: The same factor structure (same pattern of zero and non-zero loadings) holds across groups.
  2. Metric (Weak) Invariance: Factor loadings are invariant ($\boldsymbol{\Lambda}_{\text{focal}} = \boldsymbol{\Lambda}_{\text{ref}}$). Required for meaningful cross-group comparison of covariances and regression slopes.
  3. Scalar (Strong) Invariance: Item intercepts are invariant ($\boldsymbol{\nu}_{\text{focal}} = \boldsymbol{\nu}_{\text{ref}}$). Strictly required for meaningful comparison of latent factor means across groups.
  4. Strict Invariance: Unique item residual variances are invariant ($\boldsymbol{\Theta}_{\text{focal}} = \boldsymbol{\Theta}_{\text{ref}}$).

Comparing group manifest means when scalar invariance fails is scientifically invalid: observed differences reflect measurement bias rather than genuine population disparities.

§3.12 Practical Rules for Constructing Psychometric Scales

When researchers develop novel self-report psychometric inventories, adhering to established psychometric standards prevents common measurement artifacts:

  1. Item Pool Generation: Begin with a broad pool of candidate items (typically 3 to 4 times the desired final scale length) generated via comprehensive literature review, qualitative stakeholder interviews, and clinical focus groups.
  2. Cognitive Pretesting: Administer candidate items to a pilot cohort using "think-aloud" protocols to identify ambiguous wording, confusing double negatives, and culturally specific idioms.
  3. Item Difficulty and Discrimination Indexing: In binary or Likert items, examine item endorsement frequency $p$. Items with extreme difficulty ($p < 0.10$ or $p > 0.90$) exhibit low variance and provide minimal discrimination, whereas items near $p \approx 0.50$ maximize statistical information.
  4. Corrected Item-Total Correlations: Compute the point-biserial or Pearson correlation between each item and the total scale score excluding that item ($r_{i(t-i)}$). Discard any items exhibiting $r_{i(t-i)} < 0.30$.
  5. Cross-Validation of Factorial Structure: Perform Exploratory Factor Analysis (EFA) on an initial calibration sample, followed by Confirmatory Factor Analysis (CFA) on an entirely independent validation sample to verify configural and metric stability.

§3.11 Unit 3 Methodological Synthesis and Review

Measurement theory provides the rigorous quantitative bridge between unobservable scientific constructs and empirical manifest data:

  • Match statistical analysis procedures strictly to the Stevens scale level of measurement (nominal, ordinal, interval, ratio).
  • Understand the foundational decomposition of Classical Test Theory ($\sigma_X^2 = \sigma_T^2 + \sigma_E^2$) and recognize that random measurement error severely attenuates observed correlations.
  • Move beyond the limitations of Cronbach's alpha by verifying essential $\tau$-equivalence, reporting McDonald's $\omega_{\text{total}}$, and evaluating dimensionality via factor analysis.
  • Deploy the Campbell-Fiske MTMM matrix to decouple substantive construct variance from common method variance.
  • Leverage modern Item Response Theory (1PL, 2PL, 3PL) and information functions to optimize measurement precision across the trait continuum and verify measurement invariance across demographic groups.
THREE WORKED QUESTIONS

Step-by-Step Statistical Research Methods Solutions

Three original questions connect calculations, definitions and interpretation. Open each solution to follow the reasoning.

Solved Problem Example 3.1: Cronbach's Alpha and Split-Half Prophecy Calculation

A neuropsychological test assessing executive working memory consists of $k = 6$ items. A sample of $n = 500$ clinical participants yields the following item-level variance properties:

  • Sum of individual item variances: $\sum_{j=1}^6 \sigma_j^2 = 7.20$.
  • Total composite test score variance: $\sigma_X^2 = 24.00$.

Furthermore, an alternative psychometrician divides the 6-item test into two 3-item half-tests (odd vs. even items, $X_{\text{odd}}$ and $X_{\text{even}}$). The sample Pearson correlation between the two half-tests is $r_{\text{half}} = 0.65$.

  1. Compute Cronbach's alpha ($\alpha$) for the full 6-item scale.
  2. Compute the split-half reliability of the full test using the Spearman-Brown Prophecy Formula.
  3. If the test developers want to increase the test's reliability to $\rho_{\text{target}} = 0.90$ by adding parallel items, compute the required number of additional items $k_{\text{new}}$.
Compute Cronbach's Alpha

Apply the formula for Cronbach's Alpha:

$$\alpha = \frac{k}{k - 1} \left( 1 - \frac{\sum_{j=1}^k \sigma_j^2}{\sigma_X^2} \right)$$

Substitute $k = 6$, $\sum \sigma_j^2 = 7.20$, and $\sigma_X^2 = 24.00$:

$$\frac{k}{k - 1} = \frac{6}{5} = 1.20$$
$$1 - \frac{7.20}{24.00} = 1 - 0.30 = 0.70$$
$$\alpha = 1.20 \times 0.70 = \mathbf{0.840}$$

Cronbach's alpha for the 6-item working memory scale is $0.840$, indicating good internal consistency under the assumption of essential $\tau$-equivalence.

Compute Full-Test Reliability via Spearman-Brown Prophecy

The observed correlation between the two half-tests is $r_{\text{half}} = 0.65$. To estimate the reliability of the doubled (full) test ($k = 2$):

$$\rho_{\text{full}} = \frac{2 r_{\text{half}}}{1 + r_{\text{half}}} = \frac{2 \times 0.65}{1 + 0.65} = \frac{1.30}{1.65} \approx \mathbf{0.7879}$$

The split-half estimate of reliability is $0.788$. The difference between $\alpha = 0.840$ and $\rho_{\text{split}} = 0.788$ may reflect slight disparities in item difficulty or variance between odd and even halves.

Calculate the Required Test Length to Achieve Target Reliability

We wish to increase reliability from the current baseline $\rho_{11} = 0.840$ to $\rho_{kk} = 0.90$. The general Spearman-Brown formula is:

$$\rho_{kk} = \frac{k \rho_{11}}{1 + (k - 1)\rho_{11}}$$

Solve algebraically for the lengthening multiplier $k$:

$$\rho_{kk} [1 + (k - 1)\rho_{11}] = k \rho_{11}$$
$$\rho_{kk} + k \rho_{kk} \rho_{11} - \rho_{kk} \rho_{11} = k \rho_{11}$$
$$\rho_{kk} (1 - \rho_{11}) = k \rho_{11} (1 - \rho_{kk})$$
$$k = \frac{\rho_{kk} (1 - \rho_{11})}{\rho_{11} (1 - \rho_{kk})}$$

Substitute $\rho_{11} = 0.840$ and target $\rho_{kk} = 0.90$:

$$k = \frac{0.90 \times (1 - 0.840)}{0.840 \times (1 - 0.90)} = \frac{0.90 \times 0.160}{0.840 \times 0.10} = \frac{0.144}{0.084} \approx 1.7143$$

The test must be lengthened by a factor of $1.7143$. The total number of required items is:

$$N_{\text{items}} = 6 \times 1.7143 = 10.286 \implies \mathbf{11 \text{ items}}$$

The test developers must add $11 - 6 = \mathbf{5 \text{ additional parallel items}}$ to achieve a reliability exceeding $0.90$.

Solved Problem Example 3.2: Attenuation Correction in Bivariate Correlation

An epidemiologist investigates the relationship between Chronic Psychological Stress ($X$) and Inflammatory Serum Biomarker Interleukin-6 ($Y$). Because both constructs are measured with imperfect instruments:

  • Psychological Stress is measured via a self-report questionnaire with reliability $\rho_{XX'} = 0.64$.
  • Interleukin-6 is measured from a single-tube blood draw with laboratory assay reliability $\rho_{YY'} = 0.49$.
  • The sample Pearson correlation observed between the manifest scores is $r_{XY} = 0.28$.

The funding review committee argues that a correlation of $0.28$ (explaining only $r^2 = 7.84\%$ of variance) is clinically trivial and uninteresting.

  1. State Spearman's formula for correction for attenuation (disattenuated correlation).
  2. Compute the estimated true correlation $r_{T_X T_Y}$ between the latent constructs free from measurement error.
  3. Compute the percentage of variance in true latent Interleukin-6 explained by true latent Psychological Stress, and evaluate the committee's conclusion.
State Spearman's Attenuation Theorem

In 1904, Charles Spearman proved that random measurement error always biases (attenuates) observed bivariate correlations toward zero. Let $X = T_X + E_X$ and $Y = T_Y + E_Y$, with uncorrelated errors $\mathrm{Cov}(E_X, E_Y) = \mathrm{Cov}(T_X, E_Y) = \mathrm{Cov}(T_Y, E_X) = 0$. The expected observed correlation is:

$$\rho_{XY} = \frac{\mathrm{Cov}(T_X + E_X, T_Y + E_Y)}{\sigma_X \sigma_Y} = \frac{\mathrm{Cov}(T_X, T_Y)}{\sigma_X \sigma_Y}$$

Multiply and divide by $\sigma_{T_X} \sigma_{T_Y}$:

$$\rho_{XY} = \frac{\mathrm{Cov}(T_X, T_Y)}{\sigma_{T_X} \sigma_{T_Y}} \times \frac{\sigma_{T_X}}{\sigma_X} \times \frac{\sigma_{T_Y}}{\sigma_Y} = \rho_{T_X T_Y} \sqrt{\rho_{XX'}} \sqrt{\rho_{YY'}}$$

Solving for the true latent correlation yields the Disattenuated Correlation:

$$r_{T_X T_Y} = \frac{r_{XY}}{\sqrt{\rho_{XX'} \rho_{YY'}}}$$
Compute the Disattenuated Correlation

Substitute the empirical parameters:

  • Observed correlation $r_{XY} = 0.28$
  • Stress scale reliability $\rho_{XX'} = 0.64 \implies \sqrt{\rho_{XX'}} = \sqrt{0.64} = 0.80$
  • Biomarker assay reliability $\rho_{YY'} = 0.49 \implies \sqrt{\rho_{YY'}} = \sqrt{0.49} = 0.70$

Compute the denominator:

$$\sqrt{\rho_{XX'} \rho_{YY'}} = 0.80 \times 0.70 = 0.56$$

Compute the disattenuated correlation:

$$r_{T_X T_Y} = \frac{0.28}{0.56} = \mathbf{0.500}$$

The true correlation between the latent biological constructs is $0.500$—nearly double the observed manifest value of $0.28$!

Evaluate Explained Variance and Substantive Impact

Compute the true coefficient of determination:

$$R_{\text{true}}^2 = (r_{T_X T_Y})^2 = (0.500)^2 = \mathbf{0.250} \quad (25.0\%)$$

While the naive observed data suggested stress explained only $7.84\%$ of biomarker variance, in reality, true latent psychological stress accounts for $25.0\%$ of the variance in latent inflammatory response. Over $68\%$ of the apparent unexplained variance in the observed data was an artifact of unreliable measurement instrumentation ($\rho_{XX'} = 0.64$ and $\rho_{YY'} = 0.49$). The committee's dismissal was based on a fundamental measurement failure: confusing instrument noise with biological irrelevance.

Solved Problem Example 3.3: McDonald's Omega from Confirmatory Factor Analysis

A psychiatric research team conducts a Confirmatory Factor Analysis (CFA) on a $4$-item scale measuring Generalized Anxiety Disorder. The standardized factor loadings onto a single latent factor $\xi$ are estimated as:

  • Item 1 (Restlessness): $\lambda_1 = 0.70$
  • Item 2 (Fatigue): $\lambda_2 = 0.50$
  • Item 3 (Difficulty Concentrating): $\lambda_3 = 0.80$
  • Item 4 (Muscle Tension): $\lambda_4 = 0.60$

Assume a standard congeneric model where items are standardized ($\mathrm{Var}(Y_j) = 1.0$), factor variance $\mathrm{Var}(\xi) = 1.0$, and error variances satisfy $\theta_j = 1 - \lambda_j^2$.

  1. Compute the unique error variance $\theta_j$ for each of the four items.
  2. Calculate McDonald's composite reliability $\omega_{\text{total}}$.
  3. Compute Cronbach's $\alpha$ on the model-implied correlation matrix and verify that $\alpha \le \omega_{\text{total}}$.
Calculate Item Error Variances

For standardized congeneric items:

$$\theta_j = 1 - \lambda_j^2$$

Evaluate for each item:

  • $\theta_1 = 1 - (0.70)^2 = 1 - 0.49 = \mathbf{0.51}$
  • $\theta_2 = 1 - (0.50)^2 = 1 - 0.25 = \mathbf{0.75}$
  • $\theta_3 = 1 - (0.80)^2 = 1 - 0.64 = \mathbf{0.36}$
  • $\theta_4 = 1 - (0.60)^2 = 1 - 0.36 = \mathbf{0.64}$

Compute the sum of error variances:

$$\sum_{j=1}^4 \theta_j = 0.51 + 0.75 + 0.36 + 0.64 = \mathbf{2.26}$$
Calculate McDonald's Omega Total

Compute the sum of factor loadings:

$$\sum_{j=1}^4 \lambda_j = 0.70 + 0.50 + 0.80 + 0.60 = \mathbf{2.60}$$

Square the sum of loadings (latent true score variance):

$$\left( \sum_{j=1}^4 \lambda_j \right)^2 = (2.60)^2 = \mathbf{6.76}$$

Total composite score variance is:

$$\sigma_X^2 = \left( \sum_{j=1}^4 \lambda_j \right)^2 + \sum_{j=1}^4 \theta_j = 6.76 + 2.26 = \mathbf{9.02}$$

Compute McDonald's Omega:

$$\omega_{\text{total}} = \frac{\left( \sum \lambda_j \right)^2}{\sigma_X^2} = \frac{6.76}{9.02} \approx \mathbf{0.7494}$$

The composite reliability of the anxiety scale is $0.749$ (or $74.9\%$).

Compute Model-Implied Cronbach's Alpha and Compare

The model-implied correlation between item $j$ and item $l$ is $r_{jl} = \lambda_j \lambda_l$. Compute the off-diagonal covariance sum:

$$\sum_{j \ne l} \sigma_{jl} = \left( \sum_{j=1}^4 \lambda_j \right)^2 - \sum_{j=1}^4 \lambda_j^2$$
$$\sum_{j=1}^4 \lambda_j^2 = 0.49 + 0.25 + 0.64 + 0.36 = 1.74$$
$$\sum_{j \ne l} \sigma_{jl} = 6.76 - 1.74 = 5.02$$

The sum of item variances is $\sum \sigma_j^2 = 4 \times 1.0 = 4.0$. Total scale variance is $\sigma_X^2 = 4.0 + 5.02 = 9.02$. Compute Cronbach's Alpha:

$$\alpha = \frac{k}{k - 1} \left( 1 - \frac{\sum \sigma_j^2}{\sigma_X^2} \right) = \frac{4}{3} \left( 1 - \frac{4.0}{9.02} \right) = 1.3333 \times (1 - 0.44346) = 1.3333 \times 0.55654 \approx \mathbf{0.7421}$$

Observe that:

$$\alpha = 0.7421 < \omega_{\text{total}} = 0.7494$$

Because the factor loadings are heterogeneous ($\lambda \in [0.50, 0.80]$, violating $\tau$-equivalence), Cronbach's alpha strictly under-estimates the true composite reliability. McDonald's omega provides the mathematically accurate representation of scale reliability.