Email updates

Keep up to date with the latest news and content from Trials and BioMed Central.

Open Access Methodology

A statistical framework for quantifying clinical equipoise for individual cases during randomized controlled surgical trials

Nicholas R Parsons1*, Yuri Kulikov1, Alan Girling2 and Damian Griffin1

Author Affiliations

1 Warwick Medical School, University of Warwick, Coventry, CV4 7AL, UK

2 Department of Public Health, Epidemiology and Biostatistics, University of Birmingham, Birmingham, B15 2TT, UK

For all author emails, please log on.

Trials 2011, 12:258  doi:10.1186/1745-6215-12-258


The electronic version of this article is the complete one and can be found online at: http://www.trialsjournal.com/content/12/1/258


Received:25 July 2011
Accepted:13 December 2011
Published:13 December 2011

© 2011 Parsons et al; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

Randomised controlled trials are being increasingly used to evaluate new surgical interventions. There are a number of problematic methodological issues specific to surgical trials, the most important being identifying whether patients are eligible for recruitment into the trial. This is in part due to the diversity in practice patterns across institutions and the enormous range of available interventions that often leads to a low level of agreement between clinicians about both the value and the appropriate choice of intervention. We argue that a clinician should offer patients the option of recruitment into a trial, even if the clinician is not individually in a position of equipoise, if there is collective (clinical) equipoise amongst the wider clinical community about the effectiveness of a proposed intervention (the clinical equipoise principle). We show how this process can work using data collected from an ongoing trial of a surgical intervention.

Results

We describe a statistical framework for the assessment of uncertainty prior to patient recruitment to a clinical trial using a panel of expert clinical assessors and techniques for eliciting, pooling and modelling of expert opinions. The methodology is illustrated using example data from the UK Heel Fracture Trial. The statistical modelling provided results that were clear and simple to present to clinicians and showed how decisions regarding recruitment were influenced by both the collective opinion of the expert panel and the type of decision rule selected.

Conclusions

The statistical framework presented has potential to identify eligible patients and assist in the simplification of eligibility criteria which might encourage greater participation in clinical trials evaluating surgical interventions.

Keywords:
Equipoise; Randomised controlled trial; Surgery; Statistical model

1 Background

There is an increasing demand for randomised controlled trials (RCTs) in surgery to provide high quality evaluation of new interventions; we use the word intervention synonymously with treatment, procedure or surgical procedure. In a background of ever evolving and improving healthcare, differences between interventions for the same condition are often small, substantially increasing the risk of biased estimation of treatment effects in simple (uncontrolled) observational studies of the interventions [1]. The need for the kind of high level evidence provided by RCTs for surgical interventions is clear [2], although a number of methodological issues have been raised for surgical trials [1,3]. One of the most important issues being recruitment, and specifically identifying whether patients are eligible for entry into a trial.

The existing tremendous diversity in practice patterns across institutions coupled with an ever increasing range of available interventions suggests a low level of agreement between clinicians about both the value of many interventions and the appropriate choice of intervention [4]. A present or imminent controversy in the expert medical community about a choice between interventions is called clinical (or collective) equipoise. Equipoise is the point where we are equally poised in our beliefs about the potential benefits of a particular intervention [5]; i.e. is intervention A better than intervention B. Clinical equipoise is present "if there is genuine uncertainty within the expert medical community - not necessarily on the part of the individual investigator - about the preferred treatment" [5]. In many cases the only way to resolve collective uncertainty about the optimum intervention choice is to undertake a clinical trial. Individual equipoise relates to a single clinician, i.e. the position where he or she has no preference amongst a range of available treatments. It is subject to change for a host of reasons, including peer pressure, the results of potentially imperfect studies and the influence of advertising. Freedman [5] argues that global clinical equipoise should override the individual clinician's lack of equipoise. Clinicians should subsume their personal views and recruit patients into a trial, even if not individually in a position of equipoise themselves. This view is implicitly accepted by society in the form of ethics committees, which must ensure that the treatments being compared are reasonable options before trial participants are sought. Often, for a treatment that is not completely novel, this is demonstrated by the presence of clinical equipoise in an expert and/or wider medical community. Once ethics committee permission has been granted, it then becomes an individual clinician's decision whether the offer of entry into the trial is appropriate for an individual patient [6]. Unfortunately, the varied preferences expressed (which may be rational, anecdotal or irrational) between individual institutions and between individual surgeons within and between institutions often make patient recruitment to trials very challenging.

Statistically the level of individual uncertainty about the effectiveness of an intervention can be quantified by a (subjective) probability, which is assigned to a specific hypothesis and is personal and varies with an individual's knowledge and expertise. "A measure of a state of knowledge" [7] is provided by the Bayesian concept of subjective probability. The process of expert evaluation about the effectiveness of a proposed intervention in an RCT is synonymous with elicitation of a Bayesian prior; i.e. a statement of knowledge prior to performing an experiment or trial usually stated in the form of a probability density. There are a number of approaches to turning informally expressed ideas into a mathematical prior distribution, with no consensus as to the optimal method of determination for a process that is usually problem specific [8]. We choose to elicit the subjective opinion of a panel of experts as a basis for decision making regarding the eligibility of a patient for recruitment to an RCT [9]. This has the advantage of being dynamic and flexible, in the sense that it is quite feasible that opinions will change during the course of a trial, for example with the publication of related research [8], or as experience accumulates amongst clinicians as to how best to undertake a surgical procedure.

Methods for formal measurement of clinical uncertainty, as a prelude to a clinical trial have been suggested previously [10] and measures of surgeon's equipoise in the setting of surgical trials have also been reported [11]. However, we develop these ideas further, using techniques for eliciting subjective judgements before a trial [12-14] and introduce a novel framework for decision making regarding recruitment to an RCT that we hope will be easily understood by clinicians and implemented in real time during the course of a trial. It is particularly challenging recruiting patients to trials comparing operative to non-operative treatments or a standard against a new but popular well-marketed treatment. Therefore we develop a statistical framework to model clinical equipoise (Section 3), using a parametric and a nonparametric approach, for data collected from a clinical trial comparing conservative and operative treatment for displaced fractures of the calcaneus. The results of applying the models are reported in Sections 3.7-10 and we draw conclusions in Section 4.

2 Methods

Using available web design tools a method was developed to capture the opinions of clinicians in real time for individual patients (cases) in an ongoing RCT. It comprised of a virtual expert panel giving their opinion about the effectiveness of a proposed treatment for individual patients based on online clinical details; the individual assessments were then synthesized and fed back electronically to the lead clinical investigator. This process is described in greater detail below.

