Intended for healthcare professionals

CCBYNC Open access

Association between home insulation and hospital admission rates: retrospective cohort study using linked data from a national intervention programme

BMJ 2020; 371 doi: (Published 29 December 2020) Cite this as: BMJ 2020;371:m4571

Linked Editorial

Housing as a public health investment

  1. Caroline Fyfe, doctoral student1,
  2. Lucy Telfar,
  3. Barnard, senior research fellow1,
  4. Philippa Howden-Chapman, professor1,
  5. Jeroen Douwes, professor2
  1. 1He Kianga Oranga, Department of Public Health, University of Otago Wellington, PO Box 7343, Wellington South 6242, New Zealand
  2. 2Centre for Public Health Research, Massey University Wellington, New Zealand
  1. Correspondence to: C Fyfe fyfca725{at}
  • Accepted 21 October 2020


Objectives To investigate whether retrofitting insulation into homes can reduce cold associated hospital admission rates among residents and to identify whether the effect varies between different groups within the population and by type of insulation.

Design A quasi-experimental retrospective cohort study using linked datasets to evaluate a national intervention programme.

Participants 994 317 residents of 204 405 houses who received an insulation subsidy through the Energy Efficiency and Conservation Authority Warm-up New Zealand: Heat Smart retrofit programme between July 2009 and June 2014.

Main outcome measure A difference-in-difference approach was used to compare the change in hospital admissions of the study population post-insulation with the change in hospital admissions of the control population that did not receive the intervention over the same two timeframes. Relative rate ratios were used to compare the two groups.

Results 234 873 hospital admissions occurred during the study period. Hospital admission rates after the intervention increased in the intervention and control groups for all population categories and conditions with the exception of acute hospital admissions among Pacific Peoples (rate ratio 0.94, 95% confidence interval 0.90 to 0.98), asthma (0.92, 0.86 to 0.99), cardiovascular disease (0.90, 0.88 to 0.93), and ischaemic heart disease for adults older than 65 years (0.79, 0.74 to 0.84). Post-intervention increases were, however, significantly lower (11%) in the intervention group compared with the control group (relative rate ratio 0.89, 95% confidence interval 0.88 to 0.90), representing 9.26 (95% confidence interval 9.05 to 9.47) fewer hospital admissions per 1000 in the intervention population. Effects were more pronounced for respiratory disease (0.85, 0.81 to 0.90), asthma in all age groups (0.80, 0.70 to 0.90), and ischaemic heart disease in those older than 65 years (0.75, 0.66 to 0.83).

Conclusion This study showed that a national home insulation intervention was associated with reduced hospital admissions, supporting previous research, which found an improvement in self-reported health.


Excess winter mortality and morbidity are often greatest in countries with relatively mild climates, including New Zealand, probably because of poorer thermal housing standards compared with well insulated houses in colder climates.123 Cold houses are also associated with indoor dampness and mould, both consistently linked to an increased risk of respiratory symptoms, with strong evidence for asthma exacerbations and respiratory infections, and limited evidence for asthma development.456 Several intervention studies involving improved thermal insulation or heat supply, or both have shown a reduction in respiratory symptoms and improved self-reported health, but these studies were either qualitative78 or did not have sufficient statistical power to obtain significant results for hospital admissions.9 Taken together, the evidence suggests a role for indoor cold in increased mortality and morbidity, with improved insulation having the potential to mitigate these effects.10111213 As relatively few intervention studies or natural experiments of sufficient size have been conducted, however, a gap remains in quantitative evidence for the role of insulation in mitigating the health impacts of cold housing.

Most homes in New Zealand are cold by international standards and do not meet the World Health Organization’s recommended minimum indoor temperature of 18°C.1011121415 Most of the housing stock was built between 1950 and 2000 and consists of timber framed, single storey, detached houses with either a brick or a weatherboard skin.1617 Before a 1977 amendment to the building code (NZS 421P) there was no requirement to insulate homes. Since then, new homes are required to be insulated, but retrofitting insulation into existing homes has been at the discretion of the homeowner. In 2005 an estimated 60% of houses built before 1979 lacked any form of ceiling insulation, and 68% were without underfloor insulation. In addition, 37% of houses built in the 1980s and 26% built in the 1990s had levels of ceiling insulation that were insufficient to meet even the 1977 standard.16

