Association of osteopontin with kidney function and kidney failure in chronic kidney disease patients: the GCKD study.

BACKGROUND: Osteopontin (OPN), synthesized in the thick ascending limb of Henle's loop and in the distal tubule, is involved in the pathogenesis of kidney fibrosis, a hallmark of kidney failure (KF). In a cohort of chronic kidney disease (CKD) patients, we evaluated OPN's association with kidney markers and KF. METHODS: OPN was measured from baseline serum samples of German Chronic Kidney Disease study participants. Cross-sectional regression models for estimated glomerular filtration rate (eGFR) and urinary albumin-to-creatinine ratio (UACR) as well as Cox regression models for all-cause mortality and KF were evaluated to estimate the OPN effect. Additionally, the predictive ability of OPN and time-dependent population-attributable fraction were evaluated. RESULTS: Over a median follow-up of 6.5 years, 471 KF events and 629 deaths occurred among 4950 CKD patients. One-unit higher log(OPN) was associated with 5.5 mL/min/1.73 m2 lower eGFR [95% confidence interval (95% CI) -6.4 to -4.6] and 1% change in OPN with 0.7% higher UACR (estimated effect 0.7, 95% CI 0.6-0.8). Moreover, higher OPN levels were associated with a higher risk of KF [hazard ratio (HR) 1.4, 95% CI 1.2-1.7] and all-cause mortality (HR 1.5, 95% CI 1.3-1.8). After 6 years, 31% of the KF events could be attributed to higher OPN levels (95% CI 3%-56%). CONCLUSIONS: In this study, higher OPN levels were associated with kidney function markers worsening and a higher risk for adverse outcomes. A larger proportion of KF could be attributed to higher OPN levels, warranting further research on OPN with regards to its role in CKD progression and possible treatment options.

Translation and Psychometric Evaluation of the German Version of the Thirst Distress Scale for Patients With Heart Failure.

BACKGROUND: In patients with chronic heart failure, thirst can be perceived as an intensive and burdensome symptom, which may have a negative impact on patients' quality of life. To initiate thirst-relieving interventions, assessment of thirst and its related distress is essential. At the time of this study, no instrument was available to evaluate thirst distress in patients with heart failure in Germany. OBJECTIVE: The aims of this study were to translate the "Thirst Distress Scale for patients with Heart Failure" (TDS-HF) from English into German and to test validity and reliability of the scale. METHODS: The English version of the TDS-HF was translated into German. A linguistically and culturally sensitive forward-and-backward translation was performed. Psychometric evaluation included confirmatory factor analysis, reliability in terms of internal consistency, and concurrent validity. RESULTS: Eighty-four hospitalized patients (mean age, 72 ± 10 years; 29% female; mean left ventricular ejection fraction, 36% ± 12%; 62% New York Heart Association functional classes III-IV, 45% on fluid restriction) from an acute care hospital were involved in the study. The item-total correlation ranged from 0.58 to 0.78. Interitem correlations varied between 0.37 and 0.79. Internal consistency was high, with a Cronbach α of 0.89. There was a high correlation between the total score of the TDS-HF and the visual analog scale to assess thirst intensity ( r = 0.72, P ≤ .001), and a low correlation with fluid restriction ( r = 0.35, P = .002). CONCLUSIONS: The evaluation of the German TDS-HF showed satisfactory psychometric properties in this sample. The instrument is usable for further research and additional psychometric testing.

Harmonizing Heterogeneous Endpoints in Coronavirus Disease 2019 Trials Without Loss of Information.

OBJECTIVES: Many trials investigate potential effects of treatments for coronavirus disease 2019. To provide sufficient information for all involveddecision-makers (clinicians, public health authorities, and drug regulatory agencies), a multiplicity of endpoints must be considered. The objectives are to provide hands-on statistical guidelines for harmonizing heterogeneous endpoints in coronavirus disease 2019 clinical trials. DESIGN: Randomized controlled trials for patients infected with coronavirus disease 2019. SETTING: General methods that apply to any randomized controlled trial for patients infected with coronavirus disease 2019. PATIENTS: Coronavirus disease 2019 positive individuals. INTERVENTIONS: None. MEASUREMENTS AND MAIN RESULTS: We develop a multistate model that is based on hospitalization, mechanical ventilation, death, and discharge. These events are both categories of the ordinal endpoint recommended by the World Health Organization and also within the core outcome set of the Core Outcome Measures in Effectiveness Trials initiative for coronavirus disease 2019 trials. To support our choice of states in the multistate model, we also perform a brief review of registered coronavirus disease 2019 clinical trials. Based on the multistate model, we give recommendation for compact, informative illustration of time-dynamic treatment effects and explorative statistical analysis. A majority of coronavirus disease 2019 clinical trials collect information on mechanical ventilation, hospitalization, and death. Using reconstructed and real data of coronavirus disease 2019 trials, we show how a stacked probability plot provides a detailed understanding of treatment effects on the patients' course of hospital stay. It contributes to harmonizing multiple endpoints and differing lengths of follow-up both within and between trials. CONCLUSIONS: All ongoing clinical trials should include a stacked probability plot in their statistical analysis plan as descriptive analysis. While primary analysis should be on an early endpoint with appropriate capability to be a surrogate (parameter), our multistate model provides additional detailed descriptive information and links results within and between coronavirus disease 2019 trials.

