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-causemortality.

Conditional power (CP) serves as a widely utilized approach for futility monitoring in group sequential designs. However, adopting the CP methods may lead to inadequate control of the type II error rate at the desired level. In this study, we introduce a flexible beta spending function tailored to regulate the type II error rate while employing CP based on a predetermined standardized effect size for futility monitoring (a so-called CP-beta spending function). This function delineates the expenditure of type II error rate across the entirety of the trial. Unlike other existing beta spending functions, the CP-beta spending function seamlessly incorporates beta spending concept into the CP framework, facilitating precise stagewise control of the type II error rate during futility monitoring. In addition, the stopping boundaries derived from the CP-beta spending function can be calculated via integration akin to other traditional beta spending function methods. Furthermore, the proposed CP-beta spending function accommodates various thresholds on the CP-scale at different stages of the trial, ensuring its adaptability across different information time scenarios. These attributes render the CP-beta spending function competitive among other forms of beta spending functions, making it applicable to any trials in group sequential designs with straightforward implementation. Both simulation study and example from an acute ischemic stroke trial demonstrate that the proposed method accurately captures expected power, even when the initially determined sample size does not consider futility stopping, and exhibits a good performance in maintaining overall type I error rates for evidentfutility.

Many statistical models have been proposed in the literature for the analysis of longitudinal data. One may propose to model two or more correlated longitudinal processes simultaneously, with a goal of understanding their association over time. Joint modeling is then required to carefully study the association structure among the outcomes as well as drawing joint inferences about the different outcomes. In this study, we sought to model the associations among six nutrition outcomes while circumventing the computational challenge posed by their clustered and high-dimensional nature. We analyzed data from a 2 2 randomized crossover trial conducted in Kenya, to compare the effect of high-dose and low-dose iodine in household salt on systolic blood pressure (SBP) and diastolic blood pressure (DBP) in women of reproductive age and their household matching pair of school-aged children. Two additional outcomes, namely, urinary iodine concentration (UIC) in women and children were measured repeatedly to monitor the amount of iodine excreted through urine. We extended the model proposed by Mwangi etal. (2021, Communications in Statistics: Case Studies, Data Analysis and Applications, 7(3), 413-431) allowing flexible piecewise joint models for six outcomes to depend on separate random effects, which are themselves correlated. This entailed fitting 15 bivariate general linear mixed models and deriving inference for the joint model using pseudo-likelihood theory. We analyzed the outcomes separately and jointly using piecewise linear mixed-effects (PLME) model and further validated the results using current state-of-the-art Jones and Kenward methodology (JKME model) used for analyzing randomized crossover trials. The results indicate that high-dose iodine in salt significantly reduced blood pressure (BP) compared to low-dose iodine in salt. Estimates for the random effects and residual error components showed that SBP and DBP had strong positive correlation, with effect of the random slope indicating that significantly related outcomes are strongly associated in their evolution. There was a moderately strong inverse relationship between evolutions of UIC and BP both in women and children. These findings confirmed the original hypothesis that high-dose iodine salt has significant lowering effect on BP. We further sought to evaluate the performance of our proposed PLME model against the widely used JKME model, within the multivariate joint modeling framework through a simulation study mimicking a crossover design. From our findings, the multivariate joint PLME model performed exceptionally well both in estimation of random-effects matrix (G) and Hessian matrix (H), allowing satisfactory model convergence during estimation. It allowed a more complex fit to the data with both random intercepts and slopes effects compared to the multivariate joint JKME model that allowed for random intercepts only. When a hierarchical viewpoint is adopted, in the sense that outcomes are specified conditionally upon random effects, the variance-covariance matrix of the random effects must be positive definite. In some cases, additional random effects could explain much variability in the data, thus improving precision in estimation of the estimands (effect size) parameters. The key highlight in this evaluation shows that multivariate joint JKME model is a powerful tool especially while fitting mixed models with random intercepts only, in crossover design settings. Addition of random slopes may lead to model complexities in most cases, resulting in unsatisfactory model convergence during estimation. To circumvent convergence pitfalls, extention of JKME model to PLME model allows a more flexible fit to the data (generated from crossover design settings), especially in the multivariate joint modelingframework.

In this contribution, the Schemper-Henderson measure of explained variation for survival outcomes is extended to accommodate competing events (CEs) in addition to events of interest. The extension is achieved by moving from the unconditional and conditional survival functions of the original measure to unconditional and conditional cumulative incidence functions, the latter obtained, for example, from Fine and Gray models. In the absence of CEs, the original measure is obtained as a special case. We define explained variation on the population level and provide two different types of estimates. Recently, the authors have achieved a multiplicative decomposition of explained variation into degrees of necessity and degrees of sufficiency. These measures are also extended to the case of competing risks survival data. A SAS macro and an R function are provided to facilitate application. Interesting empirical properties of the measures are explored on the population level and by an extensive simulation study. Advantages of the approach are exemplified by an Austrian study of breast cancer with a high proportion of CEs.