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.

Valgus Hip With High Antetorsion Causes Pain Through Posterior Extraarticular FAI

BACKGROUND: Valgus hips with increased antetorsion present with lack of external rotation and posterior hip pain that is aggravated with hip extension and external rotation. This may be the result of posterior femoroacetabular impingement (FAI). QUESTIONS/PURPOSES: We asked whether (1) the range of motion (ROM); (2) the location of anterior and posterior bony collision zones; and (3) the prevalence of extraarticular impingement differ between valgus hips with increased antetorsion compared with normal hips and hips with idiopathic FAI. METHODS: Surface models based on CT scan reconstructions of 13 valgus hips with increased antetorsion, 22 hips with FAI, and 27 normal hips were included. Validated three-dimensional collision detection software was used to quantify the simulated hip ROM and the location of impingement on the acetabular and the femoral sides. RESULTS: Hips with coxa valga and antetorsion showed decreased extension, external rotation, and adduction, whereas internal rotation in 90° of flexion was increased. Impingement zones were more anteroinferior on the femur and posteroinferior on the acetabular (pelvic) side; and the zones were more frequently extraarticular, posterior, or to a lesser degree anterior against the inferior iliac spine. We found a higher prevalence of extraarticular impingement for valgus hips with increased antetorsion. CONCLUSIONS: Valgus hips with increased antetorsion predispose to posterior extraarticular FAI and to a lesser degree anteroinferior spine impingement. LEVEL OF EVIDENCE: Level II, prognostic study. See Guidelines for Authors for a complete description of levels of evidence.

Internal Models, Vestibular Cognition, and Mental Imagery: Conceptual Considerations

Vestibular cognition has recently gained attention. Despite numerous experimental and clinical demonstrations, it is not yet clear what vestibular cognition really is. For future research in vestibular cognition, adopting a computational approach will make it easier to explore the underlying mech- anisms. Indeed, most modeling approaches in vestibular science include a top-down or a priori component. We review recent Bayesian optimal observer models, and discuss in detail the conceptual value of prior assumptions, likelihood and posterior estimates for research in vestibular cognition. We then consider forward models in vestibular processing, which are required in order to distinguish between sensory input that is induced by active self-motion, and sensory input that is due to passive self-motion. We suggest that forward models are used not only in the service of estimating sensory states but they can also be drawn upon in an offline mode (e.g., spatial perspective transformations), in which interaction with sensory input is not desired. A computational approach to vestibular cogni- tion will help to discover connections across studies, and it will provide a more coherent framework for investigating vestibular cognition.

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.

Predicting the population impact of chlamydia screening programmes: comparative mathematical modelling study

BACKGROUND: Published individual-based, dynamic sexual network modelling studies reach different conclusions about the population impact of screening for Chlamydia trachomatis. The objective of this study was to conduct a direct comparison of the effect of organised chlamydia screening in different models. METHODS: Three models simulating population-level sexual behaviour, chlamydia transmission, screening and partner notification were used. Parameters describing a hypothetical annual opportunistic screening program in 16-24 year olds were standardised, whereas other parameters from the three original studies were retained. Model predictions of the change in chlamydia prevalence were compared under a range of scenarios. RESULTS: Initial overall chlamydia prevalence rates were similar in women but not men and there were age and sex-specific differences between models. The number of screening tests carried out was comparable in all models but there were large differences in the predicted impact of screening. After 10 years of screening, the predicted reduction in chlamydia prevalence in women aged 16-44 years ranged from 4% to 85%. Screening men and women had a greater impact than screening women alone in all models. There were marked differences between models in assumptions about treatment seeking and sexual behaviour before the start of the screening intervention. CONCLUSIONS: Future models of chlamydia transmission should be fitted to both incidence and prevalence data. This meta-modelling study provides essential information for explaining differences between published studies and increasing the utility of individual-based chlamydia transmission models for policy making.

Analysis of Phenotypic Variation in Childhood Wheezing Disorders

