A Hierarchical Bayesian Model for Measuring Motion Adaptation

Previous research has shown that motion imagery draws on the same neural circuits that are involved in perception of motion, thus leading to a motion aftereffect (Winawer et al., 2010). Imagined stimuli can induce a similar shift in participants’ psychometric functions as neural adaptation due to a perceived stimulus. However, these studies have been criticized on the grounds that they fail to exclude the possibility that the subjects might have guessed the experimental hypothesis, and behaved accordingly (Morgan et al., 2012). In particular, the authors claim that participants can adopt arbitrary response criteria, which results in similar changes of the central tendency μ of psychometric curves as those shown by Winawer et al. (2010).

Bayesian statistics in oncology: a guide for the clinical investigator

The rise of evidence-based medicine as well as important progress in statistical methods and computational power have led to a second birth of the >200-year-old Bayesian framework. The use of Bayesian techniques, in particular in the design and interpretation of clinical trials, offers several substantial advantages over the classical statistical approach. First, in contrast to classical statistics, Bayesian analysis allows a direct statement regarding the probability that a treatment was beneficial. Second, Bayesian statistics allow the researcher to incorporate any prior information in the analysis of the experimental results. Third, Bayesian methods can efficiently handle complex statistical models, which are suited for advanced clinical trial designs. Finally, Bayesian statistics encourage a thorough consideration and presentation of the assumptions underlying an analysis, which enables the reader to fully appraise the authors' conclusions. Both Bayesian and classical statistics have their respective strengths and limitations and should be viewed as being complementary to each other; we do not attempt to make a head-to-head comparison, as this is beyond the scope of the present review. Rather, the objective of the present article is to provide a nonmathematical, reader-friendly overview of the current practice of Bayesian statistics coupled with numerous intuitive examples from the field of oncology. It is hoped that this educational review will be a useful resource to the oncologist and result in a better understanding of the scope, strengths, and limitations of the Bayesian approach.

Differential healing response attributed to culprit lesions of patients with acute coronary syndromes and stable coronary artery after implantation of drug-eluting stents: An optical coherence tomography study

BACKGROUND Pathology studies have shown delayed arterial healing in culprit lesions of patients with acute coronary syndrome (ACS) compared with stable coronary artery disease (CAD) after placement of drug-eluting stents (DES). It is unknown whether similar differences exist in-vivo during long-term follow-up. Using optical coherence tomography (OCT), we assessed differences in arterial healing between patients with ACS and stable CAD five years after DES implantation. METHODS AND RESULTS A total of 88 patients comprised of 53 ACS lesions with 7864 struts and 35 stable lesions with 5298 struts were suitable for final OCT analysis five years after DES implantation. The analytical approach was based on a hierarchical Bayesian random-effects model. OCT endpoints were strut coverage, malapposition, protrusion, evaginations and cluster formation. Uncovered (1.7% vs. 0.7%, adjusted p=0.041) or protruding struts (0.50% vs. 0.13%, adjusted p=0.038) were more frequent among ACS compared with stable CAD lesions. A similar trend was observed for malapposed struts (1.33% vs. 0.45%, adj. p=0.072). Clusters of uncovered or malapposed/protruding struts were present in 34.0% of ACS and 14.1% of stable patients (adj. p=0.041). Coronary evaginations were more frequent in patients with ST-elevation myocardial infarction compared with stable CAD patients (0.16 vs. 0.13 per cross section, p=0.027). CONCLUSION Uncovered, malapposed, and protruding stent struts as well as clusters of delayed healing may be more frequent in culprit lesions of ACS compared with stable CAD patients late after DES implantation. Our observational findings suggest a differential healing response attributable to lesion characteristics of patients with ACS compared with stable CAD in-vivo.

Immune reconstitution inflammatory syndrome in patients starting antiretroviral therapy for HIV infection: a systematic review and meta-analysis

