Zinpentraxin Alfa for Idiopathic Pulmonary Fibrosis: The Randomized Phase III STARSCAPE Trial.

Rationale: A phase II trial reported clinical benefit over 28 weeks in patients with idiopathic pulmonary fibrosis (IPF) who received zinpentraxin alfa. Objectives: To investigate the efficacy and safety of zinpentraxin alfa in patients with IPF in a phase III trial. Methods: This 52-week phase III, double-blind, placebo-controlled, pivotal trial was conducted at 275 sites in 29 countries. Patients with IPF were randomized 1:1 to intravenous placebo or zinpentraxin alfa 10 mg/kg every 4 weeks. The primary endpoint was absolute change from baseline to Week 52 in FVC. Secondary endpoints included absolute change from baseline to Week 52 in percent predicted FVC and 6-minute walk distance. Safety was monitored via adverse events. Post hoc analysis of the phase II and phase III data explored changes in FVC and their impact on the efficacy results. Measurements and Main Results: Of 664 randomized patients, 333 were assigned to placebo and 331 to zinpentraxin alfa. Four of the 664 randomized patients were never administered study drug. The trial was terminated early after a prespecified futility analysis that demonstrated no treatment benefit of zinpentraxin alfa over placebo. In the final analysis, absolute change from baseline to Week 52 in FVC was similar between placebo and zinpentraxin alfa (-214.89 ml and -235.72 ml; P = 0.5420); there were no apparent treatment effects on secondary endpoints. Overall, 72.3% and 74.6% of patients receiving placebo and zinpentraxin alfa, respectively, experienced one or more adverse events. Post hoc analysis revealed that extreme FVC decline in two placebo-treated patients resulted in the clinical benefit of zinpentraxin alfa reported by phase II. Conclusions: Zinpentraxin alfa treatment did not benefit patients with IPF over placebo. Learnings from this program may help improve decision making around trials in IPF. Clinical trial registered with www.clinicaltrials.gov (NCT04552899).

Modeling Risk Factors for Intraindividual Variability: A Mixed-Effects Beta-Binomial Model Applied to Cognitive Function in Older People in the English Longitudinal Study of Ageing.

Cognitive functioning in older age profoundly impacts quality of life and health. While most research on cognition in older age has focused on mean levels, intraindividual variability (IIV) around this may have risk factors and outcomes independent of the mean value. Investigating risk factors associated with IIV has typically involved deriving a summary statistic for each person from residual error around a fitted mean. However, this ignores uncertainty in the estimates, prohibits exploring associations with time-varying factors, and is biased by floor/ceiling effects. To address this, we propose a mixed-effects location scale beta-binomial model for estimating average probability and IIV in a word recall test in the English Longitudinal Study of Ageing. After adjusting for mean performance, an analysis of 9,873 individuals across 7 (mean = 3.4) waves (2002-2015) found IIV to be greater at older ages, with lower education, in females, with more difficulties in activities of daily living, in later birth cohorts, and when interviewers recorded issues potentially affecting test performance. Our study introduces a novel method for identifying groups with greater IIV in bounded discrete outcomes. Our findings have implications for daily functioning and care, and further work is needed to identify the impact for future health outcomes.

Merging or ensembling: integrative analysis in multiple neuroimaging studies.

The aim of this paper is to systematically investigate merging and ensembling methods for spatially varying coefficient mixed effects models (SVCMEM) in order to carry out integrative learning of neuroimaging data obtained from multiple biomedical studies. The "merged" approach involves training a single learning model using a comprehensive dataset that encompasses information from all the studies. Conversely, the "ensemble" approach involves creating a weighted average of distinct learning models, each developed from an individual study. We systematically investigate the prediction accuracy of the merged and ensemble learners under the presence of different degrees of interstudy heterogeneity. Additionally, we establish asymptotic guidelines for making strategic decisions about when to employ either of these models in different scenarios, along with deriving optimal weights for the ensemble learner. To validate our theoretical results, we perform extensive simulation studies. The proposed methodology is also applied to 3 large-scale neuroimaging studies.

Longitudinal associations of plasma amino acid levels with recovery from malarial coma.