Recurrent wheezing or asthma is a common problem in children that has increased considerably in prevalence in the past few decades. The causes and underlying mechanisms are poorly understood and it is thought that a numb er of distinct diseases causing similar symptoms are involved. Due to the lack of a biologically founded classification system, children are classified according to their observed disease related features (symptoms, signs, measurements) into phenotypes. The objectives of this PhD project were a) to develop tools for analysing phenotypic variation of a disease, and b) to examine phenotypic variability of wheezing among children by applying these tools to existing epidemiological data. A combination of graphical methods (multivariate co rrespondence analysis) and statistical models (latent variables models) was used. In a first phase, a model for discrete variability (latent class model) was applied to data on symptoms and measurements from an epidemiological study to identify distinct phenotypes of wheezing. In a second phase, the modelling framework was expanded to include continuous variability (e.g. along a severity gradient) and combinations of discrete and continuo us variability (factor models and factor mixture models). The third phase focused on validating the methods using simulation studies. The main body of this thesis consists of 5 articles (3 published, 1 submitted and 1 to be submitted) including applications, methodological contributions and a review. The main findings and contributions were: 1) The application of a latent class model to epidemiological data (symptoms and physiological measurements) yielded plausible pheno types of wheezing with distinguishing characteristics that have previously been used as phenotype defining characteristics. 2) A method was proposed for including responses to conditional questions (e.g. questions on severity or triggers of wheezing are asked only to children with wheeze) in multivariate modelling.ii 3) A panel of clinicians was set up to agree on a plausible model for wheezing diseases. The model can be used to generate datasets for testing the modelling approach. 4) A critical review of methods for defining and validating phenotypes of wheeze in children was conducted. 5) The simulation studies showed that a parsimonious parameterisation of the models is required to identify the true underlying structure of the data. The developed approach can deal with some challenges of real-life cohort data such as variables of mixed mode (continuous and categorical), missing data and conditional questions. If carefully applied, the approach can be used to identify whether the underlying phenotypic variation is discrete (classes), continuous (factors) or a combination of these. These methods could help improve precision of research into causes and mechanisms and contribute to the development of a new classification of wheezing disorders in children and other diseases which are difficult to classify.

Effectiveness of non-steroidal anti-inflammatory drugs for the treatment of pain in knee and hip osteoarthritis: a network meta-analysis.

BACKGROUND Non-steroidal anti-inflammatory drugs (NSAIDs) are the backbone of osteoarthritis pain management. We aimed to assess the effectiveness of different preparations and doses of NSAIDs on osteoarthritis pain in a network meta-analysis. METHODS For this network meta-analysis, we considered randomised trials comparing any of the following interventions: NSAIDs, paracetamol, or placebo, for the treatment of osteoarthritis pain. We searched the Cochrane Central Register of Controlled Trials (CENTRAL) and the reference lists of relevant articles for trials published between Jan 1, 1980, and Feb 24, 2015, with at least 100 patients per group. The prespecified primary and secondary outcomes were pain and physical function, and were extracted in duplicate for up to seven timepoints after the start of treatment. We used an extension of multivariable Bayesian random effects models for mixed multiple treatment comparisons with a random effect at the level of trials. For the primary analysis, a random walk of first order was used to account for multiple follow-up outcome data within a trial. Preparations that used different total daily dose were considered separately in the analysis. To assess a potential dose-response relation, we used preparation-specific covariates assuming linearity on log relative dose. FINDINGS We identified 8973 manuscripts from our search, of which 74 randomised trials with a total of 58 556 patients were included in this analysis. 23 nodes concerning seven different NSAIDs or paracetamol with specific daily dose of administration or placebo were considered. All preparations, irrespective of dose, improved point estimates of pain symptoms when compared with placebo. For six interventions (diclofenac 150 mg/day, etoricoxib 30 mg/day, 60 mg/day, and 90 mg/day, and rofecoxib 25 mg/day and 50 mg/day), the probability that the difference to placebo is at or below a prespecified minimum clinically important effect for pain reduction (effect size [ES] -0·37) was at least 95%. Among maximally approved daily doses, diclofenac 150 mg/day (ES -0·57, 95% credibility interval [CrI] -0·69 to -0·46) and etoricoxib 60 mg/day (ES -0·58, -0·73 to -0·43) had the highest probability to be the best intervention, both with 100% probability to reach the minimum clinically important difference. Treatment effects increased as drug dose increased, but corresponding tests for a linear dose effect were significant only for celecoxib (p=0·030), diclofenac (p=0·031), and naproxen (p=0·026). We found no evidence that treatment effects varied over the duration of treatment. Model fit was good, and between-trial heterogeneity and inconsistency were low in all analyses. All trials were deemed to have a low risk of bias for blinding of patients. Effect estimates did not change in sensitivity analyses with two additional statistical models and accounting for methodological quality criteria in meta-regression analysis. INTERPRETATION On the basis of the available data, we see no role for single-agent paracetamol for the treatment of patients with osteoarthritis irrespective of dose. We provide sound evidence that diclofenac 150 mg/day is the most effective NSAID available at present, in terms of improving both pain and function. Nevertheless, in view of the safety profile of these drugs, physicians need to consider our results together with all known safety information when selecting the preparation and dose for individual patients. FUNDING Swiss National Science Foundation (grant number 405340-104762) and Arco Foundation, Switzerland.