In 2009, the New Zealand Energy Efficiency and Conservation Authority established an insulation subsidy programme for existing homes, Warm-up New Zealand: Heat Smart (WUNZ). This programme ran from June 2009 to July 2014 and represented a government investment of NZ$347m (£181m; $241m; €203m). As insulation levels had been revised several times since 1977, houses that received a WUNZ subsidy needed to retrofit insulation that conformed to the 2007 insulation standard (NZBC clause H1/AS1).13 This standard sets out the requirements for the thermal resistance (R value) of insulation according to the climate zone in which the house is situated (see fig 1). All owners of houses built before 2000 were eligible; however, the level of the subsidy was based on two income categories—those on a low income received a 60% subsidy, whereas those on a general income received a 33% subsidy, capped at NZ$1300 in each case. At the time the WUNZ programme operated, the average cost of insulating a home was estimated to be around NZ$3000. Some local authorities offered interest-free loans (repaid through rates) to cover the balance of the cost.18 A heater subsidy for houses with no fixed form of heating, or to replace an open fire or inefficient log burner, was initially included as part of the overall WUNZ budget. This subsidy was set at NZ$500 for households on a general income and NZ$1200 for those on a low income.18 It was phased out, ending in December 2012, after an economic evaluation of the programme concluded that insulating houses provided greater benefit than subsidising heaters.19

Fig 1
Fig 1

Map of New Zealand climate zones. R value is the thermal resistance of insulation set by the 2007 insulation standard

To assess whether retrofitting insulation was associated with reduced hospital admissions, we linked a dataset of more than 200 000 houses insulated over the five years of the WUNZ subsidy programme to health data.


The study involved a retrospective cohort of all residents in New Zealand whose houses were insulated between 1 July 2009 and 31 December 2011 through the WUNZ subsidy programme. Residents of these houses were identified by address from primary health organisation records. These records include everybody registered with a general practitioner, and they are updated every quarter. We collected data for residents of WUNZ houses for three years before and three years after the insulation intervention. A difference-in-difference approach was used to compare the before and after intervention change in hospital admissions of the study population with that of a control population over the same period. To ensure that the control group was similar to the intervention group, we included all residents whose houses were insulated through the WUNZ subsidy programme, but at a later date, between 1 January 2012 and 30 June 2014.

The intervention included ceiling insulation, underfloor insulation, and a top-up of existing ceiling insulation when it fell short of the 2007 standard. The type and amount of insulation installed depended on what already existed in the house and access to spaces requiring insulation. The programme did not extend to provision of double glazing or cavity wall insulation; both are rare in houses built before 2007. In a limited number of houses (11.6%) a clean heating device was also installed. Most of these (78.0%) were heat pumps, with a further 19.7% comprising energy efficient log burners. The remainder were made up of wood pellet burners (1.4%) and flued gas heaters (1.0%).

Data linkage

We used a WUNZ dataset as the basis for this study. The dataset contained records of all houses that received an insulation subsidy between July 2009 and June 2014, including address, date of the insulation retrofit, type and amount of insulation provided, and clean heating device installed, when applicable. We geocoded addresses using ArcGIS10.3 (a spatial analysis software package that can convert address information to a coordinate point on a map) to match with a Land Information New Zealand registered address. We then geocoded primary health organisation records for June 2006 to July 2014, matching them to the same Land Information New Zealand registered address format; and joined them to the WUNZ dataset using the address string (fig 2).

Fig 2
Fig 2

Address matching process to link WUNZ houses to residents using primary health organisation records. *A person-house pair represents a record for each person in each house, if a person lived in more than one WUNZ house, then a new record is created for each house they inhabited. WUNZ=Warm-up New Zealand: Heat Smart programme; PHO=primary health organisation; EECA=Energy Efficiency and Conservation Authority

We assigned intervention and control groups and baseline and follow-up timeframes from the joined records. We collected baseline data for everyone who lived in the house for the three years before the installation of insulation, with the comparison baseline for the control group being between three and six years before insulation was installed. Follow-up included everyone who lived in the house during the three years after the insulation had been installed or until July 2014, if that was sooner. The comparison follow-up for the control group was the three years before insulation was installed (fig 3).

Fig 3
Fig 3

Method for assignment of intervention and control groups and baseline and follow-up timeframes

Person years contributed were calculated by the number of quarter years an individual lived in a WUNZ house (person-house pair). If a person moved from one WUNZ house to another, a new person-house pair was created. A nested cohort of those with six years continual occupancy (three years before insulation and three years after insulation, or six years before insulation in the control group) was also identified.