BACKGROUND: Disordered amino acid metabolism is observed in cerebral malaria (CM). This study sought to determine whether abnormal amino acid concentrations were associated with level of consciousness in children recovering from coma. Twenty-one amino acids and coma scores were quantified longitudinally and the data were analysed for associations. METHODS: In a prospective observational study, 42 children with CM were enrolled. Amino acid levels were measured at entry and at frequent intervals thereafter and consciousness was assessed by Blantyre Coma Scores (BCS). Thirty-six healthy children served as controls for in-country normal amino acid ranges. Logistic regression was employed using a generalized linear mixed-effects model to assess associations between out-of-range amino acid levels and BCS. RESULTS: At entry 16/21 amino acid levels were out-of-range. Longitudinal analysis revealed 10/21 out-of-range amino acids were significantly associated with BCS. Elevated phenylalanine levels showed the highest association with low BCS. This finding held when out-of-normal-range data were analysed at each sampling time. CONCLUSION: Longitudinal data is provided for associations between abnormal amino acid levels and recovery from CM. Of 10 amino acids significantly associated with BCS, elevated phenylalanine may be a surrogate for impaired clearance of ether lipid mediators of inflammation and may contribute to CM pathogenesis.

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.

Classified functional mixed effects model prediction.

In nowadays biomedical research, there has been a growing demand for making accurate prediction at subject levels. In many of these situations, data are collected as longitudinal curves and display distinct individual characteristics. Thus, prediction mechanisms accommodated with functional mixed effects models (FMEM) are useful. In this paper, we developed a classified functional mixed model prediction (CFMMP) method, which adapts classified mixed model prediction (CMMP) to the framework of FMEM. Performance of CFMMP against functional regression prediction based on simulation studies and the consistency property of CFMMP estimators are explored. Real-world applications of CFMMP are illustrated using real world examples including data from the hormone research menstrual cycles and the diffusion tensor imaging.

A latent variable mixed-effects location scale model that also considers between-person differences in the autocorrelation.

In public health research an increasing number of studies is conducted in which intensive longitudinal data is collected in an experience sampling or a daily diary design. Typically, the resulting data is analyzed with a mixed-effects model or mixed-effects location scale model because they allow one to examine a host of interesting longitudinal research questions. Here, we introduce an extension of the mixed-effects location scale model in which measurement error of the observed variables is considered by a latent factor model and in which-in addition to the mean-or location-related effects-the residual variance of the latent factor and the parameters of the autoregressive process of this latent factor can differ between persons. We show how to estimate the parameters of the model with a maximum likelihood approach, whose performance is also compared with a Bayesian approach in a small simulation study. We illustrate the models using a real data example and end with a discussion in which we suggest questions for future research.

A novel block-coordinate gradient descent algorithm for simultaneous grouped selection of fixed and random effects in joint modeling.

Joint models for longitudinal and time-to-event data are receiving increasing attention owing to its capability of capturing the possible association between these two types of data. Typically, a joint model consists of a longitudinal submodel for longitudinal processes and a survival submodel for the time-to-event response, and links two submodels by common covariates that may carry both fixed and random effects. However, research gaps still remain on how to simultaneously select fixed and random effects from the two submodels under the joint modeling framework efficiently and effectively. In this article, we propose a novel block-coordinate gradient descent (BCGD) algorithm to simultaneously select multiple longitudinal covariates that may carry fixed and random effects in the joint model. Specifically, for the multiple longitudinal processes, a linear mixed effect model is adopted where random intercepts and slopes serve as essential covariates of the trajectories, and for the survival submodel, the popular proportional hazard model is employed. A penalized likelihood estimation is used to control the dimensionality of covariates in the joint model and estimate the unknown parameters, especially when estimating the covariance matrix of random effects. The proposed BCGD method can successfully capture the useful covariates of both fixed and random effects with excellent selection power, and efficiently provide a relatively accurate estimate of fixed and random effects empirically. The simulation results show excellent performance of the proposed method and support its effectiveness. The proposed BCGD method is further applied on two real data sets, and we examine the risk factors for the effects of different heart valves, differing on type of tissue, implanted in the aortic position and the risk factors for the diagnosis of primary biliary cholangitis.

Comparative performance of contact plate metod and swab method for surface microbial contamination on medical fabrics.