The Face of Power

Objective: In humans and other animals, open, expansive postures (compared to contracted postures) are evolutionary developed expressions of power and have been shown to cause neuroendocrine and behavioral changes (Carney, Cuddy, & Yap, 2010). In the present study we aimed to investigate whether power postures have a bearing on the participant’s facial appearance and whether others are able to distinguish faces after “high power posing” from faces after “low power posing”. Methods: 16 models were photographed 4-5 minutes after having adopted high and low power postures. Two different high power and two different low power postures were held for 2 minutes each. Power-posing sessions were performed on two consecutive days. High and low power photographs of each model were paired and an independent sample of 100 participants were asked to pick the more dominant and the more likeable face of each pair. Results: Photographs that were taken after adopting high power postures were chosen significantly more often as being more dominant looking. There was no preference when asked to choose the more likeable photograph (chance level). A further independent sample rated each photograph for head tilt, making it unlikely that dominance ratings were caused merely by the posture of the head. Consistently, facial width-to-height ratio did not differ between faces after high and low power posing. Conclusions: Postures associated with high power affect facial appearance, leading to a more dominant looking face. This finding may have implications for everyday life, for instance when a dominant appearance is needed.

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.

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.

Management of children with autism spectrum disorder in the dental setting: concerns, behavioural approaches and recommendations

OBJECTIVES This article reviews the present literature on the issues encountered while coping with children with autistic spectrum disorder from the dental perspective. The autistic patient profile and external factors affecting the oral health status of this patient population are discussed upon the existing body of evidence. MATERIAL AND METHODS The MEDLINE database was searched using the terms 'Autistic Disorder', 'Behaviour Control/methods', 'Child', 'Dental care for disabled', 'Education', 'Oral Health', and 'Pediatric Dentistry' to locate related articles published up to January 2013. RESULTS Most of the relevant studies indicate poor oral hygiene whereas they are inconclusive regarding the caries incidence in autistic individuals. Undergraduate dental education appears to determine the competence of dental professionals to treat developmentally disabled children and account partly for compromised access to dental care. Dental management of an autistic child requires in-depth understanding of the background of the autism and available behavioural guidance theories. The dental professional should be flexible to modify the treatment approach according to the individual patient needs.

Expectation Conditional Maximization-Based Deformable Shape Registration

This paper addresses the issue of matching statistical and non-rigid shapes, and introduces an Expectation Conditional Maximization-based deformable shape registration (ECM-DSR) algorithm. Similar to previous works, we cast the statistical and non-rigid shape registration problem into a missing data framework and handle the unknown correspondences with Gaussian Mixture Models (GMM). The registration problem is then solved by fitting the GMM centroids to the data. But unlike previous works where equal isotropic covariances are used, our new algorithm uses heteroscedastic covariances whose values are iteratively estimated from the data. A previously introduced virtual observation concept is adopted here to simplify the estimation of the registration parameters. Based on this concept, we derive closed-form solutions to estimate parameters for statistical or non-rigid shape registrations in each iteration. Our experiments conducted on synthesized and real data demonstrate that the ECM-DSR algorithm has various advantages over existing algorithms.

Predicting smoking cessation and its relapse in HIV-infected patients: the Swiss HIV Cohort Study

OBJECTIVES: The aim of the study was to assess whether prospective follow-up data within the Swiss HIV Cohort Study can be used to predict patients who stop smoking; or among smokers who stop, those who start smoking again. METHODS: We built prediction models first using clinical reasoning ('clinical models') and then by selecting from numerous candidate predictors using advanced statistical methods ('statistical models'). Our clinical models were based on literature that suggests that motivation drives smoking cessation, while dependence drives relapse in those attempting to stop. Our statistical models were based on automatic variable selection using additive logistic regression with component-wise gradient boosting. RESULTS: Of 4833 smokers, 26% stopped smoking, at least temporarily; because among those who stopped, 48% started smoking again. The predictive performance of our clinical and statistical models was modest. A basic clinical model for cessation, with patients classified into three motivational groups, was nearly as discriminatory as a constrained statistical model with just the most important predictors (the ratio of nonsmoking visits to total visits, alcohol or drug dependence, psychiatric comorbidities, recent hospitalization and age). A basic clinical model for relapse, based on the maximum number of cigarettes per day prior to stopping, was not as discriminatory as a constrained statistical model with just the ratio of nonsmoking visits to total visits. CONCLUSIONS: Predicting smoking cessation and relapse is difficult, so that simple models are nearly as discriminatory as complex ones. Patients with a history of attempting to stop and those known to have stopped recently are the best candidates for an intervention.