Patients who met the initial trial inclusion criteria were identified and approached by a member of the research team to alert them to the possibility of participating in a trial. They were then asked permission for their anonymized clinical details to be distributed among a panel of experts/clinicians for an opinion regarding the effectiveness of the proposed treatment. Clinical data from consented patients were made available on a secure website managed by eLab at the University of Warwick, and all panel experts/clinicians were alerted by email and text message (if requested) to the posting of a new patient and asked to offer their personal opinion on the likely success of the proposed treatments. The assessment scale is described in more detail for the specific example of the UK Heel Fracture Trial. Initially the system was tested in a pilot study with seven surgeons from five UK hospitals. Ten retrospective calcaneal fracture cases were selected to represent typical variability. The surgeons followed the instructions on the website with online and telephone technical support available; no specific training was given. When voting on all ten cases was completed, surgeons were asked to fill in an evaluation questionnaire. Voting on a single case never took longer than 5 minutes and the available clinical information was found sufficient and the whole process user friendly by all participating surgeons.

After the successful pilot study the system was introduced as an independent component of the UK Heel Fracture Trial, which compared conservative and operative treatment for displaced fractures of the calcaneus. The study had separate ethical approval and a consent form, in addition to the main trial. This allowed inclusion both of those patients who took part in the UK Heel Fracture Trial and those who declined, as soon as the patient met the trial eligibility criteria. To avoid interference with the clinical course, patients were asked permission to use their data at the 6 weeks follow-up clinic or later. Their anonymous clinical data including X-rays and CT images were posted to a secure website. The expert assessment panel included 12 surgeons from 9 hospitals. All surgeons were foot and ankle specialists and acted as principal investigators in their individual trial centres.

After assessing the clinical data available for a given patient, the surgeon was able to scroll down to an interactive scale, featuring bars (initially set at zero) above each of seven outcome categories indicating whether after surgical intervention the patient would get "much worse" (1), "significantly worse" (2), "a bit worse" (3), "no difference" (4), "a bit better" (5), "significantly better" (6) or "much better" (7). A left-click of the mouse and a drag allowed each outcome prognosis bar to be set to a desired percentage, which was reported numerically over the bar. Once the assessment summed to 100% (reflected in a digital window in the upper left corner of the scale) the submit button allowed the data to be sent to the trial lead for analysis. The UK Heel Fracture Trial compared operative (surgical) and non-operative (conservative) treatment. Surgical techniques are becoming widespread for calcaneal fracture, but do have associated risks, therefore it was important for the clinician to assess the improvement potential relative to the risks for this procedure. Belief, in the context we describe here, that surgery can make a patient better implies intention to do surgery, while disbelief implies intention to avoid surgical intervention, hence to choose the conservative option. The question posed to the expert panel can and should be tailored to the specific trial. For the UK Heel Fracture Trial the experts were asked to compare operative (surgical) and non-operative (conservative) treatment, which although strongly contrasting treatment options may vary in the exact detail of the constituent components. For studies with less contrasting treatment options (e.g. two types of surgery) the question to experts may simply be whether the test intervention would be better or worse for a patient, compared to a control (standard) intervention.

Table 1 shows four examples of data elicited from between 4 and 6 clinical experts, not necessarily the same individuals labelled as 1 to 6, who provided their opinions on the effectiveness of surgical compared to non-surgical intervention after fracture of the calcaneus. As expected there are clear differences in the both the locations and shapes of the individual distributions for a number of these cases and indeed a number of clear similarities for other cases. For instance, the opinions of the clinicians vary widely for case 1; clinical expert 3 is reasonably confident that the patient will improve significantly after treatment whereas for expert 4 the most likely outcome of treatment is that the condition of the patient will be unchanged. There is much clearer agreement for case 4 where all the experts expect the patient to worsen significantly after treatment. How do we use these data to decide whether a patient (case) is eligible for recruitment to a clinical trial? We propose two approaches here to model the opinions obtained from each expert clinician, a parametric model based on a Beta distribution (Section 3.2) and a nonparametric model based on estimated means and standard deviations (Section 3.3) that characterise expert opinions using concepts of belief, disbelief and uncertainty. The belief, disbelief and uncertainty are visualized using a ternary plot that displays these characteristics in a manner that allows them to be compared to decision rules that partition the opinion space. Finally, resampling methods are used to draw inferences concerning the sufficiency of evidence from the clinical experts to patient eligibility for recruitment

Table 1. Assessment of the likely effectiveness of surgical intervention after fracture of the calcaneus for four example cases and up to six clinical experts

3 Results and Discussion

3.1 Expert opinion

An opinion regarding the effectiveness of a procedure can be thought of as comprising of three distinctive aspects; belief, disbelief and uncertainty. Belief represents the tendency of an expert to expect a particular treatment to perform better than an alternative (control intervention) for a particular patient (case); i.e. the tendency for the experts to score cases in the higher end categories of the rating scale of Table 1. Conversely, the level of disbelief is equated with the tendency for an intervention to have a worse outcome as compared to a control intervention; i.e. the tendency for the experts to score cases in the lower end categories of the rating scale. The uncertainty associated with the belief and disbelief represents the spread of the data across the opinion range; i.e. all the scores might be concentrated in the central category (no difference) or be spread equally between all categories in Table 1 - we would have equal belief in these two scenarios but a maximum difference in uncertainty.

Borrowing from the notation of subjective logic [15,16], we label the belief, disbelief and uncertainty associated with an opinion for expert i as bi, di and ui, and apply the constraint that

<a onClick="popup('http://www.trialsjournal.com/content/12/1/258/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.trialsjournal.com/content/12/1/258/mathml/M1">View MathML</a>

(1)

where the triplet πi = {bi, di, ui} is described as the opinion of expert i. Intuitively it makes sense that there should be a constraint on these characteristics, as expressed in (1), as clearly when we have a maximum level of belief in a procedure we must necessarily have zero disbelief and uncertainty. Similarly, when there is a maximum level of uncertainty there clearly must be zero levels of belief and disbelief. The constraint that our levels of belief, disbelief and uncertainty must sum to unity is of course a matter of convenience, in an analogous manner to that in conventional probability where the same constraint is used. It seems reasonable, using statistical arguments, that we should scale our levels of belief and disbelief about the effectiveness of a procedure by the associated uncertainty. That is we are interested in the quantities b/u and d/u, in the same way we might want to normalize a treatment difference in an RCT by the associated standard deviation measuring the spread or uncertainty in the estimated difference to give an effect size. In order to estimate b, d and u, we need to develop a model for the clinical expert assessment data.

3.2 Parametric model

3.2.1 Assessment pooling

