Estimation of the linear mixed integrated Ornstein-Uhlenbeck model.

The linear mixed model with an added integrated Ornstein-Uhlenbeck (IOU) process (linear mixed IOU model) allows for serial correlation and estimation of the degree of derivative tracking. It is rarely used, partly due to the lack of available software. We implemented the linear mixed IOU model in Stata and using simulations we assessed the feasibility of fitting the model by restricted maximum likelihood when applied to balanced and unbalanced data. We compared different (1) optimization algorithms, (2) parameterizations of the IOU process, (3) data structures and (4) random-effects structures. Fitting the model was practical and feasible when applied to large and moderately sized balanced datasets (20,000 and 500 observations), and large unbalanced datasets with (non-informative) dropout and intermittent missingness. Analysis of a real dataset showed that the linear mixed IOU model was a better fit to the data than the standard linear mixed model (i.e. independent within-subject errors with constant variance)

NHS waiting lists and evidence of national or local failure: analysis of health service data.

OBJECTIVES: To investigate the national distribution of prolonged waiting for elective day case and inpatient surgery, and to examine associations of prolonged waiting with markers of NHS capacity, activity in the independent sector, and need. SETTING: NHS hospital trusts in England. POPULATION: People waiting for elective treatment in the specialties of general surgery; ear, nose and throat surgery; ophthalmic surgery; and trauma and orthopaedic surgery. MAIN OUTCOME MEASURE: Numbers of people waiting six months or longer (prolonged waiting). Characteristics of trusts with large numbers waiting six months or longer were examined by using logistic regression. RESULTS: The distribution of numbers of people waiting for day case or elective surgery in all the specialties examined was highly positively skewed. Between 52% and 83% of patients waiting longer than six months in the specialties studied were found in one quarter of trusts, which in turn contributed 23-45% of the national throughput specific to the specialty. In general, there was little evidence to show that capacity (measured by numbers of operating theatres, dedicated day case theatres, available beds, and bed occupancy rate) or independent sector activity were associated with prolonged waiting, although exceptions were noted for individual specialties. There was consistent evidence showing an increase in prolonged waiting, with increased numbers of anaesthetists across all specialties and with increased bed occupancy rates for ear, nose and throat surgery. Markers of greater need for health care, such as deprivation score and rate of limiting long term illness, were inversely associated with prolonged waiting. CONCLUSION: In most instances, substantial numbers of patients waiting unacceptably long periods for elective surgery were limited to a small number of hospitals. Little and inconsistent support was found for associations of prolonged waiting with markers of capacity, independent sector activity, or need in the surgical specialties examined

The association between BMI and mortality using offspring BMI as an indicator of own BMI: large intergenerational mortality study

Objectives To obtain valid estimates of the association between body mass index (BMI) and mortality by using offspring BMI as an instrumental variable for own BMI

Joint modelling rationale for chained equations.

BACKGROUND: Chained equations imputation is widely used in medical research. It uses a set of conditional models, so is more flexible than joint modelling imputation for the imputation of different types of variables (e.g. binary, ordinal or unordered categorical). However, chained equations imputation does not correspond to drawing from a joint distribution when the conditional models are incompatible. Concurrently with our work, other authors have shown the equivalence of the two imputation methods in finite samples. METHODS: Taking a different approach, we prove, in finite samples, sufficient conditions for chained equations and joint modelling to yield imputations from the same predictive distribution. Further, we apply this proof in four specific cases and conduct a simulation study which explores the consequences when the conditional models are compatible but the conditions otherwise are not satisfied. RESULTS: We provide an additional "non-informative margins" condition which, together with compatibility, is sufficient. We show that the non-informative margins condition is not satisfied, despite compatible conditional models, in a situation as simple as two continuous variables and one binary variable. Our simulation study demonstrates that as a consequence of this violation order effects can occur; that is, systematic differences depending upon the ordering of the variables in the chained equations algorithm. However, the order effects appear to be small, especially when associations between variables are weak. CONCLUSIONS: Since chained equations is typically used in medical research for datasets with different types of variables, researchers must be aware that order effects are likely to be ubiquitous, but our results suggest they may be small enough to be negligible

Association Between Risk-of-Bias Assessments and Results of Randomized Trials in Cochrane Reviews: The ROBES Meta-Epidemiologic Study.