In patients with HIV-1 infection who are starting combination antiretroviral therapy (ART), the incidence of immune reconstitution inflammatory syndrome (IRIS) is not well defined. We did a meta-analysis to establish the incidence and lethality of the syndrome in patients with a range of previously diagnosed opportunistic infections, and examined the relation between occurrence and the degree of immunodeficiency. Systematic review identified 54 cohort studies of 13 103 patients starting ART, of whom 1699 developed IRIS. We calculated pooled cumulative incidences with 95% credibility intervals (CrI) by Bayesian methods and did a random-effects metaregression to analyse the relation between CD4 cell count and incidence of IRIS. In patients with previously diagnosed AIDS-defining illnesses, IRIS developed in 37.7% (95% CrI 26.6-49.4) of those with cytomegalovirus retinitis, 19.5% (6.7-44.8) of those with cryptococcal meningitis, 15.7% (9.7-24.5) of those with tuberculosis, 16.7% (2.3-50.7) of those with progressive multifocal leukoencephalopathy, and 6.4% (1.2-24.7) of those with Kaposi's sarcoma, and 12.2% (6.8-19.6) of those with herpes zoster. 16.1% (11.1-22.9) of unselected patients starting ART developed any type of IRIS. 4.5% (2.1-8.6) of patients with any type of IRIS died, 3.2% (0.7-9.2) of those with tuberculosis-associated IRIS died, and 20.8% (5.0-52.7) of those with cryptococcal meningitis died. Metaregression analyses showed that the risk of IRIS is associated with CD4 cell count at the start of ART, with a high risk in patients with fewer than 50 cells per microL. Occurrence of IRIS might therefore be reduced by initiation of ART before immunodeficiency becomes advanced.

Mortality from road traffic accidents in Switzerland: longitudinal and spatial analyses

Road traffic accidents (RTA) are an important cause of premature death. We examined socio-demographic and geographical determinants of RTA mortality in Switzerland by linking 2000 census data to RTA mortality records 2000-2005 (ICD-10 codes V00-V99). Data from 5.5 million residents aged 18-94 years, 1744 study areas, and 1620 RTA deaths were analyzed, including 978 deaths (60.4%) in motor vehicle occupants, 254 (15.7%) in motorcyclists, 107 (6.6%) in cyclists, and 259 (16.0%) in pedestrians. Weibull survival models and Bayesian methods were used to calculate hazard ratios (HR), and standardized mortality ratios (SMR) across study areas. Adjusted HR comparing women with men ranged from 0.04 (95% CI 0.02-0.07) in motorcyclists to 0.43 (95% CI 0.32-0.56) in pedestrians. There was a u-shaped relationship with age in motor vehicle occupants and motorcyclists. In cyclists and pedestrians, mortality increased after age 55 years. Mortality was higher in individuals with primary education (HR 1.53; 95% CI 1.29-1.81), and higher in single (HR 1.24; 95% CI 1.05-1.46), widowed (HR 1.31; 95% CI 1.05-1.65) and divorced individuals (HR 1.62; 95% CI 1.33-1.97), compared to persons with tertiary education or married persons. The association with education was particularly strong for pedestrians (HR 1.87; 95% CI 1.20-2.91). RTA mortality increased with decreasing population density of study areas for motor vehicle occupants (test for trend p<0.0001) and motorcyclists (p=0.0021) but not for cyclists (p=0.39) or pedestrians (p=0.29). SMR standardized for socio-demographic and geographical variables ranged from 82 to 190. Prevention efforts should aim to reduce inequities across socio-demographic and educational groups, and across geographical areas, with interventions targeted at high-risk groups and areas, and different traffic users, including pedestrians.

Patterns of morphological changes and hybridization between sympatric whitefish morphs (Coregonus spp.) in a Swiss lake: a role for eutrophication?

Whitefish, genus Coregonus, show exceptional levels of phenotypic diversity with sympatric morphs occurring in numerous postglacial lakes in the northern hemisphere. Here, we studied the effects of human-induced eutrophication on sympatric whitefish morphs in the Swiss lake, Lake Thun. In particular, we addressed the questions whether eutrophication (i) induced hybridization between two ecologically divergent summer-spawning morphs through a loss of environmental heterogeneity, and (ii) induced rapid adaptive morphological changes through changes in the food web structure. Genetic analysis based on 11 microsatellite loci of 282 spawners revealed that the pelagic and the benthic morph represent highly distinct gene pools occurring at different relative proportions on all seven known spawning sites. Gill raker counts, a highly heritable trait, showed nearly discrete distributions for the two morphs. Multilocus genotypes characteristic of the pelagic morph had more gill rakers than genotypes characteristic of benthic morph. Using Bayesian methods, we found indications of recent but limited introgressive hybridization. Comparisons with historical gill raker data yielded median evolutionary rates of 0.24 haldanes and median selection intensities of 0.27 for this trait in both morphs for 1948-2004 suggesting rapid evolution through directional selection at this trait. However, phenotypic plasticity as an alternative explanation for this phenotypic change cannot be discarded. We hypothesize that both the temporal shifts in mean gill raker counts and the recent hybridization reflect responses to changes in the trophic state of the lake induced by pollution in the 1960s, which created novel selection pressures with respect to feeding niches and spawning site preferences.