BACKGROUND: The contact plate method is widely accepted and used in various fields where hygiene and contamination levels are crucial. Evidence regarding the applicability of the contact plate method for sampling fabric microbial contamination levels in real medical environments was limited. This study aimed to assess the applicability of the contact plate method for detecting microbial contamination on medical fabrics in a real healthcare environment, thereby providing a benchmark for fabric microbial sampling methods. METHODS: In a level three obstetrics ward of a hospital, twenty-four privacy curtains adjacent to patient beds were selected for this study. The contact plate and swab method were used to collect microbial samples from the privacy curtains on the 1st, 7th, 14th, and 28th days after they were hung. The total colony count on each privacy curtain surface was calculated, and microbial identification was performed. RESULTS: After excluding the effects of time, room type, and curtain location on the detected microbial load, the linear mixed-effects model analysis showed that contact plate method yielded lower colony counts compared to swab method (P < 0.001). However, the contact plate method isolated more microbial species than swab method (P < 0.001). 291 pathogenic strains were isolated using the contact plate method and 133 pathogenic strains were isolated via the swab method. There was no difference between the two sampling methods in the detection of gram-negative bacteria (P = 0.089). Furthermore, the microbial load on curtains in double-occupancy rooms was lower than those in triple-occupancy rooms (P = 0.021), and the microbial load on curtains near windows was lower than that near doors (P = 0.004). CONCLUSION: Contact plate method is superior to swab method in strain isolation. Swab method is more suitable for evaluating the bacterial contamination of fabrics.

A framework for longitudinal latent factor modelling of treatment response in clinical trials with applications to Psoriatic Arthritis and Rheumatoid Arthritis.

OBJECTIVE: Clinical trials involve the collection of a wealth of data, comprising multiple diverse measurements performed at baseline and follow-up visits over the course of a trial. The most common primary analysis is restricted to a single, potentially composite endpoint at one time point. While such an analytical focus promotes simple and replicable conclusions, it does not necessarily fully capture the multi-faceted effects of a drug in a complex disease setting. Therefore, to complement existing approaches, we set out here to design a longitudinal multivariate analytical framework that accepts as input an entire clinical trial database, comprising all measurements, patients, and time points across multiple trials. METHODS: Our framework composes probabilistic principal component analysis with a longitudinal linear mixed effects model, thereby enabling clinical interpretation of multivariate results, while handling data missing at random, and incorporating covariates and covariance structure in a computationally efficient and principled way. RESULTS: We illustrate our approach by applying it to four phase III clinical trials of secukinumab in Psoriatic Arthritis (PsA) and Rheumatoid Arthritis (RA). We identify three clinically plausible latent factors that collectively explain 74.5% of empirical variation in the longitudinal patient database. We estimate longitudinal trajectories of these factors, thereby enabling joint characterisation of disease progression and drug effect. We perform benchmarking experiments demonstrating our method's competitive performance at estimating average treatment effects compared to existing statistical and machine learning methods, and showing that our modular approach leads to relatively computationally efficient model fitting. CONCLUSION: Our multivariate longitudinal framework has the potential to illuminate the properties of existing composite endpoint methods, and to enable the development of novel clinical endpoints that provide enhanced and complementary perspectives on treatment response.

A meta-analysis of microbial thermal inactivation in low moisture foods.

Microbial thermal inactivation in low moisture foods is challenging due to enhanced thermal resistance of microbes and low thermal conductivity of food matrices. In this study, we leveraged the body of previous work on this topic to model key experimental features that determine microbial thermal inactivation in low moisture foods. We identified 27 studies which contained 782 mean D-values and developed linear mixed-effect models to assess the effect of microorganism type, matrix structure and composition, water activity, temperature, and inoculation and recovery methods on cell death kinetics. Intraclass correlation statistics (I2) and conditional R2 values of the linear mixed effects models were: E. coli (R2-0.91, I2-83%), fungi (R2-0.88, I2-85%), L. monocytogenes (R2-0.84, I2-75%), Salmonella (R2-0.69, I2-46%). Finally, global response surface models (RSM) were developed to further study the non-linear effect of aw and temperature on inactivation. The fit of these models varied by organisms from R2 0.88 (E. coli) to 0.35 (fungi). Further dividing the Salmonella data into individual RSM models based on matrix structure improved model fit to R2 0.90 (paste-like products) and 0.48 (powder-like products). This indicates a negative relationship between data diversity and model performance.

Modeling spatial correlation between earthquake insured losses in New Zealand: A mixed-effects analysis.

Earthquake insurance is a critical risk management strategy that contributes to improving recovery and thus greater resilience of individuals. Insurance companies construct premiums without taking into account spatial correlations between insured assets. This leads to potentially underestimating the risk, and therefore the exceedance probability curve. We here propose a mixed-effects model to estimate losses per ward that is able to account for heteroskedasticity and spatial correlation between insured losses. Given the significant impact of earthquakes in New Zealand due to its particular geographical and demographic characteristics, the government has established a public insurance company that collects information about the insured buildings and any claims lodged. We thus develop a two-level variance component model that is based on earthquake losses observed in New Zealand between 2000 and 2021. The proposed model aims at capturing the variability at both the ward and territorial authority levels and includes independent variables, such as seismic hazard indicators, the number of usual residents, and the average dwelling value in the ward. Our model is able to detect spatial correlation in the losses at the ward level thus increasing its predictive power and making it possible to assess the effect of spatially correlated claims that may be considerable on the tail of loss distribution.