Flaws in the design of randomized trials may bias intervention effect estimates and increase between-trial heterogeneity. Empirical evidence suggests that these problems are greatest for subjectively assessed outcomes. For the Risk of Bias in Evidence Synthesis (ROBES) Study, we extracted risk-of-bias judgements (for sequence generation, allocation concealment, blinding, and incomplete data) from a large collection of meta-analyses published in the Cochrane Library (issue 4; April 2011). We categorized outcome measures as mortality, other objective outcome, or subjective outcome, and we estimated associations of bias judgements with intervention effect estimates using Bayesian hierarchical models. Among 2,443 randomized trials in 228 meta-analyses, intervention effect estimates were, on average, exaggerated in trials with high or unclear (versus low) risk-of-bias judgements for sequence generation (ratio of odds ratios (ROR) = 0.91, 95% credible interval (CrI): 0.86, 0.98), allocation concealment (ROR = 0.92, 95% CrI: 0.86, 0.98), and blinding (ROR = 0.87, 95% CrI: 0.80, 0.93). In contrast to previous work, we did not observe consistently different bias for subjective outcomes compared with mortality. However, we found an increase in between-trial heterogeneity associated with lack of blinding in meta-analyses with subjective outcomes. Inconsistency in criteria for risk-of-bias judgements applied by individual reviewers is a likely limitation of routinely collected bias assessments. Inadequate randomization and lack of blinding may lead to exaggeration of intervention effect estimates in randomized trials

All-cause hospitalisation among people living with HIV according to gender, mode of HIV acquisition, ethnicity, and geographical origin in Europe and North America: findings from the ART-CC cohort collaboration

BACKGROUND: Understanding demographic disparities in hospitalisation is crucial for the identification of vulnerable populations, interventions, and resource planning. METHODS: Data were from the Antiretroviral Therapy Cohort Collaboration (ART-CC) on people living with HIV in Europe and North America, followed up between January, 2007 and December, 2020. We investigated differences in all-cause hospitalisation according to gender and mode of HIV acquisition, ethnicity, and combined geographical origin and ethnicity, in people living with HIV on modern combination antiretroviral therapy (cART). Analyses were performed separately for European and North American cohorts. Hospitalisation rates were assessed using negative binomial multilevel regression, adjusted for age, time since cART intitiaion, and calendar year. FINDINGS: Among 23 594 people living with HIV in Europe and 9612 in North America, hospitalisation rates per 100 person-years were 16·2 (95% CI 16·0-16·4) and 13·1 (12·8-13·5). Compared with gay, bisexual, and other men who have sex with men, rates were higher for heterosexual men and women, and much higher for men and women who acquired HIV through injection drug use (adjusted incidence rate ratios ranged from 1·2 to 2·5 in Europe and from 1·2 to 3·3 in North America). In both regions, individuals with geographical origin other than the region of study generally had lower hospitalisation rates compared with those with geographical origin of the study country. In North America, Indigenous people and Black or African American individuals had higher rates than White individuals (adjusted incidence rate ratios 1·9 and 1·2), whereas Asian and Hispanic people living with HIV had somewhat lower rates. In Europe there was a lower rate in Asian individuals compared with White individuals. INTERPRETATION: Substantial disparities exist in all-cause hospitalisation between demographic groups of people living with HIV in the current cART era in high-income settings, highlighting the need for targeted support. FUNDING: Royal Free Charity and the National Institute on Alcohol Abuse and Alcoholism

Oral anticoagulants for primary prevention, treatment and secondary prevention of venous thromboembolic disease, and for prevention of stroke in atrial fibrillation:systematic review, network meta-analysis and cost-effectiveness analysis