We obtained personal and health information from the New Zealand Ministry of Health for residents we identified as having lived in a WUNZ house during the study period. These were matched to each person-house linked record using the resident’s national health index number. This number is assigned to each person who receives healthcare in New Zealand. Geographical area and area level deprivation, based on the New Zealand index of deprivation 2013,20 were linked to the joined datasets using the meshblock in which the house was situated. A meshblock is the smallest geographical unit for which statistical data are available.


We excluded people aged 0 years or less or 90 years or more at baseline and when deprivation 10th or sex were unknown or unassigned, as the data quality for these records was deemed inadequate.21 Records were also excluded when there was a gap in continuous residence of more than six months, owing to uncertainty about whether the address was a permanent residence (fig 2). Participants who died during the study were also excluded.

Data analysis

Total person years and hospital admissions contributing to the study were calculated. We also analysed subgroups of the population stratified by age, sex, ethnicity, area level deprivation (10ths aggregated into fifths for the analysis), climate zone, and type of insulation installed. All data analyses were conducted using SAS 9.4.

Poisson regression analysis was used to calculate baseline and follow-up hospital admission rates (number of hospital admissions divided by number of person years), determine rate ratios (hospital admission rate at follow-up divided by hospital admission rate at baseline) for the intervention and control groups, and to calculate relative rate ratios (intervention rate ratio divided by control rate ratio). Furthermore, we undertook Poisson regression analysis for the population subgroups. The use of person years rather than individuals as a denominator; variations in sex, ethnicity, and age within households; and changes in household composition over the six years covered by the study meant it was not possible to adjust for clustering in analyses involving the total population. To determine the impact on results of clustering by household we used a separate analysis undertaken on the six year continual occupancy nested cohort.

Absolute difference in the number of hospital admissions per 1000 population was measured by subtracting the change in adjusted hospital admission rates between baseline and follow-up in the intervention group from the change in adjusted hospital admission rates between baseline and follow-up in the control group.

Analysis of hospital admissions was undertaken based on a list of conditions identified by a previous study as housing related potentially avoidable hospital admissions.22 Because WUNZ did not include interventions to address household crowding that were provided in that study, we further refined this list to include only conditions associated with low ambient temperature, based on an evaluation of literature from more recent housing studies.132324 The final list of cold associated hospital admissions selected for analysis was: acute hospital admissions (non-elective admissions: excluding maternity and injury related), respiratory disease (international classification of diseases, 10th revision (ICD-10) codes J00–J99), cardiovascular disease (I00–I99) in people aged 15–65 years and in adults older than 65 years, upper respiratory tract infection (J00–J06), lower respiratory tract infection (J09–J16, J18, J20–J22), asthma (J45–J46), chronic obstructive pulmonary disease (J40–44, J47) in adults older than 65 years, and ischaemic heart disease (I20–I25) in adults older than 65 years.

Unless otherwise stated, analyses were conducted for all ages—children younger than 15 years, adults aged 15-65 years, and adults older than 65 years. We also conducted a sensitivity analysis, excluding Christchurch City and Waimakariri, to control for potential confounding by the earthquakes that occurred in Christchurch in September 2010 and February 2011. Finally, as participants who died had higher hospital admission rates, we compared results obtained with and without those who died during the study.

Public and patient involvement

No patients were involved in setting the research question or the outcome measures, nor were they involved in developing plans for recruitment, design, or implementation of the study. No patients were asked to advise on interpretation or writing up of results. There are no plans to disseminate the results of the research to study participants or the relevant patient community.


Overall, 1 611 182 person-house pairs, representing 204 405 houses, 994 317 people, were included in the study. The nested cohort with six years continual occupancy consisted of 15.9% of people (n=157 338), 24.2% of houses (n=49 457), and 9.8% of person years (n=1 027 881).

A total of 699 970 hospital admissions were identified in the study population between July 2006 and June 2014, of which 268 344 (38.3%) occurred during baseline or follow-up. After exclusions, 234 873 hospital admissions remained for analysis, with 35.3% (n=82 958) of these in the six year continual occupancy nested cohort.

