3  Measures of Evidence

3.1 General

3.1.1 Indirect

  • If we assume suspect A did not commit the crime, his behavior in the following days is very unusual, so let’s act as if A is the perpetrator.
  • We have only two suspects. The likelihood that suspect A committed the crime is small, so we turn attention to suspect B even though we don’t know the chance that B did it.

3.1.2 Direct

  • A police officer witnessed suspect B committing the crime (no analogy in drug development)
  • Based on DNA, fingerprint, and motive, the probability that suspect B committed the crime is 0.98
  • Note: it is impossible to assign an absolute probability that the suspect is guilty without having a pre-investigation prior probability of guilt (but one can compute relative guilt ratios without this)

3.2 Frequentist

  • p-value: P(data as or more extreme than what we observed | H0 true)
  • Probability here = long-run relative frequency
  • NHST: null hypothesis significance test
  • Backwards in time/information flow
  • Doesn’t relate to clinical significance (Mark et al. (2016))
  • Requires arbitrary multiplicity adjustments because considers what could have happened
  • Don’t want to know P(batter who just got a hit is left-handed)

Analogy to Diagnostic Testing

  • Sensitivity: P(T+ | D+) [power]
  • Specificity: P(T- | D-) [1 - type I “error”]
  • Bayes’ rule: P(D+|T) = P(T|D+) P(D+) / P(T)
    =sens x prevalence / (sens x prev + (1-spec)(1-prev))
  • High sens and spec easily overcome by low prev
  • sens and spec vary with patient
  • What could have happened but didn’t important, just as need to correct p-values for multiple data looks in sequential RCTs
    • need to correct sens and spec for verification/referral bias
    • after using Bayes’ rule to get P(D|T) the corrections cancel out
  • P(D|T) is directly actionable and defines its own error risk
    • P(cancer) = 0.8 implies P(unnecessary biopsy)=0.2

Back to p-values

  • p=“degree to which the data are embarrassed by the null hypothesis” (Maxwell (2004))
  • Only can provide evidence against H0, never evidence in favor of something
  • Efficacy inferred from having much evidence against “no efficacy”
  • If set α=0.05, type I assertion probability never reduces even as n→∞
  • Likelihood school of inference: both type I and type II assertion probabilities → 0
  • Type I α perhaps useful at study design stage
  • After study completes, can only know if type I “error” was committed if true effect is exactly zero
    • but then would not have needed the study
  • Type I assertion prob is a long-run operating characteristic for sequence of studies
  • Consider sequence of p-values from them
    • type I prob α means P(p-value < α | zero effect) = α
  • Neither p-value nor α are probs of a decision error
  • p-value = P(data more extreme than ours | no effect)
    • not a false + prob for experiment at hand
    • that would require a prior

A basic difficulty for most students is the proper formulation of the alternatives H0 and H1 for any given problem and the consequent determination of the proper critical region (upper tail, lower tail, two-sided). …

Comment. Small wonder that students have trouble. They may be trying to think. …

More on the teaching of statistics. Little advancement in the teaching of statistics is possible, and little hope for statistical methods to be useful in the frightful problems that face man today, until the literature and classroom be rid of terms so deadening to scientific enquiry as null hypothesis, population (in place of frame), true value, level of significance for comparison of treatments, representative sample.

Statistical significance of B over A thus conveys no knowledge, no basis for action.—Deming (1975)

The null-hypothesis significance test treats ‘acceptance’ or ‘rejection’ of a hypothesis as though these were decisions one makes. But a hypothesis is not something, like a piece of pie offered for dessert, which can be accepted or rejected by a voluntary physical action. Acceptance or rejection of a hypothesis is a cognitive process, a degree of believing or disbelieving which, if rational, is not a matter of choice but determined solely by how likely it is, given the evidence, that the hypothesis is true. – Rozeboom (1960) quoted by EJ Wagenmakers and Q Gronau

… Another concern is that Bayesian methods do not control error rates as indicated by p values. … This concern is countered by repeated demonstrations that error rates are extremely difficult to pin down because they are based on sampling and testing intentions. — Kruschke & Liddell (2017)