The assessment of the likely effectiveness of the intervention x was scored on a discrete valued symmetric scale with descriptive terms selected to imply an even spacing between categories. For our selected example, the seven-category ordinal scale, described in Section 2, was transformed onto the interval [0,1] as follows; <a onClick="popup('http://www.trialsjournal.com/content/12/1/258/mathml/M2','MathML',630,470);return false;" target="_blank" href="http://www.trialsjournal.com/content/12/1/258/mathml/M2">View MathML</a>, <a onClick="popup('http://www.trialsjournal.com/content/12/1/258/mathml/M3','MathML',630,470);return false;" target="_blank" href="http://www.trialsjournal.com/content/12/1/258/mathml/M3">View MathML</a>, <a onClick="popup('http://www.trialsjournal.com/content/12/1/258/mathml/M4','MathML',630,470);return false;" target="_blank" href="http://www.trialsjournal.com/content/12/1/258/mathml/M4">View MathML</a>, <a onClick="popup('http://www.trialsjournal.com/content/12/1/258/mathml/M5','MathML',630,470);return false;" target="_blank" href="http://www.trialsjournal.com/content/12/1/258/mathml/M5">View MathML</a>, <a onClick="popup('http://www.trialsjournal.com/content/12/1/258/mathml/M6','MathML',630,470);return false;" target="_blank" href="http://www.trialsjournal.com/content/12/1/258/mathml/M6">View MathML</a>, <a onClick="popup('http://www.trialsjournal.com/content/12/1/258/mathml/M7','MathML',630,470);return false;" target="_blank" href="http://www.trialsjournal.com/content/12/1/258/mathml/M7">View MathML</a> and <a onClick="popup('http://www.trialsjournal.com/content/12/1/258/mathml/M8','MathML',630,470);return false;" target="_blank" href="http://www.trialsjournal.com/content/12/1/258/mathml/M8">View MathML</a>. This retains the implicit spacing of the ordinal scale and centres the new scale at the same point as the original scale. Equivalent arguments can be constructed for ordinal scales with different numbers of categories.

Let xi, where 0 ≤ xi ≤1, quantify the likely effectiveness of a procedure for individual expert i as part of a panel of n experts. The distribution of xi is assumed to follow an approximate Beta distribution (Figure 1), a continuous probability distribution defined on the interval (0,1) and parameterized by two positive parameters, denoted by α and β, that modify the shape of the distribution. The Beta distribution is widely used for modelling random probabilities, particularly in the context of Bayesian analysis [17] and has been used to describe not only variability within a population as in a conventional statistical model, but also to describe the subjective degree of belief in a Bayesian sense [8]. Expressed mathematically, the probability density function for xi is

thumbnailFigure 1. Beta distributions B(α, β) for various values of parameters α and β.

<a onClick="popup('http://www.trialsjournal.com/content/12/1/258/mathml/M9','MathML',630,470);return false;" target="_blank" href="http://www.trialsjournal.com/content/12/1/258/mathml/M9">View MathML</a>,

where Γ(.) is the gamma function and parameters αi ≥ 1 and βi ≥ 1, requiring that the distribution be unimodal or at the extreme case, when αi = βI = 1, uniform. In the surgical trial setting described here, it seems unlikely that for instance a u-shaped distribution for xi (e.g. α = 0.5 and β = 0.5) would be plausible.

The multiplicative pooled assessment [18,19] of the expert panel is obtained as

<a onClick="popup('http://www.trialsjournal.com/content/12/1/258/mathml/M10','MathML',630,470);return false;" target="_blank" href="http://www.trialsjournal.com/content/12/1/258/mathml/M10">View MathML</a>

where f0(x) follows a Beta distribution with parameters <a onClick="popup('http://www.trialsjournal.com/content/12/1/258/mathml/M11','MathML',630,470);return false;" target="_blank" href="http://www.trialsjournal.com/content/12/1/258/mathml/M11">View MathML</a> and <a onClick="popup('http://www.trialsjournal.com/content/12/1/258/mathml/M12','MathML',630,470);return false;" target="_blank" href="http://www.trialsjournal.com/content/12/1/258/mathml/M12">View MathML</a>. This provides a pooled assessment that represents the intersection of the beliefs of the expert panel [19].

3.2.2 Opinion model

In order to translate the assessments from the panel of n experts to a collective expert opinion, the measures <a onClick="popup('http://www.trialsjournal.com/content/12/1/258/mathml/M13','MathML',630,470);return false;" target="_blank" href="http://www.trialsjournal.com/content/12/1/258/mathml/M13">View MathML</a> and <a onClick="popup('http://www.trialsjournal.com/content/12/1/258/mathml/M14','MathML',630,470);return false;" target="_blank" href="http://www.trialsjournal.com/content/12/1/258/mathml/M14">View MathML</a> (Section 3.1), that characterise the pooled opinion, are related to the pooled assessment parameters <a onClick="popup('http://www.trialsjournal.com/content/12/1/258/mathml/M15','MathML',630,470);return false;" target="_blank" href="http://www.trialsjournal.com/content/12/1/258/mathml/M15">View MathML</a> and <a onClick="popup('http://www.trialsjournal.com/content/12/1/258/mathml/M16','MathML',630,470);return false;" target="_blank" href="http://www.trialsjournal.com/content/12/1/258/mathml/M16">View MathML</a>. Equating the level of belief expressed by an expert to the pooled assessments, it is clear that <a onClick="popup('http://www.trialsjournal.com/content/12/1/258/mathml/M15','MathML',630,470);return false;" target="_blank" href="http://www.trialsjournal.com/content/12/1/258/mathml/M15">View MathML</a> must be proportional to <a onClick="popup('http://www.trialsjournal.com/content/12/1/258/mathml/M13','MathML',630,470);return false;" target="_blank" href="http://www.trialsjournal.com/content/12/1/258/mathml/M13">View MathML</a>, that is a larger value of <a onClick="popup('http://www.trialsjournal.com/content/12/1/258/mathml/M15','MathML',630,470);return false;" target="_blank" href="http://www.trialsjournal.com/content/12/1/258/mathml/M15">View MathML</a> represents a greater degree of belief; at the extreme as <a onClick="popup('http://www.trialsjournal.com/content/12/1/258/mathml/M17','MathML',630,470);return false;" target="_blank" href="http://www.trialsjournal.com/content/12/1/258/mathml/M17">View MathML</a>, then <a onClick="popup('http://www.trialsjournal.com/content/12/1/258/mathml/M18','MathML',630,470);return false;" target="_blank" href="http://www.trialsjournal.com/content/12/1/258/mathml/M18">View MathML</a> and <a onClick="popup('http://www.trialsjournal.com/content/12/1/258/mathml/M19','MathML',630,470);return false;" target="_blank" href="http://www.trialsjournal.com/content/12/1/258/mathml/M19">View MathML</a>, when we have maximum belief we must have minimum uncertainty. Similarly arguments lead to <a onClick="popup('http://www.trialsjournal.com/content/12/1/258/mathml/M16','MathML',630,470);return false;" target="_blank" href="http://www.trialsjournal.com/content/12/1/258/mathml/M16">View MathML</a> being proportional to <a onClick="popup('http://www.trialsjournal.com/content/12/1/258/mathml/M14','MathML',630,470);return false;" target="_blank" href="http://www.trialsjournal.com/content/12/1/258/mathml/M14">View MathML</a>; a larger value of <a onClick="popup('http://www.trialsjournal.com/content/12/1/258/mathml/M16','MathML',630,470);return false;" target="_blank" href="http://www.trialsjournal.com/content/12/1/258/mathml/M16">View MathML</a> represents a greater degree of disbelief. Although, clearly from example (a) in Figure 1, when the pooled Beta distribution parameter estimates are at their minimum and <a onClick="popup('http://www.trialsjournal.com/content/12/1/258/mathml/M20','MathML',630,470);return false;" target="_blank" href="http://www.trialsjournal.com/content/12/1/258/mathml/M20">View MathML</a> and <a onClick="popup('http://www.trialsjournal.com/content/12/1/258/mathml/M21','MathML',630,470);return false;" target="_blank" href="http://www.trialsjournal.com/content/12/1/258/mathml/M21">View MathML</a> then there is maximum uncertainty (<a onClick="popup('http://www.trialsjournal.com/content/12/1/258/mathml/M22','MathML',630,470);return false;" target="_blank" href="http://www.trialsjournal.com/content/12/1/258/mathml/M22">View MathML</a>) and minimum belief and disbelief, <a onClick="popup('http://www.trialsjournal.com/content/12/1/258/mathml/M23','MathML',630,470);return false;" target="_blank" href="http://www.trialsjournal.com/content/12/1/258/mathml/M23">View MathML</a>. Formalizing these arguments leads to the following expressions that satisfy all these conditions

