Неправильный расчет отношения риска / относительного риска с мета-пакетом escalc() для мета-анализа
У меня быстрый (и мне жаль, если глупый) вопрос.
Я использую пакет метафора и эту функцию:
data<-escalc(measure = "RR", ai=1, bi=0, ci=43, di=443, append=TRUE)
рассчитать относительный риск.
В то время как мои ручные вычисления говорят о 11.30 (log(RR)=2.425), функция escalc выдает 8.40 (log(RR)=2.13).
Я на самом деле взял эти значения из исследовательской работы и знаю, что 11.30 - это правильный результат, но я понятия не имею, почему этот расчет дает 8.40.
Любая помощь очень приветствуется
1 ответ
Решение
Когда одна из ячеек равна 0, escalc по умолчанию применяет коррекцию непрерывности +1/2. Если вы не хотите этого, используйте add=0, Увидеть help(escalc),