ABCtoolbox: a versatile toolkit for approximate Bayesian computations

Background The estimation of demographic parameters from genetic data often requires the computation of likelihoods. However, the likelihood function is computationally intractable for many realistic evolutionary models, and the use of Bayesian inference has therefore been limited to very simple models. The situation changed recently with the advent of Approximate Bayesian Computation (ABC) algorithms allowing one to obtain parameter posterior distributions based on simulations not requiring likelihood computations. Results Here we present ABCtoolbox, a series of open source programs to perform Approximate Bayesian Computations (ABC). It implements various ABC algorithms including rejection sampling, MCMC without likelihood, a Particle-based sampler and ABC-GLM. ABCtoolbox is bundled with, but not limited to, a program that allows parameter inference in a population genetics context and the simultaneous use of different types of markers with different ploidy levels. In addition, ABCtoolbox can also interact with most simulation and summary statistics computation programs. The usability of the ABCtoolbox is demonstrated by inferring the evolutionary history of two evolutionary lineages of Microtus arvalis. Using nuclear microsatellites and mitochondrial sequence data in the same estimation procedure enabled us to infer sex-specific population sizes and migration rates and to find that males show smaller population sizes but much higher levels of migration than females. Conclusion ABCtoolbox allows a user to perform all the necessary steps of a full ABC analysis, from parameter sampling from prior distributions, data simulations, computation of summary statistics, estimation of posterior distributions, model choice, validation of the estimation procedure, and visualization of the results.

Patient-Oriented SCORAD (PO-SCORAD): a new self-assessment scale in atopic dermatitis validated in Europe

BACKGROund: Patient-oriented medicine is an emerging concept, encouraged by the World Health Organization, to greater involvement of the patient in the management of chronic diseases. The Patient-Oriented SCORing Atopic Dermatitis (PO-SCORAD) index is a self-assessment score allowing the patient to comprehensively evaluate the actual course of atopic dermatitis (AD), using subjective and objective criteria derived mainly from the SCORAD, a validated AD severity clinical assessment tool.

Computations of Vector-Valued Siegel Modular Forms

We carry out some computations of vector-valued Siegel modular forms of degree two, weight (k, 2) and level one, and highlight three experimental results: (1) we identify a rational eigenform in a three-dimensional space of cusp forms; (2) we observe that non-cuspidal eigenforms of level one are not always rational; (3) we verify a number of cases of conjectures about congruences between classical modular forms and Siegel modular forms. Our approach is based on Satoh's description of the module of vector-valued Siegel modular forms of weight (k, 2) and an explicit description of the Hecke action on Fourier expansions. (C) 2013 Elsevier Inc. All rights reserved.

Impact of complications in total ankle replacement and ankle arthrodesis analyzed with a validated outcome measurement

Major modifications in the design and techniques of total ankle replacement have challenged the perception that ankle arthrodesis is the treatment of choice for end-stage ankle arthritis. High complication and revision rates have been reported after both procedures.

Validated capillary electrophoresis assay for the simultaneous enantioselective determination of propafenone and its major metabolites in biological samples

A robust, inexpensive, and fully validated CE method for the simultaneous determination of the enantiomers of propafenone (PPF), 5-hydroxy-propafenone (5OH-PPF) and N-despropyl-propafenone (NOR-PPF) in serum and in in vitro media is described. It is based upon liquid-liquid extraction at alkaline pH followed by analysis of the reconstituted extract by CE in presence of a pH 2.0 running buffer composed of 100 mM sodium phosphate, 19% methanol, and 0.6% highly sulfated beta-CD. For each compound, the S-enantiomers are shown to migrate ahead of their antipodes, and the overall run time is about 30 min. Enantiomer levels between 25 and 1000 ng/mL provide linear calibration graphs, and the LOD for all enantiomers is between 10 and 12 ng/mL. The assay is shown to be suitable for the determination of the enantiomers of PPF and its metabolites in in vitro incubations comprising human liver microsomes or single CYP450 enzymes (SUPERSOMES). Incubations with CYP2D6 SUPERSOMES revealed, for the first time, the simultaneous formation of the enantiomers of 5OH-PPF and NOR-PPF with that enzyme. CE data can be used for the evaluation of the enzymatic N-dealkylation and hydroxylation rates.

