Outlier effects on robust joint modelling of longitudinal and survival data

Robust joint modelling is an emerging field of research. Through the advancements in electronic patient healthcare records, the popularly of joint modelling approaches has grown rapidly in recent years providing simultaneous analysis of longitudinal and survival data. This research advances previous work through the development of a novel robust joint modelling methodology for one of the most common types of standard joint models, that which links a linear mixed model with a Cox proportional hazards model. Through t-distributional assumptions, longitudinal outliers are accommodated with their detrimental impact being down weighed and thus providing more efficient and reliable estimates. The robust joint modelling technique and its major benefits are showcased through the analysis of Northern Irish end stage renal disease patients. With an ageing population and growing prevalence of chronic kidney disease within the United Kingdom, there is a pressing demand to investigate the detrimental relationship between the changing haemoglobin levels of haemodialysis patients and their survival. As outliers within the NI renal data were found to have significantly worse survival, identification of outlying individuals through robust joint modelling may aid nephrologists to improve patient's survival. A simulation study was also undertaken to explore the difference between robust and standard joint models in the presence of increasing proportions and extremity of longitudinal outliers. More efficient and reliable estimates were obtained by robust joint models with increasing contrast between the robust and standard joint models when a greater proportion of more extreme outliers are present. Through illustration of the gains in efficiency and reliability of parameters when outliers exist, the potential of robust joint modelling is evident. The research presented in this thesis highlights the benefits and stresses the need to utilise a more robust approach to joint modelling in the presence of longitudinal outliers.

HDL-C and HDL-C/ApoA-I predict long-term progression of glycemia in established type 2 diabetes

OBJECTIVE: Low HDL cholesterol (HDL-C) and small HDL particle size may directly promote hyperglycemia. We evaluated associations of HDL-C, apolipoprotein A-I (apoA-I), and HDL-C/apoA-I with insulin secretion, insulin resistance, HbA1c, and long-term glycemic deterioration, reflected by initiation of pharmacologic glucose control.

RESEARCH DESIGN AND METHODS: The 5-year Fenofibrate Intervention and Event Lowering in Diabetes (FIELD) study followed 9,795 type 2 diabetic subjects. We calculated baseline associations of fasting HDL-C, apoA-I, and HDL-C/apoA-I with HbA1c and, in those not taking exogenous insulin (n = 8,271), with estimated β-cell function (homeostasis model assessment of β-cell function [HOMA-B]) and insulin resistance (HOMA-IR). Among the 2,608 subjects prescribed lifestyle only, Cox proportional hazards analysis evaluated associations of HDL-C, apoA-I, and HDL-C/apoA-I with subsequent initiation of oral hypoglycemic agents (OHAs) or insulin.

RESULTS: Adjusted for age and sex, baseline HDL-C, apoA-I, and HDL-C/apoA-I were inversely associated with HOMA-IR (r = -0.233, -0.134, and -0.230; all P < 0.001; n = 8,271) but not related to HbA1c (all P > 0.05; n = 9,795). ApoA-I was also inversely associated with HOMA-B (r = -0.063; P = 0.002; n = 8,271) adjusted for age, sex, and HOMA-IR. Prospectively, lower baseline HDL-C and HDL-C/apoA-I levels predicted greater uptake (per 1-SD lower: hazard ratio [HR] 1.13 [CI 1.07-1.19], P < 0.001; and HR 1.16 [CI 1.10-1.23], P < 0.001, respectively) and earlier uptake (median 12.9 and 24.0 months, respectively, for quartile 1 vs. quartile 4; both P < 0.01) of OHAs and insulin, with no difference in HbA1c thresholds for initiation (P = 0.87 and P = 0.81). Controlling for HOMA-IR and triglycerides lessened both associations, but HDL-C/apoA-I remained significant.

CONCLUSIONS: HDL-C, apoA-I, and HDL-C/apoA-I were associated with concurrent insulin resistance but not HbA1c. However, lower HDL-C and HDL-C/apoA-I predicted greater and earlier need for pharmacologic glucose control.