If the design were unknown, then it is not possible to calculate a P value. … Every practicing statistician must deal with data from experiments the designs of which have been compromised. For example, clinical trials are plagued with missing data, patients lost to follow-up, patients on the wrong dosing schedule, and so forth. Practicing statisticians cannot take the unconditional perspective too seriously or they cannot do statistics! — Berry (1987)

Other Subtle Problems with p-values

  • Proof by contradiction doesn’t directly apply outside of non-fuzzy logic
  • p is not P(results as impressive as ours | H0)
    It is P(more impressive)
  • Computed under assumption of zero treatment effect; doesn’t formally entertain harm
  • Can argue that H0 is always false
  • NHST entails fixing n; many studies stop with p=0.06 when could have added 20 more subjects
    • frequentist approach to sample size re-estimation requires discounting of data from first wave of subjects
  • See also Nuzzo (2014)

The following is almost but not quite the reasoning of null hypothesis rejection:

If the null hypothesis is correct, then this datum (D) can not occur.
It has, however, occurred. Therefore the null hypothesis is false.

If this were the reasoning of H0 testing, then it would be formally correct. … But this is not the reasoning of NHST. Instead, it makes this reasoning probabilistic, as follows:

If the null hypothesis is correct, then these data are highly unlikely.
These data have occurred.
Therefore, the null hypothesis is highly unlikely.

By making it probabilistic, it becomes invalid. … the syllogism becomes formally incorrect and leads to a conclusion that is not sensible:

If a person is an American, then he is probably not a member of Congress. (TRUE, RIGHT?)
This person is a member of Congress.
Therefore, he is probably not an American. (Pollard & Richardson, 1987)
.. The illusion of attaining improbability or the Bayesian Id’s wishful thinking error …


Induction has long been a problem in the philosophy of science. Meehl (1990) attributed to the distinguished philosopher Morris Raphael Cohen the saying “All logic texts are divided into two parts. In the first part, on deductive logic, the fallacies are explained; in the second part, on inductive logic, they are committed.” — Cohen (1994)

A person is interested in a probability model. But guided by the philosophy of p-values, he asks no questions about this model, and instead asks what is the probability, given the data and some other model, which is not the model of interest, of seeing an ad hoc statistic larger than some value. (Any change in a model produces a different model.) Since there are an infinite number of models that are not the model of interest, and since there are an infinite number of statistics, the creation of p-values can go on forever. Yet none have anything to say about the model of interest.

Why? Fisher (1970) said: “Belief in null hypothesis as an accurate representation of the population sampled is confronted by a logical disjunction: Either the null is false, or the p-value has attained by chance an exceptionally low value.”

Fisher’s “logical disjunction” is evidently not one, since the either-or describes different propositions. A real disjunction can however be found: Either the null is false and we see a small p-value, or the null is true and we see a small p-value. Or just: Either the null is true or it is false and we see a small p-value. Since “Either the null is true or it is false”” is a tautology, and is therefore necessarily true, we are left with, “We see a small p-value.” The p-value casts no light on the truth or falsity of the null.

Frequentist theory claims, assuming the truth of the null, we can equally likely see any p-value whatsoever. And since we always do (see any value), all p-values are logically evidence for the null and not against it. Yet practice insists small p-value is evidence the null is (likely) false. That is because people argue: For most small p-values I have seen in the past, the null has been false; I now see a new small p-value, therefore the null hypothesis in this new problem is likely false. That argument works, but it has no place in frequentist theory (which anyway has innumerable other difficulties).

Any use of p-values in deciding model truth thus involves a fallacy or misunderstanding. This is formally proven by Briggs (2016, chap. 9), a work which I draw from to suggest a replacement for p-values, which is this. Clients ask, “What’s the probability that if I know X, Y will be true?” Instead of telling them that, we give them p-values. — Briggs (2017)