Table 1 lists the personal characteristics of the intervention and control groups. The intervention group had a lower proportion of Māori than the control group (16.9% v 20.6%) as well as a lower proportion of Pacific Peoples (6.2% v 8.6%). The intervention group also had fewer people in small area level, socioeconomic most deprived fifth (20.5% v 24.0%) and a greater number in (colder) climate zone 3 (25.1% v 18.3%). A greater number of houses in the intervention group had heating installed (18.2% v 5.0%) and a lower number had whole house (ceiling and underfloor) insulation (49.2% v 56.6%). The six year continual occupancy nested cohort was older than the total population, had particularly low representation among those aged 15 to 45 years, and had a lower proportion of Māori and Pacific Peoples.

Table 1

Personal characteristics of intervention and control groups

View this table:

Hospital admissions by condition

Hospital admission rates after the intervention increased in the intervention and control groups for all population categories and conditions, with the exception of Pacific Peoples (rate ratio 0.94, 95% confidence interval 0.90 to 0.98), asthma (0.92, 0.86 to 0.99), cardiovascular disease (0.90, 0.88 to 0.93), and ischaemic heart disease for adults older than 65 years (total population: 0.79, 0.75 to 0.84; six year continual residency nested cohort: 0.76, 0.65 to 0.90). All reductions in hospital admission rates were in the intervention group.

The overall increase in hospital admission rates was significantly reduced in the intervention group compared with control group (relative rate ratio 0.89, 95% confidence interval 0.88 to 0.90), representing 9.26 (9.05 to 9.47) fewer hospital admissions per 1000 in the intervention population. This effect was more pronounced for respiratory disease (0.85, 0.81 to 0.90) and cardiovascular disease (0.86, 0.82 to 0.90) when analysed separately (fig 4). Among specific conditions, reductions in comparative rates of hospital admission were greatest for asthma in all ages (0.80, 0.70 to 0.90), chronic obstructive pulmonary disease in adults older than 65 years (0.80, 0.70 to 0.91), and ischaemic heart disease in adults older than 65 years (0.75, 0.66 to 0.83) (fig 4). Although the effect did not reach statistical significance in the six year continual residency nested cohort for these more specific disease categories, a reduction in comparative rates of hospital admission remained for all conditions except chronic obstructive pulmonary disease.

Fig 4
Fig 4

Difference in difference in rate ratios for disease specific hospital admissions between intervention and control groups before and after insulation was retrofitted

After adjustment for clustering, using the six year continual residency nested cohort, analyses showed a small impact for all hospital admissions combined, with an increase in the relative rate ratio from 0.89 (0.87 to 0.91) to 0.91 (0.87 to 0.94); the impact increased when results were disaggregated. For all subsequent analyses related to the six year continual occupancy nested cohort, the results are presented after adjustment for clustering (table 2). Repeating analyses including participants who had died during the study did not appreciably change the results: the relative rate ratios for analyses on overall hospital admission rates without and with deaths changed from 0.89 (0.87 to 0.91) to 0.91 (0.90 to 0.93).

Table 2

Relative rate ratio for acute hospital admissions* by population subgroup

View this table:

Sociodemographic characteristics

Pacific Peoples experienced a greater comparative reduction in hospital admission rates (0.82, 0.77 to 0.88) than Māori (0.93, 0.89 to 0.97); however, when the six year continual occupancy nested cohort was examined, the effect size was similar for both: Pacific Peoples (0.86, 0.75 to 0.98) and Māori (0.86, 0.77 to 0.95).

When age groups were analysed separately, differences in relative rates of hospital admissions were small. However, owing to higher hospital admission rates among people aged 65 years or older, the absolute reduction in hospital admissions per 1000 population was far higher for this group (28.62, 27.32 to 30.57).

Little difference was found in the relative rate ratio between deprivation fifths. However, owing to increases in hospital admission rates as levels of deprivation increased, the absolute reduction in hospital admissions in the intervention group was higher in the fourth fifth of deprivation 12.61 (11.89 to 13.18) and in the fifth fifth 11.70 (11.03 to 12.23) per 1000 population compared with the first fifth 7.41 (7.07 to 7.83) per 1000 population.

Climate zone

Little difference was found in the relative rate ratios between climate zones. When a sensitivity analysis was conducted to adjust for the impact of the earthquakes in Christchurch in 2010 and 2011, the difference in effect of insulation between the total population (0.87, 0.85 to 0.90) and the six year continual residency nested cohort (0.96, 0.93 to 0.99) in climate zone 3 was no longer evident; total population (0.91, 0.86 to 0.96) and six year continual residency nested cohort (0.91, 0.84 to 0.97).

Intervention type