HER2 in high-risk rectal cancer patients treated in EXPERT-C, a randomized phase II trial of neoadjuvant capecitabine and oxaliplatin (CAPOX) and chemoradiotherapy (CRT) with or without cetuximab.

BACKGROUND: HER2 is an established therapeutic target in breast and gastric cancers. The role of HER2 in rectal cancer is unclear, as conflicting data on the prevalence of HER2 expression in this disease have been reported. We evaluated the prevalence of HER2 and its impact on the outcome of high-risk rectal cancer patients treated with neoadjuvant CAPOX and CRT±cetuximab in the EXPERT-C trial. PATIENTS AND METHODS: Eligible patients with available tumour tissue for HER2 analysis were included. HER2 expression was determined by immunohistochemistry (IHC) in pre-treatment biopsies and/or surgical specimens (score 0-3+). Immunostaining was scored according to the consensus panel recommendations on HER2 scoring for gastric cancer. Tumours with equivocal IHC result (2+) were tested for HER2 amplification by D-ISH. Tumours with IHC 3+ or D-ISH ratio ≥2.0 were classified as HER2+. The impact of HER2 on primary and secondary end points of the study was analysed. RESULTS: Of 164 eligible study patients, 104 (63%) biopsy and 114 (69%) surgical specimens were available for analysis. Only 3 of 104 (2.9%) and 3 of 114 (2.6%) were HER2+, respectively. In 77 patients with paired specimens, concordance for HER2 status was found in 74 (96%). Overall, 141 patients were assessable for HER2 and 6 out of 141 (4.3%) had HER2 overexpression and/or amplification. The median follow-up was 58.6 months. HER2 was not associated with a difference in the outcome for any of the study end points, including in the subset of 90 KRAS/BRAF wild-type patients treated±cetuximab. CONCLUSIONS: Based on the low prevalence of expression as recorded in the EXPERT-C trial, HER2 does not appear to represent a useful therapeutic target in high-risk rectal cancer. However, the role of HER2 as a potential predictive biomarker of resistance to anti-EGFR-based treatments and a therapeutic target in anti-EGFR refractory metastatic colorectal cancer (CRC) warrants further investigation. TRIAL REGISTRATION: ISRCTN Register: 99828560.

Epirubicin, oxaliplatin, and capecitabine with or without panitumumab for patients with previously untreated advanced oesophagogastric cancer (REAL3): a randomised, open-label phase 3 trial.

BACKGROUND: EGFR overexpression occurs in 27-55% of oesophagogastric adenocarcinomas, and correlates with poor prognosis. We aimed to assess addition of the anti-EGFR antibody panitumumab to epirubicin, oxaliplatin, and capecitabine (EOC) in patients with advanced oesophagogastric adenocarcinoma. METHODS: In this randomised, open-label phase 3 trial (REAL3), we enrolled patients with untreated, metastatic, or locally advanced oesophagogastric adenocarcinoma at 63 centres (tertiary referral centres, teaching hospitals, and district general hospitals) in the UK. Eligible patients were randomly allocated (1:1) to receive up to eight 21-day cycles of open-label EOC (epirubicin 50 mg/m(2) and oxaliplatin 130 mg/m(2) on day 1 and capecitabine 1250 mg/m(2) per day on days 1-21) or modified-dose EOC plus panitumumab (mEOC+P; epirubicin 50 mg/m(2) and oxaliplatin 100 mg/m(2) on day 1, capecitabine 1000 mg/m(2) per day on days 1-21, and panitumumab 9 mg/kg on day 1). Randomisation was blocked and stratified for centre region, extent of disease, and performance status. The primary endpoint was overall survival in the intention-to-treat population. We assessed safety in all patients who received at least one dose of study drug. After a preplanned independent data monitoring committee review in October, 2011, trial recruitment was halted and panitumumab withdrawn. Data for patients on treatment were censored at this timepoint. This study is registered with ClinicalTrials.gov, number NCT00824785. FINDINGS: Between June 2, 2008, and Oct 17, 2011, we enrolled 553 eligible patients. Median overall survival in 275 patients allocated EOC was 11.3 months (95% CI 9.6-13.0) compared with 8.8 months (7.7-9.8) in 278 patients allocated mEOC+P (hazard ratio [HR] 1.37, 95% CI 1.07-1.76; p=0.013). mEOC+P was associated with increased incidence of grade 3-4 diarrhoea (48 [17%] of 276 patients allocated mEOC+P vs 29 [11%] of 266 patients allocated EOC), rash (29 [11%] vs two [1%]), mucositis (14 [5%] vs none), and hypomagnesaemia (13 [5%] vs none) but reduced incidence of haematological toxicity (grade ≥ 3 neutropenia 35 [13%] vs 74 [28%]). INTERPRETATION: Addition of panitumumab to EOC chemotherapy does not increase overall survival and cannot be recommended for use in an unselected population with advanced oesophagogastric adenocarcinoma. FUNDING: Amgen, UK National Institute for Health Research Biomedical Research Centre.