<a onClick="popup('http://www.trialsjournal.com/content/12/1/258/mathml/M24','MathML',630,470);return false;" target="_blank" href="http://www.trialsjournal.com/content/12/1/258/mathml/M24">View MathML</a>

(2-3)

Solving equations (2) and (3), along with the condition that <a onClick="popup('http://www.trialsjournal.com/content/12/1/258/mathml/M25','MathML',630,470);return false;" target="_blank" href="http://www.trialsjournal.com/content/12/1/258/mathml/M25">View MathML</a> (equation 1), yields the following expressions that characterize the relationship between the triplet <a onClick="popup('http://www.trialsjournal.com/content/12/1/258/mathml/M26','MathML',630,470);return false;" target="_blank" href="http://www.trialsjournal.com/content/12/1/258/mathml/M26">View MathML</a> and the parameters <a onClick="popup('http://www.trialsjournal.com/content/12/1/258/mathml/M15','MathML',630,470);return false;" target="_blank" href="http://www.trialsjournal.com/content/12/1/258/mathml/M15">View MathML</a> and <a onClick="popup('http://www.trialsjournal.com/content/12/1/258/mathml/M16','MathML',630,470);return false;" target="_blank" href="http://www.trialsjournal.com/content/12/1/258/mathml/M16">View MathML</a>,

<a onClick="popup('http://www.trialsjournal.com/content/12/1/258/mathml/M27','MathML',630,470);return false;" target="_blank" href="http://www.trialsjournal.com/content/12/1/258/mathml/M27">View MathML</a>

(4-6)

where the triplet <a onClick="popup('http://www.trialsjournal.com/content/12/1/258/mathml/M26','MathML',630,470);return false;" target="_blank" href="http://www.trialsjournal.com/content/12/1/258/mathml/M26">View MathML</a> clearly satisfies <a onClick="popup('http://www.trialsjournal.com/content/12/1/258/mathml/M25','MathML',630,470);return false;" target="_blank" href="http://www.trialsjournal.com/content/12/1/258/mathml/M25">View MathML</a>; a more detailed derivation of equations (4)-(6) is provided elsewhere [15,16]. Thus, when <a onClick="popup('http://www.trialsjournal.com/content/12/1/258/mathml/M28','MathML',630,470);return false;" target="_blank" href="http://www.trialsjournal.com/content/12/1/258/mathml/M28">View MathML</a>, <a onClick="popup('http://www.trialsjournal.com/content/12/1/258/mathml/M29','MathML',630,470);return false;" target="_blank" href="http://www.trialsjournal.com/content/12/1/258/mathml/M29">View MathML</a> and the pooled opinion is total uncertainty (ignorance); see example (a) in Figure 1. If parameters <a onClick="popup('http://www.trialsjournal.com/content/12/1/258/mathml/M15','MathML',630,470);return false;" target="_blank" href="http://www.trialsjournal.com/content/12/1/258/mathml/M15">View MathML</a> and <a onClick="popup('http://www.trialsjournal.com/content/12/1/258/mathml/M16','MathML',630,470);return false;" target="_blank" href="http://www.trialsjournal.com/content/12/1/258/mathml/M16">View MathML</a> are greater than unity but equal, we have equal belief and disbelief; for example (b) in Figure 1 where <a onClick="popup('http://www.trialsjournal.com/content/12/1/258/mathml/M30','MathML',630,470);return false;" target="_blank" href="http://www.trialsjournal.com/content/12/1/258/mathml/M30">View MathML</a> and <a onClick="popup('http://www.trialsjournal.com/content/12/1/258/mathml/M31','MathML',630,470);return false;" target="_blank" href="http://www.trialsjournal.com/content/12/1/258/mathml/M31">View MathML</a>. As <a onClick="popup('http://www.trialsjournal.com/content/12/1/258/mathml/M15','MathML',630,470);return false;" target="_blank" href="http://www.trialsjournal.com/content/12/1/258/mathml/M15">View MathML</a> increases relative to <a onClick="popup('http://www.trialsjournal.com/content/12/1/258/mathml/M16','MathML',630,470);return false;" target="_blank" href="http://www.trialsjournal.com/content/12/1/258/mathml/M16">View MathML</a> the belief increases and the uncertainty decreases and conversely as <a onClick="popup('http://www.trialsjournal.com/content/12/1/258/mathml/M16','MathML',630,470);return false;" target="_blank" href="http://www.trialsjournal.com/content/12/1/258/mathml/M16">View MathML</a> increases relative to <a onClick="popup('http://www.trialsjournal.com/content/12/1/258/mathml/M15','MathML',630,470);return false;" target="_blank" href="http://www.trialsjournal.com/content/12/1/258/mathml/M15">View MathML</a> the disbelief increases and the uncertainty decreases; these two scenarios are illustrated in examples (d) and (c) in Figure 1, where <a onClick="popup('http://www.trialsjournal.com/content/12/1/258/mathml/M32','MathML',630,470);return false;" target="_blank" href="http://www.trialsjournal.com/content/12/1/258/mathml/M32">View MathML</a>, <a onClick="popup('http://www.trialsjournal.com/content/12/1/258/mathml/M33','MathML',630,470);return false;" target="_blank" href="http://www.trialsjournal.com/content/12/1/258/mathml/M33">View MathML</a> and <a onClick="popup('http://www.trialsjournal.com/content/12/1/258/mathml/M34','MathML',630,470);return false;" target="_blank" href="http://www.trialsjournal.com/content/12/1/258/mathml/M34">View MathML</a> and <a onClick="popup('http://www.trialsjournal.com/content/12/1/258/mathml/M35','MathML',630,470);return false;" target="_blank" href="http://www.trialsjournal.com/content/12/1/258/mathml/M35">View MathML</a>, <a onClick="popup('http://www.trialsjournal.com/content/12/1/258/mathml/M36','MathML',630,470);return false;" target="_blank" href="http://www.trialsjournal.com/content/12/1/258/mathml/M36">View MathML</a> and <a onClick="popup('http://www.trialsjournal.com/content/12/1/258/mathml/M37','MathML',630,470);return false;" target="_blank" href="http://www.trialsjournal.com/content/12/1/258/mathml/M37">View MathML</a>.