BACKGROUND: Warfarin is effective for stroke prevention in atrial fibrillation (AF), but anticoagulation is underused in clinical care. The risk of venous thromboembolic disease during hospitalisation can be reduced by low-molecular-weight heparin (LMWH): warfarin is the most frequently prescribed anticoagulant for treatment and secondary prevention of venous thromboembolism (VTE). Warfarin-related bleeding is a major reason for hospitalisation for adverse drug effects. Warfarin is cheap but therapeutic monitoring increases treatment costs. Novel oral anticoagulants (NOACs) have more rapid onset and offset of action than warfarin, and more predictable dosing requirements.OBJECTIVE: To determine the best oral anticoagulant/s for prevention of stroke in AF and for primary prevention, treatment and secondary prevention of VTE.DESIGN: Four systematic reviews, network meta-analyses (NMAs) and cost-effectiveness analyses (CEAs) of randomised controlled trials.SETTING: Hospital (VTE primary prevention and acute treatment) and primary care/anticoagulation clinics (AF and VTE secondary prevention).PARTICIPANTS: Patients eligible for anticoagulation with warfarin (stroke prevention in AF, acute treatment or secondary prevention of VTE) or LMWH (primary prevention of VTE).INTERVENTIONS: NOACs, warfarin and LMWH, together with other interventions (antiplatelet therapy, placebo) evaluated in the evidence network.MAIN OUTCOME MEASURES: Efficacy Stroke, symptomatic VTE, symptomatic deep-vein thrombosis and symptomatic pulmonary embolism. Safety Major bleeding, clinically relevant bleeding and intracranial haemorrhage. We also considered myocardial infarction and all-cause mortality and evaluated cost-effectiveness.DATA SOURCES: MEDLINE and PREMEDLINE In-Process & Other Non-Indexed Citations, EMBASE and The Cochrane Library, reference lists of published NMAs and trial registries. We searched MEDLINE and PREMEDLINE In-Process & Other Non-Indexed Citations, EMBASE and The Cochrane Library. The stroke prevention in AF review search was run on the 12 March 2014 and updated on 15 September 2014, and covered the period 2010 to September 2014. The search for the three reviews in VTE was run on the 19 March 2014, updated on 15 September 2014, and covered the period 2008 to September 2014.REVIEW METHODS: Two reviewers screened search results, extracted and checked data, and assessed risk of bias. For each outcome we conducted standard meta-analysis and NMA. We evaluated cost-effectiveness using discrete-time Markov models.RESULTS: Apixaban (Eliquis(rcledR), Bristol-Myers Squibb, USA; Pfizer, USA) [5 mg bd (twice daily)] was ranked as among the best interventions for stroke prevention in AF, and had the highest expected net benefit. Edoxaban (Lixiana(rcledR), Daiichi Sankyo, Japan) [60 mg od (once daily)] was ranked second for major bleeding and all-cause mortality. Neither the clinical effectiveness analysis nor the CEA provided strong evidence that NOACs should replace postoperative LMWH in primary prevention of VTE. For acute treatment and secondary prevention of VTE, we found little evidence that NOACs offer an efficacy advantage over warfarin, but the risk of bleeding complications was lower for some NOACs than for warfarin. For a willingness-to-pay threshold of > £5000, apixaban (5 mg bd) had the highest expected net benefit for acute treatment of VTE. Aspirin or no pharmacotherapy were likely to be the most cost-effective interventions for secondary prevention of VTE: our results suggest that it is not cost-effective to prescribe NOACs or warfarin for this indication.CONCLUSIONS: NOACs have advantages over warfarin in patients with AF, but we found no strong evidence that they should replace warfarin or LMWH in primary prevention, treatment or secondary prevention of VTE.LIMITATIONS: These relate mainly to shortfalls in the primary data: in particular, there were no head-to-head comparisons between different NOAC drugs.FUTURE WORK: Calculating the expected value of sample information to clarify whether or not it would be justifiable to fund one or more head-to-head trials.STUDY REGISTRATION: This study is registered as PROSPERO CRD42013005324, CRD42013005331 and CRD42013005330.FUNDING: The National Institute for Health Research Health Technology Assessment programme

Ongoing monitoring of data clustering in multicenter studies

