Individualised variable-interval risk-based screening for sight-threatening diabetic retinopathy: the Liverpool Risk Calculation Engine

Diabetologia, Aug 2017

Aims/hypothesis Individualised variable-interval risk-based screening offers better targeting and improved cost-effectiveness in screening for diabetic retinopathy. We developed a generalisable risk calculation engine (RCE) to assign personalised intervals linked to local population characteristics, and explored differences in assignment compared with current practice. Methods Data from 5 years of photographic screening and primary care for people with diabetes, screen negative at the first of > 1 episode, were combined in a purpose-built near-real-time warehouse. Covariates were selected from a dataset created using mixed qualitative/quantitative methods. Markov modelling predicted progression to screen-positive (referable diabetic retinopathy) against the local cohort history. Retinopathy grade informed baseline risk and multiple imputation dealt with missing data. Acceptable intervals (6, 12, 24 months) and risk threshold (2.5%) were established with patients and professional end users. Results Data were from 11,806 people with diabetes (46,525 episodes, 388 screen-positive). Covariates with sufficient predictive value were: duration of known disease, HbA1c, age, systolic BP and total cholesterol. Corrected AUC (95% CIs) were: 6 months 0.88 (0.83, 0.93), 12 months 0.90 (0.87, 0.93) and 24 months 0.91 (0.87, 0.94). Sensitivities/specificities for a 2.5% risk were: 6 months 0.61, 0.93, 12 months 0.67, 0.90 and 24 months 0.82, 0.81. Implementing individualised RCE-based intervals would reduce the proportion of people becoming screen-positive before the allocated screening date by > 50% and the number of episodes by 30%. Conclusions/interpretation The Liverpool RCE shows sufficient performance for a local introduction into practice before wider implementation, subject to external validation. This approach offers potential enhancements of screening in improved local applicability, targeting and cost-effectiveness.

A PDF file should load here. If you do not see its contents the file may be temporarily unavailable at the journal website or you do not have a PDF plug-in installed and enabled in your browser.

Alternatively, you can download the file locally and open with any standalone PDF reader:

Individualised variable-interval risk-based screening for sight-threatening diabetic retinopathy: the Liverpool Risk Calculation Engine

