R - рассчитать доверительный интервал сгруппированных данных

Question

Предположим, у вас есть набор данных с именем data следующим образом:

Gender Pneumonia_Incidence lower_CI upper_CI
Male   38000               30000    44000
Female 34000               32000    38000

Теперь я хотел бы рассчитать общую частоту пневмонии, которую можно легко сделать:

sum(data$Pneumonia_Incidence)

Но как я могу рассчитать нижнюю и верхнюю КИ для этой оценки? Допустимо ли просто суммировать значения lower_CI и upper_CI или мне нужно сделать что-то еще?

Answer 1

Как рассчитать нижний и верхний значения CI для этой оценки?

Вы не можете использовать предоставленную вами информацию. Вам необходимо знать дисперсию каждой предполагаемой распространенности. Это должно быть известно, поскольку кто-то вычислил доверительные интервалы . С их помощью вы можете получить объединенную оценку для общей дисперсии и затем рассчитать общий доверительный интервал.

Если это важно для вас, я настоятельно рекомендую вам обратиться к квалифицированному статистику или, по крайней мере, к авторитетному учебнику. Это не простая математика средней школы. Могут быть и другие проблемы, такие как вес выборки. Я, конечно, не стал бы здесь обращаться за статистической консультацией, и, вероятно, даже в этом другом месте, которое часто упоминается. Что если ваш начальник спросит вас, как вы рассчитали доверительный интервал? Вы бы сказали, что обращались к inte rnet?

Допустимо ли просто суммировать значения lower_CI и upper_CI ...

Нет. Отклонения связаны с размерами выборки. Учти это. Представьте себе две группы: одну с очень большим размером выборки и одну с очень маленькой. Группа с большим размером выборки будет иметь более узкий доверительный интервал, чем группа с небольшим размером выборки. Если вы просто добавите два интервала, вы получите общий интервал, который будет одинаково взвешен обеими группами, что на первый взгляд кажется неверным. Это предвзятая оценка.

... или мне нужно сделать что-то еще?

Проконсультироваться со статистиком. :)