XBP1s levels are implicated in the biology and outcome of myeloma mediating different clinical outcomes to thalidomide-based treatments.

Immunoglobulin production by myeloma plasma cells depends on the unfolded protein response for protein production and folding. Recent studies have highlighted the importance of IRE1alpha and X box binding protein 1 (XBP1), key members of this pathway, in normal B-plasma cell development. We have determined the gene expression levels of IRE1alpha, XBP1, XBP1UNSPLICED (XBP1u), and XBP1SPLICED (XBP1s) in a series of patients with myeloma and correlated findings with clinical outcome. We show that IRE1alpha and XBP1 are highly expressed and that patients with low XBP1s/u ratios have a significantly better overall survival. XBP1s is an independent prognostic marker and can be used with beta2 microglobulin and t(4;14) to identify a group of patients with a poor outcome. Furthermore, we show the beneficial therapeutic effects of thalidomide in patients with low XBP1s/u ratios. This study highlights the importance of XBP1 in myeloma and its significance as an independent prognostic marker and as a predictor of thalidomide response.

The effectiveness of a nurse practitioner-led pain management team in long-term care: A mixed methods study

BACKGROUND: Considering the high rates of pain as well as its under-management in long-term care (LTC) settings, research is needed to explore innovations in pain management that take into account limited resource realities. It has been suggested that nurse practitioners, working within an inter-professional model, could potentially address the under-management of pain in LTC.

OBJECTIVES: This study evaluated the effectiveness of implementing a nurse practitioner-led, inter-professional pain management team in LTC in improving (a) pain-related resident outcomes; (b) clinical practice behaviours (e.g., documentation of pain assessments, use of non-pharmacological and pharmacological interventions); and, (c) quality of pain medication prescribing practices.

METHODS: A mixed method design was used to evaluate a nurse practitioner-led pain management team, including both a quantitative and qualitative component. Using a controlled before-after study, six LTC homes were allocated to one of three groups: 1) a nurse practitioner-led pain team (full intervention); 2) nurse practitioner but no pain management team (partial intervention); or, 3) no nurse practitioner, no pain management team (control group). In total, 345 LTC residents were recruited to participate in the study; 139 residents for the full intervention group, 108 for the partial intervention group, and 98 residents for the control group. Data was collected in Canada from 2010 to 2012.

RESULTS: Implementing a nurse practitioner-led pain team in LTC significantly reduced residents' pain and improved functional status compared to usual care without access to a nurse practitioner. Positive changes in clinical practice behaviours (e.g., assessing pain, developing care plans related to pain management, documenting effectiveness of pain interventions) occurred over the intervention period for both the nurse practitioner-led pain team and nurse practitioner-only groups; these changes did not occur to the same extent, if at all, in the control group. Qualitative analysis highlighted the perceived benefits of LTC staff about having access to a nurse practitioner and benefits of the pain team, along with barriers to managing pain in LTC.