The nil hypothesis is always false. Tukey (1991) wrote that “It is foolish to ask ‘Are the effects of A and B different?’ They are always different—for some decimal place”. Schmidt (1992) … reminded researchers that, given the fact that the nill hypothesis is always false, the rate of Type I errors is 0%, not 5%, and that only Type II errors can be made. — Cohen (1994)

(The use of p-values) is a lamentable demonstration of the credulity with which modern scientists will abandon biologic wisdom in favor of any quantitative ideology that offers the specious allure of a mathematical replacement for sensible thought. — Feinstein (1977)

It is incomparably more useful to have a plausible range for the value of a parameter than to know, with whatever degree of certitude, what single value is untenable. — Oakes (1986)

Issues with Confidence Limits

  • CLs have only long-term interpretations and give false impression that all values within the interval are equally likely
  • Cannot control which interval for which you want a probability statement

I see that the 0.95 confidence interval for the mean blood pressure difference is [2,7]. But I want to know the confidence I should have in it being in the interval [0,5] and you’re telling me it can’t be computed with frequentist confidence intervals? — Wagenmakers et al. (2017)

The worry is that, when data are weak and there is strong prior information that is not being used, classical methods can give answers that are not just wrong—that’s no dealbreaker, it’s accepted in statistics that any method will occasionally give wrong answers—but clearly wrong; wrong not only just conditional on the unknown parameter but also conditional on the data. Scientifically inappropriate conclusions. That’s the meaning of ‘poor calibration.’ Even this, in some sense, should not be a problem—after all, if a method gives you a conclusion that you know is wrong, you can just set it aside, right?—but, unfortunately, many users of statistics consider to take p < 0.05 or p < 0.01 comparisons as ‘statistically significant’ and to use these as a motivation to accept their favored alternative hypothesis. This has led to such farces, in recent claims, in leading psychology journals that various small experiments have demonstrated the existence of extra-sensory perception or huge correlations between menstrual cycle and voting, and so on. — Gelman & Hennig (2017)

So what happened with the development of efficacy measures is we developed a whole new field called biostatistics. It had been sort of an orphan corner of mathematics until the Kefauver-Harris Amendments, and there had been extremely important advances in how do you study efficacy of drugs. Most of it devolves down to whether or not you’re likely to see a benefit more than chance alone would predict. But how likely and how much benefit was left for some free-floating kind of notion by the FDA. So any benefit in essence, more than any toxicity in essence, would lead to licensure. That has led to what I call “small effectology.”

Nortin Hadler, Interviewed by Tom Ashbrook On Point, WBUR radio, 2016-03-29, 15:26

p-value: the chance that someone else’s data are more extreme than mine if H0 is true, not the chance that H0 is true given my data

Aside from ignoring applicable pre-study data, the p-value is at least monotonically related to what we need. But it is not calibrated to be on a scale meant for optimum decision making.

The criterion of p < .05 says that we should be willing to tolerate a 5% false alarm rate in decisions to reject the null value. In general, frequentist decision rules are driven by a desire to limit the probability of false alarms. The probability of false alarm (i.e., the p value) is based on the set of all possible test results that might be obtained by sampling fictitious data from a particular null hypothesis in a particular way (such as with fixed sample size or for fixed duration) and examining a particular suite of tests (such as various contrasts among groups). Because of the focus on false alarm rates, frequentist practice is replete with methods for adjusting decision thresholds for different suites of intended tests. …

Bayesian decisions are not based on false alarm rates from counterfactual sampling distributions of hyopthetical data. Instead, Bayesian decisions are based on the posterior distribution from the actual data. — Kruschke & Liddell (2017)

… Neyman and Pearson outline the price that must be paid to enjoy the purported benefits of objectivity: We must abandon our ability to measure evidence, or judge truth, in an individual experiment. … Hypothesis tests are equivalent to a system of justice that is not concerned with which individual defendent is found guity or innocent (that is , “whether each separate hypothesis is true or false”) but tries instead to control the overall number of incorrect verdicts (that is, “in the long run of experience, we shall not often be wrong”). Controlling mistakes in the long run is a laudable goal, but just as our sense of justice demands that individual persons be correctly judged, scientific intuition says that we should try to draw the proper conclusions from individual studies. — Goodman (1999)