Predicting Potential Prevention Effects on Hospital Burden of Nosocomial Infections: A Multistate Modeling Approach.

OBJECTIVES: Hospital-acquired infections (HAIs) place a substantial burden on health systems. Tools are required to quantify the change in this burden as a result of a preventive intervention. We aim to estimate how much a reduction in the rate of hospital-acquired infections translates into a change in hospital mortality and length of stay. METHODS: Using multistate modelling and competing risks methodology, we created a tool to estimate the reduction in burden after the introduction of a preventive effect on the infection rate. The tool requires as inputs the patients' length of hospital stay, patients' infection information (status, time), patients' final outcome (discharged alive, dead), and a preventive effect. We demonstrated the methods on both simulated data and 3 published data sets from Germany, France, and Spain. RESULTS: A hypothetical prevention that cuts the infection rate in half would result in 21 lives and 2212 patient-days saved in French ventilator-associated pneumonia data, 61 lives and 3125 patient-days saved in Spanish nosocomial infection data, and 20 lives and 1585 patient-days saved in German nosocomial pneumonia data. CONCLUSIONS: Our tool provides a quick and easy means of acquiring an impression of the impact a preventive measure would have on the burden of an infection. The tool requires quantities routinely collected and computation can be done with a calculator. R code is provided for researchers to determine the burden in various settings with various effects. Furthermore, cost data can be used to get the financial benefit of the reduction in burden.

Quantifying and communicating the burden of COVID-19.

BACKGROUND: An essential aspect of preventing further COVID-19 outbreaks and to learn for future pandemics is the evaluation of different political strategies, which aim at reducing transmission of and mortality due to COVID-19. One important aspect in this context is the comparison of attributable mortality. METHODS: We give a comprehensive overview of six epidemiological measures that are used to quantify COVID-19 attributable mortality (p-score, standardized mortality ratio, absolute number of excess deaths, per capita rate, z-score and the population attributable fraction). RESULTS: By defining the six measures based on observed and expected deaths, we explain their relationship. Moreover, three publicly available data examples serve to illustrate the interpretational strengths and weaknesses of the various measures. Finally, we give recommendation which measures are suitable for an evaluation of public health strategies against COVID-19. The R code to reproduce the results is available as online supplementary material. CONCLUSION: The number of excess deaths should be always reported together with the population attributable fraction, the p-score or the standardized mortality ratio instead of a per capita rate. For a complete picture of COVID-19 attributable mortality, quantifying and communicating its relative burden also to a lay audience is of major importance.

Absolute mortality risk assessment of COVID-19 patients: the Khorshid COVID Cohort (KCC) study.

BACKGROUND: Already at hospital admission, clinicians require simple tools to identify hospitalized COVID-19 patients at high risk of mortality. Such tools can significantly improve resource allocation and patient management within hospitals. From the statistical point of view, extended time-to-event models are required to account for competing risks (discharge from hospital) and censoring so that active cases can also contribute to the analysis. METHODS: We used the hospital-based open Khorshid COVID Cohort (KCC) study with 630 COVID-19 patients from Isfahan, Iran. Competing risk methods are used to develop a death risk chart based on the following variables, which can simply be measured at hospital admission: sex, age, hypertension, oxygen saturation, and Charlson Comorbidity Index. The area under the receiver operator curve was used to assess accuracy concerning discrimination between patients discharged alive and dead. RESULTS: Cause-specific hazard regression models show that these baseline variables are associated with both death, and discharge hazards. The risk chart reflects the combined results of the two cause-specific hazard regression models. The proposed risk assessment method had a very good accuracy (AUC = 0.872 [CI 95%: 0.835-0.910]). CONCLUSIONS: This study aims to improve and validate a personalized mortality risk calculator based on hospitalized COVID-19 patients. The risk assessment of patient mortality provides physicians with additional guidance for making tough decisions.