3.3 Nonparametric model

An alternative nonparametric formulation for belief, disbelief and uncertainty allows a more general approach to that described in Section 3.2. Defining μi and σi as the mean and standard deviation of the assessment of the effectiveness of the intervention xi for expert i, where xi is in the range [0,1]. Then the uncertainty (ui), belief (bi) and disbelief (di) associated with an opinion for expert i can be expressed as <a onClick="popup('http://www.trialsjournal.com/content/12/1/258/mathml/M38','MathML',630,470);return false;" target="_blank" href="http://www.trialsjournal.com/content/12/1/258/mathml/M38">View MathML</a>, bi = μi(1-ui) and bi = di =(1-μi)(1-ui); as 0 ≤ xi ≤ 1, then 0 ≤ ui ≤ 1 and the measures satisfy equation (1). For example using the data from Table 1 for expert 4 from case 3, the weighted mean and standard deviation, based on the transformed seven-category ordinal scale described in Section 3.2.1 <a onClick="popup('http://www.trialsjournal.com/content/12/1/258/mathml/M39','MathML',630,470);return false;" target="_blank" href="http://www.trialsjournal.com/content/12/1/258/mathml/M39">View MathML</a> with weights given by (5,15,20,20,20,15,5), are μ = 0.5 and σ = 0.226, and so u = 0.204 and b = d = 0.398. Multiplicative pooling leads directly to estimates for the opinion triplet <a onClick="popup('http://www.trialsjournal.com/content/12/1/258/mathml/M26','MathML',630,470);return false;" target="_blank" href="http://www.trialsjournal.com/content/12/1/258/mathml/M26">View MathML</a>, with weights given by the nth root of the product of the individual expert weights, in an analogous manner to that described in Section 3.2.1 for the parametric model.

In fact the expressions for uncertainty, belief and disbelief for the Beta model in equations (4)-(6) follow directly from the above expressions for u, b and d, based on μ and σ, after some rescaling, by noting that the mean and variance of the Beta distribution are α/(α+β) and αβ/{(α+β)2 (α+β+1)} respectively.

3.4 Opinion space

As proposed by Jøsang [15], a ternary plot provides a convenient method of representing the triplet of belief, disbelief and uncertainty that constitute a pooled expert opinion. A ternary plot represents the ratios of the three variables as positions in an equilateral triangle, where each base, or side, of the triangle represents a proportion, with the point of the triangle opposite that base representing a proportion equal to one. As a proportion increases in any one sample, the point representing that sample moves from the base to the opposite point of the triangle. For instance, when <a onClick="popup('http://www.trialsjournal.com/content/12/1/258/mathml/M28','MathML',630,470);return false;" target="_blank" href="http://www.trialsjournal.com/content/12/1/258/mathml/M28">View MathML</a> (maximum uncertainty) the opinion is mapped to the apex of the equilateral triangle, whereas when <a onClick="popup('http://www.trialsjournal.com/content/12/1/258/mathml/M30','MathML',630,470);return false;" target="_blank" href="http://www.trialsjournal.com/content/12/1/258/mathml/M30">View MathML</a> there is equal belief, disbelief and uncertainty and the pooled opinion is mapped to the centre of the triangle. The cases representing greater levels of belief and greater levels of disbelief are mapped towards the right-hand and left-hand vertices of the triangle respectively.

3.5 Decision rules

In order to determine the level of equipoise that should be satisfied for a clinical trial to be considered ethical Johnson et al. [20] conducted an ethometric study to investigate how much clinical equipoise can be disturbed before potential trial subjects deem it to be unethical. A series of hypothetical clinical trial scenarios were presented to people from a broad range of societal and geographical groups within the UK. They were asked to specify the level of collective doubt between two treatment modalities that they would accept if casting a vote on an ethics committee. Johnson et al. [20] defined the 80:20 rule, that represented the split in equipoise that should be allowed for a trial to be judged to be ethical and recommended its use as an appropriate tool for deciding whether recruitment is ethically justifiable; based on their empirical evidence that less than 3% of subjects questioned thought that a trial should morally be undertaken if equipoise was beyond this point. By way of comparison, an alternative mean threshold rule might consider it ethical to recruit patients if the mean clinical effectiveness (μ), estimated as α/(α+β) for the Beta distribution, were within pre-determined limits. For instance, it might be considered ethical to recruit patients into a trial if the mean clinical effectiveness were in the range 0.4 ≤ μ ≤ 0.7.

The 80:20 and mean threshold equipoise decision rules can be mapped onto the opinion space and visualized on a ternary plot. For the Beta model (Section 3.2), the former rule can be mapped on to the ternary plot by iteratively finding solutions for Beta distribution parameters, α and β, that give estimates for the probability density function equal to 0.2 and 0.8 to the left and right of the central point on the expert rating scale, and for the latter rule by simply solving equations (4)-(6) using the constraint that μ(α+β)=α.

3.6 Hypothesis testing

The significance of the estimated pooled opinion (<a onClick="popup('http://www.trialsjournal.com/content/12/1/258/mathml/M40','MathML',630,470);return false;" target="_blank" href="http://www.trialsjournal.com/content/12/1/258/mathml/M40">View MathML</a>) is assessed using resampling. For the Beta model for Section 3.2, pooled assessment parameters <a onClick="popup('http://www.trialsjournal.com/content/12/1/258/mathml/M41','MathML',630,470);return false;" target="_blank" href="http://www.trialsjournal.com/content/12/1/258/mathml/M41">View MathML</a> and <a onClick="popup('http://www.trialsjournal.com/content/12/1/258/mathml/M42','MathML',630,470);return false;" target="_blank" href="http://www.trialsjournal.com/content/12/1/258/mathml/M42">View MathML</a> are estimated for Sm, a set of size n constructed by sampling with replacement from {1,2...n}; for example for the pooled assessment of 5 experts Sm might be {1,2,2,4,1} or {5,3,3,1,1}. This process is repeated many times by random construction of Sm to give empirical bootstrap [21] distributions <a onClick="popup('http://www.trialsjournal.com/content/12/1/258/mathml/M43','MathML',630,470);return false;" target="_blank" href="http://www.trialsjournal.com/content/12/1/258/mathml/M43">View MathML</a> and <a onClick="popup('http://www.trialsjournal.com/content/12/1/258/mathml/M44','MathML',630,470);return false;" target="_blank" href="http://www.trialsjournal.com/content/12/1/258/mathml/M44">View MathML</a>, and thereby <a onClick="popup('http://www.trialsjournal.com/content/12/1/258/mathml/M45','MathML',630,470);return false;" target="_blank" href="http://www.trialsjournal.com/content/12/1/258/mathml/M45">View MathML</a>. From this empirical distribution, a bootstrap confidence interval for <a onClick="popup('http://www.trialsjournal.com/content/12/1/258/mathml/M40','MathML',630,470);return false;" target="_blank" href="http://www.trialsjournal.com/content/12/1/258/mathml/M40">View MathML</a> is derived for the purpose of hypothesis testing. A similar resampling scheme can also be developed simply for the nonparametric model of Section 3.3.

