[R-meta] robust variance estimator in meta-analyses of rare events (proportions)

[Guido Schwarzer]

Am 11.10.18 um 16:20 schrieb Pier-Alexandre Tardif:

> Hello again,
>
> Thanks for the follow up. I have compared the two methods (inverse variance with double arcsine vs mixed-effects logistic) and I still have one question.
>
> 1) Summary proportions obtained with either model are really low and relatively similar:
> 1st model (inverse variance, double arcsine): 0.01306 [95% CI: 0.00561-0.02271]
> 2nd model (mixed-effects logistic model): 0.0167 [95% CI: 0.0099-0.0281])
>
> In terms of 'presentation', would it make sense to rescale these proportions and how could we then interpret them?

Yes, you can rescale them. Actually, R function metaprop() from meta has 
an argument 'pscale' for this:

library(meta)
print(summary(metaprop(3:4, c(5678, 1234), pscale = 10000)), digits = 1)
Number of studies combined: k = 2

                      events      95%-CI  z p-value
Fixed effect model     14.9 [7.1; 31.2] --      --
Random effects model   13.4 [2.3; 78.9] --      --

Quantifying heterogeneity:
tau^2 = 1.3581; H = 2.38; I^2 = 82.3% [25.5%; 95.8%]

Test of heterogeneity:
     Q d.f. p-value
  5.65    1  0.0175

Details on meta-analytical method:
- Inverse variance method
- DerSimonian-Laird estimator for tau^2
- Logit transformation
- Clopper-Pearson confidence interval for individual studies
- Events per 10000 observations

As you can see in the last line, results are expressed as events per 
10000 observations (for pscale = 10000).

Best wishes,
Guido

	[[alternative HTML version deleted]]