Joint analysis of duration of ventilation, length of intensive care, and mortality of COVID-19 patients: a multistate approach.

BACKGROUND: The clinical progress of patients hospitalized due to COVID-19 is often associated with severe pneumonia which may require intensive care, invasive ventilation, or extracorporeal membrane oxygenation (ECMO). The length of intensive care and the duration of these supportive therapies are clinically relevant outcomes. From the statistical perspective, these quantities are challenging to estimate due to episodes being time-dependent and potentially multiple, as well as being determined by the competing, terminal events of discharge alive and death. METHODS: We used multistate models to study COVID-19 patients' time-dependent progress and provide a statistical framework to estimate hazard rates and transition probabilities. These estimates can then be used to quantify average sojourn times of clinically important states such as intensive care and invasive ventilation. We have made two real data sets of COVID-19 patients (n = 24* and n = 53**) and the corresponding statistical code publically available. RESULTS: The expected lengths of intensive care unit (ICU) stay at day 28 for the two cohorts were 15.05* and 19.62** days, while expected durations of mechanical ventilation were 7.97* and 9.85** days. Predicted mortality stood at 51%* and 15%**. Patients mechanically ventilated at the start of the example studies had a longer expected duration of ventilation (12.25*, 14.57** days) compared to patients non-ventilated (4.34*, 1.41** days) after 28 days. Furthermore, initially ventilated patients had a higher risk of death (54%* and 20%** vs. 48%* and 6%**) after 4 weeks. These results are further illustrated in stacked probability plots for the two groups from time zero, as well as for the entire cohort which depicts the predicted proportions of the patients in each state over follow-up. CONCLUSIONS: The multistate approach gives important insights into the progress of COVID-19 patients in terms of ventilation duration, length of ICU stay, and mortality. In addition to avoiding frequent pitfalls in survival analysis, the methodology enables active cases to be analyzed by allowing for censoring. The stacked probability plots provide extensive information in a concise manner that can be easily conveyed to decision makers regarding healthcare capacities. Furthermore, clear comparisons can be made among different baseline characteristics.

The population-attributable fraction for time-dependent exposures using dynamic prediction and landmarking.

The public health impact of a harmful exposure can be quantified by the population-attributable fraction (PAF). The PAF describes the attributable risk due to an exposure and is often interpreted as the proportion of preventable cases if the exposure was extinct. Difficulties in the definition and interpretation of the PAF arise when the exposure of interest depends on time. Then, the definition of exposed and unexposed individuals is not straightforward. We propose dynamic prediction and landmarking to define and estimate a PAF in this data situation. Two estimands are discussed which are based on two hypothetical interventions that could prevent the exposure in different ways. Considering the first estimand, at each landmark the estimation problem is reduced to a time-independent setting. Then, estimation is simply performed by using a generalized-linear model accounting for the current exposure state and further (time-varying) covariates. The second estimand is based on counterfactual outcomes, estimation can be performed using pseudo-values or inverse-probability weights. The approach is explored in a simulation study and applied on two data examples. First, we study a large French database of intensive care unit patients to estimate the population-benefit of a pathogen-specific intervention that could prevent ventilator-associated pneumonia caused by the pathogen Pseudomonas aeruginosa. Moreover, we quantify the population-attributable burden of locoregional and distant recurrence in breast cancer patients.

Investigating the Impact of Early Valve Surgery on Survival in Staphylococcus aureus Infective Endocarditis Using a Marginal Structural Model Approach: Results of a Large, Prospectively Evaluated Cohort.

BACKGROUND: The impact of valve surgery on outcomes of Staphylococcus aureus infective endocarditis (SAIE) remains controversial. We tested the hypothesis that early valve surgery (EVS) improves survival by using a novel approach that allows for inclusion of major confounders in a time-dependent way. METHODS: EVS was defined as valve surgery within 60 days. Univariable and multivariable Cox regression analyses were performed. To account for treatment selection bias, we additionally used a weighted Cox model (marginal structural model) that accounts for time-dynamic imbalances between treatment groups. To address survivor bias, EVS was included as a time-dependent variable. Follow-up of patients was 1 year. RESULTS: Two hundred and three patients were included in the analysis; 50 underwent EVS. All-cause mortality at day 30 was 26%. In the conventional multivariable Cox regression model, the effect of EVS on the death hazard was 0.85 (95% confidence interval [CI], .47-1.52). Using the weighted Cox model, the death hazard rate (HR) of EVS was 0.71 (95% CI, .34-1.49). In subgroup analyses, no survival benefit was observed in patients with septic shock (HR, 0.80 [CI, .26-2.46]), in NVIE (HR, 0.76 [CI, .33-1.71]) or PVIE (HR, 1.02 [CI, .29-3.54]), or in patients with EVS within 14 days (HR, 0.97 [CI, .46-2.07]). CONCLUSIONS: Using both a conventional Cox regression model and a weighted Cox model, we did not find a survival benefit for patients who underwent EVS in our cohort. Until results of randomized controlled trials are available, EVS in SAIE should be based on individualized decisions of an experienced multidisciplinary team. CLINICAL TRIALS REGISTRATION: German Clinical Trials registry (DRKS00005045).