This resampling methodology represents the variability in opinion that might be obtained for any combination of experts in the panel, including in principle a panel composed entirely of a single expert, and as such represents the full range of possible opinions for the selected population of experts. For the relative small panel of experts in our example, exhaustive permutation resampling [21] is the preferred option, but this may be computational unrealistic for large n where bootstrapping with M = 1000 would be sufficient.

3.7 Beta distribution fitting

The outlined statistical framework is illustrated using the example data introduced in Section 2 (Table 1). We focus here on the Beta model (Section 3.2) as an exemplar, as this fits our data well and is computational slightly more complex to implement than the nonparametric method. Statistical analysis was undertaken in the statistical software package R [22]. Code to replicate the analysis presented here is available on request from the corresponding author.

The parameters of the Beta distribution were estimated for each clinical expert for the four cases shown in Table 1 using the fitdistr function available in the MASS [23] library in the statistical software package R [22]. This function estimates parameters for a range of univariate distributions, including the Beta distribution, using maximum-likelihood methods. For the four example cases introduced in Section 2 the pooled parameter estimates were <a onClick="popup('http://www.trialsjournal.com/content/12/1/258/mathml/M46','MathML',630,470);return false;" target="_blank" href="http://www.trialsjournal.com/content/12/1/258/mathml/M46">View MathML</a>, <a onClick="popup('http://www.trialsjournal.com/content/12/1/258/mathml/M47','MathML',630,470);return false;" target="_blank" href="http://www.trialsjournal.com/content/12/1/258/mathml/M47">View MathML</a>, <a onClick="popup('http://www.trialsjournal.com/content/12/1/258/mathml/M48','MathML',630,470);return false;" target="_blank" href="http://www.trialsjournal.com/content/12/1/258/mathml/M48">View MathML</a>, <a onClick="popup('http://www.trialsjournal.com/content/12/1/258/mathml/M49','MathML',630,470);return false;" target="_blank" href="http://www.trialsjournal.com/content/12/1/258/mathml/M49">View MathML</a> and <a onClick="popup('http://www.trialsjournal.com/content/12/1/258/mathml/M50','MathML',630,470);return false;" target="_blank" href="http://www.trialsjournal.com/content/12/1/258/mathml/M50">View MathML</a>, <a onClick="popup('http://www.trialsjournal.com/content/12/1/258/mathml/M51','MathML',630,470);return false;" target="_blank" href="http://www.trialsjournal.com/content/12/1/258/mathml/M51">View MathML</a>, <a onClick="popup('http://www.trialsjournal.com/content/12/1/258/mathml/M52','MathML',630,470);return false;" target="_blank" href="http://www.trialsjournal.com/content/12/1/258/mathml/M52">View MathML</a>, <a onClick="popup('http://www.trialsjournal.com/content/12/1/258/mathml/M53','MathML',630,470);return false;" target="_blank" href="http://www.trialsjournal.com/content/12/1/258/mathml/M53">View MathML</a>. The fitted distributions for each clinical expert and pooled estimates are shown in Figure 2.

thumbnailFigure 2. Fitted Beta distributions for each clinical expert (---) and pooled estimates (-) for each case.

3.8 Opinions

The pooled parameter estimates from the Beta distribution fitting for the four example cases were used to estimates the belief, disbelief and uncertainty using equations (4)-(6); this gave the following estimates, <a onClick="popup('http://www.trialsjournal.com/content/12/1/258/mathml/M54','MathML',630,470);return false;" target="_blank" href="http://www.trialsjournal.com/content/12/1/258/mathml/M54">View MathML</a>, <a onClick="popup('http://www.trialsjournal.com/content/12/1/258/mathml/M55','MathML',630,470);return false;" target="_blank" href="http://www.trialsjournal.com/content/12/1/258/mathml/M55">View MathML</a>, <a onClick="popup('http://www.trialsjournal.com/content/12/1/258/mathml/M56','MathML',630,470);return false;" target="_blank" href="http://www.trialsjournal.com/content/12/1/258/mathml/M56">View MathML</a>, <a onClick="popup('http://www.trialsjournal.com/content/12/1/258/mathml/M57','MathML',630,470);return false;" target="_blank" href="http://www.trialsjournal.com/content/12/1/258/mathml/M57">View MathML</a>, <a onClick="popup('http://www.trialsjournal.com/content/12/1/258/mathml/M58','MathML',630,470);return false;" target="_blank" href="http://www.trialsjournal.com/content/12/1/258/mathml/M58">View MathML</a>, <a onClick="popup('http://www.trialsjournal.com/content/12/1/258/mathml/M59','MathML',630,470);return false;" target="_blank" href="http://www.trialsjournal.com/content/12/1/258/mathml/M59">View MathML</a>, <a onClick="popup('http://www.trialsjournal.com/content/12/1/258/mathml/M60','MathML',630,470);return false;" target="_blank" href="http://www.trialsjournal.com/content/12/1/258/mathml/M60">View MathML</a>, <a onClick="popup('http://www.trialsjournal.com/content/12/1/258/mathml/M61','MathML',630,470);return false;" target="_blank" href="http://www.trialsjournal.com/content/12/1/258/mathml/M61">View MathML</a> and <a onClick="popup('http://www.trialsjournal.com/content/12/1/258/mathml/M62','MathML',630,470);return false;" target="_blank" href="http://www.trialsjournal.com/content/12/1/258/mathml/M62">View MathML</a>, <a onClick="popup('http://www.trialsjournal.com/content/12/1/258/mathml/M63','MathML',630,470);return false;" target="_blank" href="http://www.trialsjournal.com/content/12/1/258/mathml/M63">View MathML</a>, <a onClick="popup('http://www.trialsjournal.com/content/12/1/258/mathml/M64','MathML',630,470);return false;" target="_blank" href="http://www.trialsjournal.com/content/12/1/258/mathml/M64">View MathML</a>, <a onClick="popup('http://www.trialsjournal.com/content/12/1/258/mathml/M65','MathML',630,470);return false;" target="_blank" href="http://www.trialsjournal.com/content/12/1/258/mathml/M65">View MathML</a>. Inspection of Figure 2, indicates that there appears to be significant belief for case 2 that the patient will improve after treatment (surgery) and conversely significant disbelief in the effectiveness of the treatment for case for case 4; this is reflected in the large (> 0.6) estimates of b and d for cases 2 and 4 respectively. Also, there is significant uncertainty, seen by the flatness of the curves in Figure 2(c), in the collective opinions of the experts for case 3; this is apparent in the large level of uncertainty for this case, relative to the other cases.

3.9 Decision rules

In order to determine whether an opinion provides sufficient evidence for eligibility for recruitment to a clinical trial, we must first define a decision rule. Here we focus on two rules, the 80:20 [20] and the mean threshold rules; although the procedures described here are equally applicable to many more rules that could potentially be defined. The 80:20 and mean threshold rules partition the opinion space, visualized by the ternary plot, into regions that determine whether the patient can or cannot ethically be recruited to a trial.