3.2.1 Computing p-values Using Simulation

  • Simulations help one to grasp theory
  • One-sample problem for a normal mean μ
  • Single-arm study, μ > 0 denotes efficacy
  • Assume σ = 1
  • n=30, true μ = 0.3
  • Simulate 100,000 studies
  • Also compute P(result approx. as impressive as observed | μ = 0)
    within window of width 0.1
Code
n <- 30
set.seed(1)
y     <- rnorm(n, 0.3, sd=1)            # generate data
ybar  <- mean(y)                        # observed mean
ucl   <- ybar + qnorm(0.95) / sqrt(30)  # upper C.L.
# Run 100,000 studies and compute their sample means:
repeated.ybar <- rnorm(100000, 0, sd=sqrt(1/30))
# TRUE/FALSE variables are converted to 1/0 when taking the mean
# This is an easy way to compute a proportion
p     <- mean(repeated.ybar >= ybar)
pa    <- mean(repeated.ybar >= ybar & repeated.ybar <= ybar + 0.1)
repeated.ucl <- repeated.ybar + qnorm(0.95) / sqrt(30)
cover <- mean(repeated.ucl >= 0)
cat('Observed mean           : ', round(ybar, 3), '\n',
    'Upper 0.95 1-sided CL   : ', round(ucl, 3), '\n',
    'One-sided p-value       : ', round(p, 4), '\n',
    'Exact p-value           : ', round(1 - pnorm(ybar, 0, sd=1/sqrt(30)),
                                  4), '\n',
    'Confidence coverage     : ', round(cover, 4), '\n',
    'P(Approx. as impressive): ', round(pa, 4), '\n',
    sep='')
Observed mean           : 0.382
Upper 0.95 1-sided CL   : 0.683
One-sided p-value       : 0.0186
Exact p-value           : 0.0181
Confidence coverage     : 0.949
P(Approx. as impressive): 0.0146
  • Modify simulation for 2 data looks, with stopping rule
  • First look after n=15 with nominal type I error 0.05
  • Stop study if mean exceeds corresponding cutoff
  • Otherwise use mean of n=30 as final estimate and basis for test
  • Compute two ‘nominal’ p-values and then the actual p-value under the stopping rule
  • Assumptions exposed
    • intended look is actually carried out
    • look is ignored if p > 0.05
Code
set.seed(1)
# Make first look
y1 <- rnorm(n / 2, 0.3, sd=1)
ybar1 <- mean(y1)
# Make second look
y2 <- rnorm(n / 2, 0.3, sd=1)
ybar2 <- mean(c(y1, y2))      # combine to get n=30
ybar.at.stop  <- ifelse(ybar1 * sqrt(15) >= qnorm(0.95), ybar1, ybar2)
ybar.at.stopb <- ifelse(ybar1 * sqrt(15) >= qnorm(0.975),ybar1, ybar2)
# Run 100,000 studies.  For each get mean with n=15 and 30 and apply the same stopping rule
repeated.ybar1 <- rnorm(100000, 0, sd=sqrt(1/15))
# Compute overall mean with n=30:
repeated.ybar2 <- (repeated.ybar1 + rnorm(100000, 0, sd=sqrt(1/15))) / 2
# Compute estimate of mean at last look, using 2 stopping rules
repeated.ybar  <- ifelse(repeated.ybar1 * sqrt(15) >= qnorm(0.95),
                         repeated.ybar1, repeated.ybar2)
repeated.ybarb <- ifelse(repeated.ybar1 * sqrt(15) >= qnorm(0.975),
                         repeated.ybar1, repeated.ybar2)
pval1          <- mean(repeated.ybar1 >= ybar1)  # ordinary p-value n=15
pval2          <- mean(repeated.ybar2 >= ybar2)  # ordinary p-value n=30
# P-value accounting for multiple looks using alpha=0.05
pval           <- mean(repeated.ybar  >= ybar.at.stop)
# Same using alpha=0.025
pvalb          <- mean(repeated.ybarb >= ybar.at.stopb)