Widespread Signals of Convergent Adaptation to High Altitude in Asia and America

Living at high altitude is one of the most difficult challenges that humans had to cope with during their evolution. Whereas several genomic studies have revealed some of the genetic bases of adaptations in Tibetan, Andean, and Ethiopian populations, relatively little evidence of convergent evolution to altitude in different continents has accumulated. This lack of evidence can be due to truly different evolutionary responses, but it can also be due to the low power of former studies that have mainly focused on populations from a single geographical region or performed separate analyses on multiple pairs of populations to avoid problems linked to shared histories between some populations. We introduce here a hierarchical Bayesian method to detect local adaptation that can deal with complex demographic histories. Our method can identify selection occurring at different scales, as well as convergent adaptation in different regions. We apply our approach to the analysis of a large SNP data set from low- and high-altitude human populations from America and Asia. The simultaneous analysis of these two geographic areas allows us to identify several candidate genome regions for altitudinal selection, and we show that convergent evolution among continents has been quite common. In addition to identifying several genes and biological processes involved in high-altitude adaptation, we identify two specific biological pathways that could have evolved in both continents to counter toxic effects induced by hypoxia.

An empirical comparison of methods for meta-analysis of diagnostic accuracy showed hierarchical models are necessary

OBJECTIVE: Meta-analysis of studies of the accuracy of diagnostic tests currently uses a variety of methods. Statistically rigorous hierarchical models require expertise and sophisticated software. We assessed whether any of the simpler methods can in practice give adequately accurate and reliable results. STUDY DESIGN AND SETTING: We reviewed six methods for meta-analysis of diagnostic accuracy: four simple commonly used methods (simple pooling, separate random-effects meta-analyses of sensitivity and specificity, separate meta-analyses of positive and negative likelihood ratios, and the Littenberg-Moses summary receiver operating characteristic [ROC] curve) and two more statistically rigorous approaches using hierarchical models (bivariate random-effects meta-analysis and hierarchical summary ROC curve analysis). We applied the methods to data from a sample of eight systematic reviews chosen to illustrate a variety of patterns of results. RESULTS: In each meta-analysis, there was substantial heterogeneity between the results of different studies. Simple pooling of results gave misleading summary estimates of sensitivity and specificity in some meta-analyses, and the Littenberg-Moses method produced summary ROC curves that diverged from those produced by more rigorous methods in some situations. CONCLUSION: The closely related hierarchical summary ROC curve or bivariate models should be used as the standard method for meta-analysis of diagnostic accuracy.

Fully Automatic Segmentation of Brain Tumor Images using Support Vector Machine Classiffication in Combination with Hierarchical Conditional Random Field Regularization

Delineating brain tumor boundaries from magnetic resonance images is an essential task for the analysis of brain cancer. We propose a fully automatic method for brain tissue segmentation, which combines Support Vector Machine classification using multispectral intensities and textures with subsequent hierarchical regularization based on Conditional Random Fields. The CRF regularization introduces spatial constraints to the powerful SVM classification, which assumes voxels to be independent from their neighbors. The approach first separates healthy and tumor tissue before both regions are subclassified into cerebrospinal fluid, white matter, gray matter and necrotic, active, edema region respectively in a novel hierarchical way. The hierarchical approach adds robustness and speed by allowing to apply different levels of regularization at different stages. The method is fast and tailored to standard clinical acquisition protocols. It was assessed on 10 multispectral patient datasets with results outperforming previous methods in terms of segmentation detail and computation times.

Hierarchical accumulation of RyR post-translational modifications drives disease progression in dystrophic cardiomyopathy

AIMS:Duchenne muscular dystrophy (DMD) is a muscle disease with serious cardiac complications. Changes in Ca(2+) homeostasis and oxidative stress were recently associated with cardiac deterioration, but the cellular pathophysiological mechanisms remain elusive. We investigated whether the activity of ryanodine receptor (RyR) Ca(2+) release channels is affected, whether changes in function are cause or consequence and which post-translational modifications drive disease progression. METHODS AND RESULTS:Electrophysiological, imaging, and biochemical techniques were used to study RyRs in cardiomyocytes from mdx mice, an animal model of DMD. Young mdx mice show no changes in cardiac performance, but do so after ∼8 months. Nevertheless, myocytes from mdx pups exhibited exaggerated Ca(2+) responses to mechanical stress and 'hypersensitive' excitation-contraction coupling, hallmarks of increased RyR Ca(2+) sensitivity. Both were normalized by antioxidants, inhibitors of NAD(P)H oxidase and CaMKII, but not by NO synthases and PKA antagonists. Sarcoplasmic reticulum Ca(2+) load and leak were unchanged in young mdx mice. However, by the age of 4-5 months and in senescence, leak was increased and load was reduced, indicating disease progression. By this age, all pharmacological interventions listed above normalized Ca(2+) signals and corrected changes in ECC, Ca(2+) load, and leak. CONCLUSION:Our findings suggest that increased RyR Ca(2+) sensitivity precedes and presumably drives the progression of dystrophic cardiomyopathy, with oxidative stress initiating its development. RyR oxidation followed by phosphorylation, first by CaMKII and later by PKA, synergistically contributes to cardiac deterioration.