The most effective intervention types were whole house (ceiling and underfloor) insulation (relative rate ratio 0.88, 0.85 to 0.90) and ceiling insulation alone (0.88, 0.84 to 0.92). Installing heating alongside insulation did not seem to increase the effectiveness of the intervention (0.89, 0.84 to 0.95) (table 2). When underfloor and top-up of existing insulation were installed without other interventions, the effect also failed to reach statistical significance.


This study found a statistically significant reduction in hospital admission rates associated with home insulation when the intervention group was compared with a control group of residents over the same two timeframes who did not receive the intervention. This finding was consistent with previously reported improvements in self-reported health after insulation interventions.789 The reduction in hospital admission rates was more pronounced for specific respiratory conditions. An earlier intervention study in New Zealand also found improvements in parent reported wheeze (adjusted odds ratio 0.51, 0.32 to 0.81) and dry cough (0.50, 0.31 to 0.82) among children with asthma, and a non-statistically significant reduction in hospital admissions for respiratory conditions.25

Vulnerable populations

When population vulnerability was considered—for people who had remained in the same house throughout the study (six year continual occupancy nested cohort), Māori and Pacific Peoples particularly benefitted from having insulation retrofitted. This was also true for Pacific Peoples in the total population. These findings are important as Māori and Pacific Peoples have poorer health outcomes than European and other ethnicities in New Zealand. For example, hospital admission rates in 2013 for all respiratory disease in Māori were twice those for the European and other population (relative risk 2.1), with rates for Pacific Peoples more than two and a half times higher (relative risk 2.6).26

Across deprivation fifths, no pattern to relative improvement in health outcomes was observed. It has been suggested that the size and quality of dwellings are more strongly correlated with the thermal efficiency of housing than is neighbourhood deprivation.27 An analysis of the pathway between housing improvements and health conditions in Glasgow, Scotland concluded that although energy efficiency measures could help manage symptoms, they did not prevent the occurrence of disease.28 In addition, another study cautioned that evaluation of energy efficiency measures provided only a snapshot of historical exposure and that improvements could not mitigate health impacts of cold and damp related exposures built up over a lifetime.8 However, as hospital admission rates increased with increasing levels of deprivation leading to greater absolute reductions in hospital admissions in the highest deprivation fifths, targeting more deprived communities could still provide an overall greater impact on health than untargeted interventions.29

The effect of home insulation was similar across age groups. The greatest reduction in hospital admissions per 1000 population were, however, shown in people aged 65 years or older, because of higher hospital admission rates in this age group. We also found the most pronounced effect for ischaemic heart disease related hospital admissions in those older than 65 years. Older people (≥65 years) have a greater vulnerability to cold housing because of comorbidities and more time spent in the home environment.2730 An analysis of the first year of WUNZ found that older people were over-represented in the study population (compared with the 2006 census), linked to their ability to access capital under the WUNZ subsidy model, and targeting by service providers.31

Intervention type

Some variation was found in effect on hospital admission rates by type of insulation installed. Whole house (ceiling and underfloor) insulation, ceiling insulation, and installation of a heater alongside insulation, had similar statistically significant results for acute hospital admissions in the total population. This finding was consistent with a study in England, which recorded the lowest indoor temperatures among houses without ceiling insulation,27 and a study in Glasgow, Scotland, which identified increased odds of recovery from a circulatory condition when central heating was part of a housing improvement programme.28 The lack of increased effectiveness identified when a heater was included in the intervention reflected findings of an earlier WUNZ study 13; however, current results should be treated with caution given the early phasing out of the heater subsidy. The primary aim of the WUNZ heater subsidy was to provide a clean heating source to replace open fires and inefficient log burners. It was targeted towards reducing outdoor air pollution and was acknowledged at the time as having the potential to provide less warmth than some of the solid fuel burners being replaced.32 In addition, the intervention would only have been effective if households used the new heaters. Several studies found that New Zealanders’ use of home heating was relatively low and often restricted to the main living area.1033 This was particularly true of low income households, who paid the highest proportion of their income on electricity and were more likely to use prepayment metres, something that allowed them to use self-disconnection as a method of managing their energy budgets.34

Limitations of this study

Our study has several limitations. Linking existing datasets involved a quasi-experimental study method that can increase risk of bias due to non-randomisation and an inability to match intervention and control houses.33 This risk was managed by using the difference-in-difference model and sourcing the intervention and control populations from the same programme.3536