Peer sexual harassment among 10-year-olds: Roles, genders, classroom occurrence, and associations with emotional problems.

This is the first study examining peer sexual harassment among 10-year-olds (N = 985), studying how being a victim, perpetrator, or witness relates to emotional problems, and how these associations are moderated by gender and class occurrence of sexual harassment. Results showed that 45% of the participants reported victimization, 17% perpetration, and 60% witnessing sexual harassment, with vast overlaps between roles. Victimization and witnessing were related to more emotional problems. Victimized girls reported more emotional problems than boys, but girls who perpetrated reported fewer emotional problems than boys. Associations between peer sexual harassment and emotional problems varied across classrooms. Our findings highlight the occurrence of peer sexual harassment in younger ages, emphasizing an ecological perspective when addressing it in school.

Permutation Tests for Assessing Potential Non-Linear Associations between Treatment Use and Multivariate Clinical Outcomes.

In many psychometric applications, the relationship between the mean of an outcome and a quantitative covariate is too complex to be described by simple parametric functions; instead, flexible nonlinear relationships can be incorporated using penalized splines. Penalized splines can be conveniently represented as a linear mixed effects model (LMM), where the coefficients of the spline basis functions are random effects. The LMM representation of penalized splines makes the extension to multivariate outcomes relatively straightforward. In the LMM, no effect of the quantitative covariate on the outcome corresponds to the null hypothesis that a fixed effect and a variance component are both zero. Under the null, the usual asymptotic chi-square distribution of the likelihood ratio test for the variance component does not hold. Therefore, we propose three permutation tests for the likelihood ratio test statistic: one based on permuting the quantitative covariate, the other two based on permuting residuals. We compare via simulation the Type I error rate and power of the three permutation tests obtained from joint models for multiple outcomes, as well as a commonly used parametric test. The tests are illustrated using data from a stimulant use disorder psychosocial clinical trial.

Linear Mixed-Effects Models for Dependent Data: Power and Accuracy in Parameter Estimation.

Linear mixed-effects models have been increasingly used to analyze dependent data in psychological research. Despite their many advantages over ANOVA, critical issues in their analyses remain. Due to increasing random effects and model complexity, estimation computation is demanding, and convergence becomes challenging. Applied users need help choosing appropriate methods to estimate random effects. The present Monte Carlo simulation study investigated the impacts when the restricted maximum likelihood (REML) and Bayesian estimation models were misspecified in the estimation. We also compared the performance of Akaike information criterion (AIC) and deviance information criterion (DIC) in model selection. Results showed that models neglecting the existing random effects had inflated Type I errors, unacceptable coverage, and inaccurate R-squared measures of fixed and random effects variation. Furthermore, models with redundant random effects had convergence problems, lower statistical power, and inaccurate R-squared measures for Bayesian estimation. The convergence problem is more severe for REML, while reduced power and inaccurate R-squared measures were more severe for Bayesian estimation. Notably, DIC was better than AIC in identifying the true models (especially for models including person random intercept only), improving convergence rates, and providing more accurate effect size estimates, despite AIC having higher power than DIC with 10 items and the most complicated true model.

Secure Messaging Use Among Patients with Depression: An Analysis Using Real-World Data.

Background: Although depression is one of the most common mental health disorders outpacing other diseases and conditions, poor access to care and limited resources leave many untreated. Secure messaging (SM) offers patients an online means to bridge this gap by communicating nonurgent medical questions. We focused on self-care health management behaviors and delved into SM initiation as the initial act of engagement and SM exchanges as continuous engagement patterns. This study examined whether those with depression might be using SM more than those without depression. Methods: Patient portal data were obtained from a large academic medical center's electronic health records spanning 5 years, from January 2018 to December 2022. We organized and analyzed SM initiations and exchanges using the linear mixed-effects modeling technique. Results: Our predictors correlated with SM initiations, accounting for 25.1% of variance explained. In parallel, 24.9% of SM exchanges were attributable to these predictors. Overall, our predictors demonstrate stronger associations with SM exchanges. Discussion: We examined patients with and without depression across 2,629 zip codes over five years. Our findings reveal that the predictors affecting SM initiations and exchanges are multifaceted, with certain predictors enhancing its utilization and others impeding it. Conclusions: SM telehealth service provided support to patients with mental health needs to a greater extent than those without. By increasing access, fostering better communication, and efficiently allocating resources, telehealth services not only encourage patients to begin using SM but also promote sustained interaction through ongoing SM exchanges.

