Modeling multilevel sleep transitional data via Poisson log-linear multilevel models

This paper proposes Poisson log-linear multilevel models to investigate population variability in sleep state transition rates. We specifically propose a Bayesian Poisson regression model that is more flexible, scalable to larger studies, and easily fit than other attempts in the literature. We further use hierarchical random effects to account for pairings of individuals and repeated measures within those individuals, as comparing diseased to non-diseased subjects while minimizing bias is of epidemiologic importance. We estimate essentially non-parametric piecewise constant hazards and smooth them, and allow for time varying covariates and segment of the night comparisons. The Bayesian Poisson regression is justified through a re-derivation of a classical algebraic likelihood equivalence of Poisson regression with a log(time) offset and survival regression assuming piecewise constant hazards. This relationship allows us to synthesize two methods currently used to analyze sleep transition phenomena: stratified multi-state proportional hazards models and log-linear models with GEE for transition counts. An example data set from the Sleep Heart Health Study is analyzed

Modeling Multilevel Sleep Transitional Data Via Poisson Log-Linear Multilevel Models

This paper proposes Poisson log-linear multilevel models to investigate population variability in sleep state transition rates. We specifically propose a Bayesian Poisson regression model that is more flexible, scalable to larger studies, and easily fit than other attempts in the literature. We further use hierarchical random effects to account for pairings of individuals and repeated measures within those individuals, as comparing diseased to non-diseased subjects while minimizing bias is of epidemiologic importance. We estimate essentially non-parametric piecewise constant hazards and smooth them, and allow for time varying covariates and segment of the night comparisons. The Bayesian Poisson regression is justified through a re-derivation of a classical algebraic likelihood equivalence of Poisson regression with a log(time) offset and survival regression assuming piecewise constant hazards. This relationship allows us to synthesize two methods currently used to analyze sleep transition phenomena: stratified multi-state proportional hazards models and log-linear models with GEE for transition counts. An example data set from the Sleep Heart Health Study is analyzed

A unified approach to modeling multivariate binary data using copulas over partitions

Many seemingly disparate approaches for marginal modeling have been developed in recent years. We demonstrate that many current approaches for marginal modeling of correlated binary outcomes produce likelihoods that are equivalent to the proposed copula-based models herein. These general copula models of underlying latent threshold random variables yield likelihood based models for marginal fixed effects estimation and interpretation in the analysis of correlated binary data. Moreover, we propose a nomenclature and set of model relationships that substantially elucidates the complex area of marginalized models for binary data. A diverse collection of didactic mathematical and numerical examples are given to illustrate concepts

RESTRICTED LIKELIHOOD RATIO TESTS FOR FUNCTIONAL EFFECTS IN THE FUNCTIONAL LINEAR MODEL

The goal of our article is to provide a transparent, robust, and computationally feasible statistical approach for testing in the context of scalar-on-function linear regression models. In particular, we are interested in testing for the necessity of functional effects against standard linear models. Our methods are motivated by and applied to a large longitudinal study involving diffusion tensor imaging of intracranial white matter tracts in a susceptible cohort. In the context of this study, we conduct hypothesis tests that are motivated by anatomical knowledge and which support recent findings regarding the relationship between cognitive impairment and white matter demyelination. R-code and data are provided to reproduce the application

MODELING SLEEP FRAGMENTATION IN POPULATIONS OF SLEEP HYPNOGRAMS

We introduce methods for the analysis of large populations of sleep architectures (hypnograms) that respect the 5-state 20-transition-type structure defined by the American Academy of Sleep Medicine. By applying these methods to the hypnograms of 5598 subjects from the Sleep Heart Health Study we: 1) provide the firrst analysis of sleep hypnogram data of such size and complexity in a community cohort with a 4-level comorbidity; 2) compare 5-state 20-transition-type sleep to 3-state 6-transition-type sleep for a check of feasibility and information-loss; 3) extend current approaches to multivariate survival data analysis to populations of time-to-transition processes; and 4) provide scalable solutions for data analyses required by the case study. This allows us to provide detailed new insights into the association between sleep apnea and sleep architecture. Supporting R as well as SAS code and data are included in the online supplementary materials

Mixed effect Poisson log-linear models for clinical and epidemiological sleep hypnogram data

Bayesian Poisson log-linear multilevel models scalable to epidemiological studies are proposed to investigate population variability in sleep state transition rates. Hierarchical random effects are used to account for pairings of individuals and repeated measures within those individuals, as comparing diseased to non-diseased subjects while minimizing bias is of importance. Essentially, non-parametric piecewise constant hazards are estimated and smoothed, allowing for time-varying covariates and segment of the night comparisons. The Bayesian Poisson regression is justified through a re-derivation of a classical algebraic likelihood equivalence of Poisson regression with a log(time) offset and survival regression assuming exponentially distributed survival times. Such re-derivation allows synthesis of two methods currently used to analyze sleep transition phenomena: stratified multi-state proportional hazards models and log-linear models with GEE for transition counts. An example data set from the Sleep Heart Health Study is analyzed. This paper is a revamping of Modeling multilevel sleep transitional data via Poisson log-linear multilevel models available at: http://www.bepress.com/jhubiostat/paper212

LASAGNA PLOTS: A SAUCY ALTERNATIVE TO SPAGHETTI PLOTS

Longitudinal repeated measures data has often been visualized with spaghetti plots for continuous out- comes. For large datasets, this often leads to over-plotting and consequential obscuring of trends in the data. This is primarily due to overlapping of trajectories. Here, we suggest a framework called lasagna plot ting that constrains the subject-specific trajectories to prevent overlapping and utilizes gradients of color to depict the outcome. Dynamic sorting and visualization is demonstrated as an exploratory data analysis tool. Supplemental material in the form of sample R code additional illustrated examples are available online

Pregnancy outcomes in a malaria-exposed Malian cohort of women of child-bearing age

In Sub-Saharan Africa, malaria continues to be associated with adverse pregnancy outcomes including stillbirth, early neonatal death, preterm delivery, and low birth weight. Current preventive measures are insufficient and new interventions are urgently needed. However, before such interventions can be tested in pregnant women, background information on pregnancy outcomes in this target population must be collected. We conducted an observational study in Ouélessébougou, Mali, a malaria-endemic area where first antenatal visit commonly occurs during the second trimester of pregnancy, hindering calculation of miscarriage rate in the population. To accurately determine the rate of miscarriage, 799 non-pregnant women of child-bearing age were enrolled and surveyed via monthly follow up visits that included pregnancy tests. Out of 505 women that completed the study, 364 became pregnant and 358 pregnancies were analyzed: 43 (12%) resulted in miscarriage, 28 (65.1%) occurred during the first trimester of pregnancy. We also determined rates of stillbirth, neonatal death, preterm delivery, and small for gestational age. The results showed high rate of miscarriage during the first trimester and established a basis to evaluate new interventions to prevent pregnancy malaria. This survey design enabled identification of first trimester miscarriages that are often missed by studies conducted in antenatal clinics.Clinical trial registration[https://clinicaltrials.gov/], identifier [NCT0297 4608]