As a result of inaccuracies in address data resulting in records being assigned to an incorrect Land Information New Zealand address string, the study was at risk of both type I and type II errors. In addition, patients who did not update their general practitioner within three months of a change of address and residents who were not registered with a general practice would not have been identified through the data linkage process described and would therefore not have been included in the study. Some evidence suggested that those on a low income who rent and recent migrant populations were more likely to not be registered with a general practice, owing to higher levels of mobility and lower levels of contact with health services.21 This could have reduced the strength of the relation between improvements in home insulation and health outcomes across deprivation fifths.

The ending of the heater subsidy in December 2012 limited interpretation of results for the subgroup of households that received heaters alongside insulation. Because of the way the cohort was split between intervention and control groups, using a mid-point of December 2011, most of those who received a heater subsidy were in the intervention group, further complicating comparisons between the two groups.

Although there was no selection or prioritisation process for obtaining a subsidy, the motivation of householders to apply for one might cause unobserved differences between the intervention and control groups. Were the information available and the potential study population large enough, matching by key indicators would have further added to the internal validity of the study.36 We used the six year continual occupancy nested cohort as a comparator, to control for differences between baseline and follow-up and also to determine whether clustering by household influenced the results. These indicated small differences in the relative rate ratio for larger groups but did show greater differences when the datasets were disaggregated, resulting in fewer statistically significant findings for specific disease categories.


The study showed that the retrofitting of home insulation was associated with a reduction in frequency of hospital admissions. The size of the reduction was more pronounced for respiratory disease. Future studies should examine in more detail the differences in the effectiveness of insulation among different subgroups within the population, to better target interventions to those who would benefit most from them.

The impact of the intervention on cold associated hospital admissions varied with the type of insulation retrofitted, and no additional benefit was indicated from combining insulation and a heater. This might have been due to the low numbers in this group, particularly among the control population. However, concerns about the energy costs of running heaters and low expectations of warmth among residents might have led to under-utilisation of the heaters provided. Further study should focus on whether improvements in energy efficiency result in long term behavioural changes, and on the impact these changes can have on cold associated ill health.

What is already known on this topic

  • Excess winter mortality and morbidity are often greater in countries with relatively mild climates than in countries with colder climates

  • This paradox has been linked to poorer thermal efficiency of housing in temperate countries leading to colder indoor temperatures and creating an environment more susceptible to damp and mould

  • Several intervention studies have shown that improving thermal efficiency can lead to (most often self-reported) improvement in symptoms of cold associated disease

What this study adds

  • Retrofitting home insulation was associated with a reduction in the frequency of hospital admissions

  • The magnitude of the reduction was more pronounced for respiratory disease


We thank the following people and organisations for their support: David Briggs (Imperial College London) for advice on study design; Stephen Haslett (Massey University) for overview and advice on statistical analysis; James Stanley (University of Otago Wellington) and Marine Corbin and Mathu Shanthakumar (Massey University) for statistical advice and support with SAS software; and the Energy Efficiency and Conservation Authority and New Zealand Ministry of Health for technical advice and provision of the datasets.


  • Contributors: CF conducted the study under the guidance of LTB, PHC, and JD. All authors read, edited, and approved the final manuscript. The corresponding author attests that all listed authors meet authorship criteria and that no others meeting the criteria have been omitted.

  • Funding: No external funding was sought, or provided, for this study.

  • Competing interests: All authors have completed the ICMJE uniform disclosure form at and declare: no support from any organisation for the submitted; no financial relationships with any organisations that might have an interest in the submitted work in the previous three years; no other relationships or activities that could appear to have influenced the submitted work.

  • Ethical approval: The study was approved by the New Zealand Ministry of Health, Health and Disability Ethics Committee (15CEN8AM01).

  • Data sharing: No additional data available.

  • The lead author (CF) affirms that the manuscript is an honest, accurate, and transparent account of the study being reported; that no important aspects of the study have been omitted; any discrepancies from the study as planned (and, if relevant, registered) have been explained.

  • Dissemination to participants and related patient and public communities: It is not possible to disseminate results of the study to individuals as all records have been anonymised as part of the data linkage process. Results can be disseminated to the organisations that provided the original datasets, on request.

  • Provenance and peer review: Not commissioned; externally peer reviewed.

  • Publisher’s note: Published maps are provided without any warranty of any kind, either express or implied. BMJ remains neutral with regard to jurisdictional claims in published maps.

This is an Open Access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited and the use is non-commercial. See: