Evaluation of the proportional hazards assumption and covariate adjustment methods in comparative surgical observational studies with time-to-event endpoints.

INTRODUCTION: Comparative studies on surgical treatments with time-to-event endpoints have provided substantial evidence for clinical practice, but the accurate use of survival data analysis and the control of confounding bias remain big challenges. METHODS: This was a survey of surgical studies with survival outcomes published in four general medical journals and five general surgical journals in 2021. The two most concerned statistical issues were evaluated, including confounding control by propensity score analysis (PSA) or multivariable analysis and testing of proportional hazards (PH) assumption in Cox model. RESULTS: A total of 74 studies were included, comprising 63 observational studies and 11 randomized controlled trials. Among the observational studies, the proportion of studies utilizing PSA in surgical oncology and non-oncology studies was similar (40.9 % versus 36.8 %, P = 0.762). However, the former reported a significantly lower proportion of PH assumption assessments compared to the latter (13.6 % versus 42.1 %, P = 0.020). Twenty-five observational studies (25/63) used PSA methods, but two-thirds of them (17/25) showed unclear balance of baseline data after PSA. And the proportion of PH assumption testing after PSA was slightly lower than that before PSA, but the difference was not statistically significant (24.0 % versus 28.0 %, P = 0.317). Comprehensive suggestions were given on confounding control in survival analysis and alternative resolutions for non-compliance with PH assumption. CONCLUSION: This study highlights suboptimal reporting of PH assumption evaluation in observational surgical studies both before and after PSA. Efforts and consensus are needed with respect to the underlying assumptions of statistical methods.

Statistical Analysis of Self-Reported Health Conditions in Cohort Studies: Handling of Missing Onset Age.

OBJECTIVE: This paper discusses methodological challenges in epidemiological association analysis of a time-to-event outcome and hypothesized risk factors, where age/time at the onset of the outcome may be missing in some cases, a condition commonly encountered when the outcome is self-reported. STUDY DESIGN AND SETTING: A cohort study with long-term follow-up for outcome ascer- tainment such as the Childhood Cancer Survivor Study (CCSS), a large cohort study of 5-year survivors of childhood cancer diagnosed in 1970-1999 in which occurrences and age at onset of various chronic health conditions (CHCs) are self-reported in surveys. Simple methods for handling missing onset age and their potential bias in the exposure-outcome association infer- ence are discussed. The interval-censored method is discussed as a remedy for handling this problem. The finite sample performance of these approaches is compared through Monte Carlo simulations. Examples from the CCSS include four CHCs (diabetes, myocardial infarction, osteoporosis/osteopenia, and growth hormone deficiency). RESULTS: The interval-censored method is usable in practice using the standard statisti- cal software. The simulation study showed that the regression coefficient estimates from the 'Interval censored' method consistently displayed reduced bias and, in most cases, smaller stan- dard deviations, resulting in smaller mean square errors, compared to those from the simple approaches, regardless of the proportion of subjects with an event of interest, the proportion of missing onset age, and the sample size. CONCLUSION: The interval-censored method is a statistically valid and practical approach to the association analysis of self-reported time-to-event data when onset age may be missing. While the simpler approaches that force such data into complete data may enable the standard analytic methods to be applicable, there is considerable loss in both accuracy and precision relative to the interval-censored method.

A Shared-Frailty Spatial Scan Statistic Model for Time-to-Event Data.

Spatial scan statistics are well-known methods widely used to detect spatial clusters of events. Furthermore, several spatial scan statistics models have been applied to the spatial analysis of time-to-event data. However, these models do not take account of potential correlations between the observations of individuals within the same spatial unit or potential spatial dependence between spatial units. To overcome this problem, we have developed a scan statistic based on a Cox model with shared frailty and that takes account of the spatial dependence between spatial units. In simulation studies, we found that (i) conventional models of spatial scan statistics for time-to-event data fail to maintain the type I error in the presence of a correlation between the observations of individuals within the same spatial unit and (ii) our model performed well in the presence of such correlation and spatial dependence. We have applied our method to epidemiological data and the detection of spatial clusters of mortality in patients with end-stage renal disease in northern France.