The population-attributable fraction for time-dependent exposures and competing risks-A discussion on estimands.

The population-attributable fraction (PAF) quantifies the public health impact of a harmful exposure. Despite being a measure of significant importance, an estimand accommodating complicated time-to-event data is not clearly defined. We discuss current estimands of the PAF used to quantify the public health impact of an internal time-dependent exposure for data subject to competing outcomes. To overcome some limitations, we proposed a novel estimand that is based on dynamic prediction by landmarking. In a profound simulation study, we discuss interpretation and performance of the various estimands and their estimators. The methods are applied to a large French database to estimate the health impact of ventilator-associated pneumonia for patients in intensive care.

The Impact of Early Adequate Treatment on Extubation and Discharge Alive of Patients With Pseudomonas aeruginosa-Related Ventilator-Associated Pneumonia.

OBJECTIVES: We aim to examine the effect of early adequate treatment in comparison with inadequate or delayed treatment on being extubated or discharged alive over time, in patients with Pseudomonas aeruginosa-related ventilator-associated pneumonia. DESIGN: Retrospective analyses of a prospective observational multicenter cohort study. SETTING: ICU. PATIENTS: Patients of the French prospective database (OUTCOMEREA) were included if they acquired a ventilator-associated pneumonia due to P. aeruginosa between 1997 and 2014 and were mechanically ventilated for more than 48 hours. INTERVENTIONS: Early adequate treatment in comparison with inadequate or delayed adequate treatment. MEASUREMENTS AND MAIN RESULTS: Multistate models were applied to estimate the time-dependent probability of being extubated or discharged alive, and separate Cox regression analyses were used to assess the treatment effect on all important events that influence the outcome of interest. A propensity score-adjusted innovative regression technique was used for a combined and comprehensive patient-relevant summary effect measure. No evidence was found for a difference between adequate and inadequate or delayed treatment on being extubated or discharged alive. However, for all patients, the probability of being extubated or discharged alive remains low and does not exceed 50% even 40 days after a P. aeruginosa-related ventilator-associated pneumonia. CONCLUSIONS: Early adequate treatment does not seem to be associated with an improved prognosis. Its potential benefit requires further investigation in larger observational studies.

Bias due to censoring of deaths when calculating extra length of stay for patients acquiring a hospital infection.

BACKGROUND: In many studies the information of patients who are dying in the hospital is censored when examining the change in length of hospital stay (cLOS) due to hospital-acquired infections (HIs). While appropriate estimators of cLOS are available in literature, the existence of the bias due to censoring of deaths was neither mentioned nor discussed by the according authors. METHODS: Using multi-state models, we systematically evaluate the bias when estimating cLOS in such a way. We first evaluate the bias in a mathematically closed form assuming a setting with constant hazards. To estimate the cLOS due to HIs non-parametrically, we relax the assumption of constant hazards and consider a time-inhomogeneous Markov model. RESULTS: In our analytical evaluation we are able to discuss challenging effects of the bias on cLOS. These are in regard to direct and indirect differential mortality. Moreover, we can make statements about the magnitude and direction of the bias. For real-world relevance, we illustrate the bias on a publicly available prospective cohort study on hospital-acquired pneumonia in intensive-care. CONCLUSION: Based on our findings, we can conclude that censoring the death cases in the hospital and considering only patients discharged alive should be avoided when estimating cLOS. Moreover, we found that the closed mathematical form can be used to describe the bias for settings with constant hazards.

A case-cohort approach for multi-state models in hospital epidemiology.