Background: Multicenter study designs have several advantages, but the possibility of non-random measurement error resulting from procedural differences between the centers is a special concern. While it is possible to address and correct for some measurement error through statistical analysis, proactive data monitoring is essential to ensure high-quality data collection. Methods: In this article, we describe quality assurance efforts aimed at reducing the effect of measurement error in a recent follow-up of a large cluster-randomized controlled trial through periodic evaluation of intraclass correlation coefficients (ICCs) for continuous measurements. An ICC of 0 indicates the variance in the data is not due to variation between the centers, and thus the data are not clustered by center. Results: Through our review of early data downloads, we identified several outcomes (including sitting height, waist circumference, and systolic blood pressure) with higher than expected ICC values. Further investigation revealed variations in the procedures used by pediatricians to measure these outcomes. We addressed these procedural inconsistencies through written clarification of the protocol and refresher training workshops with the pediatricians. Further data monitoring at subsequent downloads showed that these efforts had a beneficial effect on data quality (sitting height ICC decreased from 0.92 to 0.03, waist circumference from 0.10 to 0.07, and systolic blood pressure from 0.16 to 0.12). Conclusions: We describe a simple but formal mechanism for identifying ongoing problems during data collection. The calculation of the ICC can easily be programmed and the mechanism has wide applicability, not just to cluster randomized controlled trials but to any study with multiple centers or with multiple observers

The duration of protection of school-aged BCG vaccination in England: a population-based case-control study.

Background: Evidence of protection from childhood Bacillus Calmette-Guerin (BCG) against tuberculosis (TB) in adulthood, when most transmission occurs, is important for TB control and resource allocation. Methods: We conducted a population-based case-control study of protection by BCG given to children aged 12-13 years against tuberculosis occurring 10-29 years later. We recruited UK-born White subjects with tuberculosis and randomly sampled White community controls. Hazard ratios and 95% confidence intervals (CIs) were estimated using case-cohort Cox regression, adjusting for potential confounding factors, including socio-economic status, smoking, drug use, prison and homelessness. Vaccine effectiveness (VE = 1 - hazard ratio) was assessed at successive intervals more than 10 years following vaccination. Results: We obtained 677 cases and 1170 controls after a 65% response rate in both groups. Confounding by deprivation, education and lifestyle factors was slight 10-20 years after vaccination, and more evident after 20 years. VE 10-15 years after vaccination was 51% (95% CI 21, 69%) and 57% (CI 33, 72%) at 15-20 years. Subsequently, BCG protection appeared to wane; 20-25 years VE = 25% (CI -14%, 51%) and 25-29 years VE = 1% (CI -84%, 47%). Based on multiple imputation of missing data (in 17% subjects), VE estimated in the same intervals after vaccination were similar [56% (CI 33, 72%), 57% (CI 36, 71%), 25% (-10, 48%), 21% (-39, 55%)]. Conclusions: School-aged BCG vaccination offered moderate protection against tuberculosis for at least 20 years, which is longer than previously thought. This has implications for assessing the cost-effectiveness of BCG vaccination and when evaluating new TB vaccines

Heterogeneity in outcomes of treated HIV-positive patients in Europe and North America: relation with patient and cohort characteristics

Background HIV cohort collaborations, which pool data from diverse patient cohorts, have provided key insights into outcomes of antiretroviral therapy (ART). However, the extent of, and reasons for, between-cohort heterogeneity in rates of AIDS and mortality are unclear. Methods We obtained data on adult HIV-positive patients who started ART from 1998 without a previous AIDS diagnosis from 17 cohorts in North America and Europe. Patients were followed up from 1 month to 2 years after starting ART. We examined between-cohort heterogeneity in crude and adjusted (age, sex, HIV transmission risk, year, CD4 count and HIV-1 RNA at start of ART) rates of AIDS and mortality using random-effects meta-analysis and meta-regression. Results During 61 520 person-years, 754/38 706 (1.9%) patients died and 1890 (4.9%) progressed to AIDS. Between-cohort variance in mortality rates was reduced from 0.84 to 0.24 (0.73 to 0.28 for AIDS rates) after adjustment for patient characteristics. Adjusted mortality rates were inversely associated with cohorts' estimated completeness of death ascertainment [excellent: 96-100%, good: 90-95%, average: 75-89%; mortality rate ratio 0.66 (95% confidence interval 0.46-0.94) per category]. Mortality rate ratios comparing Europe with North America were 0.42 (0.31-0.57) before and 0.47 (0.30-0.73) after adjusting for completeness of ascertainment. Conclusions Heterogeneity between settings in outcomes of HIV treatment has implications for collaborative analyses, policy and clinical care. Estimated mortality rates may require adjustment for completeness of ascertainment. Higher mortality rate in North American, compared with European, cohorts was not fully explained by completeness of ascertainment and may be because of the inclusion of more socially marginalized patients with higher mortality ris