Logistic random effects regression models: a comparison of statistical packages for binary and ordinal outcomes

<p>Abstract</p> <p>Background</p> <p>Logistic random effects models are a popular tool to analyze multilevel also called hierarchical data with a binary or ordinal outcome. Here, we aim to compare different statistical software implementations of these models.</p> <p>Methods</p> <p>We used individual patient data from 8509 patients in 231 centers with moderate and severe Traumatic Brain Injury (TBI) enrolled in eight Randomized Controlled Trials (RCTs) and three observational studies. We fitted logistic random effects regression models with the 5-point Glasgow Outcome Scale (GOS) as outcome, both dichotomized as well as ordinal, with center and/or trial as random effects, and as covariates age, motor score, pupil reactivity or trial. We then compared the implementations of frequentist and Bayesian methods to estimate the fixed and random effects. Frequentist approaches included R (lme4), Stata (GLLAMM), SAS (GLIMMIX and NLMIXED), MLwiN ([R]IGLS) and MIXOR, Bayesian approaches included WinBUGS, MLwiN (MCMC), R package MCMCglmm and SAS experimental procedure MCMC.</p> <p>Three data sets (the full data set and two sub-datasets) were analysed using basically two logistic random effects models with either one random effect for the center or two random effects for center and trial. For the ordinal outcome in the full data set also a proportional odds model with a random center effect was fitted.</p> <p>Results</p> <p>The packages gave similar parameter estimates for both the fixed and random effects and for the binary (and ordinal) models for the main study and when based on a relatively large number of level-1 (patient level) data compared to the number of level-2 (hospital level) data. However, when based on relatively sparse data set, i.e. when the numbers of level-1 and level-2 data units were about the same, the frequentist and Bayesian approaches showed somewhat different results. The software implementations differ considerably in flexibility, computation time, and usability. There are also differences in the availability of additional tools for model evaluation, such as diagnostic plots. The experimental SAS (version 9.2) procedure MCMC appeared to be inefficient.</p> <p>Conclusions</p> <p>On relatively large data sets, the different software implementations of logistic random effects regression models produced similar results. Thus, for a large data set there seems to be no explicit preference (of course if there is no preference from a philosophical point of view) for either a frequentist or Bayesian approach (if based on vague priors). The choice for a particular implementation may largely depend on the desired flexibility, and the usability of the package. For small data sets the random effects variances are difficult to estimate. In the frequentist approaches the MLE of this variance was often estimated zero with a standard error that is either zero or could not be determined, while for Bayesian methods the estimates could depend on the chosen "non-informative" prior of the variance parameter. The starting value for the variance parameter may be also critical for the convergence of the Markov chain.</p

Genome-wide copy number variation study in anorectal malformations

Anorectal malformations (ARMs, congenital obstruction of the anal opening) are among the most common birth defects requiring surgical treatment (2-5/10 000 live-births) and carry significant chronic morbidity. ARMs present either as isolated or as part of the phenotypic spectrum of some chromosomal abnormalities or monogenic syndromes. The etiology is unknown. To assess the genetic contribution to ARMs, we investigated single-nucleotide polymorphisms and copy number variations (CNVs) at genome-wide scale. A total of 363 Han Chinese sporadic ARM patients and 4006 Han Chinese controls were included. Overall, we detected a 1.3-fold significant excess of rare CNVs in patients. Stratification of patients by presence/absence of other congenital anomalies showed that while syndromic ARM patients carried significantly longer rare duplications than controls (P = 0.049), non-syndromic patients were enriched with both rare deletions and duplications when compared with controls (P = 0.00031). Twelve chromosomal aberrations and 114 rare CNVs were observed in patients but not in 868 controls nor 11 943 healthy individuals from the Database of Genomic Variants. Importantly, these aberrations were observed in isolated ARM patients. Gene-based analysis revealed 79 genes interfered by CNVs in patients only. In particular, we identified a de novo DKK4 duplication. DKK4 is a member of the WNT signaling pathway which is involved in the development of the anorectal region. In mice, Wnt disruption results in ARMs. Our data suggest a role for rare CNVs not only in syndromic but also in isolated ARM patients and provide a list of plausible candidate genes for the disorder.postprin

Methodological considerations in injury burden of disease studies across Europe: a systematic literature review

Background Calculating the disease burden due to injury is complex, as it requires many methodological choices. Until now, an overview of the methodological design choices that have been made in burden of disease (BoD) studies in injury populations is not available. The aim of this systematic literature review was to identify existing injury BoD studies undertaken across Europe and to comprehensively review the methodological design choices and assumption parameters that have been made to calculate years of life lost (YLL) and years lived with disability (YLD) in these studies. Methods We searched EMBASE, MEDLINE, Cochrane Central, Google Scholar, and Web of Science, and the grey literature supplemented by handsearching, for BoD studies. We included injury BoD studies that quantified the BoD expressed in YLL, YLD, and disability-adjusted life years (DALY) in countries within the European Region between early-1990 and mid-2021. Results We retrieved 2,914 results of which 48 performed an injury-specific BoD assessment. Single-country independent and Global Burden of Disease (GBD)-linked injury BoD studies were performed in 11 European countries. Approximately 79% of injury BoD studies reported the BoD by external cause-of-injury. Most independent studies used the incidence-based approach to calculate YLDs. About half of the injury disease burden studies applied disability weights (DWs) developed by the GBD study. Almost all independent injury studies have determined YLL using national life tables. Conclusions Considerable methodological variation across independent injury BoD assessments was observed; differences were mainly apparent in the design choices and assumption parameters towards injury YLD calculations, implementation of DWs, and the choice of life table for YLL calculations. Development and use of guidelines for performing and reporting of injury BoD studies is crucial to enhance transparency and comparability of injury BoD estimates across Europe and beyond

Walkway lengths for measuring walking speed in stroke rehabilitation

Author name used in this publication: Shamay S. M. Ng2011-2012 > Academic research: refereed > Publication in refereed journalVersion of RecordPublishe

Differential effects of alpha-adrenoceptor agonists on relaxation in the mesenteric artery and aorta of the rat

Poster session 7: Adrenoceptors / Adrenergic Nerve : poster no. P7-