Methods for comparing center‐specific survival outcomes using direct standardization

The evaluation of center‐specific outcomes is often through survival analysis methods. Such evaluations must account for differences in the distribution of patient characteristics across centers. In the context of censored event times, it is also important that the measure chosen to evaluate centers not be influenced by imbalances in the center‐specific censoring distributions. The practice of using center indicators in a hazard regression model is often invalid, inconvenient, or undesirable to carry out. We propose a semiparametric version of the standardized rate ratio (SRR) useful for the evaluation of centers with respect to a right‐censored event time. The SRR for center j can be interpreted as the ratio of the expected number of deaths in the total population (if the total population were in fact subject to the center j mortality hazard) to the observed number of events. The proposed measure is not affected by differences in center‐specific covariate or censoring distributions. Asymptotic properties of the proposed estimators are derived, with finite‐sample properties examined through simulation studies. The proposed methods are applied to national kidney transplant data. Copyright © 2014 John Wiley & Sons, Ltd.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/106887/1/sim6089.pd

Bayesian modeling and inference for clinical trials with partial retrieved data following dropout

In randomized clinical trials, it is common that patients may stop taking their assigned treatments and then switch to a standard treatment (standard of care available to the patient) but not the treatments under investigation. While the availability of limited retrieved data on patients who switch to standard treatment, called off-protocol data, could be highly valuable in assessing the associated treatment effect with the experimental therapy, it leads to a complex data structure requiring the development of models that link the information of per-protocol data with the off-protocol data. In this paper, we develop a novel Bayesian method to jointly model longitudinal treatment measurements under various dropout scenarios. Specifically, we propose a multivariate normal mixed-effects model for repeated measurements from the assigned treatments and the standard treatment, a multivariate logistic regression model for those stopping the assigned treatments, logistic regression models for those starting a standard treatment off protocol, and a conditional multivariate logistic regression model for completely withdrawing from the study. We assume that withdrawing from the study is non-ignorable but intermittent missingness is assumed to be at random. Various properties of the proposed model are examined. An efficient Markov chain Monte Carlo sampling algorithm is developed. A real data set from a clinical trial is analyzed in detail via the proposed method

Conflict diagnostics in directed acyclic graphs, with applications in bayesian evidence synthesis