Competing Risks Analysis of Kidney Transplant Waitlist Outcomes: Two Important Statistical Perspectives.

Modern competing risks analysis has 2 primary goals in clinical epidemiology as follows: (i) to maximize the clinician's knowledge of etiologic associations existing between potential predictor variables and various cause-specific outcomes via cause-specific hazard models, and (ii) to maximize the clinician's knowledge of noteworthy differences existing in cause-specific patient risk via cause-specific subdistribution hazard models (cumulative incidence functions [CIFs]). A perfect application exists in analyzing the following 4 distinct outcomes after listing for a deceased donor kidney transplant (DDKT): (i) receiving a DDKT, (ii) receiving a living donor kidney transplant (LDKT), (iii) waitlist removal due to patient mortality or a deteriorating medical condition, and (iv) waitlist removal due to other reasons. It is important to realize that obtaining a complete understanding of subdistribution hazard ratios (HRs) is simply not possible without first having knowledge of the multivariable relationships existing between the potential predictor variables and the cause-specific hazards (perspective #1), because the cause-specific hazards form the "building blocks" of CIFs. In addition, though we believe that a worthy and practical alternative to estimating the median waiting-time-to DDKT is to ask, "what is the conditional probability of the patient receiving a DDKT, given that he or she would not previously experience one of the competing events (known as the cause-specific conditional failure probability)," only an appropriate estimator of this conditional type of cumulative incidence should be used (perspective #2). One suggested estimator, the well-known "one minus Kaplan-Meier" approach (censoring competing events), simply does not represent any probability in the presence of competing risks and will almost always produce biased estimates (thus, it should never be used).

Modeling Recovery Housing Retention and Program Outcomes by Justice Involvement among Residents in Virginia, USA: An Observational Study.

Living in recovery housing can improve addiction recovery and desistance outcomes. This study examined whether retention in recovery housing and types of discharge outcomes (completed, "neutral," and "negative" outcomes) differed for clients with recent criminal legal system (CLS) involvement. Using data from 101 recovery residences certified by the Virginia Association of Recovery Residences based on 1,978 individuals completing the REC-CAP assessment, competing risk analyses (cumulative incidence function, restricted mean survival time, and restricted mean time lost) followed by the marginalization of effects were implemented to examine program outcomes at final discharge. Residents with recent CLS involvement were more likely to be discharged for positive reasons (successful completion of their goals) and premature/negative reasons (e.g., disciplinary releases) than for neutral reasons. Findings indicate that retention for 6-18 months is essential to establish and maintain positive discharge outcomes, and interventions should be developed to enhance retention in recovery residents with recent justice involvement.

Time to prerandomization seizure count design sufficiently assessed the safety and tolerability of perampanel for the treatment of primary generalized tonic-clonic seizures.

OBJECTIVE: Static assignment of participants in randomized clinical trials to placebo or ineffective treatment confers risk from continued seizures. An alternative trial design of time to exceed prerandomization monthly seizure count (T-PSC) has replicated the efficacy conclusions of traditionally designed trials, with shorter exposure to placebo and ineffective treatment. Trials aim to evaluate efficacy as well as safety and tolerability; therefore, we evaluated whether this T-PSC design also could replicate the trial's safety and tolerability conclusions. METHODS: We retrospectively applied the T-PSC design to analyze treatment-emergent adverse events (TEAEs) from a blinded, placebo-controlled trial of perampanel for primary generalized tonic-clonic seizures (NCT01393743). The safety analysis set consisted of 81 and 82 participants randomized to perampanel and placebo arms, respectively. We evaluated the incidences of TEAEs, treatment-related TEAEs, serious TEAEs, and TEAEs of special interest that occurred before T-PSC relative to those observed during the full-length trial. RESULTS: Of the 67 and 59 participants who experienced TEAEs in the perampanel and placebo arms during full-length trial, 66 (99%) and 54 (92%) participants experienced TEAEs with onset occurring before T-PSC, respectively. When limited to treatment-related TEAEs, 55 of 56 (98%) and 32 of 37 (86%) participants reported treatment-related TEAEs that occurred before T-PSC in the perampanel and placebo arms, respectively. There were more TEAEs after T-PSC with placebo as compared to perampanel (Fisher exact odds ratio = 8.6, p = .035), which resulted in overestimation of the difference in TEAE rate. There was a numerical reduction in serious TEAEs (3/13 occurred after T-PSC, one in placebo and two in perampanel). SIGNIFICANCE: Almost all TEAEs occurred before T-PSC. More treatment-related TEAEs occurred after T-PSC for participants randomized to placebo than perampanel, which may be due to either a shorter T-PSC or delayed time to TEAE for placebo.

Bridging the gap between two-stage and joint models: The case of tumor growth inhibition and overall survival models.

Many clinical trials generate both longitudinal biomarker and time-to-event data. We might be interested in their relationship, as in the case of tumor size and overall survival in oncology drug development. Many well-established methods exist for analyzing such data either sequentially (two-stage models) or simultaneously (joint models). Two-stage modeling (2stgM) has been challenged (i) for not acknowledging that biomarkers are endogenous covariable to the survival submodel and (ii) for not propagating the uncertainty of the longitudinal biomarker submodel to the survival submodel. On the other hand, joint modeling (JM), which properly circumvents both problems, has been criticized for being time-consuming, and difficult to use in practice. In this paper, we explore a third approach, referred to as a novel two-stage modeling (N2stgM). This strategy reduces the model complexity without compromising the parameter estimate accuracy. The three approaches (2stgM, JM, and N2stgM) are formulated, and a Bayesian framework is considered for their implementation. Both real and simulated data were used to analyze the performance of such approaches. In all scenarios, our proposal estimated the parameters approximately as JM but without being computationally expensive, while 2stgM produced biased results.

Lost in the plot: missing visual elements in Kaplan-Meier plots of phase III oncology trials.

Missing visual elements (MVE) in Kaplan-Meier (KM) curves can misrepresent data, preclude curve reconstruction, and hamper transparency. This study evaluated KM plots of phase III oncology trials. MVE were defined as an incomplete y-axis range or missing number at risk table in a KM curve. Surrogate endpoint KM curves were additionally evaluated for complete interpretability, defined by (1) reporting the number of censored patients and (2) correspondence of the disease assessment interval with the number at risk interval. Among 641 trials enrolling 518 235 patients, 116 trials (18%) had MVE in KM curves. Industry sponsorship, larger trials, and more recently published trials were correlated with lower odds of MVE. Only 3% of trials (15 of 574) published surrogate endpoint KM plots with complete interpretability. Improvements in the quality of KM curves of phase III oncology trials, particularly for surrogate endpoints, are needed for greater interpretability, reproducibility, and transparency in oncology research.

Genetic and early life factors influence on time-to-multiple sclerosis diagnosis: A UK Biobank study.

BACKGROUND: Previous investigations into multiple sclerosis (MS) risk factors predominantly relied on retrospective studies, which do not consider different follow-up times and assume a constant risk effect throughout lifetime. OBJECTIVE: We aimed to evaluate the impact of genetic and early life factors on MS diagnosis by employing a time-to-event analysis in a prospective cohort. METHODS: We used the UK Biobank data, considering the observation period from birth up to 31 December 2022. We considered genetic risk, using a multiple sclerosis polygenic risk score (MS-PRS), and various early life factors. Tobacco smoking and infectious mononucleosis diagnosis were also considered as time-varying variables along the follow-up. Using a Cox proportional hazards model, we examined the associations between these factors and MS diagnosis instantaneous risk. RESULTS: We analyzed 345,027 participants, of which 1669 had an MS diagnosis. Our analysis revealed age-dependent effects for sex (females vs males) and higher MS-PRS, with greater hazard ratios observed in young adults. CONCLUSION: The age-dependent effects suggest that retrospective studies could have underestimated sex and genetic variants' risk roles during younger ages. Therefore, we emphasize the importance of a time-to-event approach using longitudinal data to better characterize age-dependent risk effects.

Familywise error for multiple time-to-event endpoints in a group sequential design.

We investigate the familywise error rate (FWER) for time-to-event endpoints evaluated using a group sequential design with a hierarchical testing procedure for secondary endpoints. We show that, in this setup, the correlation between the log-rank test statistics at interim and at end of study is not congruent with the canonical correlation derived for normal-distributed endpoints. We show, both theoretically and by simulation, that the correlation also depends on the level of censoring, the hazard rates of the endpoints, and the hazard ratio. To optimize operating characteristics in this complex scenario, we propose a simulation-based method to assess the FWER which, better than the alpha-spending approach, can inform the choice of critical values for testing secondary endpoints.

Evaluating prognostic biomarkers for survival outcomes subject to informative censoring.

Prognostic biomarkers for survival outcomes are widely used in clinical research and practice. Such biomarkers are often evaluated using a C-index as well as quantities based on time-dependent receiver operating characteristic curves. Existing methods for their evaluation generally assume that censoring is uninformative in the sense that the censoring time is independent of the failure time with or without conditioning on the biomarker under evaluation. With focus on the C-index and the area under a particular receiver operating characteristic curve, we describe and compare three estimation methods that account for informative censoring based on observed baseline covariates. Two of them are straightforward extensions of existing plug-in and inverse probability weighting methods for uninformative censoring. By appealing to semiparametric theory, we also develop a doubly robust, locally efficient method that is more robust than the plug-in and inverse probability weighting methods and typically more efficient than the inverse probability weighting method. The methods are evaluated and compared in a simulation study, and applied to real data from studies of breast cancer and heart failure.

Bayesian survival analysis with INLA.

This tutorial shows how various Bayesian survival models can be fitted using the integrated nested Laplace approximation in a clear, legible, and comprehensible manner using the INLA and INLAjoint R-packages. Such models include accelerated failure time, proportional hazards, mixture cure, competing risks, multi-state, frailty, and joint models of longitudinal and survival data, originally presented in the article "Bayesian survival analysis with BUGS." In addition, we illustrate the implementation of a new joint model for a longitudinal semicontinuous marker, recurrent events, and a terminal event. Our proposal aims to provide the reader with syntax examples for implementing survival models using a fast and accurate approximate Bayesian inferential approach.

Predictability of buprenorphine-naloxone treatment retention: A multi-site analysis combining electronic health records and machine learning.

BACKGROUND AND AIMS: Opioid use disorder (OUD) and opioid dependence lead to significant morbidity and mortality, yet treatment retention, crucial for the effectiveness of medications like buprenorphine-naloxone, remains unpredictable. Our objective was to determine the predictability of 6-month retention in buprenorphine-naloxone treatment using electronic health record (EHR) data from diverse clinical settings and to identify key predictors. DESIGN: This retrospective observational study developed and validated machine learning-based clinical risk prediction models using EHR data. SETTING AND CASES: Data were sourced from Stanford University's healthcare system and Holmusk's NeuroBlu database, reflecting a wide range of healthcare settings. The study analyzed 1800 Stanford and 7957 NeuroBlu treatment encounters from 2008 to 2023 and from 2003 to 2023, respectively. MEASUREMENTS: Predict continuous prescription of buprenorphine-naloxone for at least 6 months, without a gap of more than 30 days. The performance of machine learning prediction models was assessed by area under receiver operating characteristic (ROC-AUC) analysis as well as precision, recall and calibration. To further validate our approach's clinical applicability, we conducted two secondary analyses: a time-to-event analysis on a single site to estimate the duration of buprenorphine-naloxone treatment continuity evaluated by the C-index and a comparative evaluation against predictions made by three human clinical experts. FINDINGS: Attrition rates at 6 months were 58% (NeuroBlu) and 61% (Stanford). Prediction models trained and internally validated on NeuroBlu data achieved ROC-AUCs up to 75.8 (95% confidence interval [CI] = 73.6-78.0). Addiction medicine specialists' predictions show a ROC-AUC of 67.8 (95% CI = 50.4-85.2). Time-to-event analysis on Stanford data indicated a median treatment retention time of 65 days, with random survival forest model achieving an average C-index of 65.9. The top predictor of treatment retention identified included the diagnosis of opioid dependence. CONCLUSIONS: US patients with opioid use disorder or opioid dependence treated with buprenorphine-naloxone prescriptions appear to have a high (∼60%) treatment attrition by 6 months. Machine learning models trained on diverse electronic health record datasets appear to be able to predict treatment continuity with accuracy comparable to that of clinical experts.

Associations of Urine and Plasma Metabolites with Kidney Failure and Death in a CKD Cohort.

RATIONALE & OBJECTIVE: Biomarkers that enable better identification of persons with chronic kidney disease (CKD) who are at higher risk for disease progression and adverse events are needed. This study sought to identify urine and plasma metabolites associated with progression of kidney disease. STUDY DESIGN: Prospective metabolome-wide association study. SETTING & PARTICIPANTS: Persons with CKD enrolled in the German CKD Study (GCKD) with metabolite measurements; with external validation within the Atherosclerosis Risk in Communities Study. EXPOSURES: 1,513 urine and 1,416 plasma metabolites (Metabolon, Inc.) measured at study entry using untargeted mass spectrometry. OUTCOMES: Main endpoints were kidney failure (KF), and a composite endpoint of KF, eGFR <15 mL/min/1.73m2, or 40% decline in eGFR (CKE). Death from any cause was a secondary endpoint. After a median of 6.5 years follow-up, 500 persons experienced KF, 1,083 experienced CKE and 680 died. ANALYTICAL APPROACH: Time-to-event analyses using multivariable proportional hazard regression models in a discovery-replication design, with external validation. RESULTS: 5,088 GCKD participants were included in analyses of urine metabolites and 5,144 in analyses of plasma metabolites. Among 182 unique metabolites, 30 were significantly associated with KF, 49 with CKE, and 163 with death. The strongest association with KF was observed for plasma hydroxyasparagine (hazard ratio: 1.95, 95% confidence interval: 1.68-2.25). An unnamed metabolite measured in plasma and urine was significantly associated with KF, CKE, and death. External validation of the identified associations of metabolites with KF or CKE revealed direction-consistency for 88% of observed associations. Selected associations of 18 metabolites with study outcomes have not been previously reported. LIMITATIONS: Use of observational data and semi-quantitative metabolite measurements at a single time point. CONCLUSIONS: The observed associations between metabolites and KF, CKE or death in persons with CKD confirmed previously reported findings and also revealed several associations not previously described. These findings warrant confirmatory research in other study cohorts.

A nonparametric proportional risk model to assess a treatment effect in time-to-event data.

Time-to-event analysis often relies on prior parametric assumptions, or, if a semiparametric approach is chosen, Cox's model. This is inherently tied to the assumption of proportional hazards, with the analysis potentially invalidated if this assumption is not fulfilled. In addition, most interpretations focus on the hazard ratio, that is often misinterpreted as the relative risk (RR), the ratio of the cumulative distribution functions. In this paper, we introduce an alternative to current methodology for assessing a treatment effect in a two-group situation, not relying on the proportional hazards assumption but assuming proportional risks. Precisely, we propose a new nonparametric model to directly estimate the RR of two groups to experience an event under the assumption that the risk ratio is constant over time. In addition to this relative measure, our model allows for calculating the number needed to treat as an absolute measure, providing the possibility of an easy and holistic interpretation of the data. We demonstrate the validity of the approach by means of a simulation study and present an application to data from a large randomized controlled trial investigating the effect of dapagliflozin on all-cause mortality.

A Comparison of Additional Benefit Assessment Methods for Time-to-Event Endpoints Using Hazard Ratio Point Estimates or Confidence Interval Limits by Means of a Simulation Study.

BACKGROUND: For time-to-event endpoints, three additional benefit assessment methods have been developed aiming at an unbiased knowledge about the magnitude of clinical benefit of newly approved treatments. The American Society of Clinical Oncology (ASCO) defines a continuous score using the hazard ratio point estimate (HR-PE). The European Society for Medical Oncology (ESMO) and the German Institute for Quality and Efficiency in Health Care (IQWiG) developed methods with an ordinal outcome using lower and upper limits of the 95% HR confidence interval (HR-CI), respectively. We describe all three frameworks for additional benefit assessment aiming at a fair comparison across different stakeholders. Furthermore, we determine which ASCO score is consistent with which ESMO/IQWiG category. METHODS: In a comprehensive simulation study with different failure time distributions and treatment effects, we compare all methods using Spearman's correlation and descriptive measures. For determination of ASCO values consistent with categories of ESMO/IQWiG, maximizing weighted Cohen's Kappa approach was used. RESULTS: Our research depicts a high positive relationship between ASCO/IQWiG and a low positive relationship between ASCO/ESMO. An ASCO score smaller than 17, 17 to 20, 20 to 24, and greater than 24 corresponds to ESMO categories. Using ASCO values of 21 and 38 as cutoffs represents IQWiG categories. LIMITATIONS: We investigated the statistical aspects of the methods and hence implemented slightly reduced versions of all methods. CONCLUSIONS: IQWiG and ASCO are more conservative than ESMO, which often awards the maximal category independent of the true effect and is at risk of overcompensating with various failure time distributions. ASCO has similar characteristics as IQWiG. Delayed treatment effects and underpowered/overpowered studies influence all methods in some degree. Nevertheless, ESMO is the most liberal one. HIGHLIGHTS: For the additional benefit assessment, the American Society of Clinical Oncology (ASCO) uses the hazard ratio point estimate (HR-PE) for their continuous score. In contrast, the European Society for Medical Oncology (ESMO) and the German Institute for Quality and Efficiency in Health Care (IQWiG) use the lower and upper 95% HR confidence interval (HR-CI) to specific thresholds, respectively. ESMO generously assigns maximal scores, while IQWiG is more conservative.This research provides the first comparison between IQWiG and ASCO and describes all three frameworks for additional benefit assessment aiming for a fair comparison across different stakeholders. Furthermore, thresholds for ASCO consistent with ESMO and IQWiG categories are determined, enabling a comparison of the methods in practice in a fair manner.IQWiG and ASCO are the more conservative methods, while ESMO awards high percentages of maximal categories, especially with various failure time distributions. ASCO has similar characteristics as IQWiG. Delayed treatment effects and under/-overpowered studies influence all methods. Nevertheless, ESMO is the most liberal one. An ASCO score smaller than 17, 17 to 20, 20 to 24, and greater than 24 correspond to the categories of ESMO. Using ASCO values of 21 and 38 as cutoffs represents categories of IQWiG.

Bayesian pairwise meta-analysis of time-to-event outcomes in the presence of non-proportional hazards: A simulation study of flexible parametric, piecewise exponential and fractional polynomial models.

BACKGROUND: Traditionally, meta-analysis of time-to-event outcomes reports a single pooled hazard ratio assuming proportional hazards (PH). For health technology assessment evaluations, hazard ratios are frequently extrapolated across a lifetime horizon. However, when treatment effects vary over time, an assumption of PH is not always valid. The Royston-Parmar (RP), piecewise exponential (PE), and fractional polynomial (FP) models can accommodate non-PH and provide plausible extrapolations of survival curves beyond observed data. METHODS: Simulation study to assess and compare the performance of RP, PE, and FP models in a Bayesian framework estimating restricted mean survival time difference (RMSTD) at 50 years from a pairwise meta-analysis with evidence of non-PH. Individual patient data were generated from a mixture Weibull distribution. Twelve scenarios were considered varying the amount of follow-up data, number of trials in a meta-analysis, non-PH interaction coefficient, and prior distributions. Performance was assessed through bias and mean squared error. Models were applied to a metastatic breast cancer example. RESULTS: FP models performed best when the non-PH interaction coefficient was 0.2. RP models performed best in scenarios with complete follow-up data. PE models performed well on average across all scenarios. In the metastatic breast cancer example, RMSTD at 50-years ranged from -14.6 to 8.48 months. CONCLUSIONS: Synthesis of time-to-event outcomes and estimation of RMSTD in the presence of non-PH can be challenging and computationally intensive. Different approaches make different assumptions regarding extrapolation and sensitivity analyses varying key assumptions are essential to check the robustness of conclusions to different assumptions for the underlying survival function.

Time-to-Event Endpoints in Imaging Biomarker Studies.

Time-to-event endpoints are widely used as measures of patients' well-being and indicators of prognosis. In imaging-based biomarker studies, there are increasingly more studies that focus on examining imaging biomarkers' prognostic or predictive utilities on those endpoints, whether in a trial or an observational study setting. In this educational review article, we briefly introduce some basic concepts of time-to-event endpoints and point out potential pitfalls in the context of imaging biomarker research in hope of improving radiologists' understanding of related subjects. Besides, we have included some review and discussions on the benefits of using time-to-event endpoints and considerations on selecting overall survival or progression-free survival for primary analysis. LEVEL OF EVIDENCE: 5 TECHNICAL EFFICACY: Stage 3.

Exploring factors influencing time from dispatch to unit availability according to the transport decision in the pre-hospital setting: an exploratory study.

BACKGROUND: Efficient resource distribution is important. Despite extensive research on response timings within ambulance services, nuances of time from unit dispatch to becoming available still need to be explored. This study aimed to identify the determinants of the duration between ambulance dispatch and readiness to respond to the next case according to the patients' transport decisions. METHODS: Time from ambulance dispatch to availability (TDA) analysis according to the patients' transport decision (Transport versus Non-Transport) was conducted using R-Studio™ for a data set of 93,712 emergency calls managed by a Middle Eastern ambulance service from January to May 2023. Log-transformed Hazard Ratios (HR) were examined across diverse parameters. A Cox regression model was utilised to determine the influence of variables on TDA. Kaplan-Meier curves discerned potential variances in the time elapsed for both cohorts based on demographics and clinical indicators. A competing risk analysis assessed the probabilities of distinct outcomes occurring. RESULTS: The median duration of elapsed TDA was 173 min for the transported patients and 73 min for those not transported. The HR unveiled Significant associations in various demographic variables. The Kaplan-Meier curves revealed variances in TDA across different nationalities and age categories. In the competing risk analysis, the 'Not Transported' group demonstrated a higher incidence of prolonged TDA than the 'Transported' group at specified time points. CONCLUSIONS: Exploring TDA offers a novel perspective on ambulance services' efficiency. Though promising, the findings necessitate further exploration across diverse settings, ensuring broader applicability. Future research should consider a comprehensive range of variables to fully harness the utility of this period as a metric for healthcare excellence.

[Special issue PRO] A demonstration of estimands and sensitivity analyses for time-to-deterioration of patient reported outcomes.

In oncology trials, health-related quality of life (HRQoL), specifically patient-reported symptom burden and functional status, can support the interpretation of survival endpoints, such as progression-free survival. However, applying time-to-event endpoints to patient-reported outcomes (PRO) data is challenging. For example, in time-to-deterioration analyses clinical events such as disease progression are common in many settings and are often handled through censoring the patient at the time of occurrence; however, disease progression and HRQoL are often related leading to informative censoring. Special consideration to the definition of events and intercurrent events (ICEs) is necessary. In this work, we demonstrate time-to-deterioration of PRO estimands and sensitivity analyses to answer research questions using composite, hypothetical, and treatment policy strategies applied to a single endpoint of disease-related symptoms. Multiple imputation methods under both the missing-at-random and missing-not-at-random assumptions are used as sensitivity analyses of primary estimands. Hazard ratios ranged from 0.52 to 0.66 over all the estimands and sensitivity analyses modeling a robust treatment effect favoring the treatment in time to disease symptom deterioration or death. Differences in the estimands include how people who experience disease progression or discontinue the randomized treatment due to AEs are accounted for in the analysis. We use the estimand framework to define interpretable and principled approaches for different time-to-deterioration research questions and provide practical recommendations. Reporting the proportions of patient events and patient censoring by reason helps understand the mechanisms that drive the results, allowing for optimal interpretation.