The division lines between the regions for the 80:20 rule were determined iteratively (using an interval search method) by finding estimates of the Beta distribution parameters α and β that exactly divided the probability density 80% and 20% around equipoise, and projecting these estimates into the opinion space using equations (4)-(6). This process was achieved using an implementation of the uniroot function in R [22]. After discussion with the clinical experts it became clear that the point of equipoise for the assessment scale described in Section 2 for the 80:20 rule was not located centrally but was in fact located at the division between the 'No difference' and the 'A bit better' categories. That is, because surgery was seen to be an active intervention for a condition that required treatment, the point of equipoise was located slightly to the right of the centre point of the assessment scale; which for our definition of the assessment scale is at <a onClick="popup('http://www.trialsjournal.com/content/12/1/258/mathml/M66','MathML',630,470);return false;" target="_blank" href="http://www.trialsjournal.com/content/12/1/258/mathml/M66">View MathML</a> rather than at <a onClick="popup('http://www.trialsjournal.com/content/12/1/258/mathml/M67','MathML',630,470);return false;" target="_blank" href="http://www.trialsjournal.com/content/12/1/258/mathml/M67">View MathML</a> on the interval (0, 1). The asymmetry that this implies for the 80:20 decision rule is clear in Figure 3. The mean threshold rule divided the opinion space into three distinct regions μ < 0.4, 0.4 ≤ μ ≤ 0.7 and μ < 0.7 characterised by the thresholds 0.4 and 0.7 for the mean, that determined whether the intervention was likely to be effective. The divisions between regions were mapped onto the opinion space by solving equations (4)-(6) using the constraint that μ(α+β) = α. For instance for μ = 0.7 and u = 0, then b = 0.7 and d = 0.3 and when d = 0 then <a onClick="popup('http://www.trialsjournal.com/content/12/1/258/mathml/M68','MathML',630,470);return false;" target="_blank" href="http://www.trialsjournal.com/content/12/1/258/mathml/M68">View MathML</a> and <a onClick="popup('http://www.trialsjournal.com/content/12/1/258/mathml/M69','MathML',630,470);return false;" target="_blank" href="http://www.trialsjournal.com/content/12/1/258/mathml/M69">View MathML</a>; these points define the intersections between the upper division boundary with the lower and right edges of the ternary plot in Figure 3.

thumbnailFigure 3. Estimated triplets for all permutations of opinions with the 80:20 (---) and mean (--) decision rules. The best estimate of collective opinion is given by the large symbol (•).

3.10 Hypothesis testing

The exhaustive permutation test described in Section 3.6 was applied to each of the test cases. This gave 462, 126, 35 and 126 combinations of opinions for the four cases that used respectively 6, 5, 4 and 5 expert clinical assessors. The belief, disbelief and uncertainty for all the combinations of opinion were estimated for each of the four cases and plotted along with the decision rules in Figure 3.

The 'cloud' of points for each case represents the variability due to the range of opinions expressed by the expert assessors. Where there were considerable differences of opinion, for instance for case 1, there was a much wider spread of points than where there was overall agreement amongst the experts about the likelihood of success of the intervention, for instance for case 2 or 4. It is instructive to look at one particular opinion triplet to more fully understand the meaning of the ternary plots.

For case 1, the opinion triplet π = {0.712,0.211,0.077} located towards the lower right hand vertex of the ternary plot has very high belief and low uncertainty. This is the opinion associated with six replicates of the assessment of clinical expert 3 for case 1 (see Table 1), who had a strong belief that the patient would get significantly better after treatment. If this expert assessor were indeed representative of the wider population of experts, then it would certainly be unethical for the patient to be recruited to the trial and consequentially the opinion for this potential scenario is located to the right of the 80:20 and mean threshold decision rules.

Labelling the regions to the lower right and lower left of the plots to the right and left of the 80:20 and mean threshold decision rule partition curves as the 'belief' and 'disbelief' regions, allows us to count the number of opinions falling within these regions for each case and rule; see Table 2. Defining the null hypothesis to be that a case should not be recruited to the trial, Table 2 provides evidence for this hypothesis and suggests appropriate p-values based on the 80:20 rule for the four cases to be 0.026 (i.e. 12/462), 0.333, 0.000 and 1.000 and based on the mean threshold rule to be 0.011 (i.e. 5/462), 0.000, 0.029 and 1.000. Testing at the 5% level (two-sided) indicates that for the 80:20 rule cases 1 and 3 would be eligible for recruitment and for the mean threshold rule cases 1, 2 and 3 would be eligible for recruitment. For this decision making process to have some validity, the decision rule and the significance level would clearly need to be stated before data collection was undertaken.

Table 2. Opinion counts by case, decision region and rule, and the total number of opinion combinations available for the exhaustive permutation test

4 Conclusions

We describe a statistical framework for the assessment of clinical uncertainty, as a prelude to a clinical trial and demonstrate, using data from the UK Heel Fracture Trial, how expert opinions can be pooled, modelled and presented on a ternary plot that represents an opinion space. Individual cases can then be assessed in relation to decision rules mapped onto the opinion space, providing clear and rapid decisions regarding trial eligibility. The methodology has potential to identify eligible patients and assist in the simplification of eligibility criteria which might encourage greater participation in clinical trials.

Methods for the assessment of clinical uncertainty, as a prelude to a clinical trial, have been suggested previously [10,11]. However, the methodology described here is the first attempt at a structured statistical framework to undertake this type of analysis. Beta distributions were fitted to assessments of the likely effectiveness of an intervention elicited from a virtual panel of experts and pooled using methods familiar to exponents of determining expert probabilities [19]. Opinions were expressed using previously suggested [15] definitions of belief, disbelief and uncertainty that we believe fully characterised the clinical expert assessments. Our analysis restricted the choice of Beta distributions for modelling to unimodal forms (α ≥ 1 and β ≥ 1). This was not a concern for the examples described here or indeed more widely for other data we have explored in the setting of surgical trials. However, it is in principle possible in other applications that the most likely assessment of clinical effectiveness of an intervention is that a patient would either get much better or much worse with any other outcome being extremely unlikely. In this setting belief, disbelief and uncertainty as expressed in equations (4)-(6) would not be defined. For the data presented here the Beta model proved to be the most informative, however where this is not the case the nonparametric methods described, based on estimated means and standard deviations, provide useful alternatives for any distribution on the interval [0,1]. Although the examples described here all use seven point likert type scales for elicitation, the statistical framework introduced would work equally well with any type of ordered categorical assessment scale.

Expert opinions are pooled here using multiplicative methods [19], as we felt that this best represented clinical equipoise [24] and the views of the experts consulted for the example data; i.e. that all experts opinions were 'correct' and the pooling should represent the consensus based on the intersection of beliefs. However, our view is pragmatic and we see no reason why additive pooling could not be used in preference to multiplicative pooling, particularly if it was felt that the latter method was giving too much weight to the assessment of one or more 'over-confident' individual experts.