Complex stochastic models represented by directed acyclic graphs (DAGs) are increasingly employed to synthesise multiple, imperfect and disparate sources of evidence, to estimate quantities that are difficult to measure directly. The various data sources are dependent on shared parameters and hence have the potential to conflict with each other, as well as with the model. In a Bayesian framework, the model consists of three components: the prior distribution, the assumed form of the likelihood and structural assumptions. Any of these components may be incompatible with the observed data. The detection and quantification of such conflict and of data sources that are inconsistent with each other is therefore a crucial component of the model criticism process. We first review Bayesian model criticism, with a focus on conflict detection, before describing a general diagnostic for detecting and quantifying conflict between the evidence in different partitions of a DAG. The diagnostic is a p-value based on splitting the information contributing to inference about a "separator" node or group of nodes into two independent groups and testing whether the two groups result in the same inference about the separator node(s). We illustrate the method with three comprehensive examples: an evidence synthesis to estimate HIV prevalence; an evidence synthesis to estimate influenza case-severity; and a hierarchical growth model for rat weights.This work was supported by the Medical Research Council [Unit Programme Numbers U105260566 and U105260557

A semi-parametric approach to estimate risk functions associated with multi-dimensional exposure profiles: application to smoking and lung cancer

A common characteristic of environmental epidemiology is the multi-dimensional aspect of exposure patterns, frequently reduced to a cumulative exposure for simplicity of analysis. By adopting a flexible Bayesian clustering approach, we explore the risk function linking exposure history to disease. This approach is applied here to study the relationship between different smoking characteristics and lung cancer in the framework of a population based case control study

Variable Selection in Covariate Dependent Random Partition Models: an Application to Urinary Tract Infection

Lower urinary tract symptoms can indicate the presence of urinary tract infection (UTI), a condition that if it becomes chronic requires expensive and time consuming care as well as leading to reduced quality of life. Detecting the presence and gravity of an infection from the earliest symptoms is then highly valuable. Typically, white blood cell (WBC) count measured in a sample of urine is used to assess UTI. We consider clinical data from 1341 patients in their first visit in which UTI (i.e. WBC ≥ 1) is diagnosed. In addition, for each patient, a clinical profile of 34 symptoms was recorded. In this paper, we propose a Bayesian nonparametric regression model based on the Dirichlet process prior aimed at providing the clinicians with a meaningful clustering of the patients based on both the WBC (response variable) and possible patterns within the symptoms profiles (covariates). This is achieved by assuming a probability model for the symptoms as well as for the response variable. To identify the symptoms most associated to UTI, we specify a spike and slab base measure for the regression coefficients: this induces dependence of symptoms selection on cluster assignment. Posterior inference is performed through Markov Chain Monte Carlo methods

A hierarchical Bayesian approach to multiple testing in disease mapping

We propose a Bayesian approach to multiple testing in disease mapping. This study was motivated by a real example regarding the mortality rate for lung cancer, males, in the Tuscan region (Italy). The data are relative to the period 1995–1999 for 287 municipalities. We develop a tri-level hierarchical Bayesian model to estimate for each area the posterior classification probability that is the posterior probability that the municipality belongs to the set of non-divergent areas. We show also the connections of our model with the false discovery rate approach. Posterior classification probabilities are used to explore areas at divergent risk from the reference while controlling for multiple testing. We consider both the Poisson-Gamma and the Besag, York and Mollié model to account for extra Poisson variability in our Bayesian formulation. Posterior inference on classification probabilities is highly dependent on the choice of the prior. We perform a sensitivity analysis and suggest how to rely on subject-specific information to derive informative a priori distributions. Hierarchical Bayesian models provide a sensible way to model classification probabilities in the context of disease mapping

Using funnel plots in public health surveillance

<p>Abstract</p> <p>Background</p> <p>Public health surveillance is often concerned with the analysis of health outcomes over small areas. Funnel plots have been proposed as a useful tool for assessing and visualizing surveillance data, but their full utility has not been appreciated (for example, in the incorporation and interpretation of risk factors).</p> <p>Methods</p> <p>We investigate a way to simultaneously focus funnel plot analyses on direct policy implications while visually incorporating model fit and the effects of risk factors. Health survey data representing modifiable and nonmodifiable risk factors are used in an analysis of 2007 small area motor vehicle mortality rates in Alberta, Canada.</p> <p>Results</p> <p>Small area variations in motor vehicle mortality in Alberta were well explained by the suite of modifiable and nonmodifiable risk factors. Funnel plots of raw rates and of risk adjusted rates lead to different conclusions; the analysis process highlights opportunities for intervention as risk factors are incorporated into the model. Maps based on funnel plot methods identify areas worthy of further investigation.</p> <p>Conclusions</p> <p>Funnel plots provide a useful tool to explore small area data and to routinely incorporate covariate relationships in surveillance analyses. The exploratory process has at each step a direct and useful policy-related result. Dealing thoughtfully with statistical overdispersion is a cornerstone to fully understanding funnel plots.</p

The natural history of <i>Chlamydia trachomatis </i>infection in women:a multi-parameter evidence synthesis

Background and objectives: The evidence base supporting the National Chlamydia Screening Programme, initiated in 2003, has been questioned repeatedly, with little consensus on modelling assumptions, parameter values or evidence sources to be used in cost-effectiveness analyses. The purpose of this project was to assemble all available evidence on the prevalence and incidence of Chlamydia trachomatis (CT) in the UK and its sequelae, pelvic inflammatory disease (PID), ectopic pregnancy (EP) and tubal factor infertility (TFI) to review the evidence base in its entirety, assess its consistency and, if possible, arrive at a coherent set of estimates consistent with all the evidence. Methods: Evidence was identified using ‘high-yield’ strategies. Bayesian Multi-Parameter Evidence Synthesis models were constructed for separate subparts of the clinical and population epidemiology of CT. Where possible, different types of data sources were statistically combined to derive coherent estimates. Where evidence was inconsistent, evidence sources were re-interpreted and new estimates derived on a post-hoc basis. Results: An internally coherent set of estimates was generated, consistent with a multifaceted evidence base, fertility surveys and routine UK statistics on PID and EP. Among the key findings were that the risk of PID (symptomatic or asymptomatic) following an untreated CT infection is 17.1% [95% credible interval (CrI) 6% to 29%] and the risk of salpingitis is 7.3% (95% CrI 2.2% to 14.0%). In women aged 16–24 years, screened at annual intervals, at best, 61% (95% CrI 55% to 67%) of CT-related PID and 22% (95% CrI 7% to 43%) of all PID could be directly prevented. For women aged 16–44 years, the proportions of PID, EP and TFI that are attributable to CT are estimated to be 20% (95% CrI 6% to 38%), 4.9% (95% CrI 1.2% to 12%) and 29% (95% CrI 9% to 56%), respectively. The prevalence of TFI in the UK in women at the end of their reproductive lives is 1.1%: this is consistent with all PID carrying a relatively high risk of reproductive damage, whether diagnosed or not. Every 1000 CT infections in women aged 16–44 years, on average, gives rise to approximately 171 episodes of PID and 73 of salpingitis, 2.0 EPs and 5.1 women with TFI at age 44 years. Conclusions and research recommendations: The study establishes a set of interpretations of the major studies and study designs, under which a coherent set of estimates can be generated. CT is a significant cause of PID and TFI. CT screening is of benefit to the individual, but detection and treatment of incident infection may be more beneficial. Women with lower abdominal pain need better advice on when to seek early medical attention to avoid risk of reproductive damage. The study provides new insights into the reproductive risks of PID and the role of CT. Further research is required on the proportions of PID, EP and TFI attributable to CT to confirm predictions made in this report, and to improve the precision of key estimates. The cost-effectiveness of screening should be re-evaluated using the findings of this report. Funding: The Medical Research Council grant G0801947