SEN support from the start of school and its impact on unplanned hospital utilisation in children with cleft lip and palate: a demonstration target trial emulation protocol using ECHILD ========================================================================================================================================================================================= * Vincent Grigori Nguyen * Anna Zylbergsztejn * Katie Harron * Tamsin Ford * Kristine Black-Hawkins * Kate Boddy * Johnny Downs * Martin Doyle * Matthew Lilliman * Jacob Matthews * Stuart Logan * Jugnoo Rahi * Ruth Gilbert * Lorraine Dearden * Bianca De Stavola ## Abstract Special Educational Needs (SEN) provision for school children provides extra support and reasonable adjustments for children and young people with additional educational, behavioural or health needs to ensure equal education opportunities; for example those born with a healthcare need such as cleft lip and palate may be provided SEN to aid with challenges in communications. However, there is limited knowledge of whether SEN provisions impact academic or health outcomes in such a population and conducting a randomised controlled trial to establish this evidence is not plausible. In lieu of randomised controlled trials, target trial emulation methods can be used in attempt to answer causal questions using observational data whilst reducing confounding and other biases likely to arise with such data. The Education and Child Health Insights from Linked Data (ECHILD) dataset could be used as part of trial emulation methods to understand the impact of SEN provisions on academic and healthcare outcomes. ECHILD is the first dataset to hold longitudinal school, health and social care data on all pupils in England, obtained by linking the National Pupil Database (NPD) with Hospital Episode Statistics (HES). In this protocol, we describe how the ECHILD dataset could be used to explore and conduct a target trial emulation to evaluate whether children who were born with cleft lip and palate would have different unplanned hospital utilisation if they received or did not receive SEN provisions by Year 1 (specifically by January in their second year of school) when they are aged 5 or 6. **Methods** Focussing on the population of children who are identified as having been born with cleft lip and palate, an intervention of varying levels of SEN provision (including no SEN provision) by January of the second year of school, and an outcome of unplanned hospital utilisation, we apply a trial emulation design to reduce confounding when using observational data to investigate the causal impact of SEN on unplanned hospital admissions. Our target population is children born 2001-2014 who had a recording of cleft lip and palate in HES and who started their second year of primary school (Year 1) in a state school between 2006 and 2019; children with a first recording of cleft lip and palate after Year 1 were excluded (these were pupils who likely immigrated to England after birth). We intend to use a time window of SEN provision assignment between the start of school (reception) and by the January school census in Year 1. Using target trial emulation, we aim to estimate the average treatment effect of SEN provision on the number of unplanned hospital admissions (including admissions to accident and emergency) between the January school census in Year 1 and Year 6 (the end of primary school, when children are 10-11 years old). **Ethics and dissemination** Permissions to use linked, de-identified data from Hospital Episode Statistics and the National Public Database were granted by DfE (DR200604.02B) and NHS Digital (DARS-NIC-381972). Ethical approval for the ECHILD project was granted by the National Research Ethics Service (17/LO/1494), NHS Health Research Authority Research Ethics Committee (20/EE/0180) and UCL Great Ormond Street Institute of Child Health’s Joint Research and Development Office (20PE06). Stakeholders (academics, clinicians, educators and child/young people advocacy groups) will consistently be consulted to refine populations, interventions and outcomes of studies that use the ECHILD dataset to conduct target trial emulation. Scientific, lay and policy briefings will be produced to inform public health policy through partners in the Department of Education and the Department of Health and Social Care. ## Introduction Special Educational Need (SEN) provisions aim to deliver reasonable adjustments in children and young people (CYP) in a school environment who need additional health, educational, or behavioural support, including children with complex health needs or learning disability. SEN provides support to those in need using a variety of facilities including support from a teaching assistant, assistance with communication, special learning programmes and support with physical needs. SEN provisions are divided into two levels, including: SEN Support (this is often known as Action, Action Plus or non-Statemented SEN) and Education and Health Care Plan (EHCP, this is previously known as Statemented SEN). SEN Support is organised at a school or college level and provides access to children and young people in need of SEN with support that may include teaching assistants who aid in communications, special learning programmes and supporting physical needs. An EHCP is organised by local authorities for children and young people who require further adjustments and often require additional resources (compared to SEN Support) to aid in education, health and social care needs. Currently, there is limited research on the impact of SEN on academic and healthcare outcomes in populations who are in need of SEN. To establish the causal effect of SEN on outcomes, randomised controlled trials would have to be conducted, however such study designs are not always feasible due to the human, time, financial and ethical costs associated. As SEN support is universally available in primary schools in England, conducting a randomised controlled trial would be unfeasible and possibly unethical at least for certain groups of children. In lieu of randomised controlled trials, observational studies provides a pragmatic, data-driven, observational study alternative when trials are not possible. One major challenge with using observational data when compared to data collected from randomised controlled trials is the risk of confounding, particularly, confounding by indication where assignment to treatment is not random and is often related to the severity of a medical condition. However, attentive study design can mitigate such biases in observational data by emulating the protocol of an equivalent randomised controlled trial (Hernán and Robins, 2016). An example of dataset that can used to evaluate the causal effect of SEN on healthcare outcomes using trial emulation is the ECHILD dataset, ([https://www.adruk.org/our-work/browse-all-projects/echild-linking-childrens-health-and-education-data-for-england-142/](https://www.adruk.org/our-work/browse-all-projects/echild-linking-childrens-health-and-education-data-for-england-142/)). The ECHILD dataset is the first dataset in England to link academic data with secondary care hospital data for all pupils, and can be used to investigate the associations between health, education and social care (Mc Grath-Lone et al., 2021). The ECHILD dataset links data from the National Pupil Database (NPD) and Hospital Episode Statistics (HES) and currently includes children and young people in England who were born between 1 September 1995 and 31 August 2020. Therefore, the ECHILD dataset provides the opportunity to conduct observational studies with long-term follow-up; currently follow-up is up to age 25 years (from birth in 1995 until hospitalisations in 2020). With data on clinical conditions, social care status, hospital visits, academic attainment and SEN provisions in school, the ECHILD dataset enables investigation of associations between a variety of exposures and outcomes when adjusting for different confounders in populations at risk. Examples of population at risk that can be phenotyped in hospital data include children with major congenital anomalies (e.g. Down Syndrome), cerebral palsy, learning disabilities, epilepsy, diabetes, and premature birth. Exposures in the education data include provision of support for SEN, Free School Meal status and measures of deprivation. Outcomes in ECHILD potentially include both educational and health related outcomes such as (specific types of) hospitalisations, national examination results (at multiple key stages) and exclusions from school. In this study protocol, we describe how we aim to use the ECHILD dataset and design an appropriate study that emulates a target trial that would address the question of interest. As an initial demonstration of the ECHILD dataset, we focus on the impact of receiving different types of SEN provision (including none) on unplanned hospital utilisation in children who were born with cleft lip and palate (Bell et al., 2016). We address how to evaluate whether an intervention recorded in educational data can impact outcomes recorded in health data. Incidences of cleft lip and palate impact 900 newborns in England yearly and impact communications (hearing and speech), dental health (Gallagher and Collett, 2019) and psychosocial health. Cleft lip and palate are associated with lower academic attainment (Fitzsimons et al., 2018) and have been linked to a three-fold increase in hospitalisations in Australia when compared to those without cleft lip and palate for all ages (Bell et al., 2016). Prior observational studies have suggested that extra support when starting school may be beneficial to children with cleft lip and palate for academic outcomes (Fitzsimons et al., 2018). Whilst Fitzsimons et al., 2018, has shown that children with cleft lip and palate are more likely to receive SEN provisions, particularly the types of SEN that concentrate on speech, language and communication support, to our knowledge, no study has compared the impact of extra help in terms of SEN using a control group of children with cleft lip and palate who have less or no support for any outcome. ## Methods ### Study Design and Setting The study will be a retrospective observational study based on data from the ECHILD dataset. The ECHILD database has been described previously (Mc Grath-Lone et al., 2021). Briefly, the ECHILD dataset consists of pseudo-anonymised linked data from the National Pupil Dataset and Hospital Episode Statistics for children who were born between 1 September 1995 and 31 August 2020 in England. ### Dataset and Linkage The data source we will use is the ECHILD database, a pseudo-anonymised dataset that links the NPD with HES. A full description of the ECHILD dataset has previously been described (Mc Grath-Lone et al., 2021). In brief, the ECHILD’s extract of the National Pupil Database contains data from academic terms (Summer, Autumn and Spring) between 2006 and 2020 and contains details on school, local authority, age, gender, ethnicity, first language, socioeconomic status, free school meal status, absence related data, social care/children in need related data and SEN status. The National Pupil Database is considered near universal as it covers all state schools in England from Key Stage 1 to Key Stage 4. The ECHILD’s extract of Hospital Episode Statistics contains details on admitted patient care, outpatient appointments, accident and emergency utilisation, and critical care between 1997 until 2021, It contains details on age, gender, ethnicity, clinical information recorded during hospital admissions, including details of diagnoses, and operations. Hospital Episode Statistics is considered near complete as it covers 99% of public hospital activity in England (Herbert et al., 2017). HES records since 1998 are also linked to ONS Mortality data covering information on causes and timing of deaths. The linkage coverage periods are described in (Mc Grath-Lone et al., 2021) and ECHILD has been shown to have a high linkage rate between NPD and HES (Libuy et al., 2021a) of 95%. Such high linkage rates are attributable through a two-stage linkage process (Libuy et al., 2021b). ### Population Our population consists of children born with a cleft lip and palate and followed from Year 1 of school (the second year of formal school attendance) between 2006 and 2019 (i.e. born between 2001 and 2014). To identify pupils who started Year 1 between 2006 and 2019, the earliest recording of “1” from the NCActualYear (National Curriculum Actual Year) variable in the NPD dataset will be used. To identify pupils with cleft lip and palate, ICD10 codes will be applied to HES diagnoses and OPCS codes will be applied to operations (Appendix 1 - Table 1 - Table 3); for each pupil, the earliest recorded date in HES would be considered the “diagnosis” date. We will remove pupils whose first recording of cleft lip and palate in HES is after their first year in school to evaluate the incident use of SEN; however as cleft lip and palate are often identified at or before birth, we suspect this would encompass a small number of pupils, for example, those born outside of England. To categorise cleft severity, we will use the following categories: cleft lip only, cleft palate only, unilateral cleft lip and palate, and bilateral cleft lip and palate(algorithm/method in Appendix 1 Table 2). ### Intervention and Exposure Variables The intervention assessed is the level of SEN (including none) recorded by Year 1 of school (age 5 to 6); whilst SEN support status can change throughout a child/young person’s educational journey, our implementation of trial emulation focusses on an intention-to-treat analysis (ITT), that is it analyses the assignment of treatment and not whether treatment was consistently adhered to (or provided). Due to the potential needs based on the severity of cleft lip and palate, we aim to classify the exposure variable as “levels of SEN” as opposed to a binary outcome (i.e. SEN vs no SEN).This is because in the event that the majority of Bilateral Cleft Lip + Palate (or more severe types of cleft lip and palate) contains small numbers in terms of “no SEN”, we could compare “SEN Support” against “EHCP”. For this reason, our exposure variable would be classified as: No SEN, SEN Support and EHCP by Year 1 and use the appropriate level as reference (i.e. SEN Support would be used as the “control” group if all children had received some form of SEN). To establish SEN status at Year 1, we will use the January (Autumn) census in Year 1 of school due to funding being calculated using these censuses. See Table 1 for a list of variables describing SEN in the NPD. In Figure 1 we demonstrate how we will identify the population and each subpopulation (based upon cleft severity), and how they will be classified according to the exposure variables (levels of SEN). View this table: [Table 1:](http://medrxiv.org/content/early/2022/04/05/2022.04.01.22273280/T1) Table 1: List of variables recording SEN in the National Pupil Database (NPD) ![Figure 1:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/04/05/2022.04.01.22273280/F1.medium.gif) [Figure 1:](http://medrxiv.org/content/early/2022/04/05/2022.04.01.22273280/F1) Figure 1: Flowchart of how the proposed population would be derived. Our population would consist of those who are identified in HES with cleft lip and palate before Year 1 and, who are linkable to NPD. To allow for heterogeneity of effects we aim to sub-classify our population based upon cleft lip and palate severity. SEN numbers will be based upon SEN recorded by Year 1. No SEN represents either never receiving SEN or receiving SEN after Year 1. ### Follow-up Our population will be followed-up from the time exposure is assessed (January Census of Year 1) till the end of primary school (Year 6) or end of study whichever occurred first. Children and young people may be considered lost to follow-up if they were off rolled (no longer in education), transferred outside of an English school that received government funding or died (based upon linked Office of National Statistics mortality linkage). Whilst the ECHILD data can be used to follow up children beyond the end of primary school (and for some individuals, beyond secondary school), we will limit our follow-up period in this protocol to the end of Year 6 for two reasons: 1) changing to secondary (or middle) school will re-evaluate the need for SEN, and many pupils may no longer be provided SEN, and 2) the time between the assignment of SEN considered here and the outcome may be too long if beyond Year 6, with the outcome affected by many intermediate factors (pupils starting Year 1 in 2006, would have 15 years of follow-up in HES). ### Outcome Variables The primary outcome variable is unplanned hospital utilisation. To identify unplanned hospital admissions in Admitted Patient Care, we will use the admission method variable in HES (admimeth) (Appendix 2 for the case definition from HES Admitted Patient Care). For hospital utilisation that did not require an overnight admission, we will use the HES Accident and Emergency dataset to account for non-admitted unplanned hospital utilisation (Harron et al., 2018). We aim to combine Admitted Patient Care and “Accident and Emergency” datasets to create a timeline of unplanned hospitalisation between the Autumn census in Year 1 and the end of Year 6. To measure unplanned hospitalisations, we aim to use two methods 1) the number of unplanned hospital admissions and 2) the number of unplanned days in hospital from study entry until study exit. When calculating the number of admissions using HES, we plan to combine admissions where the discharge date of one admission (disdate) was within one day of the next admission’s admission date (admidate). ### Covariates To account for non-random SEN assignment to this population of children, we will use information on several covariates that are known or suspected to influence SEN assignment and the outcome of interest. These include socio-demographic, health-related and education-related variables referring to the child. Specifically, these included: gender, Income Deprivation Affecting Children Index (IDACI) quintiles for lower super output area of residence, free-school meal status, comorbidities, birth month, clinical treatments for cleft lip and palate prior to SEN level assignment, birth characteristics such as gestational age and birth weight and school and local authority level data. An additional school-level variable we aim to include is the type of school (mainstream or special). We will also include the interaction between type of school and level of SEN to account for the likely heterogeneity in the effect of SEN across schools. The distribution of these potential confounders by population subgroup will be examined (see table 2 for an outline) and directed acyclic graphs representing the assumed relationships among these variables, SEN exposure, and the outcome of interest will be drawn to identify the variables that will be controlled for within the analysis using DAGitty ver 3.0. View this table: [Table 2:](http://medrxiv.org/content/early/2022/04/05/2022.04.01.22273280/T2) Table 2: Demographics Table of the Distribution of the sociodemographic, educational and health characteristics amongst children with different severity of Cleft. The table will include N and row percentages (when applicable) ### Bias To mitigate the likely confounding bias affecting observational data, we will adopt the Target Trial Emulation (TTE) framework. TTE enables observational data to be mapped to a hypothetical target experimental trial counterpart by creating the specification of an ideal (pragmatic) trial and using this as a basis to shape the observational study design. TTE consists of (1) defining the specifications of a hypothetical target experimental trial of the causal question of interest (including the corresponding effect), (2) emulating the specifications of the ideal target trial using observational data and (3) estimating the effects of interest using the emulated trial data. The first component of TTE involves defining an inclusion/exclusion criterion on entry, a treatment strategy (including time), follow-up frequency and modality, outcome measures, effects of interest (estimands) and the estimation method. Using the second component of TTE, observational data would be wrangled to emulate the distribution of the data if it were to have been gathered in the target trial. Finally, the third component of TTE requires the adjustment for confounding and is dependent on treatment assignment to aid in deciding the analytical methodology. In Table 3 we describe the target (ideal) trial one would design to investigate the causal effect of SEN provision (by January of the second year of school) in cleft lip and palate children and the equivalent emulated trial to be generated from ECHILD. View this table: [Table 3:](http://medrxiv.org/content/early/2022/04/05/2022.04.01.22273280/T3) Table 3: Trial Emulation Framework to estimate the causal effect of SEN by Year 1 on unplanned hospitalisations in children with cleft lip and palate ### Analysis #### Explorative analyses To explore confounding by indication in the causal relationship between receiving (different levels of) SEN and the outcome(s) of interest in the target population, possibly stratified by subtype, we will use the ECHILD dataset and conduct exploratory analysis to understand 1) whether the proportion and type of SEN provision differs between severity of cleft lip and palate and 2) whether such proportions change between 2006 and 2019, separately by category of cleft lip and palate. Such analysis would allow us to understand whether the severity of the cleft lip and palate is related to receiving SEN, whether certain severities of cleft lip and palate would have enough exposed/unexposed participants and whether public policy changes through time have impacted these proportions. Figures 2 and 3 are examples using simulated data of possible representations of how cross-sectional and cumulative frequencies of SEN provision and type vary over time, separated by cleft lip and palate type. ![Figure 2:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/04/05/2022.04.01.22273280/F2.medium.gif) [Figure 2:](http://medrxiv.org/content/early/2022/04/05/2022.04.01.22273280/F2) Figure 2: Example figure of how the single ECHILD cohort is segmented into multiple sequential cohorts including diagnosis period, exposure window and follow up period. ![Figure 2:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/04/05/2022.04.01.22273280/F3.medium.gif) [Figure 2:](http://medrxiv.org/content/early/2022/04/05/2022.04.01.22273280/F3) Figure 2: Simulated data showing different proportions of types of SEN per reception cohort, per cleft severity level ![Figure 3:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/04/05/2022.04.01.22273280/F4.medium.gif) [Figure 3:](http://medrxiv.org/content/early/2022/04/05/2022.04.01.22273280/F4) Figure 3: Simulated data showing the cumulative rates of SEN for each cohort and cleft severity over time #### Causal analyses Once explanatory analyses have been conducted, alternative estimation methods will be employed to address confounding bias due to non-random assignment to SEN. We will compare results obtained assuming no-unmeasured confounding (that is, we have data on all the relevant confounders) and assuming instead that we have an instrumental variable (if there is for example variation in provision by local authorities). The first group of methods will include: outcome-based models such as g-computation and propensity score methods (Inverse probability weighting and matching) using the covariates previously mentioned. The second group of methods will only be possible if suitable instruments for SEN provision are identified, for example if there are policy changes in provision that are implemented at different times across local authorities and will use instrumental variable-based estimators. As we aim to model the number of admissions and the days spent in hospital, we will use Poisson (or negative binomial) outcome models with the logarithm of follow-up time (expressed as number of terms) as offset. The likely clustering of pupils within local authority and school will be addressed either by fitting mixed effects models or by using robust inference (or both). #### Missing Data Based upon the proportion of missingness in the data and the mechanisms of missingness, we will first use complementary non-missing data points in HES and NPD prior to data imputation; for example, using Sex variable from HES to complement missing values in the NPD variable Gender. We aim to use multiple imputation using chained equations to generate multiple possible (for example, ten alternative) values for each missing data point (Azur et al., 2011). In the event we use both multiple imputation and matching, we will conduct imputation first and then match in each alternative dataset (Leyrat et al., 2017) #### Sensitivity Analyses To account for uncertainty in the recording of observational data such as measurement errors, we aim to conduct sensitivity analyses. First, we will conduct a sensitivity analysis to mitigate against a delayed recording in SEN provisions, by expanding the exposure window to the first term in Year Two as part of the Autumn census. Second, we will conduct sensitivity analyses to account for the complexity of congenital abnormalities in addition to cleft lip and palate, stratifying by (instead of adjusting for) congenital abnormalities in children and young people with cleft lip and palate. See Figure 4 for the alternative population derivation and Appendix 4 for the socio-economic, educational and healthcare distribution table for each cleft severity and birth anomaly interaction. ![Figure 4:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/04/05/2022.04.01.22273280/F5.medium.gif) [Figure 4:](http://medrxiv.org/content/early/2022/04/05/2022.04.01.22273280/F5) Figure 4: Alternative study design that stratifies for additional complexities in contrast to adjusting for comorbidities. Abnor* = abnormalities ## Stakeholder involvement Prior to developing this protocol, two independent meetings were conducted with stakeholders (parents, pupils, teachers, etc) to understand which medical conditions are of interest and at which entry timepoints are important for child development. The first meeting was with the Department for Education’s national young SEN advisory group (FLARE) on the 18th of September 2021 and the second with the Young Persons Advisory Group for research at Great Ormond Street Hospital on the 27th November 2021. This engagement identified that school entry is an important key milestone when SEN provisions are required. Therefore, in the proposed study, we have used school start as our entry point and will generate further target trials based upon further patient engagement. ## Data Availability Aggregate results from the ECHILD dataset will be pre-printed, revised and published. Individual record-level data with personal identifiers removed is currently hosted on the Office for National Statistics Secure Research Services data-sharing service. We are grateful to the Office for National Statistics (ONS) for providing the trusted research environment for the ECHILD Database. ONS agrees that the figures and descriptions of results in the attached document may be published. This does not imply ONS acceptance of the validity of the methods used to obtain these figures, or of any analysis of the results. The ECHILD Database uses data from the Department for Education (DfE). The DfE does not accept responsibility for any inferences or conclusions derived by the authors. This work uses data provided by patients and collected by the National Health Service as part of their care and support. Source data can also be accessed by researchers by applying to NHS Digital. ## Ethics and dissemination Permissions to use linked, de-identified data from Hospital Episode Statistics and the National Public Database were granted by DfE (DR200604.02B) and NHS Digital (DARS-NIC-381972). Ethical approval for the ECHILD project was granted by the National Research Ethics Service (17/LO/1494), NHS Health Research Authority Research Ethics Committee (20/EE/0180) and UCL Great Ormond Street Institute of Child Health’s Joint Research and Development Office (20PE06). ## Data Sharing and Access Aggregate results from the ECHILD dataset will be pre-printed, revised and published. Individual record-level data with personal identifiers removed is currently hosted on the Office for National Statistics Secure Research Service’s data-sharing service. We are grateful to the Office for National Statistics (ONS) for providing the trusted research environment for the ECHILD Database. ONS agrees that the figures and descriptions of results in the attached document may be published. This does not imply ONS’ acceptance of the validity of the methods used to obtain these figures, or of any analysis of the results. The ECHILD Database uses data from the Department for Education (DfE). The DfE does not accept responsibility for any inferences or conclusions derived by the authors. This work uses data provided by patients and collected by the National Health Service as part of their care and support. Source data can also be accessed by researchers by applying to NHS Digital. ## Funding This work is/was supported by ADR UK (Administrative Data Research UK), an Economic and Social Research Council (part of UK Research and Innovation) programme (ES/V000977/1). ## Appendix View this table: [Appendix 1](http://medrxiv.org/content/early/2022/04/05/2022.04.01.22273280/T4) Appendix 1 Phenotypes (ICD10 - OPCS) codes for cleft lip and palate View this table: [Appendix 2:](http://medrxiv.org/content/early/2022/04/05/2022.04.01.22273280/T5) Appendix 2: Determining unplanned admissions in hospital episode statistics admitted patient care View this table: [Appendix 3:](http://medrxiv.org/content/early/2022/04/05/2022.04.01.22273280/T6) Appendix 3: Variables related to receiving SEN and to unplanned hospitalisations View this table: [Appendix 4](http://medrxiv.org/content/early/2022/04/05/2022.04.01.22273280/T7) Appendix 4 For each cleft severity level (Cleft Lip only, Cleft Palate only, Unilateral Cleft Lip and Palate, Bilateral Cleft Lip and Palate) and comorbidity (no comorbidities, at last one comorbidity) interaction, we will generate a separate distribution table. The table will include N and row percentages * Received April 1, 2022. * Revision received April 1, 2022. * Accepted April 5, 2022. * © 2022, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution 4.0 International), CC BY 4.0, as described at [http://creativecommons.org/licenses/by/4.0/](http://creativecommons.org/licenses/by/4.0/) ## References 1. Alterman, N., Johnson, S., Carson, C., Petrou, S., Rivero-Arias, O., Kurinczuk, J.J., Macfarlane, A., Boyle, E., Quigley, M.A., 2021. Gestational age at birth and child special educational needs: a UK representative birth cohort study. Arch. Dis. Child. 106, 842–848. [https://doi.org/10.1136/archdischild-2020-320213](https://doi.org/10.1136/archdischild-2020-320213) [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTI6ImFyY2hkaXNjaGlsZCI7czo1OiJyZXNpZCI7czo5OiIxMDYvOS84NDIiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMi8wNC8wNS8yMDIyLjA0LjAxLjIyMjczMjgwLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 2. Azur, M.J., Stuart, E.A., Frangakis, C., Leaf, P.J., 2011. Multiple imputation by chained equations: what is it and how does it work? Int. J. Methods Psychiatr. Res. 20, 40–49. [https://doi.org/10.1002/mpr.329](https://doi.org/10.1002/mpr.329) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/mpr.329&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21499542&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F04%2F05%2F2022.04.01.22273280.atom) 3. Bell, J., Nassar, N., Turner, R., Bower, C., Gillett, D., McBain, W., Raynes-Greenow, C., 2016. Hospitalisations up to adulthood for children born with orofacial clefts: Hospitalisations for people with clefts. J. Paediatr. Child Health 52, 441–448. [https://doi.org/10.1111/jpc.13024](https://doi.org/10.1111/jpc.13024) 4. Campbell, T., 2021. Special Educational Needs and Disabilities within the English primary school system: What can disproportionalities by season of birth contribute to understanding processes behind attributions and (lack of) provisions? London Scool of Economics and Political Sciences. 5. Fitzsimons, K.J., Copley, L.P., Setakis, E., Charman, S.C., Deacon, S.A., Dearden, L., van der Meulen, J.H., 2018. Early academic achievement in children with isolated clefts: a population-based study in England. Arch. Dis. Child. 103, 356–362. [https://doi.org/10.1136/archdischild-2017-313777](https://doi.org/10.1136/archdischild-2017-313777) [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTI6ImFyY2hkaXNjaGlsZCI7czo1OiJyZXNpZCI7czo5OiIxMDMvNC8zNTYiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMi8wNC8wNS8yMDIyLjA0LjAxLjIyMjczMjgwLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 6. Gallagher, E.R., Collett, B.R., 2019. Neurodevelopmental and Academic Outcomes in Children With Orofacial Clefts: A Systematic Review. Pediatrics 144, e20184027. [https://doi.org/10.1542/peds.2018-4027](https://doi.org/10.1542/peds.2018-4027) [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTA6InBlZGlhdHJpY3MiO3M6NToicmVzaWQiO3M6MTU6IjE0NC8xL2UyMDE4NDAyNyI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIyLzA0LzA1LzIwMjIuMDQuMDEuMjIyNzMyODAuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 7. Harron, K., Gilbert, R., Cromwell, D., Oddie, S., Guttmann, A., van der Meulen, J., 2018. International comparison of emergency hospital use for infants: data linkage cohort study in Canada and England. BMJ Qual. Saf. 27, 31–39. [https://doi.org/10.1136/bmjqs-2016-006253](https://doi.org/10.1136/bmjqs-2016-006253) [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoicWhjIjtzOjU6InJlc2lkIjtzOjc6IjI3LzEvMzEiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMi8wNC8wNS8yMDIyLjA0LjAxLjIyMjczMjgwLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 8. Harron, K., Gilbert, R., Cromwell, D., van der Meulen, J., 2016. Linking Data for Mothers and Babies in De-Identified Electronic Health Data. PLOS ONE 11, e0164667. [https://doi.org/10.1371/journal.pone.0164667](https://doi.org/10.1371/journal.pone.0164667) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F04%2F05%2F2022.04.01.22273280.atom) 9. Harron, K., Gilbert, R., Fagg, J., Guttmann, A., van der Meulen, J., 2021. Associations between pre-pregnancy psychosocial risk factors and infant outcomes: a population-based cohort study in England. Lancet Public Health 6, e97–e105. [https://doi.org/10.1016/S2468-2667(20)30210-3](https://doi.org/10.1016/S2468-2667(20)30210-3) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F04%2F05%2F2022.04.01.22273280.atom) 10. Herbert, A., Wijlaars, L., Zylbersztejn, A., Cromwell, D., Hardelid, P., 2017. Data Resource Profile: Hospital Episode Statistics Admitted Patient Care (HES APC). Int. J. Epidemiol. 46, 1093–1093i. [https://doi.org/10.1093/ije/dyx015](https://doi.org/10.1093/ije/dyx015) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/ije/dyx015&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28338941&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F04%2F05%2F2022.04.01.22273280.atom) 11. Hernán, M.A., Robins, J.M., 2016. Using Big Data to Emulate a Target Trial When a Randomized Trial Is Not Available: Table 1. Am. J. Epidemiol. 183, 758–764. [https://doi.org/10.1093/aje/kwv254](https://doi.org/10.1093/aje/kwv254) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/aje/kwv254&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26994063&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F04%2F05%2F2022.04.01.22273280.atom) 12. Hosozawa, M., Cable, N., Kelly, Y., Sacker, A., 2021. Gestational age on trajectories of social competence difficulties into adolescence. Arch. Dis. Child. 106, 1075–1080. [https://doi.org/10.1136/archdischild-2020-321317](https://doi.org/10.1136/archdischild-2020-321317) [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTI6ImFyY2hkaXNjaGlsZCI7czo1OiJyZXNpZCI7czoxMToiMTA2LzExLzEwNzUiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMi8wNC8wNS8yMDIyLjA0LjAxLjIyMjczMjgwLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 13. Hutchinson, J., 2021. Identifying Pupils with Special Educational Needs and Disabilities. 14. Jerrim, J., 2021. Measuring Disadvantage. 15. Jutte, D.P., Roos, N.P., Brownell, M.D., Briggs, G., MacWilliam, L., Roos, L.L., 2010. The Ripples of Adolescent Motherhood: Social, Educational, and Medical Outcomes for Children of Teen and Prior Teen Mothers. Acad. Pediatr. 10, 293–301. [https://doi.org/10.1016/j.acap.2010.06.008](https://doi.org/10.1016/j.acap.2010.06.008) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.acap.2010.06.008&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20674531&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F04%2F05%2F2022.04.01.22273280.atom) 16. Leyrat, C., Seaman, S.R., White, I.R., Douglas, I., Smeeth, L., Kim, J., Resche-Rigon, M., Carpenter, J.R., Williamson, E.J., 2017. Propensity score analysis with partially observed covariates: How should multiple imputation be used? Stat. Methods Med. Res. 0962280217713032. [https://doi.org/10.1177/0962280217713032](https://doi.org/10.1177/0962280217713032) 17. Libuy, N., Harron, K., Gilbert, R., Caulton, R., Cameron, E., Blackburn, R., 2021a. Linking education and hospital data in England: linkage process and quality. Int. J. Popul. Data Sci. 6. [https://doi.org/10.23889/ijpds.v6i1.1671](https://doi.org/10.23889/ijpds.v6i1.1671) 18. Libuy, N., Harron, K., Gilbert, R., Caulton, R., Cameron, E., Blackburn, R., 2021b. Linking education and hospital data in England: linkage process and quality. Int. J. Popul. Data Sci. 6. [https://doi.org/10.23889/ijpds.v6i1.1671](https://doi.org/10.23889/ijpds.v6i1.1671) 19. Mackay, D., Smith, G., Dobbie, R., Cooper, S.-A., Pell, J., 2013. Obstetric factors and different causes of special educational need: retrospective cohort study of 407 503 schoolchildren. BJOG Int. J. Obstet. Gynaecol. 120, 297–308. [https://doi.org/10.1111/1471-0528.12071](https://doi.org/10.1111/1471-0528.12071) 20. Mc Grath-Lone, L., Libuy, N., Harron, K., Jay, M.A., Wijlaars, L., Etoori, D., Lilliman, M., Gilbert, R., Blackburn, R., 2021. Data Resource Profile: The Education and Child Health Insights from Linked Data (ECHILD) Database. Int. J. Epidemiol. dyab149. [https://doi.org/10.1093/ije/dyab149](https://doi.org/10.1093/ije/dyab149) 21. Persson, M., Opdahl, S., Risnes, K., Gross, R., Kajantie, E., Reichenberg, A., Gissler, M., Sandin, S., 2020. Gestational age and the risk of autism spectrum disorder in Sweden, Finland, and Norway: A cohort study. PLOS Med. 17, e1003207. [https://doi.org/10.1371/journal.pmed.1003207](https://doi.org/10.1371/journal.pmed.1003207) 22. Strand, S., Lindorff, A., 2018. Ethnic disproportionality in the identification of Special Educational Needs (SEN) in England: Extent, causes and consequences. 23. Zylbersztejn, A., Gilbert, R., Hardelid, P., 2020. Developing a national birth cohort for child health research using a hospital admissions database in England: The impact of changes to data collection practices. PLOS ONE 15, e0243843. [https://doi.org/10.1371/journal.pone.0243843](https://doi.org/10.1371/journal.pone.0243843)