DISTRIBUTED REPRODUCIBLE RESEARCH USING CACHED COMPUTATIONS

The ability to make scientific findings reproducible is increasingly important in areas where substantive results are the product of complex statistical computations. Reproducibility can allow others to verify the published findings and conduct alternate analyses of the same data. A question that arises naturally is how can one conduct and distribute reproducible research? This question is relevant from the point of view of both the authors who want to make their research reproducible and readers who want to reproduce relevant findings reported in the scientific literature. We present a framework in which reproducible research can be conducted and distributed via cached computations and describe specific tools for both authors and readers. As a prototype implementation we introduce three software packages written in the R language. The cacheSweave and stashR packages together provide tools for caching computational results in a key-value style database which can be published to a public repository for readers to download. The SRPM package provides tools for generating and interacting with "shared reproducibility packages" (SRPs) which can facilitate the distribution of the data and code. As a case study we demonstrate the use of the toolkit on a national study of air pollution exposure and mortality.

Radiographic analysis of femoroacetabular impingement with Hip2Norm-reliable and validated

The purpose of this study was to validate the accuracy, consistency, and reproducibility/reliability of a new method for correction of pelvic tilt and rotation of radiographic hip parameters for pincer type of femoroacetabular impingement on an anteroposterior pelvic radiograph. Thirty cadaver hips and 100 randomized, blinded AP pelvic radiographs were used for investigation. To detect the software accuracy, the calculated femoral head coverage and classic hip parameters determined with our software were compared to reference measurements based on CT scans or conventional radiographs in a neutral orientation as gold standard. To investigate software consistency, differences among the different parameters for each cadaver pelvis were calculated when reckoned back from a random to the neutral orientation. Intra- and interobserver comparisons were used to analyze the reliability and reproducibility of all parameters. All but two parameters showed a good-to-very good accuracy with the reference measurements. No relevant systematic errors were detected in the Bland-Altman analysis. Software consistency was good-to-very good for all parameters. A good-to-very good reliability and reproducibility was found for a substantial number of the evaluated radiographic acetabular parameters. The software appears to be an accurate, consistent, reliable, and reproducible method for analysis of acetabular pathomorphologies.

Synthesis of computations and experiments for obtaining pulsatile gas flow rates from dynamic pressure difference measurements across an orifice-plate meter

The use of conventional orifice-plate meter is typically restricted to measurements of steady flows. This study proposes a new and effective computational-experimental approach for measuring the time-varying (but steady-in-the-mean) nature of turbulent pulsatile gas flows. Low Mach number (effectively constant density) steady-in-the-mean gas flows with large amplitude fluctuations (whose highest significant frequency is characterized by the value fF) are termed pulsatile if the fluctuations have a direct correlation with the time-varying signature of the imposed dynamic pressure difference and, furthermore, they have fluctuation amplitudes that are significantly larger than those associated with turbulence or random acoustic wave signatures. The experimental aspect of the proposed calibration approach is based on use of Coriolis-meters (whose oscillating arm frequency fcoriolis >> fF) which are capable of effectively measuring the mean flow rate of the pulsatile flows. Together with the experimental measurements of the mean mass flow rate of these pulsatile flows, the computational approach presented here is shown to be effective in converting the dynamic pressure difference signal into the desired dynamic flow rate signal. The proposed approach is reliable because the time-varying flow rate predictions obtained for two different orifice-plate meters exhibit the approximately same qualitative, dominant features of the pulsatile flow.

A validated assay by liquid chromatography-tandem mass spectrometry for the simultaneous quantification of elvitegravir and rilpivirine in HIV positive patients