CONCLUSIONS: The findings from this study showed that implementing a nurse practitioner-led pain team can significantly improve resident pain and functional status as well as clinical practice behaviours of LTC staff. LTC homes should employ a nurse practitioner, ideally located onsite as opposed to an offsite consultative role, to enhance inter-professional collaboration and facilitate more consistent and timely access to pain management.

Crack propagation in non-homogenous materials: Evaluation of mixed-mode SIFs, T-stress and kinking angle using a variant EFG method

A new variant of the Element-Free Galerkin (EFG) method, that combines the diffraction method, to characterize the crack tip solution, and the Heaviside enrichment function for representing discontinuity due to a crack, has been used to model crack propagation through non-homogenous materials. In the case of interface crack propagation, the kink angle is predicted by applying the maximum tangential principal stress (MTPS) criterion in conjunction with consideration of the energy release rate (ERR). The MTPS criterion is applied to the crack tip stress field described by both the stress intensity factor (SIF) and the T-stress, which are extracted using the interaction integral method. The proposed EFG method has been developed and applied for 2D case studies involving a crack in an orthotropic material, crack along an interface and a crack terminating at a bi-material interface, under mechanical or thermal loading; this is done to demonstrate the advantages and efficiency of the proposed methodology. The computed SIFs, T-stress and the predicted interface crack kink angles are compared with existing results in the literature and are found to be in good agreement. An example of crack growth through a particle-reinforced composite materials, which may involve crack meandering around the particle, is reported.

Bayesian Dependence Tests for Continuous, Binary and Mixed Continuous-Binary Variables

Tests for dependence of continuous, discrete and mixed continuous-discrete variables are ubiquitous in science. The goal of this paper is to derive Bayesian alternatives to frequentist null hypothesis significance tests for dependence. In particular, we will present three Bayesian tests for dependence of binary, continuous and mixed variables. These tests are nonparametric and based on the Dirichlet Process, which allows us to use the same prior model for all of them. Therefore, the tests are “consistent” among each other, in the sense that the probabilities that variables are dependent computed with these tests are commensurable across the different types of variables being tested. By means of simulations with artificial data, we show the effectiveness of the new tests.

Investigation of the effects of flow conditions at rotor inlet on mixed flow turbine performance for automotive applications

Current trends in the automotive industry have placed increased importance on engine downsizing for passenger vehicles. Engine downsizing often results in reduced power output and turbochargers have been relied upon to restore the power output and maintain drivability. As improved power output is required across a wide range of engine operating conditions, it is necessary for the turbocharger to operate effectively at both design and off-design conditions. One off-design condition of considerable importance for turbocharger turbines is low velocity ratio operation, which refers to the combination of high exhaust gas velocity and low turbine rotational speed. Conventional radial flow turbines are constrained to achieve peak efficiency at the relatively high velocity ratio of 0.7, due the requirement to maintain a zero inlet blade angle for structural reasons. Several methods exist to potentially shift turbine peak efficiency to lower velocity ratios. One method is to utilize a mixed flow turbine as an alternative to a radial flow turbine. In addition to radial and circumferential components, the flow entering a mixed flow turbine also has an axial component. This allows the flow to experience a non-zero inlet blade angle, potentially shifting peak efficiency to a lower velocity ratio when compared to an equivalent radial flow turbine.
This study examined the effects of varying the flow conditions at the inlet to a mixed flow turbine and evaluated the subsequent impact on performance. The primary parameters examined were average inlet flow angle, the spanwise distribution of flow angle across the inlet and inlet flow cone angle. The results have indicated that the inlet flow angle significantly influenced the degree of reaction across the rotor and the turbine efficiency. The rotor studied was a custom in-house design based on a state-of-the-art radial flow turbine design. A numerical approach was used as the basis for this investigation and the numerical model has been validated against experimental data obtained from the cold flow turbine test rig at Queen’s University Belfast. The results of the study have provided a useful insight into how the flow conditions at rotor inlet influence the performance of a mixed flow turbine.