Individualised variable-interval risk-based screening for sight-threatening diabetic retinopathy: the Liverpool Risk Calculation Engine Antonio Eleuteri 0 1 2 4 5 6 Anthony C. Fisher 0 1 2 4 5 6 Deborah M. Broadbent 0 1 2 4 5 6 Marta García-Fiñana 0 1 2 4 5 6 Christopher P. Cheyne 0 1 2 4 5 6 Amu Wang 0 1 2 4 5 6 Irene M. Stratton 0 1 2 4 5 6 Mark Gabbay 0 1 2 4 5 6 Daniel Seddon 0 1 2 4 5 6 Simon P. Harding 0 1 2 4 5 6 0 Department of Biostatistics, Institute of Translational Medicine, University of Liverpool , Liverpool , UK 1 St Paul's Eye Unit, Royal Liverpool University Hospital , Liverpool , UK 2 Department of Eye and Vision Science, Institute of Ageing and Chronic Disease, University of Liverpool , William Henry Duncan Building, 6, West Derby Street, Liverpool L7 8TX , UK 3 Simon P. Harding 4 Public Health England, Cheshire and Merseyside Screening and Immunisation Team , Liverpool , UK 5 Department of Health Services Research, University of Liverpool , Liverpool , UK 6 Gloucestershire Retinal Research Group, Cheltenham General Hospital , Cheltenham , UK Aims/hypothesis Individualised variable-interval risk-based screening offers better targeting and improved costeffectiveness in screening for diabetic retinopathy. We developed a generalisable risk calculation engine (RCE) to assign personalised intervals linked to local population characteristics, and explored differences in assignment compared with current practice. Methods Data from 5 years of photographic screening and primary care for people with diabetes, screen negative at the first of > 1 episode, were combined in a purpose-built nearreal-time warehouse. Covariates were selected from a dataset created using mixed qualitative/quantitative methods. Markov modelling predicted progression to screen-positive (referable diabetic retinopathy) against the local cohort history. Retinopathy grade informed baseline risk and multiple imputation dealt with missing data. Acceptable intervals (6, 12, A list of ISDR Study Group members is included in the electronic supplementary material (ESM). - Department of Medical Physics and Clinical Engineering, Royal Liverpool University Hospital, Liverpool, UK 5 6 7 proach offers potential enhancements of screening in improved local applicability, targeting and cost-effectiveness. Abbreviations AICc Corrected Akaike’s information criterion LDESP Liverpool Diabetes Eye Screening Programme RCE Risk calculation engine STDR Sight-threatening diabetic retinopathy Introduction Systematic screening for sight-threatening diabetic retinopathy (STDR) has been introduced in several European countries and regionally throughout the world, and has been a major driver of improved detection and early treatment. As a doubling of the global prevalence of diabetes mellitus is expected by 2030 [1], with over 10% having STDR [2], there is an urgent need to improve the cost-effectiveness of screening. While current recommendations are for annual screening intervals in most locations [3], there has been a recent move to recommend biennial screening for people with no retinopathy [4–7], including in one systematic review [8], and this was recently endorsed by the UK National Screening Committee [9]. Screening at 3-yearly intervals has been introduced in Sweden, based on data from one programme [10], and is supported as being cost-effective in a recent UK modelling study [11]. However, concerns about the safety and acceptability of extended intervals have held back adoption [12, 13]. Risk engines have been developed in recent years, including in diabetes mellitus for risk of CHD [14], and one has been proposed for diabetic retinopathy [15]. For widespread uptake, reliable flows of data need to be established and designs need to be applicable across a range of populations and health settings. As part of a programme of research to improve the targeting and cost-effectiveness of screening, we developed a generalisable personalised screening method to allow variable intervals for people with diabetes at high and low risk of developing STDR. We developed and internally validated a risk calculation engine (RCE) to estimate risk of progression to screen-positive or referable diabetic retinopathy and assign individualised screening intervals. We calculated improvement in allocation of screening interval to estimate the effect on number of screen episodes. Methods Dataset Data from established digital photographic screening (OptoMize, EMIS Group, Leeds, UK) and primary care systems (EMISweb, EMIS Group) were combined in a purpose-built data warehouse. The local ethics committee approved an opt-out approach to consent (13/NW/0196) and the research was conducted in accordance with the Declaration of Helsinki 2008. Data were collected for all people recorded in primary care as having diabetes mellitus attending the Liverpool Diabetes Eye Screening Programme (LDESP) from the systems used for routine service, anonymised and compiled before transmission to the warehouse. A set of candidate covariates was selected for the model using patient expert panels and a literature review of known risk factors (see electronic supplementary material [ESM] Methods and ESM Table 1). An RCE development dataset was extracted from the data warehouse containing covariates with ≥ 80% completeness in people with diabetes who were screen negative (non-referable retinopathy) at the first of at least two episodes that occurred in a 5 year sample period. Disease duration was defined as duration of known diabetes mellitus (first recorded date of diabetes or measure of HbA1c in primary care [11]) and assigned at the first screening episode. Values of clinical risk factors prior and nearest to the screen episode date were used. Screen-positive (the primary outcome) was defined as the presence of any of: multiple blot haemorrhages, venous beading, intra-retinal microvascular abnormalities, new vessels, pre-retinal/vitreous haemorrhage, tractional retinal detachment, exudates within 1 disc diameter (1500 μm) of the foveal centre, group of exudates within the macula or blot haemorrhages within 1 disc diameter of the foveal centre with vision worse than 6/12. Model description We selected a continuous-time Markov process to allow for a set of individuals to move independently, or transition, between states over time [16]. The patient state at each time point was defined by level of retinopathy, including separation by one or both eye involvement after Stratton et al [17] (Fig. 1:) (1) no diabetic retinopathy detected; (2) non-referable diabetic retinopathy in one eye only; (3) non-referable diabetic retinopathy in both eyes; and (4) referable diabetic retinopathy (screen-positive for at least one eye). Only one baseline screening event was used. The risks, or intensities, for each transition were entered into the model within a probability matrix containing Weibull transition intensities estimated from the data [18–20]. A detailed description is provided in the ESM Methods. The data in the RCE development dataset is an example of panel data [16, 18] where information on an individual’s disease is sampled at time points not typically coincident with the change in disease state. This interval-censoring problem is illustrated in ESM Fig. 1 and required special methods. Missing clinical data were handled using multiple imputation [21] (ESM Methods) repeated ten times to properly account for variability due to unknown values. Model fitting and covariate selection Covariates meeting the above criteria were ranked using Wald statistics. A set of nested models were built to estimate corrected Akaike’s information criterion (AICc). This method combines estimation (i.e. maximum likelihood) and model selection under a unified framework [16, 22, 23]. AIC was corrected to adjust for the number of covariates (see ESM Methods); this is a method which aims to reduce the effect of overfitting by applying a penalty for model complexity. The model with the smallest AICc was chosen to give the best fit to the data. Patient expert group A patient involvement group was developed through local and national patient groups and local advertisements with a mix of backgrounds, sex and diabetes type. The group developed their knowledge of disease and patient pathways and the assessment of risk over several tailored sessions. At the end of these, they expressed that they had sufficient knowledge to give considered input into the study design. Acceptability and feasibility were considered for the application of the RCE output for a range of risk thresholds and alternative screen intervals. Data validation and model checking We checked the development dataset using random samples of event vectors which were independently checked manually and programmatically. The model was checked for influence of outliers, regression and distributional assumptions, and Pearson-type goodnessof-fit and corrected C-index were calculated. Bootstrapping (to estimate the optimism of validation measures) and fourfold cross-validation were used for internal validation (see ESM Methods). Further internal validation was conducted using a geographical split based on the deprivation index [24], to assess whether the performances of the model were unduly affected by extremes of prevalence rates of positive screening events. Areas under the receiver operating curves were calculated as overall indicators of sensitivities and specificities. Implementation The effect of a set of risk thresholds (5%, 2.5%, 1%) on screening-interval allocation was investigated Referable DR (at least one eye) (4) using the fourfold validation sets described above, and a final threshold selected in discussion with the patient expert group. The proportion of screen-positive events that occurred before the allocated interval was calculated for each risk threshold. Overall numbers of screening episodes were calculated over a 2 year period and compared with an annual programme. Predictions in a clinical environment (using the fitted model) for new observations with missing covariates were obtained by a simple imputation strategy: we replaced the missing values of each covariate with a 75th percentile value estimated from full data at first screening to give a ‘worst case’ prediction. A small sample of cases assigned by the RCE to 6, 12 and 24 months rescreen interval were independently checked against patient records for clinical credibility. Results The data repository contained 2.48 × 1010 data fields across 302 covariates. Data extracted into the RCE development dataset were from 11,806 people with diabetes actively attending the LDESP between 20 Feb 2009 and 4 Feb 2014 (46,525 episodes, 388 screen-positive events). Data flow is shown in ESM Fig. 2. Covariates that met the entry criteria are listed in Table 1. Those that gave the best fit (total rescaled AICc = 0) and were included in the final model were: disease duration, HbA1c, age at diagnosis, systolic BP and total cholesterol. Although the retinopathy stage is not technically a covariate, it is included in the table to show the improvement in predictive power when covariates are added to the model. The risk model is summarised in the following three equations. The first expression gives the hazard rates (or intensities, or ‘risks’) of going from one state to another, for each of the six transitions shown in Fig. 1, each one dependent on the full set of covariates and the full local population history: logλij ¼ βi0j þ βiAjgeDAgeD þ βiDj iseaseDDiseaseD þ ⋯ þ βiSjPSBP ð1Þ where i , j = {1, 2, 3, 4} and βiCj is the model parameter for covariate C (or baseline intensity when C = 0); AgeD is age at diagnosis, DiseaseD is disease duration. Wald statistic Rescaled AICc % Explained likelihood From the hazard rates function (1), a transition intensity matrix is derived for the four states (each row sums to 0): 0 −λ11 λ21 Q ¼ @BB 0 0 ESM Fig. 3 shows the two baseline transition intensities to the screen-positive state in the rightmost column of Eq. (2). Probabilities of transition occurring at a specific time t are obtained by using the following equation: ð Þ ¼ exp Qt0:9 P t where the exp() operator is the matrix exponential and α = 0.9 is the estimated Weibull shape parameter. The shape parameter 0.9 reflects the clinical observation that a person is more likely to move between disease states earlier rather than later. Data and model checking The pseudo-likelihood ratio p value for the summary residuals vs time was 0.04, suggesting linearity to hold between 0 and 2 years, with a possible lack of fit beyond 2 years (ESM Fig. 4). Although the p value was below 0.05, there is not enough evidence of a lack of fit because of the small number of events relative to the model complexity [20, 22]. Cox–Snell residuals are shown in ESM Fig. 5: the calibration curve was close to the theoretical optimal calibration, showing that the model tended to give slightly pessimistic predictions of failure. The Pearson-type statistic for the Liverpool RCE model was 0.57, denoting not enough evidence to reject the null hypothesis of good fit. CrossTransition AgeD DiseaseD HbA1c Chol SBP validation showed only very small effects, i.e. the training and test performance measures were essentially the same. Fitting the model to the most deprived 65% of individuals produced only very small changes in risk allocation of the non-deprived group. The patient expert panel identified acceptable screen intervals of 6, 12 and 24 months and risks of 1% and 2.5% as acceptable risks of missing screen-positive disease at any future screen episode. Exploration of the effect of different risk thresholds on allocation to the three different screen intervals using our four-way cross-validation is shown in Table 4. As the risk threshold decreased, the proportion of incorrect screen-interval allocations decreased for screen positives (overestimation) and increased (underestimation) for screen negatives. For an annual interval, the overestimation was 50.2% and underestimation 100%; when compared with this interval, the number of individuals correctly allocated was greater for the model with any of the three thresholds. The number of intervals which were either over- or underestimated was also lower for the model at all three thresholds. Table 5 shows the detailed comparison between the overall numbers of screening episodes in an annual programme and the Liverpool RCE model. For all three risk thresholds, there was a reduction in the overall number of screening episodes required (summarised in Table 4). The research team and patient expert panel considered that a 2.5% criterion showed a satisfactory distribution across the three screening intervals and a reasonable reduction in episodes, and this was selected for implementation. Using the 2.5% threshold, the corrected C-index for the model was 0.687 and corrected AUCs (with 95% CIs) were 0.88 (0.83, 0.93) at 6 months, 0.90 (0.87, 0.93) at 12 months a Proportion of individuals allocated to 6, 12 and 24 month screening using 5%, 2.5% and 1% risk thresholds b Distribution of screen-positive and -negative events depending on risk threshold is also shown, followed by the proportion of reduction of visits compared with annual screening c Excluded data from individuals with insufficient follow-up data Screening Risk threshold 5% 5.5% 4.4% 90.1% 65.6% 34.4% 90.4% 9.6% 39.5% 2.5% 10.7% 8.6% 80.7% 78.0% 22.0% 79.8% 20.2% 29.7% 1% 26.3% 9.3% 64.4% 91.1% 8.9% 62.2% 37.8% 5.9% 0% 100% 0% 49.8% 50.2% 0% 100% Reference value Difference from standard allocation (%) Breakdown of numbers of screening episodes required over a 2 year period from each validation set (and combined) if the cut-off was 5%, 2.5% or 1% and the midpoint of positive-screen event interval was truly representative of when that positive event actually occurred, along with a comparison with the currently used standard annual screening allocation a Number of screening episodes at 6, 12 and 24 months is calculated as 4×, 2× and 1× the number of people allocated to those screening episodes, respectively, to represent the number of screenings they would attend over a 2 year period b Standard number of screening episodes based on all individuals given annual screening (i.e. 2× number of individuals) VS, validation set and 0.91 (0.87, 0.94) at 24 months. The four-way random data split gave sensitivities and specificities for a risk threshold of 2.5% at 6, 12 and 24 months, respectively: 6 months 0.61, 0.93; 12 months 0.67, 0.90; and 24 months 0.82, 0.81. Clinical review of sampled cases (n = 18) indicated that allocations to individualised screen intervals appeared reasonable. Discussion We have developed and tested an RCE in which an individual’s risk can be predicted from contemporaneous routinely collected clinical data, referenced to the clinical histories of the local population, using covariates of local relevance. The risk can be reassessed at each screening episode as new clinical information is acquired. The Markov approach we have used allows a dynamic model of the retinopathy history to be built. In a sense, the model ‘compresses’ the information about time evolution. The Markov property can be summarised by the phrase ‘The future is predicted from the past through the present’, and is particularly appropriate to our setting. The strengths of our model include our approach to tackling the data in routine screening. Retinopathy data in screening is interval censored [16] in that the event seems as if it has happened when it is detected. This may lead to biased estimates, as it ‘seems’ like the disease developed later than it actually did. Unlike other ‘classic’ model types, including the Cox model, the Markov approach can internally handle this interval censoring. In addition, it predicts the probabilities of transition for all disease states. ‘Real life’ data from routine clinical practice inevitably introduces missingness and recording errors. We embedded a model for multiple imputation of missing covariates, which was required to allow our RCE to run effectively. Potential limitations of our RCE relate to model design and some of the covariates. We did not adjust for misclassification of retinopathy during grading. This could be addressed by adding a misclassification model, but at the cost of substantially more observations and computational complexity. Some covariates were not informative in the Liverpool setting. Ethnic diversity is low and the prevalence of abnormal eGFRs <60 ml min−1 1.73 m−2 was only 14.5%. Other covariates such as social deprivation score may be worth adding. ‘Type of diabetes’ may not be accurately recorded in primary care and the increased use of insulin in type 2 diabetes makes ‘insulin usage’ an unreliable criterion. We used date of first HbA1c test to improve data on ‘duration of diabetes’, helpful especially in people with long durations, but less reliable since the introduction of HbA1c as a primary screening test. The model consistently showed good levels of prediction for the 2.5% risk threshold. The numbers of screen-positive cases with overestimated screening dates and screen-negative cases with underestimated screening dates were reduced. The majority of people were correctly allocated (78% of screen positives, 80% of screen negatives), with a reasonable allocation of (approximately) 10%:10%:80% across the 6, 12 and 24 month intervals. The number of patients who had the screen event before the allocated screening date was reduced by more than half and the overall number of screening episodes was reduced by 30%. We included a strongly embedded local patient group, which allowed us to develop an appropriate preliminary covariate list and acceptable screen intervals and risk threshold. This group developed expertise over a series of meetings and provided substantial input into design and implementation. Strong patient and professional involvement is very valuable in study design and delivery. Our RCE development process is suitable for a wide range of geographical locations and populations with a minimum prerequisite of a centrally maintained disease register with adequate historical data. Revision/addition of covariates can be accommodated based on the strength they add to a locally developed model. For example, higher prevalence of poor diabetes control or renal disease may strengthen the effect of HbA1c or eGFR. Alternative intervals including extension beyond 24 months could be developed subject to acceptability. Local populations may select alternative risk thresholds depending on the perception of risk. We give the key steps to developing and building such a system in the text box. Key steps to developing a locally relevant risk calculation engine RCE development dataset Local stakeholder engagement: Approvals and systems for regular data transfers Patient and professional groups for covariate selection First iteration RCE development dataset Data exploration and verification Review processing protocols and limit criteria for local relevance Lock risk engine development dataset Multiple imputation Covariate selection Select preliminary Markov model using all covariates (model fitting) Assign informedness to covariates Review with local patient and professional groups Finalise covariates and fix model structure RCE model Build test model Agree choice of screen intervals Run diagnostics and validation Final model version (revise if required) RCE model Secure domain for implementation model Determine frequency of updating data (recommend 2 monthly) Re-tune model (recommend 3 yearly) The use of near-real-time data and a model developed from local data in our approach is novel. Aspelund et al developed a risk-estimating model in Iceland [15]. They used a proportional hazards Weibull model informed by local retinopathy data between 1994 and 1997 and risks for covariates estimated from data published in the 1990s. ROC analysis showed a fair performance, with 59% fewer visits than annual screening. Van der Heijden et al tested this model in an up-to-date prospective cohort of people with type 2 diabetes [25]. Of a total of 8303, 3319 met the eligibility criteria, with a mean of 53 months follow-up. Discriminatory ability was good (C-statistic 0.83), but 67 of 76 people (88.2%) who developed STDR developed it after the time predicted by the model. This overestimation of risk highlights the weakness of using historical data. Hippisley-Cox and Coupland recently developed equations to predict 10 year rates of amputation and blindness using similar methods to us [26]. They studied routinely collected general practice and hospital episode data from 454,575 people with diabetes. A web-based 10 year calculator using Cox’s proportional hazards models was developed. They reported comparable C-statistics (≥ 0.73) and conducted external validation using 357 practices that used a different database. The principal limitation of this large study was the lack of validation of the diagnosis of blindness. Risk engines have been developed in other diseases including coronary heart disease, stroke and lipid therapy. The UK Prospective Diabetes Study developed a risk engine for predicting coronary heart disease [14], now in its second version (UKPDS Outcomes Model 2). We included clinical risk factors in our model. It has recently been suggested that retinopathy data are sufficient to develop a risk stratification to extend screening intervals for people at low risk [27]. This may prove to be a reasonable and pragmatic approach. We had to overcome significant challenges in developing a near-real-time data flow; this may be too difficult in some populations. However, we determined that including clinical data would aid acceptance amongst the professional community, offer better prospects for generalisability and allow inclusion of more frequent screening for high-risk individuals. Our view is supported by our own data [28] and those of others [29], and also by our patient expert group. We do recognise that, as yet, estimates of resource requirements for the effective introduction of our type of RCE are not available. External validation of models is required before general implementation [30]. However, validation methods for an approach such as ours are not well developed. An RCE comprises two principal components: (1) the dataset containing a set of covariates and the outcome of interest; and (2) the mathematical model applied to the data in the dataset. The application to a population is specific to that population. In addition to the interval censoring described above, screening data are also not proportional. This makes problematic the use of widely accepted statistics for assessing effectiveness of diagnostic tools based on Kaplan–Meier methods. An approach to validation was developed, taking these constraints into account, comprising dataset validation, model checking, internal validation (including data splitting, bootstrapping, C-index) and estimation of sensitivities/ specificities at specified intervals, all recognised internal validation methods [30]. An implementation phase will include model updating (temporal validation and model tuning) and the opportunity for comparative cross population (external) validation to correct for potential overperformance [31]. We believe that the Liverpool RCE is feasible, reliable, safe and acceptable to patients. Implementation of our RCE into routine clinical practice would offer potentially significant transfer of resources into targeting high-risk and hard-to-reach groups and improved cost-effectiveness. Based on the internal validations we have performed, it shows sufficient performance for a local introduction. However, wider implementation will require an external validation process and testing of safety and acceptability in an RCT setting [31]. Investment in IT systems will be required for implementation in large-scale health systems, such as the NHS, and to support further validation. Acknowledgements The authors are grateful to the Individualised Screening for Diabetic Retinopathy (ISDR) Patient and Public Involvement Group for essential input into design and review; to the Liverpool Care Commissioning Group for data extraction and transfer; and to the Liverpool Local Medical Committee and local general practitioners for support with establishing patient lists and consent. Data availability The Liverpool RCE Development Dataset generated and analysed during this study is not publicly available because of restrictions on data sharing and commercialisation. A fully anonymised dataset is available from the corresponding author on reasonable request. Funding This manuscript presents independent research funded by the National Institute for Health Research (NIHR; RP-PG-1210-12016). The views expressed are those of the authors, not those of the UK National Health Service, NIHR or Department of Health. MGF is part funded by NIHR Collaboration for Leadership in Applied Health Research and Care North West Coast (NIHR CLAHRC NWC). Duality of interest The authors declare that there is no duality of interest associated with this manuscript. Contribution statement All authors met ICMJE requirements, making (1) substantial contributions to conception and design, acquisition of data, or analysis and interpretation of data; (2) drafting the article or revising it critically for important intellectual content; and (3) giving final approval of the version to be published. SPH wrote the drafts of the manuscript; AE wrote the technical sections and prepared the figures. SPH is responsible for the integrity of the work as a whole. AE is guarantor for the model development and analysis. Open Access This article is distributed under the terms of the Creative C o m m o n s A t t r i b u t i o n 4 . 0 I n t e r n a t i o n a l L i c e n s e ( h t t p : / /, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 1. 15. Wild S , Roglic G , Green A , Sicree R , King H ( 2004 ) Global prevalence of diabetes: estimates for the year 2000 and projections for 2030 . Diabetes Care 227 : 1047 -1053 Yau JWY , Rogers SL , Kawasaki R et al ( 2012 ) Global prevalence and major risk factors of diabetic retinopathy . Diabetes Care 35 : 556 - 564 Scanlon PH ( 2008 ) The English national screening programme for sight-threatening diabetic retinopathy . J Med Screen 15 : 1 -4 Chalk D , Pitt M , Vaidya B , Stein K ( 2012 ) Can the retinal screening interval be safely increased to 2 years for type 2 diabetic patients without retinopathy? Diabetes Care 35 : 1663 -1668 Olafsdottir E , Stefansson E ( 2007 ) Biennial eye screening in patients with diabetes without retinopathy: 10-year experience . Br J Ophthalmol 91 : 1599 -1601 Looker HC , Nyangoma SO , Cromie DT et al ( 2013 ) Predicted impact of extending the screening interval for diabetic retinopathy: the Scottish Diabetic Retinopathy Screening programme . Aspelund T , Þórisdóttir O , Ólafsdottir E et al ( 2011 ) Individual risk assessment and information technology to optimise screening frequency for diabetic retinopathy . Diabetologia 54 : 2525 - 2532 Papoulis A ( 1991 ) Probability, random variables, and stochastic processes, 3rd edn . McGraw-Hill , New York Stratton IM, Aldington SJ , Taylor DJ , Adler AI , Scanlon PS ( 2013 ) A simple risk stratification for time to development of sightthreatening diabetic retinopathy . Diabetes Care 36 : 580 -585 Kalbfleisch JD , Lawless JF ( 1985 ) The analysis of panel data under a Markov assumption . J Am Stat Assoc 80 : 863 - 871 Kay R ( 1986 ) A Markov model for analysing cancer markers and disease states in survival studies . Biometrics 42 : 855 -865 Titman AC , Sharples LD ( 2010 ) Model diagnostics for multi-state models . Stat Methods Med Res 19 : 621 -651 van Buuren S ( 2012 ) Flexible imputation of missing data . Chapman & Hall/CRC Press, Boca Raton Harrell FE ( 2001 ) Regression modelling strategies: with applications to linear models, logistic regression, and survival analysis .

This is a preview of a remote PDF:

Antonio Eleuteri, Anthony C. Fisher, Deborah M. Broadbent, Marta García-Fiñana, Christopher P. Cheyne, Amu Wang, Irene M. Stratton, Mark Gabbay, Daniel Seddon, Simon P. Harding, for the Individualised Screening for Diabetic Retinopathy (ISDR) Study Group. Individualised variable-interval risk-based screening for sight-threatening diabetic retinopathy: the Liverpool Risk Calculation Engine, Diabetologia, 2017, 1-9, DOI: 10.1007/s00125-017-4386-0