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

Photoperiod Regulates Lean Mass Accretion, but Not Adiposity, in Growing F344 Rats Fed a High Fat Diet

yesIn this study the effects of photoperiod and diet, and their interaction, were examined for their effects on growth and body composition in juvenile F344 rats over a 4-week period. On long (16L:8D), relative to short (8L:16D), photoperiod food intake and growth rate were increased, but percentage adiposity remained constant (ca 3-4%). On a high fat diet (HFD), containing 22.8% fat (45% energy as fat), food intake was reduced, but energy intake increased on both photoperiods. This led to a small increase in adiposity (up to 10%) without overt change in body weight. These changes were also reflected in plasma leptin and lipid levels. Importantly while both lean and adipose tissue were strongly regulated by photoperiod on a chow diet, this regulation was lost for adipose, but not lean tissue, on HFD. This implies that a primary effect of photoperiod is the regulation of growth and lean mass accretion. Consistent with this both hypothalamic GHRH gene expression and serum IGF-1 levels were photoperiod dependent. As for other animals and humans, there was evidence of central hyposomatotropism in response to obesity, as GHRH gene expression was suppressed by the HFD. Gene expression of hypothalamic AgRP and CRH, but not NPY nor POMC, accorded with the energy balance status on long and short photoperiod. However, there was a general dissociation between plasma leptin levels and expression of these hypothalamic energy balance genes. Similarly there was no interaction between the HFD and photoperiod at the level of the genes involved in thyroid hormone metabolism (Dio2, Dio3, TSHβ or NMU), which are important mediators of the photoperiodic response. These data suggest that photoperiod and HFD influence body weight and body composition through independent mechanisms but in each case the role of the hypothalamic energy balance genes is not predictable based on their known function.Scottish Government (Rural and Environment Science and Analytical Services Division, http://www.scotland.gov.uk/), AWR LR LMT PJM and the BBSRC, (http://www.bbsrc.ac.uk/home/home.aspx, grant BB/K001043/1), AWR GH PJ

Amide-controlled, one-pot synthesis of tri-substituted purines generates structural diversity and analogues with trypanocidal activity

Anovel one-pot synthesis of tri-substituted purines and the discovery of purine analogues with trypanocidal activity are reported. The reaction is initiated by a metal-free oxidative coupling of primary alkoxides and diaminopyrimidines with Schiff base formation and subsequent annulation in the presence of large N,N-dimethylamides (e.g.N,N-dimethylpropanamide or larger). This synthetic route is in competition with a reaction previously-reported by our group1, allowing the generation of a combinatorial library of tri-substituted purines by the simple modification of the amide and the alkoxide employed. Among the variety of structures generated, two purine analogues displayed trypanocidal activity against the protozoan parasite Trypanosoma brucei with IC50 , 5 mM, being each of those compounds obtained through each of the synthetic pathways.J.J.D.M. thanks Spanish Ministerio de Economı´a y Competitividad for a Ramon y Cajal Fellowship. A.U.B. thanks MRC IGMM for an academic fellowship. This work was partially supported by Grant SAF2011-30528 to J.A.G.S.