Fully Automatic Segmentation of AP Pelvis X-rays via Random Forest Regression and Hierarchical Sparse Shape Composition

Knowledge of landmarks and contours in anteroposterior (AP) pelvis X-rays is invaluable for computer aided diagnosis, hip surgery planning and image-guided interventions. This paper presents a fully automatic and robust approach for landmarking and segmentation of both pelvis and femur in a conventional AP X-ray. Our approach is based on random forest regression and hierarchical sparse shape composition. Experiments conducted on 436 clinical AP pelvis x-rays show that our approach achieves an average point-to-curve error around 1.3 mm for femur and 2.2 mm for pelvis, both with success rates around 98%. Compared to existing methods, our approach exhibits better performance in both the robustness and the accuracy.

Therapeutic drug monitoring of once daily aminoglycoside dosing: comparison of two methods and investigation of the optimal blood sampling strategy.

PURPOSE Therapeutic drug monitoring of patients receiving once daily aminoglycoside therapy can be performed using pharmacokinetic (PK) formulas or Bayesian calculations. While these methods produced comparable results, their performance has never been checked against full PK profiles. We performed a PK study in order to compare both methods and to determine the best time-points to estimate AUC0-24 and peak concentrations (C max). METHODS We obtained full PK profiles in 14 patients receiving a once daily aminoglycoside therapy. PK parameters were calculated with PKSolver using non-compartmental methods. The calculated PK parameters were then compared with parameters estimated using an algorithm based on two serum concentrations (two-point method) or the software TCIWorks (Bayesian method). RESULTS For tobramycin and gentamicin, AUC0-24 and C max could be reliably estimated using a first serum concentration obtained at 1 h and a second one between 8 and 10 h after start of the infusion. The two-point and the Bayesian method produced similar results. For amikacin, AUC0-24 could reliably be estimated by both methods. C max was underestimated by 10-20% by the two-point method and by up to 30% with a large variation by the Bayesian method. CONCLUSIONS The ideal time-points for therapeutic drug monitoring of once daily administered aminoglycosides are 1 h after start of a 30-min infusion for the first time-point and 8-10 h after start of the infusion for the second time-point. Duration of the infusion and accurate registration of the time-points of blood drawing are essential for obtaining precise predictions.

Consensus on the Need for a Hierarchical List of Patient-reported Pain Outcomes for Metaanalyses of Knee Osteoarthritis Trials: An OMERACT Objective.

OBJECTIVE Although protocol registration for systematic reviews is still not mandatory, reviewers should be strongly encouraged to register the protocol to identify the methodological approach, including all outcomes of interest. This will minimize the likelihood of biased decisions in reviews, such as selective outcome reporting. A group of international experts convened to address issues regarding the need to develop hierarchical lists of outcome measurement instruments for a particular outcome for metaanalyses. METHODS Multiple outcome measurement instruments exist to measure the same outcome. Metaanalysis of knee osteoarthritis (OA) trials, and the assessment of pain as an outcome, was used as an exemplar to assess how Outcome Measures in Rheumatology (OMERACT), the Cochrane Collaboration, and other international initiatives might contribute in this area. The meeting began with formal presentations of background topics, empirical evidence from the literature, and a brief introduction to 2 existing hierarchical lists of pain outcome measurement instruments recommended for metaanalyses of knee OA trials. RESULTS After discussions, most participants agreed that there is a need to develop a methodology for generation of hierarchical lists of outcome measurement instruments to guide metaanalyses. Tools that could be used to steer development of such a prioritized list are the COSMIN checklist (Consensus-based Standards for the selection of health status Measurement Instruments) and the OMERACT Filter 2.0. CONCLUSION We list meta-epidemiological research agenda items that address the frequency of reported outcomes in trials, as well as methodologies to assess the best measurement properties (i.e., truth, discrimination, and feasibility).