cat('Sample mean at first look     :', round(ybar1, 3), '\n',
    'Sample mean at end            :', round(ybar2, 3), '\n',
    'Nominal p-value at n=15       :', round(pval1, 4), '\n',
    'Nominal p-value at n=30       :', round(pval2, 4), '\n',
    'p-value accounting for looks  :', round(pval,  4), '\n',
    'p-value " " with alpha=0.025  :', round(pvalb, 4), '\n', sep='')
Sample mean at first look     :0.401
Sample mean at end            :0.382
Nominal p-value at n=15       :0.0604
Nominal p-value at n=30       :0.018
p-value accounting for looks  :0.0585
p-value " " with alpha=0.025  :0.0367
  • True p-value accounting for 2 looks > simple p-value
  • True p-value smaller if used instead α = 0.025
  • Sampling distribution for the sample mean under our 0.05 stopping rule has a discontinuity at the value of the mean that would cause early stopping
  • Difficult to derive and use the true sampling distribution under multiple looks; most statisticians pretend only 1 look done
  • Bayesian inference doesn’t concern itself with sampling distributions
Code
spar(bty='l')
hist(repeated.ybar, nclass=100, xlab=expression(bar(Y)), main='')
abline(v=qnorm(0.95) / sqrt(15), col=gray(.85))

Sampling distribution of final estimate of the mean in a two-stage sequential single arm trial, under the null hypothesis

3.3 Bayesian

This form (probability of unknown given what is known) has enormous benefits. It is in plain language; specialized training is not needed to grasp model statements … Everything is put in terms of observables. The model is also made prominent, in the sense that it is plain there is a specific probability model with definite assumptions in use, and thus it is clear that answers will be different if a different model or different assumptions about that model are used … — Briggs (2017)

  • Relative changes in evidence are functions only of data
  • No absolute truths
  • Final evidence quantified on an absolute scale given pre-data anchor
  • Bayes’ theorem: movement of prior belief to current belief
  • Full conditioning on observables
  • Does not condition on unknowables
  • Probability statements forward in time/information so have meaning out of context
  • Multiple looks/stopping rule not relevant
  • Helpful in understanding how a Bayesian might cheat:
    • change prior after seeing data
    • hiding data, e.g. P(efficacy) = 0.95, enroll more subjects, P(efficacy) = 0.93, report previous look as final

Example Frequentist Result

  • Difference in mean SBP between treatments A, B = 6mmHg
  • p=0.01, 0.95 CL [3,9]
  • An event (?) (of ≥ 6 mmHg) of low probability has been witnessed if A=B

Corresponding Bayesian Result

  • Use a normal prior satisfying
    • pre-study chance of worsening BP is 0.5
    • pre-study chance of a large (≥ 10mmHg) improvement in BP is 0.1
  • Posterior mean BP reduction 5mmHg
  • 0.95 credible interval: [2.5, 8]
  • P(reduction in BP) = 0.97
  • P(reduction ≥ 2mmHg) = 0.9
  • P(similarity) = P(|difference| ≤ 2mmHg) easy to compute

Updating of Posterior in Sequential Trials

  • Coin flipping, prior is beta(10,10) favoring fairness
  • 100 tosses, update posterior every 10 tosses
  • After n tosses with y heads, posterior is beta(Y + 10, n - y + 10)
  • Repeat for beta(5,5) and beta(20,20) priors
  • Click on the legend to hide or display results for the 3 priors
    Double-click to show only one set
