Knitting Statistics Notes: Creating Teaching Materials with Data Analysis Code and Results Embedded

R is a widely used, free, open source software package for data analysis and graphics. Along with core functionality, there are many additional packages to enable specialist usage. One of these packages knitr' allows the integration of R and LaTeX code. This means that code, output and narrative are all included in one source document. Output is generated automatically from the code. If any changes are made to the data, all output is updated automatically. From a teaching perspective, this is particularly useful if you want to create notes with various graphs, statistics etc and you also want to provide the code that shows the student how to do this analysis in R. It is possible to hide either the code or the output. For example, if you are making introduction to statistics worksheets for non-statistics students, you might not want to include R code. If you are including exercises for students, you might not want to include the output and instead have them answer some questions. You can also set these options globally then by changing one option, create two versions of the document - one with solutions, one without, just from one source document

Myocardial infarction and stroke subsequent to urinary tract infection (MISSOURI): protocol for a self-controlled case series using linked electronic health records

Introduction: There is increasing interest in the relationship between acute infections and acute cardiovascular events. Most previous research has focused on understanding whether the risk of acute cardiovascular events increases following a respiratory tract infection. The relationship between urinary tract infections (UTIs) and acute cardiovascular events is less well studied. Therefore, the aim of this study is to determine whether there is a causal relationship between UTI and acute myocardial infarction (MI) or stroke. Methods and analysis: We will undertake a self-controlled case series study using linked anonymised general practice, hospital admission and microbiology data held within the Secure Anonymised Information Linkage (SAIL) Databank. Self-controlled case series is a relatively novel study design where individuals act as their own controls, thereby inherently controlling for time-invariant confounders. Only individuals who experience an exposure and outcome of interest are included. We will identify individuals in the SAIL Databank who have a hospital admission record for acute MI or stroke during the study period of 2010–2020. Individuals will need to be aged 30–100 during the study period and be Welsh residents for inclusion. UTI will be identified using general practice, microbiology and hospital admissions data. We will calculate the incidence of MI and stroke in predefined risk periods following an UTI and in ‘baseline’ periods (without UTI exposure) and use conditional Poisson regression models to derive incidence rate ratios. Ethics and dissemination: Data access, research permissions and approvals have been obtained from the SAIL independent Information Governance Review Panel, project number 0972. Findings will be disseminated through conferences, blogs, social media threads and peer-reviewed journals. Results will be of interest internationally to primary and secondary care clinicians who manage UTIs and may inform future clinical trials of preventative therapy

A Four-Way Comparison of Cardiac Function with Normobaric Normoxia, Normobaric Hypoxia, Hypobaric Hypoxia and Genuine High Altitude.

