Meta Analysis
Meta Analysis
Alessio Crippa
Department of Medical Epidemiology and Biostatistics
Karolinska Institutet
Systematic review Fixed-effect model Random-effects model Heterogeneity Model choice Sensitivity analysis
Systematic review
Fixed-effect model
Random-effects model
Heterogeneity
Model choice
Sensitivity analysis
Crippa Alessio
1
Systematic review Fixed-effect model Random-effects model Heterogeneity Model choice Sensitivity analysis
Systematic review
Crippa Alessio
2
Systematic review Fixed-effect model Random-effects model Heterogeneity Model choice Sensitivity analysis
Systematic review
Crippa Alessio
3
Systematic review Fixed-effect model Random-effects model Heterogeneity Model choice Sensitivity analysis
Crippa Alessio
4
Systematic review Fixed-effect model Random-effects model Heterogeneity Model choice Sensitivity analysis
Crippa Alessio
5
Systematic review Fixed-effect model Random-effects model Heterogeneity Model choice Sensitivity analysis
Crippa Alessio
6
Systematic review Fixed-effect model Random-effects model Heterogeneity Model choice Sensitivity analysis
The choice between the two models is critical. It doesn’t affect only
the computations, but also the goal of the analyses and the
interpretations of the results.
Crippa Alessio
7
Systematic review Fixed-effect model Random-effects model Heterogeneity Model choice Sensitivity analysis
Fixed-effect model
Crippa Alessio
8
Systematic review Fixed-effect model Random-effects model Heterogeneity Model choice Sensitivity analysis
Mean score
Crippa Alessio
9
Systematic review Fixed-effect model Random-effects model Heterogeneity Model choice Sensitivity analysis
Fixed-effect model
All the studies share a common (true and unknown) effect size,
denoted by θ. Effect sizes might differ because of random variability.
Yi = θ + i
Crippa Alessio
10
Systematic review Fixed-effect model Random-effects model Heterogeneity Model choice Sensitivity analysis
Study E
Study D
Study C
Study B
Study A
True effect
Crippa Alessio
11
Systematic review Fixed-effect model Random-effects model Heterogeneity Model choice Sensitivity analysis
1 1
Wi = =
VYi SE(Yi )2
VYi is the within-study variance for study, e.g. an estimate of the
(square) standard error of the observed mean score in the i-th study.
Crippa Alessio
12
Systematic review Fixed-effect model Random-effects model Heterogeneity Model choice Sensitivity analysis
Point estimation
Pk
Wi Yi
M = Pi=1
k
i=1 Wi
the sum of the products Wi Yi (effect size multiplied by weight)
divided by the sum of the weights.
Crippa Alessio
13
Systematic review Fixed-effect model Random-effects model Heterogeneity Model choice Sensitivity analysis
1
VM = Pk
i=1 Wi
p
SEM = VM
M ± 1.96 × SEM
Crippa Alessio
14
Systematic review Fixed-effect model Random-effects model Heterogeneity Model choice Sensitivity analysis
Hand calculations
study y se v w wy std_w
1 Study A 101.3 1.532 2.347 0.426 43.2 0.111
2 Study B 99.6 1.511 2.283 0.438 43.6 0.114
3 Study C 101.3 1.528 2.336 0.428 43.3 0.111
4 Study D 99.8 1.419 2.013 0.497 49.6 0.129
5 Study E 100.9 0.696 0.485 2.063 208.2 0.536
Crippa Alessio
15
Systematic review Fixed-effect model Random-effects model Heterogeneity Model choice Sensitivity analysis
101.301 · 0.426 + 99.556 · 0.438 + 101.257 · 0.428 + 99.775 · 0.497 + 100.906 · 2.063
M= =
0.426 + 0.438 + 0.428 + 0.497 + 2.063
1 1
VM = = = 0.26
0.426 + 0.438 + 0.428 + 0.497 + 2.063 3.852
√
M ± 1.96 × SEM = 100.7 ± 1.96 · 0.26 = (99.7, 101.7)
Crippa Alessio
16
Systematic review Fixed-effect model Random-effects model Heterogeneity Model choice Sensitivity analysis
Hypothesis testing
A Z -value test can be used for testing the null hypothesis that the
true efftc is equal to zero H0 : θ = 0 vs H1 : θ 6= 0
M
Z=
SEM
A two-sided p-value is obtained from the cumulative standard
normal function
Crippa Alessio
17
Systematic review Fixed-effect model Random-effects model Heterogeneity Model choice Sensitivity analysis
Crippa Alessio
18
Systematic review Fixed-effect model Random-effects model Heterogeneity Model choice Sensitivity analysis
Random-effects model
Crippa Alessio
19
Systematic review Fixed-effect model Random-effects model Heterogeneity Model choice Sensitivity analysis
90 95 100 105
Mean score
Crippa Alessio
20
Systematic review Fixed-effect model Random-effects model Heterogeneity Model choice Sensitivity analysis
The assumption of a common effect size for all the studies is often
implausible.
Crippa Alessio
21
Systematic review Fixed-effect model Random-effects model Heterogeneity Model choice Sensitivity analysis
Between-studies variation
Study E
Study D
ζC
Study C
Study B
Study A
True effect
96 100 104
Mean score
Crippa Alessio
22
Systematic review Fixed-effect model Random-effects model Heterogeneity Model choice Sensitivity analysis
Within-study variation
Study E
Study D
ζC
Study C εC
Study B
Study A
True effect
96 100 104
Mean score
Crippa Alessio
23
Systematic review Fixed-effect model Random-effects model Heterogeneity Model choice Sensitivity analysis
Random-effects model
Yi = µ + ζi + i
Crippa Alessio
24
Systematic review Fixed-effect model Random-effects model Heterogeneity Model choice Sensitivity analysis
Key points
Crippa Alessio
25
Systematic review Fixed-effect model Random-effects model Heterogeneity Model choice Sensitivity analysis
Study E
Study D
Study C
Study B
Study A
True effect
95 100 105
Crippa Alessio
26
Systematic review Fixed-effect model Random-effects model Heterogeneity Model choice Sensitivity analysis
Crippa Alessio
27
Systematic review Fixed-effect model Random-effects model Heterogeneity Model choice Sensitivity analysis
Q − (k − 1)
T 2 = max Pk ,0
Pk Wi2
i=1 Wi − i=1
P k
i=1
Wi
P 2
k k
X i=1 W Y
i i
Q= Wi Yi2 − Pk
i=1 i=1 Wi
Crippa Alessio
28
Systematic review Fixed-effect model Random-effects model Heterogeneity Model choice Sensitivity analysis
Point estimation
Pk
∗ Wi∗ Yi
M = Pi=1
k ∗
i=1 Wi
where the weights Wi∗ incorporate the between-studies variance:
1 1
Wi∗ = =
VY∗i VYi + T 2
Crippa Alessio
29
Systematic review Fixed-effect model Random-effects model Heterogeneity Model choice Sensitivity analysis
1
VM ∗ = Pk ∗
i=1 Wi
95% confidence interval for the summary effect can be calculated
using the normal approximation
∗
M ± 1.96 × SEM
A test for the the null hypothesis that the mean effect µ is equal to
0 can be obtained as in the fixed-effect meta-analysis.
Crippa Alessio
30
Systematic review Fixed-effect model Random-effects model Heterogeneity Model choice Sensitivity analysis
Hand calculations
study y se v w wy y2 wy2 w*
1 Study A 94.8 1.403 1.967 0.508 48.2 8993 4572 0.221
2 Study B 98.7 1.304 1.701 0.588 58.0 9741 5727 0.234
3 Study C 101.1 1.489 2.218 0.451 45.6 10212 4603 0.209
4 Study D 100.0 1.349 1.820 0.550 55.0 9998 5494 0.228
5 Study E 98.2 0.715 0.511 1.955 191.9 9635 18841 0.325
Crippa Alessio
31
Systematic review Fixed-effect model Random-effects model Heterogeneity Model choice Sensitivity analysis
k
X
Wi Yi2 =8993.249 · 0.508 + 9740.95 · 0.588 + 10211.553 · 0.451+
i=1
398.682
Q = 39237.17 −
4.05
11.272 − (5 − 1)
T2 = 4.93
= 2.566
4.05 − 4.05
Crippa Alessio
32
Systematic review Fixed-effect model Random-effects model Heterogeneity Model choice Sensitivity analysis
Distributional assumptions
i ∼ N(0, VYi ) or θi ∼ N(θ, VYi )
I the moment-based estimator of τ 2 does not require additional
assumption about the distribution of ζi ;
I the pooled effect M ∗ and its SE are also independent on a
distributional assumption for the random-effects;
I confidence intervals relies on approximating distributions or
central limit theorem (k >> 0).
ζi ∼ N(0, τ 2 )
I alternative estimator of τ 2 ;
I facilitates computation of confidence intervals;
I allows to give predictions for new studies.
Crippa Alessio
33
Systematic review Fixed-effect model Random-effects model Heterogeneity Model choice Sensitivity analysis
Crippa Alessio
34
Systematic review Fixed-effect model Random-effects model Heterogeneity Model choice Sensitivity analysis
Heterogeneity
Crippa Alessio
35
Systematic review Fixed-effect model Random-effects model Heterogeneity Model choice Sensitivity analysis
Crippa Alessio
36
Systematic review Fixed-effect model Random-effects model Heterogeneity Model choice Sensitivity analysis
Crippa Alessio
37
Systematic review Fixed-effect model Random-effects model Heterogeneity Model choice Sensitivity analysis
Crippa Alessio
38
Systematic review Fixed-effect model Random-effects model Heterogeneity Model choice Sensitivity analysis
Crippa Alessio
39
Systematic review Fixed-effect model Random-effects model Heterogeneity Model choice Sensitivity analysis
The Q statistic
k k
Yi − M 2
X X
2
Q= Wi (Yi − M) =
i=1 i=1
SYi
Crippa Alessio
40
Systematic review Fixed-effect model Random-effects model Heterogeneity Model choice Sensitivity analysis
Crippa Alessio
41
Systematic review Fixed-effect model Random-effects model Heterogeneity Model choice Sensitivity analysis
Crippa Alessio
42
Systematic review Fixed-effect model Random-effects model Heterogeneity Model choice Sensitivity analysis
Crippa Alessio
43
Systematic review Fixed-effect model Random-effects model Heterogeneity Model choice Sensitivity analysis
Estimating τ 2
τ 2 is the variance of the true effect sizes. It is in the same metric as
the effect and reflects the absolute amount of variation.
Der Simonian
andLaird proposed the estimator:
2 Q−df
T = max 0; C
Tau
Crippa Alessio
45
Systematic review Fixed-effect model Random-effects model Heterogeneity Model choice Sensitivity analysis
A measure of heterogeneity
Higgins et al (2002) proposed a statistic I 2 which quantifies the
amount of total variability attributed to between study
heterogeneity.
Q − df Variancebetween
I2 = × 100% = × 100%
Q Variancetotal
Advantages of I 2 are:
I intuitive interpretation, it is a percentage
I it is easy to calculate from previously published meta-analysis
I the interpretation is not dependent on the choice of measure of
association
Crippa Alessio
46
Systematic review Fixed-effect model Random-effects model Heterogeneity Model choice Sensitivity analysis
Crippa Alessio
47
Systematic review Fixed-effect model Random-effects model Heterogeneity Model choice Sensitivity analysis
Recap
Crippa Alessio
48
Systematic review Fixed-effect model Random-effects model Heterogeneity Model choice Sensitivity analysis
Crippa Alessio
49
Systematic review Fixed-effect model Random-effects model Heterogeneity Model choice Sensitivity analysis
Exploring Heterogeneity
I Exploring Heterogeneity
I do separate meta-analyses on subgroups of studies according to
study characteristics that may influence the pooled association
(gender, study design, geographical area, year of publication);
I compare the combined results with analogue to an ANOVA
model.
I Meta-regression
I similar to standard regression;
I the pooled effect size is modeled as a linear function of one or
more explanatory variables
Note: Both approaches may have low power, they require several
studies.
Crippa Alessio
50
Systematic review Fixed-effect model Random-effects model Heterogeneity Model choice Sensitivity analysis
Crippa Alessio
51
Systematic review Fixed-effect model Random-effects model Heterogeneity Model choice Sensitivity analysis
Model choice
Crippa Alessio
52
Systematic review Fixed-effect model Random-effects model Heterogeneity Model choice Sensitivity analysis
Crippa Alessio
53
Systematic review Fixed-effect model Random-effects model Heterogeneity Model choice Sensitivity analysis
It also 2
includes τ , so the weights will be smaller, and similarly their
sum Var(M) = Pk 1 .
i=1
Wi
Crippa Alessio
54
Systematic review Fixed-effect model Random-effects model Heterogeneity Model choice Sensitivity analysis
Differences in precision
Let us assume VYi = σ 2 and ni = n ∀i
1 1 σ2 σ2
VM = Pk = Pk ni
= Pk =
i=1 Wi i=1 σ 2 i=1 ni
kn
σ2 T 2
VM ∗ = +
kn k
I standard deviation of the true effect sizes (higher T 2 , less
precise);
I number of studies (higher k, more precise).
Crippa Alessio
55
Systematic review Fixed-effect model Random-effects model Heterogeneity Model choice Sensitivity analysis
Crippa Alessio
56
Systematic review Fixed-effect model Random-effects model Heterogeneity Model choice Sensitivity analysis
Fixed-effect model
- there is good reason to believe that all the studies are functionally
identical;
- the goal is to compute the common effect size (not be generalized
beyond the observed population).
Random-effects model
- the studies have enough in common but are not identical;
- the goal of the analysis is usually to generalize to a range of
scenarios.
Crippa Alessio
57
Systematic review Fixed-effect model Random-effects model Heterogeneity Model choice Sensitivity analysis
Sensitivity analysis
Crippa Alessio
58
Systematic review Fixed-effect model Random-effects model Heterogeneity Model choice Sensitivity analysis
Crippa Alessio
59
Systematic review Fixed-effect model Random-effects model Heterogeneity Model choice Sensitivity analysis
Publication Bias
Precision increases as study size increases.
Results from small studies will be more scattered than larger studies.
Crippa Alessio
60
Systematic review Fixed-effect model Random-effects model Heterogeneity Model choice Sensitivity analysis
Funnel plot
0
Standard Error
0.367
0.735
Egger’s test
Inspection of the funnel plot might be subjective. Egger (BMJ,
1997) and colleagues proposed a simple test to detect asymmetry in
the funnel plot:
One can view funnel plots as displaying evidence for “small study
effects” in general rather than publication bias in particular.
Crippa Alessio
63
Systematic review Fixed-effect model Random-effects model Heterogeneity Model choice Sensitivity analysis
Trim-and-fill analysis
Crippa Alessio
64
Systematic review Fixed-effect model Random-effects model Heterogeneity Model choice Sensitivity analysis
0
Standard Error
0.367
0.735
Cumulative meta-analysis
It allows
I the study of trends in efficacy;
I to determine when a new treatment appears to be significantly
effective or deleterious.
Crippa Alessio
66