Code
require(plotly)
x <- seq(0, 1, length=200)
p <- plot_ly(width=800, height=500)
for(ab in c(10, 5, 20)) {
    set.seed(1) 
    alpha <- beta <- ab
    lg  <- paste0('α=', alpha, ' β=', beta)
    vis <- if(ab == 10) TRUE else 'legendonly'
    Y   <- 0  
    # Plot beta distribution density function
    p <- p %>% add_lines(x=~x, y=~y, hoverinfo='none', visible=vis,
        line=list(color='blue'),
        name=paste0('Prior:', lg), legendgroup=lg,
        data=data.frame(x=x, y=dbeta(x, alpha, beta)))

    for(N in seq(10, 100, by=10)) {
        Y <- rbinom(1, 10, 0.5)  # 10 new tosses
        # Posterior distribution updated
        alpha <- alpha + Y
        beta  <- beta + 10 - Y
        p <- p %>% add_lines(x=~x, y=~y, hoverinfo='none', visible=vis,
            name=paste0('N=', N, ' ', lg), legendgroup=lg,
            line=list(color=if(N < 100) 'black' else 'red',
                opacity=if(N < 100) .95 - N / 120 else 1,
                width=N * 2 / 100),
            showlegend=FALSE,
            data=data.frame(x=x, y=dbeta(x, alpha, beta)))
      }
    }
p %>% layout(shapes=list(type='line', x0=.5, x1=.5, y0=0, y1=9, opacity=.2),
        xaxis=list(title='θ'), yaxis=list(title=''))

Prior distribution (blue) and posterior distributions as the trials progress (darkness of lines increases). The final posterior at N=100 is in red, when there were 53 heads tossed.

3.3.1 Alternative Take on the Prior

  • Can also think of the prior as augmenting or reducing the effective sample size (Wiesenfarth & Calderazzo (2019))
  • Consider a skeptical prior
  • Effectively ignoring some of the sample
  • When the data are normally distributed with mean \(\mu\) and variance \(\sigma^2\), \(\overline{Y}\) is normal with mean \(\mu\) and variance \(\frac{\sigma^{2}}{n}\)
  • When the prior is also Gaussian but with mean \(\mu_{0}\) and variance \(\sigma^{2}_{0}\), the posterior distribution of \(\mu\) is normal with the following variance and mean, respectively, if we let the precision of the prior \(\tau = \frac{1}{\sigma^{2}_{0}}\):
    • \(\sigma^{2}_{1} = (\tau + \frac{n}{\sigma^{2}})^{-1}\)
    • \(\mu_{1} = \sigma^{2}_{1} (\mu_{0} \tau + \frac{n \overline{Y}}{\sigma^{2}})\)
  • Consider special case where the data variance \(\sigma^2\) is 1 and the prior mean is zero; then the posterior variance and mean of \(\mu\) are:
    • \(\sigma^{2}_{1} = (n + \tau)^{-1}\)
    • \(\mu_{1} = \frac{n \overline{Y}}{n + \tau}\)
  • From the formulas for posterior mean and variance, the effect of the prior with variance \(\sigma^{2}_{0} = \frac{1}{\tau}\) compared to no discounting (flat prior; \(\tau=0\)) is \(\tau\) observations in a certain sense
  • But we need to see how the difference in posterior mean and variance combine to change the posterior probability and to recognize that the amount of discounting is sample size-dependent
  • For a sample size \(m\), posterior \(P(\mu > 0) = \Phi(\frac{m\overline{Y}}{\sqrt{m + \tau}})\) where \(\Phi\) is the standard normal CDF
  • Compare to no discounting (\(\tau=0\)) with a different sample size \(n\): \(P(\mu > 0) = \Phi(\sqrt{n}\overline{Y})\)
  • What is the sample size \(n\) for an undiscounted analysis giving same \(P(\mu > 0)\) as the discounted one?
  • Set \(\frac{m}{\sqrt{m + \tau}} = \sqrt{n}\) so
    \(m = \frac{n + \sqrt{n^{2} + 4n\tau}}{2}\)
  • Increase in sample size needed to overcome skepticism: \(n - m\)
  • In figure below the prior prob(\(\mu > 1\)) is also shown. This is \(1 - \Phi(\sqrt{\tau})\)
Code
spar(bty='l')
z <- list()
n <- seq(1, 100, by=2)
for(tau in c(0, 2, 5, 10))
  z[[paste0('τ=', tau, '     P(μ>1)=',
            format(1 - pnorm(sqrt(tau)), digits=3, scientific=1))]] <-
    list(x=n, y=0.5 * (n + sqrt(n^2 + 4 * n * tau)) - n)

