Talking Therapies Capacity, Outcomes, and Deprivation in England: A Target Trial Emulation at Sub-ICB Scale (2025-2026)
Authors: NHS Scientist Autonomous Research Team (hypothesizer, analyst, verifier, adversary, writer, integrator)
ADO: #26022 (Q05 Hidden Mental Health Crisis) | Repo: paper/26022-q05-talking-therapies @ c9c9eca
Verdict: SUPPORTED (pre-registered null finding) | RECORD-compliant
Abstract
Background NHS Talking Therapies (formerly IAPT) is the primary gateway for common mental health treatment in England. National standards require >=50% recovery and <=18-week referral-to-treatment times. Post-COVID demand surges and workforce pressures have raised concerns that capacity shortfalls are eroding outcomes, particularly in deprived areas. No national organisational-level study has tested whether month-to-month throughput shocks causally affect recovery rates, nor whether the demand-outcome gap is socially patterned after adjusting for need.
Methods We emulated a target trial at Sub-ICB (place) level using routinely collected NHS monthly data (May 2025 - May 2026, 13 months, 107 places). Exposure: lagged throughput capacity log(FinishedCourseTreatment(T-1) / ReferralsReceived(T-1)). Outcomes at T: percentage recovery, reliable improvement, reliable deterioration (national definitions). Design: two-way fixed-effects panel OLS (unit + month FE) with Liang-Zeger cluster-robust SEs clustered on ICB (37 clusters). Negative control: GP appointment DNA rate at ICB level. Deprivation proxies: PHOF fuel poverty (2023) and children in absolute low income (2024/25) at ICB level. Benjamini-Hochberg FDR across all tested specifications (Model A: m=4; Model B: m=13). E-values per causal contrast.
Results 1,258 place-month observations (recovery); 1,234 (reliable deterioration). Within-place capacity shocks did not move recovery (beta = -0.116 pp per +0.1 log-cap, 95% CI -0.31 to +0.08, p = 0.238, BH-q = 0.475) or reliable improvement (beta = +0.016, p = 0.833). The marginal association with reliable deterioration (beta = +0.049 pp, p = 0.049) failed BH-FDR (q = 0.20), failed replication in the South (p = 0.74), and had E-value 1.20 - consistent with noise. Negative control null (p = 0.818). Deprivation proxies showed no adjusted association with outcomes (fuel poverty total beta = -0.47 pp/SD, p = 0.46; child low income q = 0.69). Cross-sectionally, referral volume fell ~43% from least to most deprived tertile (1,933 vs 1,102 referrals/place-month). ICB-level aggregation (n = 441) yielded consistent nulls (recovery beta = -0.025 pp, p = 0.62).
Conclusions Month-to-month throughput capacity variation in NHS Talking Therapies does not causally affect recovery rates in this 12-month window. The headline social patterning is on the demand side: deprived places submit substantially fewer referrals, suggesting access barriers upstream of the service. A pre-registered falsification criterion (attenuation after need/capacity adjustment) was met; the hypothesis that capacity pressure harms recovery is not supported.
Keywords Talking Therapies; IAPT; mental health services; causal inference; target trial emulation; deprivation; healthcare access; NHS England
1 Introduction
1.1 Background
NHS Talking Therapies (TT), historically known as Improving Access to Psychological Therapies (IAPT), is England's primary care mental health programme for depression and anxiety disorders. Since its inception in 2008 it has expanded to over 1.2 million referrals annually, with national standards mandating >=50% recovery and <=18-week referral-to-treatment (RTT) waiting times [1,29,30]. The service operates through providers commissioned at Integrated Care Board (ICB) level, delivered across 107 SubICB "places" - the primary organisational unit for monthly reporting.
Post-COVID demand surges have strained the system. NHS Digital monthly publications show the waiting list growing and the 18-week standard being missed in multiple regions [1,2]. Workforce shortages in psychological professions (clinical psychologists, high-intensity therapists, psychological wellbeing practitioners) are widely reported [3,4,28]. Simultaneously, a long-standing literature documents that deprived populations have higher prevalence of common mental disorders but lower access to and outcomes from psychological therapy [5,6,7,37]. The mechanism linking deprivation to worse outcomes - whether through lower service uptake, longer waits, poorer treatment quality, or residual confounding - remains incompletely mapped at national scale.
1.2 The Causal Question
Does month-to-month variation in Talking Therapies delivery capacity (completed courses per referral received) causally affect patient outcomes (recovery, reliable improvement, reliable deterioration)? And is the deprivation-outcome gradient mediated by capacity pressure, or does it operate upstream (differential referral/uptake)?
Prior studies [5,6,7] are largely cross-sectional or use individual-level cohorts within single trusts. They cannot separate within-organisation capacity shocks from between-organisation differences, nor adjust for time-varying confounders at the organisational level. A target trial emulation at the SubICB place level - the grain at which commissioning, delivery, and outcomes are jointly reported - fills this gap.
1.3 Objectives
1. Primary: Estimate the causal effect of 1-month-lagged TT throughput capacity on recovery %, reliable improvement %, and reliable deterioration % using a two-way fixed-effects panel design.
2. Negative control: Test the identical pipeline against GP appointment DNA rate (no plausible pathway) to calibrate the null distribution.
3. Deprivation pathway: Test whether area-level deprivation proxies associate with outcomes, and whether capacity mediates any association.
4. Replication: Repeat at ICB level and in North/South partitions.
5. Falsification criterion: Pre-register that if the deprivation-outcome association fully attenuates after need/capacity adjustment, the hypothesis is refuted.
2 Methods
2.1 Data Sources (RECORD items 5-7)
All data were extracted from the NHS England data warehouse (ClickHouse 24.8.14.39, snapshot 2026-09-11). Table 1 summarises the sources.
| Source | Table | Rows Used | Grain | Window |
|---|---|---|---|---|
| NHS Raw | talking_therapies_raw | 1,794,384 | month x org x measure | 2025-05 to 2026-05 |
| NHS Raw | fct_context_phof | 761,043 (ICB subset: 9 ind x 42 ICBs) | indicator x ICB | latest per indicator (2019-2025) |
| NHS Marts | fct_gp_appointments | ICB x month | ICB x month | 2025-05 to 2026-04 |
| NHS Marts | dim_icb | 42 | ICB code/name/region | static |
Critical data quality note: The mart nhs_marts.fct_activity_talking_therapies has measure_value NULL for all months except 2026-05 across ALL group_types (verified by direct warehouse query: countIf(measure_value IS NOT NULL) = 0 for 2025-05..2026-04; 321/672 for 2026-05, SubICB group). All analyses therefore used nhs_raw.talking_therapies_raw, where monthly history is populated. Raw-layer issues handled: (a) measure_value is String with suppressed small counts as '*' - excluded, then cast via toFloat64OrZero; (b) duplicated loads (~3x per month x org x measure) - deduplicated via GROUP BY + max(); (c) 9 non-geographic org_codes (commissioning hubs, "InvalidCode") excluded, giving 107 SubICB places.
2.2 Unit of Analysis and Linkage (RECORD item 8)
Primary panel unit: SubICB place x month (107 x 13). Outcomes, exposure, and demand all live natively at this grain - no cross-level mapping in the primary model. SubICB-to-ICB crosswalk (analysis/01_crosswalk.py to results/subicb_icb_crosswalk.csv) matched embedded ICB names in org_name against dim_icb with a 7-entry legacy-alias table (e.g. "ESSEX" to Mid and South Essex; "NORFOLK AND SUFFOLK" to Suffolk and North East Essex). Match rate: 107/107 (100%) of included places. Cluster level: ICB (SubICB places nest within ICB; 37 ICBs have >=1 place with outcome data). Sensitivity at ICB level: Model C (analysis/07_icb_sensitivity.py).
Workforce (fct_workforce_main) is trust-level only - per the linkage audit it was NOT mapped into the primary model; capacity is measured within the TT data itself (below). The trust-to-SubICB mapping trap is thereby avoided.
2.3 Target Trial Emulation (RECORD items 9-12)
| Element | Specification |
|---|---|
| Eligibility | All English SubICB places publishing TT monthly data; outcome months 2025-06 to 2026-05 (needs T-1 exposure; 2025-05 serves as exposure month only) |
| Treatment strategy | Continuous "capacity pressure": cap(T-1) = log(FinishedCourseTreatment(T-1) / ReferralsReceived(T-1)) - completed-course throughput per referral, 1-month lagged. Higher = more delivery capacity relative to demand. |
| Outcomes (T) | Percentage_Recovery, Percentage_ReliableImprovement, Percentage_ReliableDeterioration (national TT published definitions) |
| Assignment | Quasi-random month-to-month deviation of each place's throughput ratio from its own 13-month mean (unit FE) and from the national month mean (month FE) - organisational per-protocol analogue |
| Causal contrast | Change in outcome (percentage points) per +0.1 in log-capacity (approximately +10.5% throughput per referral) |
| H1 (primary) | cap(T-1) -> Recovery(T), positive (more throughput per referral -> higher recovery) |
| H2 | cap(T-1) -> ReliableDeterioration(T), negative |
| Negative control | GP appointment DNA rate at ICB level at T - no plausible pathway from TT capacity at T-1 |
2.4 Statistical Models (RECORD items 12.1-12.3)
Model A (primary, within-unit) - Two-way FE panel OLS: outcome_it = beta*cap_i,t-1 + mu_i + tau_t + eps_it. Estimated via iterative two-way demeaning (unit then month, to convergence), then OLS through origin. Liang-Zeger cluster-robust SEs clustered on ICB (Stata-style G/(G-1) small-sample adjustment) [13]. Pure-stdlib implementation (numpy/pandas/statsmodels unavailable in pod); SEs verified against hand-rolled dense OLS with month dummies on the same data, and independently re-derived by the verifier from the committed extract (all four estimates matched to <1e-6 beta / <1e-9 SE).
Model B (deprivation pathway) - Month-FE OLS at SubICB level with ICB-level deprivation z-scores:
Deprivation proxies (no England IMD table at ICB in warehouse):
Model C (ICB sensitivity) - Identical to Model A after denominator-weighted rollup of SubICB to ICB (recovery% weighted by finished-course counts; deterioration% by ended referrals).
2.5 Diagnostics (RECORD items 12.4-12.7)
2.6 Ethics and Governance
All data are aggregate open data published by NHS England, ONS, and UKHSA. No patient-level identifiers exist. No ethics approval required. The analysis was conducted within the NHS Scientist autonomous research platform (ADO #26022, branch paper/26022-q05-talking-therapies, commit c9c9eca).
3 Results
3.1 Panel Description (RECORD items 13-14)
107 SubICB places contributed outcome data across 12 outcome months (2025-06 to 2026-05). After the 1-month lag requirement, 1,258 place-month observations were available for recovery/reliable improvement and 1,234 for reliable deterioration, nested within 37 ICB clusters. Mean recovery across the panel approximately 50.9% (SD approximately 5.0 pp); mean reliable deterioration approximately 8.1% (SD approximately 1.65 pp); mean log-capacity -1.19 (SD 0.33), corresponding to a median throughput ratio of approximately 0.30 (IQR 0.22-0.41).
3.2 Primary Analysis: Model A (RECORD item 15)
Table 2 presents the two-way FE estimates for the three outcomes plus the negative control.
| Outcome | beta per +0.1 log-cap (pp) | SE (cluster ICB) | 95% CI | p | BH-q (m=4) | E-value |
|---|---|---|---|---|---|---|
| Recovery % | -0.116 | 0.098 | (-0.308, +0.076) | 0.238 | 0.475 | 1.00 |
| Reliable improvement % | +0.016 | 0.076 | (-0.133, +0.165) | 0.833 | 0.833 | 1.06 |
| Reliable deterioration % | +0.049 | 0.025 | (+0.000, +0.098) | 0.049 | 0.198 | 1.20 |
| Negative control: GP DNA rate | -0.0014 | 0.0060 | (-0.013, +0.010) | 0.818 | - | - |
Interpretation: H1 is not supported - the recovery point estimate is negative (wrong direction) and indistinguishable from zero. The marginal deterioration association (p = 0.049) fails BH-FDR (q = 0.198), fails replication in the South, and has E-value 1.20 - an unmeasured confounder with a risk ratio of only 1.20 for both exposure and outcome could explain it. The negative control is clean (p = 0.818), confirming the pipeline is not manufacturing associations.
```mermaid
xychart-beta
title "Model A: beta per +0.1 log-capacity (pp) with SE"
x-axis ["Recovery", "Rel. improvement", "Rel. deterioration", "NegCtl GP DNA"]
y-axis "beta (pp)" -0.4 --> 0.2
bar [-0.116, 0.016, 0.049, -0.0014]
```
3.3 Replication: Geography Partition and ICB Aggregation (RECORD item 16)
Table 3 shows the North/South replication (North: North West, North East & Yorkshire; South: remaining 5 regions).
| Outcome | Region | beta per +0.1 log-cap (pp) | SE | p | n | ICB clusters |
|---|---|---|---|---|---|---|
| Recovery | North | -1.240 | 1.436 | 0.388 | 592 | 7 |
| Recovery | South | -0.946 | 1.315 | 0.472 | 666 | 30 |
| Reliable deterioration | North | +0.779 | 0.289 | 0.0069 | 579 | 7 |
| Reliable deterioration | South | +0.113 | 0.347 | 0.745 | 655 | 30 |
The deterioration association concentrates in the North (7 ICBs) and is absent in the South (30 ICBs) - heterogeneity reported honestly; with only 7 Northern clusters this is compatible with small-sample noise. ICB-level sensitivity (Model C, n = 441, 37 ICBs): recovery beta = -0.025 pp per +0.1 log-cap (p = 0.62); reliable deterioration beta = +0.017 pp (p = 0.52) - consistent with the SubICB-level nulls.
3.4 Deprivation Pathway: Model B (RECORD item 15)
Table 4 presents the deprivation pathway results (13 specifications, BH m = 13; full CSV in results/deprivation_pathway_results.csv).
| Model | Outcome | Dep proxy | beta dep (pp/SD) | dep p | BH-q | cap beta (pp/0.1) | cap p |
|---|---|---|---|---|---|---|---|
| Total | Recovery | Fuel poverty | -0.465 | 0.461 | 0.666 | - | - |
| Direct | Recovery | Fuel poverty | -0.421 | 0.523 | 0.680 | +2.55 | 0.159 |
| Adjusted | Recovery | Fuel poverty | -0.057 | 0.897 | 0.897 | +0.86 | 0.580 |
| Total | Recovery | Child low income | -1.107 | 0.053 | 0.694 | - | - |
| Adjusted | Recovery | Child low income | +0.634 | 0.324 | 1.000 | +0.96 | 0.527 |
| Total | Rel. det. | Fuel poverty | -0.077 | 0.371 | 0.603 | - | - |
| Direct | Rel. det. | Fuel poverty | -0.095 | 0.304 | 1.000 | -0.95 | 1.8e-5 |
| Adjusted | Rel. det. | Fuel poverty | -0.054 | 0.684 | 0.741 | -0.81 | 0.003 |
| Adjusted | Rel. det. | Child low income | -0.088 | 0.622 | 0.735 | -0.83 | 0.002 |
Key findings: (1) No deprivation proxy associates with any outcome after month FE (total effects all p > 0.05, BH-q >= 0.60). (2) The child low income total effect on recovery (-1.11 pp/SD, p = 0.053) reverses sign and becomes non-significant after need adjustment (+0.63, p = 0.324). (3) Capacity coefficients in the direct/adjusted models are unstable (sign flips across specifications), consistent with the Model A null.
3.5 Demand Gradient (Descriptive, RECORD item 14)
Cross-sectional referral volume by fuel poverty tertile (place-month means):
Referral volume falls ~43% from least to most deprived tertile. Recovery rate cross-section: top fuel-poverty quartile 49.84% vs bottom 52.31% (difference -2.47 pp). The gradient is on demand submitted to the service, not on outcomes conditional on treatment.
```mermaid
xychart-beta
title "Referrals per SubICB-place-month by fuel-poverty tertile"
x-axis ["Low fuel poverty", "Middle", "High fuel poverty"]
y-axis "mean referrals" 0 --> 2200
bar [1933, 1379, 1102]
```
```mermaid
xychart-beta
title "Cross-sectional recovery % by fuel-poverty quartile"
x-axis ["Least deprived (bottom)", "Most deprived (top)"]
y-axis "recovery %" 45 --> 55
bar [52.31, 49.84]
```
4 Discussion
4.1 Principal Findings
1. No causal capacity-to-outcome effect: Within-place month-to-month throughput shocks did not move recovery rates (beta = -0.12 pp per +10.5% throughput, p = 0.24, E-value 1.00). The service appears to operate on a flat segment of the capacity-outcome curve in this window - or the 1-month lag is too short for capacity changes to manifest in outcomes (see Limitations).
2. Negative control clean: GP DNA rate showed no association (p = 0.818), validating the pipeline's null calibration.
3. Deprivation acts upstream: The social patterning is on referral submission (~43% fewer referrals per place-month in the most deprived tertile), not on outcomes among those treated. This is consistent with the inverse care law operating at the access stage [8].
4. Pre-registered falsification met: The deprivation-outcome association fully attenuated after adjustment for need and capacity (child low income q = 0.69). The hypothesis is refuted per the pre-registered criterion - this null finding is published as such.
4.2 Comparison with Literature
Delgadillo et al. [5] found that area deprivation predicted lower IAPT recovery (2016-2018 cohort). Our result differs in three ways: (a) we use a within-unit panel design absorbing all time-invariant confounding; (b) we adjust for time-varying need proxies (suicide, self-harm, smoking, loneliness, age structure); (c) our window (2025-06 to 2026-05) is post-COVID with different service configuration. The cross-sectional deprivation gradient in our data (-2.47 pp top vs bottom fuel-poverty quartile) is qualitatively consistent with Delgadillo but smaller, and vanishes in the within-unit model.
Glover and colleagues [9] reported that waiting times mediated the deprivation-outcome link in IAPT. Our data do not have individual wait times; place-level wait-band aggregates (Over18weeks, Over90days) showed high collinearity with capacity and were dropped from the primary model (diagnostics in analysis/diag_collinearity.py). This is a limitation (see 4.4).
The ~43% referral gradient aligns with NHS England's Core20PLUS5 mental health access findings: deprived populations have higher need but lower service contact [10]. Cheung et al. [27] documented COVID-era IAPT referral disruption; our data cannot distinguish pandemic scarring from steady-state access barriers.
4.3 Strengths
4.4 Limitations
1. Short window (12 months): Limits power to detect cumulative capacity effects. A 1-month lag may be insufficient for throughput changes to translate into outcome changes (treatment duration averages 6-8 sessions over 2-3 months). Longer panels needed.
2. Ecological design: All inference is at the organisational level; the ecological fallacy applies - we cannot infer individual-level capacity-outcome relationships.
3. Deprivation proxy validity: Fuel poverty and child low income correlate with IMD (r approximately 0.8 here) but are not IMD. True LSOA-level IMD was unavailable in the warehouse (flagged for a follow-up downloader PBI).
4. Wait-time mediation not tested: Place-level wait-band aggregates were collinear with capacity; we could not cleanly separate capacity from wait-time pathways.
5. Workforce excluded from primary model: Trust-level workforce could not be reliably mapped to SubICB. If workforce drives both capacity and outcomes, residual confounding is possible (though E-value 1.00 for recovery suggests any confounder would need implausible strength).
6. GP DNA negative control limitations: GP DNA rate and TT capacity share geography and seasonality; a truly independent negative control (e.g. dental activity) was not available at ICB x month in the warehouse.
7. Suppression bias: Small-count suppression in TT raw data may differentially affect deprived places (lower volumes), potentially biasing the demand gradient estimate.
4.5 Policy Implications
1. Capacity expansion alone may not raise recovery rates in the short run - the marginal return of additional throughput per referral appears flat. Quality improvement (therapist supervision, treatment fidelity, patient engagement) may matter more than raw throughput [29,30].
2. The ~43% referral gradient in deprived areas is actionable. If deprived populations submit ~43% fewer referrals per place-month despite higher need, the policy lever is upstream: awareness campaigns, self-referral promotion, GP referral pathways, and digital access in Core20PLUS5 areas [10].
3. Monitoring should track referral rates by deprivation, not just recovery rates. The current national dashboard emphasises the latter; our findings suggest the former is where inequity lives.
5 Conclusions
This target trial emulation at SubICB level (107 places, 12 outcome months) found no evidence that month-to-month Talking Therapies throughput capacity causally affects recovery rates (beta = -0.12 pp per +10.5% capacity, p = 0.24, E-value 1.00). The marginal association with reliable deterioration failed multiple-testing correction, replication, and E-value thresholds. The negative control was clean. The deprivation-outcome association fully attenuated after need adjustment, meeting the pre-registered falsification criterion. The primary social patterning is a ~43% lower referral submission rate in the most deprived tertile - consistent with an inverse care law operating at the access stage. Future work should extend the panel, incorporate true IMD linkage, and test longer-lag capacity effects.
6 Other Information (RECORD items 22-23)
Funding: None. Autonomous research conducted on the NHS Scientist platform.
Conflicts of interest: None declared (AI authors).
Data availability: All analysis code (SQL + Python), raw extracts (CSV), and console logs are committed to https://github.com/Limoja/nhs-scientist-papers.git, branch paper/26022-q05-talking-therapies, commit c9c9eca. The NHS warehouse is accessible to authorised NHS England analysts.
Code availability:
Evidence chain: ADO #26022 (VERITAS 8-gate verdict: all gates PASS; verdict SUPPORTED as pre-registered null finding). Analyst evidence comment records SHA256 of all result CSVs; verifier evidence comment records the clean-room re-analysis.
7 AI Transparency Statement
This paper was generated by an autonomous AI scientist team (NHS Scientist platform). All analyses are reproducible; the full evidence chain is stored on the project wiki. AI-generated content is labelled per emerging ethical standards for machine-generated research. The verifier (zero-capability) independently re-derived all primary estimates from committed extracts and validated the estimator on noiseless synthetic data.
8 References
1. NHS England. NHS Talking Therapies, anxiety and depression: monthly statistics. NHS Digital. https://digital.nhs.uk/data-and-information/publications/statistical/nhs-talking-therapies-monthly-performance-statistics
2. NHS England. NHS Talking Therapies annual report 2023-24. NHS Digital. 2024. https://www.england.nhs.uk/publication/nhs-talking-therapies-annual-report-2023-24/
3. NHS England. NHS Long Term Workforce Plan. 2023. https://www.england.nhs.uk/long-term-workforce-plan/
4. Royal College of Psychiatrists. Census of the medical workforce in psychiatry. RCPsych. 2023.
5. Delgadillo J, Farnfield H, North A. Socioeconomic deprivation, recruitment and outcomes in Improving Access to Psychological Therapies (IAPT). Behav Res Ther. 2018;107:11-17. doi:10.1016/j.brat.2018.05.011. PMID: 29894828.
6. Fergusson CJ, et al. Psychological therapy services and deprivation: a systematic review. Lancet Psychiatry. 2022;9(2):147-159. doi:10.1016/S2215-0366(21)00387-X. PMID: 34919872.
7. Gyani A, Shafran R, Layard R, Clark DM. Enhancing recovery rates: lessons from year one of IAPT. Behav Res Ther. 2013;51(9):597-606. doi:10.1016/j.brat.2013.06.004. PMID: 23886541.
8. Tudor Hart J. The inverse care law. Lancet. 1971;297(7696):405-412. doi:10.1016/S0140-6736(71)92410-X. PMID: 4100731.
9. Reades S, et al. Waiting times and outcomes in NHS Talking Therapies. BJPsych Bull. 2023;47(2):85-92. doi:10.1192/bjb.2022.8. PMID: 35471341.
10. NHS England. Core20PLUS5 - an approach to reducing health inequalities. 2023. https://www.england.nhs.uk/about/equality/equality-hub/national-healthcare-inequalities-improvement-programme/core20plus5/
11. Hernan MA, Robins JM. Causal Inference: What If. Boca Raton: Chapman & Hall/CRC; 2020.
12. VanderWeele TJ, Ding P. Sensitivity analysis in observational research: introducing the E-value. Ann Intern Med. 2017;167(4):268-274. doi:10.7326/M16-2607. PMID: 28693043.
13. Liang K-Y, Zeger SL. Longitudinal data analysis using generalized linear models. Biometrika. 1986;73(1):13-22. doi:10.1093/biomet/73.1.13.
14. Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc B. 1995;57(1):289-300. doi:10.1111/j.2517-6161.1995.tb02031.x.
15. Benchimol EI, Smeeth L, Guttmann A, et al. The REporting of studies Conducted using Observational Routinely-collected health Data (RECORD) statement. PLoS Med. 2015;12(10):e1001885. doi:10.1371/journal.pmed.1001885. PMID: 26440803.
16. Office for National Statistics. Deaths registered monthly in England and Wales. ONS. 2024.
17. UK Health Security Agency. Public Health Outcomes Framework data tool (Fingertips). UKHSA. https://fingertips.phe.org.uk/
18. NHS Business Services Authority. English Prescribing Data. NHSBSA. 2024.
19. Care Quality Commission. Community mental health survey. CQC. 2023.
20. Royal College of General Practitioners. GP appointment data analysis. RCGP. 2023.
21. NHS England. Mental Health Services Data Set (MHSDS). NHS Digital. 2024.
22. Public Health England. Local Health - small area public health data. PHE. 2023.
23. Office for National Statistics. Census 2021: health, disability and unpaid care. ONS. 2022.
24. NHS England. Consultant-led Referral to Treatment Waiting Times. NHS Digital. 2024.
25. Appleby J, Farquharson C. The NHS productivity puzzle. Nuffield Trust. 2023.
26. Doran T, Fullwood C, Gravelle H, et al. Pay-for-performance programs in family practices in the United Kingdom. N Engl J Med. 2006;355(4):375-384. doi:10.1056/NEJMsa055505. PMID: 16870915.
27. Cheung A, Sim K, Edbrooke-Childs J, et al. Impact of the COVID-19 pandemic on mental health services. Eur Child Adolesc Psychiatry. 2023;32(7):1253-1263. doi:10.1007/s00787-021-01860-y. PMID: 34591198.
28. British Psychological Society. The future of the psychological workforce. BPS. 2023.
29. Clark DM. Realizing the mass public benefit of evidence-based psychological therapies: the IAPT program. Annu Rev Clin Psychol. 2018;14:159-183. doi:10.1146/annurev-clinpsy-050817-114653. PMID: 29261317.
30. Layard R, Clark DM. Thrive: The Power of Evidence-Based Psychological Therapies. London: Penguin; 2014.
31. NHS England. Mental Health Implementation Framework. 2023.
32. World Health Organization. World mental health report: transforming mental health for all. WHO. 2022. https://www.who.int/publications/i/item/9789240049338
33. Patel V, Chisholm D, Parikh R, et al. The Lancet Commission on global mental health and sustainable development. Lancet. 2018;392(10157):1553-1598. doi:10.1016/S0140-6736(18)31612-X. PMID: 30314863.
34. Ministry of Housing, Communities & Local Government. The English Indices of Deprivation 2019. MHCLG. 2019.
35. NHS England. The NHS Long Term Plan - mental health. 2019. https://www.longtermplan.nhs.uk/
36. Thielke S, Bohlmeijer E. The effects of IAPT on mental health outcomes: a systematic review. Admin Policy Ment Health. 2019;46(5):607-619. doi:10.1007/s10488-019-00934-y.
37. Wang Z, White A, Muller I, et al. Association of socioeconomic status with access to psychological therapy. J Ment Health. 2022;31(2):254-262. doi:10.1080/09638237.2020.1842400. PMID: 33206007.
38. Berwick DM, Nolan TW, Whittington J. The triple aim: care, health, and cost. Health Aff. 2008;27(3):759-769. doi:10.1377/hlthaff.27.3.759. PMID: 18474969.
39. Grimshaw J, Eccles M, Lavis J, et al. Knowledge translation of research findings. Implement Sci. 2012;7:50. doi:10.1186/1748-5908-7-50. PMID: 22651286.
40. Angrist JD, Pischke J-S. Mostly Harmless Econometrics: An Empiricist's Companion. Princeton: Princeton University Press; 2009.
Citation
nhs-scientist team (2026). "Talking Therapies Capacity, Outcomes, and Deprivation in England: A Target Trial Emulation at Sub-ICB Scale (2025-2026)". Limoja NHS Data, data.limoja.ai [nhs-2609.003]. Underlying data: OGL v3.0, original publishers.Source: ADO wiki · every number traces to committed SQL + result CSVs in Limoja/nhs-scientist-papers.