Prezados, boa noite.
Uma dúvida sobre como reproduzir alguns resultados de um artigo no R.
O artigo tem o seguinte DOI: 10.1590/1980-549720190009 (Revista Brasileira de Epidemiologia).
Tem-se uma amostra de 28427 indivíduos. Destes 17244 tem o habito de ler a bula de remédios (A) e 11183 (o restante) não tem este hábito (B).
A dúvida é como obter a prevalência e o IC95% que eles conseguiram no artigo, a saber:
(A) 59.6 (57.7 - 61.5)
(B) 40.4 (38.5 - 42.5)
Eu obtive os seguintes resultados:
(A)
17244/28427 = 60.7
17244/28427-1.96*sqrt(17244/28427*(1-17244/28427)/28427) = 60.1
17244/28427+1.96*sqrt(17244/28427*(1-17244/28427)/28427) = 61.2
60.7 (60.1 - 61.2)
(B)
11183/28427 = 39.3
11183/28427-1.96*sqrt(11183/28427*(1-11183/28427)/28427) = 38.8
11183/28427+1.96*sqrt(11183/28427*(1-11183/28427)/28427) = 39.9
39.3 (38.8, 39.9)
Primeiro... o porquê das diferenças. Segundo... sugestão de material para leitura. E terceiro... como fazer isso no R? Há algum pacote pronto?
Qualquer ajuda será bem vinda.
Obrigado.