Objective Bayesian meta-analysis for sparse discrete data

Bayesian heterogeneity in a meta-analysis with two studies and binary data

The meta-analysis of two trials is valuable in many practical situations, such as studies of rare and/or orphan diseases focussed on a single intervention. In this context, additional concerns, like small sample size and/or heterogeneity in the results obtained, might make standard frequentist and Bayesian techniques inappropriate. In a meta-analysis, moreover, the presence of between-sample heterogeneity adds model uncertainty, which must be taken into consideration when drawing inferences. We suggest that the most appropriate way to measure this heterogeneity is by clustering the samples and then determining the posterior probability of the cluster models. The meta-inference is obtained as a mixture of all the meta-inferences for the cluster models, where the mixing distribution is the posterior model probability. We present a simple two-component form of Bayesian model averaging that is unaffected by characteristics such as small study size or zero-cell counts, and which is capable of incorporating uncertainties into the estimation process. Illustrative examples are given and analysed, using real sparse binomial data.

Bayesian Meta-Analysis for Binary Data and Prior Distribution on Models

In meta-analysis, the structure of the between-sample heterogeneity plays a crucial role in estimating the meta-parameter. A Bayesian meta-analysis for binary data has recently been proposed that measures this heterogeneity by clustering the samples and then determining the posterior probability of the cluster models through model selection. The meta-parameter is then estimated using Bayesian model averaging techniques. Although an objective Bayesian meta-analysis is proposed for each type of heterogeneity, we concentrate the attention of this paper on priors over the models. We consider four alternative priors which are motivated by reasonable but different assumptions. A frequentist validation with simulated data has been carried out to analyze the properties of each prior distribution for a set of different number of studies and sample sizes. The results show the importance of choosing an adequate model prior as the posterior probabilities for the models are very sensitive to it. The hierarchical Poisson prior and the hierarchical uniform prior show a good performance when the real model is the homogeneity, or when the sample sizes are high enough. However, the uniform prior can detect the true model when it is an intermediate model (neither homogeneity nor heterogeneity) even for small sample sizes and few studies. An illustrative example with real data is also given, showing the sensitivity of the estimation of the meta-parameter to the model prior.

Meta-Analysis with Few Studies and Binary Data: A Bayesian Model Averaging Approach

In meta-analysis, the existence of between-sample heterogeneity introduces model uncertainty, which must be incorporated into the inference. We argue that an alternative way to measure this heterogeneity is by clustering the samples and then determining the posterior probability of the cluster models. The meta-inference is obtained as a mixture of all the meta-inferences for the cluster models, where the mixing distribution is the posterior model probabilities. When there are few studies, the number of cluster configurations is manageable, and the meta-inferences can be drawn with BMA techniques. Although this topic has been relatively neglected in the meta-analysis literature, the inference thus obtained accurately reflects the cluster structure of the samples used. In this paper, illustrative examples are given and analysed, using real binary data.

Serum Calcium Levels and Parkinson's Disease: A Mendelian Randomization Study

BACKGROUND

Though increasing epidemiological studies have evaluated the correlation between serum calcium contents and Parkinson's disease (PD), the results are inconsistent. At present, whether there is a causal association between serum calcium content and PD remains undetermined.

OBJECTIVE AND METHODS

This study was designed to explore the relationship between increased serum calcium contents and PD risk. In this present study, a Mendelian randomization trial was carried out using a large-scale serum calcium genome-wide association study (GWAS) dataset (N = 61,079, Europeans) and a large-scale PD GWAS dataset (N = 8,477, Europeans including 4,238 PD patients and 4,239 controls). Here, a total of four Mendelian randomization methods comprising weighted median, inverse-variance weighted meta-analysis (IVW), MR-Egger, and MR-PRESSO were used.

RESULTS

Our data concluded that genetically higher serum calcium contents were not significantly related to PD.

CONCLUSION