There has been considerable debate as to whether different modalities of simulated hypoxia induce similar cardiac responses.This was a prospective observational study of 14 healthy subjects aged 22-35 years. Echocardiography was performed at rest and at 15 and 120 minutes following two hours exercise under normobaric normoxia (NN) and under similar PiO2 following genuine high altitude (GHA) at 3,375m, normobaric hypoxia (NH) and hypobaric hypoxia (HH) to simulate the equivalent hypoxic stimulus to GHA.All 14 subjects completed the experiment at GHA, 11 at NN, 12 under NH, and 6 under HH. The four groups were similar in age, sex and baseline demographics. At baseline rest right ventricular (RV) systolic pressure (RVSP, p = 0.0002), pulmonary vascular resistance (p = 0.0002) and acute mountain sickness (AMS) scores were higher and the SpO2 lower (p<0.0001) among all three hypoxic groups (GHA, NH and HH) compared with NN. At both 15 minutes and 120 minutes post exercise, AMS scores, Cardiac output, septal S', lateral S', tricuspid S' and A' velocities and RVSP were higher and SpO2 lower with all forms of hypoxia compared with NN. On post-test analysis, among the three hypoxia groups, SpO2 was lower at baseline and 15 minutes post exercise with GHA (89.3±3.4% and 89.3±2.2%) and HH (89.0±3.1 and (89.8±5.0) compared with NH (92.9±1.7 and 93.6±2.5%). The RV Myocardial Performance (Tei) Index and RVSP were significantly higher with HH than NH at 15 and 120 minutes post exercise respectively and tricuspid A' was higher with GHA compared with NH at 15 minutes post exercise.GHA, NH and HH produce similar cardiac adaptations over short duration rest despite lower SpO2 levels with GHA and HH compared with NH. Notable differences emerge following exercise in SpO2, RVSP and RV cardiac function

New genetic loci link adipose and insulin biology to body fat distribution.

Body fat distribution is a heritable trait and a well-established predictor of adverse metabolic outcomes, independent of overall adiposity. To increase our understanding of the genetic basis of body fat distribution and its molecular links to cardiometabolic traits, here we conduct genome-wide association meta-analyses of traits related to waist and hip circumferences in up to 224,459 individuals. We identify 49 loci (33 new) associated with waist-to-hip ratio adjusted for body mass index (BMI), and an additional 19 loci newly associated with related waist and hip circumference measures (P < 5 × 10(-8)). In total, 20 of the 49 waist-to-hip ratio adjusted for BMI loci show significant sexual dimorphism, 19 of which display a stronger effect in women. The identified loci were enriched for genes expressed in adipose tissue and for putative regulatory elements in adipocytes. Pathway analyses implicated adipogenesis, angiogenesis, transcriptional regulation and insulin resistance as processes affecting fat distribution, providing insight into potential pathophysiological mechanisms

Effect of Hydrocortisone on Mortality and Organ Support in Patients With Severe COVID-19: The REMAP-CAP COVID-19 Corticosteroid Domain Randomized Clinical Trial.

Importance: Evidence regarding corticosteroid use for severe coronavirus disease 2019 (COVID-19) is limited. Objective: To determine whether hydrocortisone improves outcome for patients with severe COVID-19. Design, Setting, and Participants: An ongoing adaptive platform trial testing multiple interventions within multiple therapeutic domains, for example, antiviral agents, corticosteroids, or immunoglobulin. Between March 9 and June 17, 2020, 614 adult patients with suspected or confirmed COVID-19 were enrolled and randomized within at least 1 domain following admission to an intensive care unit (ICU) for respiratory or cardiovascular organ support at 121 sites in 8 countries. Of these, 403 were randomized to open-label interventions within the corticosteroid domain. The domain was halted after results from another trial were released. Follow-up ended August 12, 2020. Interventions: The corticosteroid domain randomized participants to a fixed 7-day course of intravenous hydrocortisone (50 mg or 100 mg every 6 hours) (n = 143), a shock-dependent course (50 mg every 6 hours when shock was clinically evident) (n = 152), or no hydrocortisone (n = 108). Main Outcomes and Measures: The primary end point was organ support-free days (days alive and free of ICU-based respiratory or cardiovascular support) within 21 days, where patients who died were assigned -1 day. The primary analysis was a bayesian cumulative logistic model that included all patients enrolled with severe COVID-19, adjusting for age, sex, site, region, time, assignment to interventions within other domains, and domain and intervention eligibility. Superiority was defined as the posterior probability of an odds ratio greater than 1 (threshold for trial conclusion of superiority >99%). Results: After excluding 19 participants who withdrew consent, there were 384 patients (mean age, 60 years; 29% female) randomized to the fixed-dose (n = 137), shock-dependent (n = 146), and no (n = 101) hydrocortisone groups; 379 (99%) completed the study and were included in the analysis. The mean age for the 3 groups ranged between 59.5 and 60.4 years; most patients were male (range, 70.6%-71.5%); mean body mass index ranged between 29.7 and 30.9; and patients receiving mechanical ventilation ranged between 50.0% and 63.5%. For the fixed-dose, shock-dependent, and no hydrocortisone groups, respectively, the median organ support-free days were 0 (IQR, -1 to 15), 0 (IQR, -1 to 13), and 0 (-1 to 11) days (composed of 30%, 26%, and 33% mortality rates and 11.5, 9.5, and 6 median organ support-free days among survivors). The median adjusted odds ratio and bayesian probability of superiority were 1.43 (95% credible interval, 0.91-2.27) and 93% for fixed-dose hydrocortisone, respectively, and were 1.22 (95% credible interval, 0.76-1.94) and 80% for shock-dependent hydrocortisone compared with no hydrocortisone. Serious adverse events were reported in 4 (3%), 5 (3%), and 1 (1%) patients in the fixed-dose, shock-dependent, and no hydrocortisone groups, respectively. Conclusions and Relevance: Among patients with severe COVID-19, treatment with a 7-day fixed-dose course of hydrocortisone or shock-dependent dosing of hydrocortisone, compared with no hydrocortisone, resulted in 93% and 80% probabilities of superiority with regard to the odds of improvement in organ support-free days within 21 days. However, the trial was stopped early and no treatment strategy met prespecified criteria for statistical superiority, precluding definitive conclusions. Trial Registration: ClinicalTrials.gov Identifier: NCT02735707

Effect of remote ischaemic conditioning on clinical outcomes in patients with acute myocardial infarction (CONDI-2/ERIC-PPCI): a single-blind randomised controlled trial.

BACKGROUND: Remote ischaemic conditioning with transient ischaemia and reperfusion applied to the arm has been shown to reduce myocardial infarct size in patients with ST-elevation myocardial infarction (STEMI) undergoing primary percutaneous coronary intervention (PPCI). We investigated whether remote ischaemic conditioning could reduce the incidence of cardiac death and hospitalisation for heart failure at 12 months. METHODS: We did an international investigator-initiated, prospective, single-blind, randomised controlled trial (CONDI-2/ERIC-PPCI) at 33 centres across the UK, Denmark, Spain, and Serbia. Patients (age >18 years) with suspected STEMI and who were eligible for PPCI were randomly allocated (1:1, stratified by centre with a permuted block method) to receive standard treatment (including a sham simulated remote ischaemic conditioning intervention at UK sites only) or remote ischaemic conditioning treatment (intermittent ischaemia and reperfusion applied to the arm through four cycles of 5-min inflation and 5-min deflation of an automated cuff device) before PPCI. Investigators responsible for data collection and outcome assessment were masked to treatment allocation. The primary combined endpoint was cardiac death or hospitalisation for heart failure at 12 months in the intention-to-treat population. This trial is registered with ClinicalTrials.gov (NCT02342522) and is completed. FINDINGS: Between Nov 6, 2013, and March 31, 2018, 5401 patients were randomly allocated to either the control group (n=2701) or the remote ischaemic conditioning group (n=2700). After exclusion of patients upon hospital arrival or loss to follow-up, 2569 patients in the control group and 2546 in the intervention group were included in the intention-to-treat analysis. At 12 months post-PPCI, the Kaplan-Meier-estimated frequencies of cardiac death or hospitalisation for heart failure (the primary endpoint) were 220 (8·6%) patients in the control group and 239 (9·4%) in the remote ischaemic conditioning group (hazard ratio 1·10 [95% CI 0·91-1·32], p=0·32 for intervention versus control). No important unexpected adverse events or side effects of remote ischaemic conditioning were observed. INTERPRETATION: Remote ischaemic conditioning does not improve clinical outcomes (cardiac death or hospitalisation for heart failure) at 12 months in patients with STEMI undergoing PPCI. FUNDING: British Heart Foundation, University College London Hospitals/University College London Biomedical Research Centre, Danish Innovation Foundation, Novo Nordisk Foundation, TrygFonden

Effect of angiotensin-converting enzyme inhibitor and angiotensin receptor blocker initiation on organ support-free days in patients hospitalized with COVID-19

IMPORTANCE Overactivation of the renin-angiotensin system (RAS) may contribute to poor clinical outcomes in patients with COVID-19. Objective To determine whether angiotensin-converting enzyme (ACE) inhibitor or angiotensin receptor blocker (ARB) initiation improves outcomes in patients hospitalized for COVID-19. DESIGN, SETTING, AND PARTICIPANTS In an ongoing, adaptive platform randomized clinical trial, 721 critically ill and 58 non–critically ill hospitalized adults were randomized to receive an RAS inhibitor or control between March 16, 2021, and February 25, 2022, at 69 sites in 7 countries (final follow-up on June 1, 2022). INTERVENTIONS Patients were randomized to receive open-label initiation of an ACE inhibitor (n = 257), ARB (n = 248), ARB in combination with DMX-200 (a chemokine receptor-2 inhibitor; n = 10), or no RAS inhibitor (control; n = 264) for up to 10 days. MAIN OUTCOMES AND MEASURES The primary outcome was organ support–free days, a composite of hospital survival and days alive without cardiovascular or respiratory organ support through 21 days. The primary analysis was a bayesian cumulative logistic model. Odds ratios (ORs) greater than 1 represent improved outcomes. RESULTS On February 25, 2022, enrollment was discontinued due to safety concerns. Among 679 critically ill patients with available primary outcome data, the median age was 56 years and 239 participants (35.2%) were women. Median (IQR) organ support–free days among critically ill patients was 10 (–1 to 16) in the ACE inhibitor group (n = 231), 8 (–1 to 17) in the ARB group (n = 217), and 12 (0 to 17) in the control group (n = 231) (median adjusted odds ratios of 0.77 [95% bayesian credible interval, 0.58-1.06] for improvement for ACE inhibitor and 0.76 [95% credible interval, 0.56-1.05] for ARB compared with control). The posterior probabilities that ACE inhibitors and ARBs worsened organ support–free days compared with control were 94.9% and 95.4%, respectively. Hospital survival occurred in 166 of 231 critically ill participants (71.9%) in the ACE inhibitor group, 152 of 217 (70.0%) in the ARB group, and 182 of 231 (78.8%) in the control group (posterior probabilities that ACE inhibitor and ARB worsened hospital survival compared with control were 95.3% and 98.1%, respectively). CONCLUSIONS AND RELEVANCE In this trial, among critically ill adults with COVID-19, initiation of an ACE inhibitor or ARB did not improve, and likely worsened, clinical outcomes. TRIAL REGISTRATION ClinicalTrials.gov Identifier: NCT0273570

Spatial analysis of health effects of large industrial incinerators in England, 1998-2008:a study using matched case-control areas

To assess whether residential proximity to industrial incinerators in England is associated with increased risk of cancer incidence and mortality

Giardiasis in North West England:faecal specimen requesting rates by GP practice

SUMMARY Many cases of giardiasis in the UK are undiagnosed and among other things, diagnosis is dependent upon the readiness of GPs to request a specimen. The aim of this study is to assess the rate of specimens requested per GP practice in Central Lancashire, to examine the differences between GP practices and to estimate the pattern of unexplained spatial variation in the practice rate of specimens after adjustment for deprivation. To achieve this, we fitted a set of binomial and Poisson regression models, with random effects for GP practice. Our analysis suggests that there were differences in the rate of specimens by GP practices (P < 0·001) for a single year, but no difference in the proportion of positive tests per specimen submitted or in the rate of positive specimens per practice population. There was a difference in the cumulative rate of positive specimens per practice population over a 9-year period (P < 0·001). Neither the specimen rate per practice for a single year nor the cumulative rate of positive specimens over multiple years demonstrated significant spatial correlation. Hence, spatial variation in the incidence of giardiasis is unlikely to be confounded by variation in GP rate of specimens