Я заинтересован в объединении оценок выживаемости сердечно-сосудистых заболеваний в возрасте 70 лет. Поскольку это оценки выживаемости, они варьируются от 0 до 1 (аналогично пропорции).У меня есть 5 исследований, и резюме оценок, 95% CI, n и SE показаны ниже.Каждая строка представляет собой исследование.
> dat
Estimate Lower Upper n SE
1 0.55 0.40 0.71 100 1.479592
2 0.23 0.15 0.35 300 2.562728
3 0.54 0.44 0.66 200 2.092459
4 0.59 0.30 0.75 400 2.959184
5 0.88 0.67 0.98 40 0.935776
dat <- structure(list(Estimate = c(0.55, 0.23, 0.54, 0.59, 0.88), Lower = c(0.4,
0.15, 0.44, 0.3, 0.67), Upper = c(0.71, 0.35, 0.66, 0.75, 0.98
), n = c(100, 300, 200, 400, 40), SE = c(1.47959183673469, 2.56272823568864,
2.09245884228672, 2.95918367346939, 0.935776042294725)), .Names = c("Estimate",
"Lower", "Upper", "n", "SE"), row.names = c(NA, -5L), class = "data.frame")
rma в пакете metafor позволяет мне вводить sei (стандартная ошибка), которая инкапсулирует информацию, которую я получаю из 95% КИ на основеисследования, но полученный доверительный интервал не ограничен между 0 и 1 (т.е. CI составляет -0,51, 1,77).Как я могу ограничить это?т.е. убедитесь, что rma обрабатывает мои dat$Estimate как значения между 0 и 1.
> rma(dat$Estimate, dat$SE, method = "DL")
Random-Effects Model (k = 5; tau^2 estimator: DL)
tau^2 (estimated amount of total heterogeneity): 0 (SE = 1.2621)
tau (square root of estimated tau^2 value): 0
I^2 (total heterogeneity / total variability): 0.00%
H^2 (total variability / sampling variability): 1.00
Test for Heterogeneity:
Q(df = 4) = 0.1380, p-val = 0.9977
Model Results:
estimate se zval pval ci.lb ci.ub
0.6302 0.5822 1.0824 0.2791 -0.5109 1.7712
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1