Because of the large variability in the pharmacokinetics of anti-HIV drugs, therapeutic drug monitoring in patients may contribute to optimize the overall efficacy and safety of antiretroviral therapy. An LC–MS/MS method for the simultaneous assay in plasma of the novel antiretroviral agents rilpivirine (RPV) and elvitegravir (EVG) has been developed to that endeavor. Plasma samples (100 μL) extraction is performed by protein precipitation with acetonitrile, and the supernatant is subsequently diluted 1:1 with 20-mM ammonium acetate/MeOH 50:50. After reverse-phase chromatography, quantification of RPV and EVG, using matrix-matched calibration samples, is performed by electrospray ionization–triple quadrupole mass spectrometry by selected reaction monitoring detection using the positive mode. The stable isotopic-labeled compounds RPV-13C6 and EVG-D6 were used as internal standards. The method was validated according to FDA recommendations, including assessment of extraction yield, matrix effects variability (<6.4%), as well as EVG and RPV short and long-term stability in plasma. Calibration curves were validated over the clinically relevant concentrations ranging from 5 to 2500 ng/ml for RPV and from 50 to 5000 ng/ml for EVG. The method is precise (inter-day CV%: 3–6.3%) and accurate (3.8–7.2%). Plasma samples were found to be stable (<15%) in all considered conditions (RT/48 h, +4°C/48 h, −20°C/3 months and 60°C/1 h). Selected metabolite profiles analysis in patients' samples revealed the presence of EVG glucuronide, that was well separated from parent EVG, allowing to exclude potential interferences through the in-source dissociation of glucuronide to parent drug. This new, rapid and robust LCMS/MS assay for the simultaneous quantification of plasma concentrations of these two major new anti-HIV drugs EVG and RPV offers an efficient analytical tool for clinical pharmacokinetics studies and routine therapeutic drug monitoring service.

Rib fractures at postmortem computed tomography (PMCT) validated against the autopsy

To evaluate the sensitivity of postmortem computed tomography (PMCT) in rib fracture detection validated against autopsy. Fifty-one forensic cases underwent a postmortem CT prior to forensic autopsy. Two image readers (radiologist and forensic pathologist) assessed high resolution CT data sets for rib fractures. Correct recognition rates (CRR), sensitivity and specificity values were calculated over all observations as well as individually for every rib and region. Additionally, for partial rib fractures the sensitivity of autopsy was calculated vice versa. 3876 entries in each study protocol (autopsy, PMCT radiologist and PMCT forensic pathologist) were investigated. A total of 690 fractures (autopsy), 491 (PMCT and radiologist) and 559 (PMCT and forensic pathologist) were detected. The CRR was 0.85. Sensitivity and specificity of PMCT for rib fracture detection were 0.63 (0.58 radiologist, 0.68 forensic pathologist) and 0.97 (both readers 0.97), respectively. Low CRR and sensitivity values were obtained for antero-lateral fractures. Partial rib fractures were better detected by PMCT. PMCT has a rather low sensitivity for rib fracture detection when validated against autopsy and indicates that clinical CT may also demonstrate a reasonable number of false negatives. Partial rib fractures often remain undetected at autopsy.

An experimentally validated finite element method for augmented vertebral bodies

Background Finite element models of augmented vertebral bodies require a realistic modelling of the cement infiltrated region. Most methods published so far used idealized cement shapes or oversimplified material models for the augmented region. In this study, an improved, anatomy-specific, homogenized finite element method was developed and validated to predict the apparent as well as the local mechanical behavior of augmented vertebral bodies. Methods Forty-nine human vertebral body sections were prepared by removing the cortical endplates and scanned with high-resolution peripheral quantitative CT before and after injection of a standard and a low-modulus bone cement. Forty-one specimens were tested in compression to measure stiffness, strength and contact pressure distributions between specimens and loading-plates. From the remaining eight, fourteen cylindrical specimens were extracted from the augmented region and tested in compression to obtain material properties. Anatomy-specific finite element models were generated from the CT data. The models featured element-specific, density-fabric-based material properties, damage accumulation, real cement distributions and experimentally determined material properties for the augmented region. Apparent stiffness and strength as well as contact pressure distributions at the loading plates were compared between simulations and experiments. Findings The finite element models were able to predict apparent stiffness (R2 > 0.86) and apparent strength (R2 > 0.92) very well. Also, the numerically obtained pressure distributions were in reasonable quantitative (R2 > 0.48) and qualitative agreement with the experiments. Interpretation The proposed finite element models have proven to be an accurate tool for studying the apparent as well as the local mechanical behavior of augmented vertebral bodies.