We propose an extensive framework for additive regression models for correlated functional responses, allowing for multiple partially nested or crossed functional random effects with flexible correlation structures for, for example, spatial, temporal, or longitudinal functional data. Additionally, our framework includes linear and nonlinear effects of functional and scalar covariates that may vary smoothly over the index of the functional response. It accommodates densely or sparsely observed functional responses and predictors which may be observed with additional error and includes both spline-based and functional principal component-based terms. Estimation and inference in this framework is based on standard additive mixed models, allowing us to take advantage of established methods and robust, flexible algorithms. We provide easy-to-use open source software in the pffr() function for the R package refund. Simulations show that the proposed method recovers relevant effects reliably, handles small sample sizes well, and also scales to larger datasets. Applications with spatially and longitudinally observed functional data demonstrate the flexibility in modeling and interpretability of results of our approach.
Background and Purpose—: Ischemic stroke of undetermined cause is a major health issue because of its high frequency and clinical relevance. Histopathologic analysis of human thrombi, retrieved from stroke patients with large-vessel occlusion during mechanical thrombectomy, may provide information about underlying pathologies. This study examines the relationship between stroke causes and histological clot composition to identify specific patterns that might help to distinguish causes of cryptogenic stroke. Methods—: Thrombi of 145 consecutive stroke patients with large-vessel occlusion were collected during intracranial mechanical recanalization. The hematoxylin and eosin–stained specimens were quantitatively analyzed in terms of the relative fractions of the main constituents (red and white blood cells and fibrin/platelets). These data, along with additional clinical and interventional parameters, were compared for different stroke subtypes, as defined by the international Trial of Org 10172 in Acute Stroke Treatment criteria. Results—: The composition of thrombi from cardioembolic and noncardioembolic stroke patients differed significantly for all main thrombus components. Cardioembolic thrombi had higher proportions of fibrin/platelets ( P =0.009), less erythrocytes ( P =0.003), and more leucocytes ( P =0.035) than noncardioembolic thrombi. Cryptogenic strokes showed strong overlap with cardioembolic strokes but not with noncardioembolic strokes, in terms of both thrombus histology and interventional and clinical outcome parameters. Conclusions—: Quantitative evaluation of thrombus composition may help to distinguish between different stroke causes. Our findings support the notion that the majority of cryptogenic strokes are cardioembolic.
1 March 2008 in Computational Statistics & Data Analysis
Summary: Distributed lag non-linear models (DLNMs) are a modelling tool for describing potentially non-linear and delayed dependencies. Here, we illustrate an extension of the DLNM framework through the use of penalized splines within generalized additive models (GAM). This extension offers built-in model selection procedures and the possibility of accommodating assumptions on the shape of the lag structure through specific penalties. In addition, this framework includes, as special cases, simpler models previously proposed for linear relationships (DLMs). Alternative versions of penalized DLNMs are compared with each other and with the standard unpenalized version in a simulation study. Results show that this penalized extension to the DLNM class provides greater flexibility and improved inferential properties. The framework exploits recent theoretical developments of GAMs and is implemented using efficient routines within freely available software. Real-data applications are illustrated through two reproducible examples in time series and survival analysis.