Analysing the determinants and consequences of hospital-acquired infections involves the evaluation of large cohorts. Infected patients in the cohort are often rare for specific pathogens, because most of the patients admitted to the hospital are discharged or die without such an infection. Death and discharge are competing events to acquiring an infection, because these individuals are no longer at risk of getting a hospital-acquired infection. Therefore, the data is best analysed with an extended survival model - the extended illness-death model. A common problem in cohort studies is the costly collection of covariate values. In order to provide efficient use of data from infected as well as uninfected patients, we propose a tailored case-cohort approach for the extended illness-death model. The basic idea of the case-cohort design is to only use a random sample of the full cohort, referred to as subcohort, and all cases, namely the infected patients. Thus, covariate values are only obtained for a small part of the full cohort. The method is based on existing and established methods and is used to perform regression analysis in adapted Cox proportional hazards models. We propose estimation of all cause-specific cumulative hazards and transition probabilities in an extended illness-death model based on case-cohort sampling. As an example, we apply the methodology to infection with a specific pathogen using a large cohort from Spanish hospital data. The obtained results of the case-cohort design are compared with the results in the full cohort to investigate the performance of the proposed method. Copyright © 2016 John Wiley & Sons, Ltd.

Basic parametric analysis for a multi-state model in hospital epidemiology.

BACKGROUND: The extended illness-death model is a useful tool to study the risks and consequences of hospital-acquired infections (HAIs). The statistical quantities of interest are the transition-specific hazard rates and the transition probabilities as well as attributable mortality (AM) and the population-attributable fraction (PAF). In the most general case calculation of these expressions is mathematically complex. METHODS: When assuming time-constant hazards calculation of the quantities of interest is facilitated. In this situation the transition probabilities can be expressed in closed mathematical forms. The estimators for AM and PAF can be easily derived from these forms. RESULTS: In this paper, we show how to explicitly calculate all the transition probabilities of an extended-illness model with constant hazards. Using a parametric model to estimate the time-constant transition specific hazard rates of a data example, the transition probabilities, AM and PAF can be directly calculated. With a publicly available data example, we show how the approach provides first insights into principle time-dynamics and data structure. CONCLUSION: Assuming constant hazards facilitates the understanding of multi-state processes. Even in a non-constant hazards setting, the approach is a helpful first step for a comprehensive investigation of complex data.

The population-attributable fraction for time-to-event data.

BACKGROUND: Even though the population-attributable fraction (PAF) is a well-established metric, it is often incorrectly estimated or interpreted not only in clinical application, but also in statistical research articles. The risk of bias is especially high in more complex time-to-event data settings. METHODS: We explain how the PAF can be defined, identified and estimated in time-to-event settings with competing risks and time-dependent exposures. By using multi-state methodology and inverse probability weighting, we demonstrate how to reduce or completely avoid severe types of biases including competing risks bias, immortal time bias and confounding due to both baseline and time-varying patient characteristics. RESULTS: The method is exemplarily applied to a real data set. Moreover, we estimate the number of deaths that were attributable to ventilator-associated pneumonia in France in the year 2016. The example demonstrates how, under certain simplifying assumptions, PAF estimates can be extrapolated to a target population of interest. CONCLUSIONS: Defining and estimating the PAF in advanced time-to-event settings within a framework that unifies causal and multi-state modelling enables to tackle common sources of bias and allows straightforward implementation with standard software packages.

Impact of adherence to individual quality-of-care indicators on the prognosis of bloodstream infection due to Staphylococcus aureus: a prospective observational multicentre cohort.

OBJECTIVES: To analyse the adherence and impact of quality-of-care indicators (QCIs) in the management of Staphylococcus aureus bloodstream infection in a prospective and multicentre cohort. METHODS: Analysis of the prospective, multicentre international S. Aureus Collaboration cohort of S. Aureus bloodstream infection cases observed between January 2013 and April 2015. Multivariable analysis was performed to evaluate the impact of adherence to QCIs on 90-day mortality. RESULTS: A total of 1784 cases were included. Overall, 90-day mortality was 29.9% and mean follow-up period was 118 days. Adherence was 67% (n = 1180/1762) for follow-up blood cultures, 31% (n = 416/1342) for early focus control, 77.6% (n = 546/704) for performance of echocardiography, 75.5% (n = 1348/1784) for adequacy of targeted antimicrobial therapy, 88.6% (n = 851/960) for adequacy of treatment duration in non-complicated bloodstream infections and 61.2% (n = 366/598) in complicated bloodstream infections. Full bundle adherence was 18.4% (n = 328/1784). After controlling for immortal time bias and potential confounders, focus control (adjusted hazard ratio = 0.76; 95% CI, 0.59-0.99; p 0.038) and adequate targeted antimicrobial therapy (adjusted hazard ratio = 0.75; 95% CI, 0.61-0.91; p 0.004) were associated with low 90-day mortality. DISCUSSION: Adherence to QCIs in S. Aureus bloodstream infection did not reach expected rates. Apart from the benefits of application as a bundle, focus control and adequate targeted therapy were independently associated with low mortality.