labcurve(z, pl=TRUE, xlab='Sample Size With No Skepticism',
         ylab='Extra Subjects Needed Due to Skepticism', adj=1)

Effect of discounting by a skeptical prior with mean zero and precision τ: the increase needed in the sample size in order to achieve the same posterior probability of μ > 0 as with the flat (non-informative) prior. τ=10 corresponds to a very skeptical prior, giving almost no chance to a large μ (μ > 1).

  • Optimistic prior: effectively adds observations

3.4 Contrasting Frequentist and Bayesian Evidence and Errors

  • Let E=true unknown efficacy measure
    • what is generating the difference in effects in the data
    • difference in true means, log odds ratio, log hazard ratio, etc.
    • E=0: no treatment effect
    • E>0: benefit of new treatment
  • Frequentist:
    • attempt to show data implausible if E=0
    • no probability statement about E; E is either 0 or nonzero
  • Bayesian:
    • probability statement about E using posterior (“current”) probs
    • E almost always thought of as continuous (P(E=0) = 0)
    • P(E > c | data)
    • c=0: get evidence for any efficacy
    • c>0: get evidence for efficacy > some amount
    • There are no “errors”
    • Errors can only be made by decision makers when actions constrained to all-or-nothing

3.4.1 Frequentist vs Bayes: Study Design

Frequentist

  • Design study to have α=0.05 β=0.1
  • Once data available, these no longer relevant since they apply to sequences of other trials, not this trial
  • α depends on intentions, β on a single value of E
  • Can also design for precision: solve for n such that 0.95 CL expected to have width w

Bayesian

  • Choose prior for E allowing for uncertainty in true effect
  • Design study to have prob ≥ 0.9 of achieving P(E > c) > 0.95 or to achieve credible interval width w

3.4.2 Type of Errors

  • Need to be careful about the use of the term ‘error’, as α is not the probability of making an error
  • α is a conditional probability of making an assertion of an effect when any such assertion is by definition wrong
  • It is a trigger/alarm probability
  • α is conditional on H0 and ignores all data

Frequentist

  • Type I assertion probability α: P(declare efficacy when E=0) = P(test stat > threshold when E=0)
  • Type II assertion probability: P(fail to declare efficacy when E=c for some particular arbitrary c)
  • α never drops no matter how large is n

Bayesian

  • P(E > c | data)
  • If act as if efficacious, P(error) is 1 - PP
  • If act as if ineffective, P(error) is PP

3.4.3 Example: p=0.03

Frequentist

  • Conclude efficacy
  • This is either right or wrong; no prob associated with true unknown E
  • Exact interpretation: if E=0 and one ran an infinite sequence of identical trials, one would see an observed E ≥ that observed 0.03 of the time

Bayesian

  • PP is its own error probability

3.4.4 Example: p=0.2

Frequentist

  • Can’t conlude E=0 but fail to have evidence for E≠0
  • No measure for P(E=0) available

Bayesian

  • Simple PP of no effect or harm: P(E < 0)

3.4.5 Clinical Significance

Frequentist

  • With large n, trivial effect can yield p < 0.05

Bayesian

  • Compute PP that true effect more than trivial

3.4.6 p=0.04, 5 other trials “negative”

Frequentist

  • No way to take the other 5 trials into account other that using non-quantitative subjective arguments

Bayesian

  • Skepticism about efficacy for current trial already captured in the prior
  • Or other trials could be used to form a prior, or Bayesian hierarchical model

3.5 Problems Caused by Use of Arbitrary Thresholds

  • Thresholding in general is arbitrary and detrimental to asessing totality of evidence
    • true for both frequentist and Bayesian
  • Leads to false confidence
    • once known that evidence measure below or above threshold, stakeholders act as if no uncertainty (Goodman (1999), Altman & Bland (1995), Greenwald et al. (1996))
  • Example honest sentence that is harder to take out of context:
    Treatment B probably (0.94) resulted in lower BP and was probably (0.78) safer than treatment A