In conclusion, we provided genetic evidence that there was no direct causal relationship between serum calcium contents and PD. Hence, calcium supplementation may not result in reduced PD risk.

Prior distributions for variance parameters in a sparse-event meta-analysis of a few small trials

In rare diseases, typically only a small number of patients are available for a randomized clinical trial. Nevertheless, it is not uncommon that more than one study is performed to evaluate a (new) treatment. Scarcity of available evidence makes it particularly valuable to pool the data in a meta-analysis. When the primary outcome is binary, the small sample sizes increase the chance of observing zero events. The frequentist random-effects model is known to induce bias and to result in improper interval estimation of the overall treatment effect in a meta-analysis with zero events. Bayesian hierarchical modeling could be a promising alternative. Bayesian models are known for being sensitive to the choice of prior distributions for between-study variance (heterogeneity) in sparse settings. In a rare disease setting, only limited data will be available to base the prior on, therefore, robustness of estimation is desirable. We performed an extensive and diverse simulation study, aiming to provide practitioners with advice on the choice of a sufficiently robust prior distribution shape for the heterogeneity parameter. Our results show that priors that place some concentrated mass on small τ values but do not restrict the density for example, the Uniform(-10, 10) heterogeneity prior on the log(τ² ) scale, show robust 95% coverage combined with less overestimation of the overall treatment effect, across varying degrees of heterogeneity. We illustrate the results with meta-analyzes of a few small trials.

Bayesian meta-analysis: The role of the between-sample heterogeneity

The random effect approach for meta-analysis was motivated by a lack of consistent assessment of homogeneity of treatment effect before pooling. The random effect model assumes that the distribution of the treatment effect is fully heterogenous across the experiments. However, other models arising by grouping some of the experiments are plausible. We illustrate on simulated binary experiments that the fully heterogenous model gives a poor meta-inference when fully heterogeneity is not the true model and that the knowledge of the true cluster model considerably improves the inference. We propose the use of a Bayesian model selection procedure for estimating the true cluster model, and Bayesian model averaging to incorporate into the meta-analysis the clustering estimation. A well-known meta-analysis for six major multicentre trials to assess the efficacy of a given dose of aspirin in post-myocardial infarction patients is reanalysed.

An objective Bayesian analysis of a crossover design via model selection and model averaging

Inference about the treatment effect in a crossover design has received much attention over time owing to the uncertainty in the existence of the carryover effect and its impact on the estimation of the treatment effect. Adding to this uncertainty is that the existence of the carryover effect and its size may depend on the presence of the treatment effect and its size. We consider estimation and testing hypothesis about the treatment effect in a two-period crossover design, assuming normally distributed response variable, and use an objective Bayesian approach to test the hypothesis about the treatment effect and to estimate its size when it exists while accounting for the uncertainty about the presence of the carryover effect as well as the treatment and period effects. We evaluate and compare the performance of the proposed approach with a standard frequentist approach using simulated data, and real data. Copyright © 2016 John Wiley & Sons, Ltd.

A Bayesian sensitivity study of risk difference in the meta-analysis of binary outcomes from sparse data

In most cases, including those of discrete random variables, statistical meta-analysis is carried out using the normal random effect model. The authors argue that normal approximation does not always properly reflect the underlying uncertainty of the original discrete data. Furthermore, in the presence of rare events the results from this approximation can be very poor. This review proposes a Bayesian meta-analysis to address binary outcomes from sparse data and also introduces a simple way to examine the sensitivity of the quantities of interest in the meta-analysis with respect to the structure dependence selected. The findings suggest that for binary outcomes data it is possible to develop a Bayesian procedure, which can be directly applied to sparse data without ad hoc corrections. By choosing a suitable class of linking distributions, the authors found that a Bayesian robustness study can be easily implemented. For illustrative purposes, an example with real data is analyzed using the proposed Bayesian meta-analysis for binomial sparse data.