Antagonistic Effects of Assortative Mating on the Evolution of Phenotypic Plasticity along Environmental Gradients.

AbstractPrevious theory has shown that assortative mating for plastic traits can maintain genetic divergence across environmental gradients despite high gene flow. Yet these models did not examine how assortative mating affects the evolution of plasticity. We here describe patterns of genetic variation across elevation for plasticity in a trait under assortative mating, using multiple-year observations of budburst date in a common garden of sessile oaks. Despite high gene flow, we found significant spatial genetic divergence for the intercept, but not for the slope, of reaction norms to temperature. We then used individual-based simulations, where both the slope and the intercept of the reaction norm evolve, to examine how assortative mating affects the evolution of plasticity, varying the intensity and distance of gene flow. Our model predicts the evolution of either suboptimal plasticity (reaction norms with a slope shallower than optimal) or hyperplasticity (slopes steeper than optimal) in the presence of assortative mating when optimal plasticity would evolve under random mating. Furthermore, a cogradient pattern of genetic divergence for the intercept of the reaction norm (where plastic and genetic effects are in the same direction) always evolves in simulations with assortative mating, consistent with our observations in the studied oak populations.

Efficient model-based bioequivalence testing.

The classical approach to analyze pharmacokinetic (PK) data in bioequivalence studies aiming to compare two different formulations is to perform noncompartmental analysis (NCA) followed by two one-sided tests (TOST). In this regard, the PK parameters area under the curve (AUC) and $C_{\max}$ are obtained for both treatment groups and their geometric mean ratios are considered. According to current guidelines by the U.S. Food and Drug Administration and the European Medicines Agency, the formulations are declared to be sufficiently similar if the $90\%$ confidence interval for these ratios falls between $0.8$ and $1.25 $. As NCA is not a reliable approach in case of sparse designs, a model-based alternative has already been proposed for the estimation of $\rm AUC$ and $C_{\max}$ using nonlinear mixed effects models. Here we propose another, more powerful test than the TOST and demonstrate its superiority through a simulation study both for NCA and model-based approaches. For products with high variability on PK parameters, this method appears to have closer type I errors to the conventionally accepted significance level of $0.05$, suggesting its potential use in situations where conventional bioequivalence analysis is not applicable.

Individual participant data meta-analysis with mixed-effects transformation models.

One-stage meta-analysis of individual participant data (IPD) poses several statistical and computational challenges. For time-to-event outcomes, the approach requires the estimation of complicated nonlinear mixed-effects models that are flexible enough to realistically capture the most important characteristics of the IPD. We present a model class that incorporates general normally distributed random effects into linear transformation models. We discuss extensions to model between-study heterogeneity in baseline risks and covariate effects and also relax the assumption of proportional hazards. Within the proposed framework, data with arbitrary random censoring patterns can be handled. The accompanying $\textsf{R}$ package tramME utilizes the Laplace approximation and automatic differentiation to perform efficient maximum likelihood estimation and inference in mixed-effects transformation models. We compare several variants of our model to predict the survival of patients with chronic obstructive pulmonary disease using a large data set of prognostic studies. Finally, a simulation study is presented that verifies the correctness of the implementation and highlights its efficiency compared to an alternative approach.

SCEBE: an efficient and scalable algorithm for genome-wide association studies on longitudinal outcomes with mixed-effects modeling.

Genome-wide association studies (GWAS) using longitudinal phenotypes collected over time is appealing due to the improvement of power. However, computation burden has been a challenge because of the complex algorithms for modeling the longitudinal data. Approximation methods based on empirical Bayesian estimates (EBEs) from mixed-effects modeling have been developed to expedite the analysis. However, our analysis demonstrated that bias in both association test and estimation for the existing EBE-based methods remains an issue. We propose an incredibly fast and unbiased method (simultaneous correction for EBE, SCEBE) that can correct the bias in the naive EBE approach and provide unbiased P-values and estimates of effect size. Through application to Alzheimer's Disease Neuroimaging Initiative data with 6 414 695 single nucleotide polymorphisms, we demonstrated that SCEBE can efficiently perform large-scale GWAS with longitudinal outcomes, providing nearly 10 000 times improvement of computational efficiency and shortening the computation time from months to minutes. The SCEBE package and the example datasets are available at https://github.com/Myuan2019/SCEBE.