12
Stata Technical Bulletin
STB-57
2. Let I be the current step number. Remove fcθ l'' values from the right end of the original Yi and estimate Δ^ based on
this trimmed set of n — fc∩ -1^ values: {Yi,..., Y ɔ/-1 ;}. Construct the next set of centered values
v v , , п—Kq 'j
γV=Yi-∆W, i = l,...,n
and estimate feθŋ using the chosen estimator for k0 applied to the set of values Yil∖
3. Increment I and repeat step 2 until an iteration L where = fcθb-1∖ Assign this common value to be the estimated
value fc0. Note that in this iteration it will also be true that ʌ(ɪə = Δ^l-1∖
4. Augment (that is, “fill”) the dataset Y with the ко imputed symmetric values
y∕ = 2∆W-yn.j+1, j = l,...,k0
and imputed standard errors
σ*=σn-j+ι, j = l,...,k0
Estimate the “trimmed and filled” value of Δ using the chosen meta-analysis method applied to the full augmented dataset
{y1,...,yn,y1*,...,κ≈}.
Conceptually, this algorithm starts with the observed data, iteratively trims (that is, removes) extreme positive studies from
the dataset until the remaining studies do not show detectable deviation from symmetry, fills (that is, imputes into the original
dataset) studies that are left-side mirrored reflections (about the center of the trimmed data) of the trimmed studies and, finally,
repeats the meta-analysis on the filled dataset to get “trimmed and filled” estimates. Each filled study is assigned the same
standard error as the trimmed study it reflects in order to maintain symmetry within the filled dataset.
Example
The method is illustrated with an example from the literature that examines the association between Chlamydia trachomatis
and oral contraceptive use derived from 29 case-control studies (Cottingham and Hunter 1992). Analysis of these data with the
publication bias tests of Begg and Mazumdar (p = 0.115) and Egger et al. (p = 0.016), as provided in metabias, suggests
that publication bias may affect the data. To examine the potential impact of publication bias on the interpretation of the data,
metatrim is invoked as follows:
. metatrim Iogor varlogor, reffeet funnel var
The random-effects model and display of the optional funnel graph are requested via options ref feet and funnel. Option var
is required because the data were provided as log-odds ratios and variances. By default, the linear estimator, L0, is used to
estimate fc0, as no other estimator was requested. metatrim provides the following output:
Note: option "var" I Pooled Method I Est |
specified. 987. CI |
Asymptotic |
No. of | |||
Fixed I Random I |
0.655 0.716 |
0.871 0.898 |
0.738 0.837 |
18.389 11.694 |
0.000 0.000 |
— 29 |
Test for heterogeneity: Q= 37.034 on 28 degrees of freedom (p= 0.118)
Moment-based estimate of between studies variance = 0.021
Trimming estimator: Linear
Meta-analysis type: Random-effects model
iteration I |
estimate |
Tn |
# to trim |
diff |
----------+- 1 I |
0.716 |
288 |
S |
— 438 |
2 I |
0.673 |
308 |
6 |
40 |
3 I |
0.660 |
313 |
7 |
16 |
4 I |
0.646 |
320 |
7 |
14 |
8 I |
0.646 |
320 |
7 |
0 |
Filled
Meta-analysis
I Method I |
Pooled Est |
987. |
CI Upper |
Asymptotic |
No. of | |
Fixed I Random I |
0.624 0.6SS |
0.842 0.831 |
0.708 0.779 |
14.969 10.374 |
0.000 0.000 |
— 36 |
Test for heterogeneity: Q= 49.412 on 35 degrees of freedom (p= 0.054)
Moment-based estimate of between studies variance = 0.031