We have presented significance tests at the 5% level to assess whether a patient might ethically be recruited to a trial. Our selection of this level for the tests was somewhat arbitrary and clearly this could be set, prior to analysis, at a higher or lower level for a different application or a less formal procedure adopted if necessary. The 80:20 rule [20], which is based on some empirical evidence, was selected as a standard for decision making regarding recruitment. The alternative mean threshold rule, as well as being intuitively reasonable, was suggested in part to encourage some debate as to what form the decision rule should take for different cases and in various settings. This is clearly an area that requires additional research.

The focus of this paper has been on developing tools for improving recruitment to trials. For those patients deemed eligible for recruitment who decide to enter an RCT, it would seem natural to use the expert evidence elicited through this process as a clinical prior, based on subjective opinion, in a formal Bayesian analysis [14].

The methodological framework discussed here has provided additional insight that would otherwise have not been available for the heel fracture trial. Although, clearly this methodology will need to be assessed in future studies to identify whether it can actually deliver improvement in trial recruitment rates. The methodological framework we describe is currently limited to two-arm trials, although we see no reason why this could not be extended to more than two treatment groups. The opinion pooling we describe is appropriate for situations where individual expert opinions may differ to a moderate or large extent, but it is not at all clear that pooling opinions where for instance experts have totally opposing views (100% belief or disbelief in treatment effectiveness) would be appropriate, as the pooled opinion would in reality represent no individual expert's opinion. Therefore we would recommend the methodology be limited to only those scenarios of the former rather than the latter type. Although we have focussed on surgical trials, we would expect the methodology described here to be applicable to any RCT where recruitment was problematic. The methodology also has clear application in pilot studies where feasibility is being assessed and also potentially as a support tool for inclusive trials where patients are allowed to select an intervention as well as being randomised in a conventional manner [25].

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

DG, chief investigator for the UK Heel Fracture Trial, and YK developed the original concepts for the study. NP developed the statistical methodology, with input from YK and AG, and wrote the first draft of the paper. All authors contributed to the paper during development and read and approved the final version of the manuscript.

Acknowledgements

The UK Heel Fracture Trial was funded by Arthritis Research UK (UK HeFT Trial Grant 15964). AG acknowledges support of his contribution to this work through the MATCH programme (EPSRC Grant GR/S29874/01) and YK acknowledges the support of Warwick Medical School (University of Warwick) for PhD sponsorship.

References

  1. Boutron I, Ravaud P, Nizaed R: The design and assessment of prospective randomised, controlled trials in orthopaedic surgery.

    Journal of Bone and Joint Surgery (British) 2007, 89:858-863. Publisher Full Text OpenURL

  2. Horton R: Surgical research or comic opera: questions, but few answers.

    Lancet 1996, 347:984-985. PubMed Abstract | Publisher Full Text OpenURL

  3. McCulloch P, Taylor I, Sasako M, Lovett B, Griffin D: Randomised trials in surgery: problems and possible solutions.

    BMJ 2002, 324:1448-1451. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  4. Halpern SD: Evidence-based equipoise and research responsiveness.

    American Journal of Bioethics 2006, 6:1-4. PubMed Abstract | Publisher Full Text OpenURL

  5. Freedman B: Equipoise and the ethics of clinical research.

    New England Journal of Medicine 1987, 317:141-145. PubMed Abstract | Publisher Full Text OpenURL

  6. Chard JA, Lilford RJ: The use of equipoise in clinical trials.

    Social Science and Medicine 1998, 47:891-898. PubMed Abstract | Publisher Full Text OpenURL

  7. Jaynes ET: Probability theory: the logic of science. Cambridge: Cambridge University Press; 2003. OpenURL

  8. Spiegelhalter DJ, Abrams KR, Myles JP: Bayesian approaches to clinical trials and health-care evaluation. New York: Wiley; 2004. OpenURL

  9. Kandane JB, Wolfson LJ: Experiences in elicitation.

    The Statistician 1998, 47:3-19. OpenURL

  10. Lilford R: Formal measurement of clinical uncertainty: prelude to a trial in perinatal medicine.

    BMJ 1994, 308:111-112. PubMed Abstract | PubMed Central Full Text OpenURL

  11. Young J, Harrison J, White G, May J, Solomon M: Developing measures of surgeons' equipoise to assess the feasibility of randomized controlled trials in vascular surgery.

    Surgery 2004, 136:1070-1076. PubMed Abstract | Publisher Full Text OpenURL

  12. Parmar MKB, Griffiths GO, Spieglelhalter DJ, Souhami RL, Altman DG, van der Scheuren E: Monitoring of large randomised clinical trials: a new approach with Bayesian methods.

    Lancet 2001, 358:375-381. PubMed Abstract | Publisher Full Text OpenURL

  13. Parmar MKB, Spiegelhalter DJ, Freedman LS: The chart trials: Bayesian design and monitoring in practice.

    Statistics in Medicine 1994, 13:1297-1312. PubMed Abstract | Publisher Full Text OpenURL

  14. Spiegelhalter DJ, Freedman LS, Parmar MKB: Bayesian approaches to randomized trials.

    Journal of the Royal Statistical Society, Series A 1994, 157:357-416. OpenURL

  15. Jøsang A: Artificial reasoning with subjective logic.

    In Proceedings of the Second Australian Workshop on Commonsense Reasoning: December 1997; Perth Edited by Nayak A, Pagnucco M. OpenURL

  16. Jøsang A: A logic for uncertain probabilities.

    International Journal of Uncertainty, Fuzziness and Knowledge-Based Systems 2001, 9:279-311. OpenURL

  17. Evans M, Hastings N, Peacock B: Beta Distribution. In Statistical Distributions. 3rd edition. New York: Wiley; 2000:34-42. OpenURL

  18. Genest C, Zidek JV: Combining probability distributions: a critique and an annotated bibliography.

    Statistical Science 1986, 1:114-135. Publisher Full Text OpenURL

  19. O'Hagan A, Buck CE, Daneshkhah A, Eiser JR, Garthwaite PH, Jenkinson DJ, Oakley JE, Rakow T: Uncertain Judgements: Eliciting Expert Probabilities. Chichester: John Wiley and Sons; 2006. OpenURL

  20. Johnson N, Lilford RJ, Brazier W: At what level of collective equipoise does a clinical trial become ethical?

    Journal of Medical Ethics 1991, 17:30-34. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  21. Efron B, Tibshirani R: An introduction to the bootstrap. Boca Raton, Florida: Chapman and Hall/CRC; 1993. OpenURL

  22. R Developmental Core Team: R: A language and environment for statistical computing. [http://www.R-project.org] webcite

    R Foundation for Statistical Computing, Vienna, Austria; 2010.

  23. Venables WN, Ripley BD: Modern applied statistics in s-plus. 4th edition. New York: Springer; 2002. OpenURL

  24. Lilford R, Jackson J: Equipoise and the ethics of randomization.

    Journal of the Royal Society of Medicine 1995, 88:552-559. PubMed Abstract | PubMed Central Full Text OpenURL

  25. Torgerson D, Sibbald B: Understanding controlled trials: What is a patient preference trial?

    BMJ 1998, 316:360. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL