Eqüidade e provisão de serviços públicos odontológicos no estado do Paraná

OBJETIVO: Analisar a associação entre indicadores socioeconômicos, de provisão de serviços públicos odontológicos e de alocação de recursos financeiros em saúde, e identificar se o sentido das associações ocorre em favor da eqüidade vertical. MÉTODOS: Foi realizado um estudo transversal de abordagem ecológica, com dados do Ministério da Saúde referentes a 399 municípios do estado do Paraná, no período de 1998 a 2005. A condição socioeconômica foi aferida por meio do Índice de Desenvolvimento Humano dos municípios, além de indicadores de renda, educação e saneamento básico, os quais foram obtidos nas bases de dados do Instituto Brasileiro de Geografia e Estatística. Os dados foram submetidos a testes estatísticos não-paramétricos: coeficiente de correlação de Spearman, Friedman e Mann-Whitney. RESULTADOS: Houve tendência redistributiva dos recursos federais transferidos aos municípios para o custeio da atenção básica, intensificada a partir do lançamento da Estratégia Saúde da Família. Observou-se expansão das ações de saúde bucal no período analisado, bem como tendência pró-eqüidade na oferta e utilização dos serviços odontológicos em atenção básica. CONCLUSÕES: Houve tendência redistributiva, ou pró-eqüidade, na provisão de serviços odontológicos no estado do Paraná, com maior provisão per capita de recursos ou serviços para municípios com piores indicadores socioeconômicos. Esta tendência se mostrou compatível com as diretrizes programáticas recentes do Ministério da Saúde.

Validação e reprodutibilidade da Escala de Evaluación de Insatisfación Corporal para Adolescentes

OBJETIVO:Validar escala de insatisfação corporal para adolescentes. MÉTODOS: Participaram do estudo 386 adolescentes, de ambos os sexos, entre dez e 17 anos de idade, de uma escola particular de ensino fundamental e médio, de São Bernardo do Campo, SP, em 2006. Foram realizadas tradução e adaptação cultural da "Escala de Evaluación de Insatisfación Corporal para Adolescentes" para o português. Foram avaliadas consistência interna por meio do coeficiente alfa de Cronbach, análise fatorial pelo método Varimax e validade discriminante pelas diferenças entre médias de estado nutricional, utilizando-se o teste de Kruskal-Wallis. Na validação concorrente, calculou-se o coeficiente de correlação de Spearman entre a escala e o índice de massa corporal, a razão circunferência quadril e a circunferência da cintura. Para reprodutibilidade, foram utilizados o teste de Wilcoxon, o coeficiente de correlação intra-classe. RESULTADOS: A escala traduzida não apresentou discordâncias significativas com a original. A escala apresentou consistência interna satisfatória para todos os subgrupos estudados (fases inicial e intermediária de adolescência, ambos os sexos) e foi capaz de discriminar os adolescentes segundo o estado nutricional. Na análise concorrente, as três medidas corporais foram correlacionadas, exceto adolescentes do sexo masculino em fase inicial, e sua reprodutibilidade foi confirmada. CONCLUSÕES: A Escala de Avaliação da Insatisfação Corporal para Adolescentes está traduzida e adaptada para o português e apresentou resultados satisfatórios, sendo recomendada para avaliação do aspecto atitudinal da imagem corporal de adolescentes.

Validade e confiabilidade da versão brasileira do Índice de Capacidade para o Trabalho

OBJETIVO: Avaliar a validade e a confiabilidade da versão brasileira de índice de capacidade para o trabalho. MÉTODOS:Estudo transversal com amostra de 475 trabalhadores de empresa do setor elétrico no estado de São Paulo (dez municípios em Campinas e região), realizado em 2005. Foram avaliados os seguintes aspectos da versão brasileira do Índice de Capacidade para o Trabalho: validade de construto, por meio de análise fatorial confirmatória e da capacidade discriminante; validade de critério, correlacionado o escore do índice com medidas de saúde auto-referidas; e confiabilidade, por meio da análise da consistência interna utilizando o coeficiente alfa de Cronbach. RESULTADOS: A análise fatorial indicou três fatores do construto capacidade para o trabalho: questões relativas aos "recursos mentais" (20,6% da variância), à autopercepção da capacidade para o trabalho (18,9% da variância) e à presença de doenças e limitações decorrentes do estado de saúde (18,4% da variância). O índice discriminou os trabalhadores segundo nível de absenteísmo, identificando média estatisticamente significativa (p<0,001) entre aqueles com absenteísmo elevado (37,2 pontos) e baixo (42,3 pontos). A análise de critério mostrou correlação do índice com todas as dimensões do estado de saúde analisadas (p<0,0001). O índice apresentou boa confiabilidade com coeficiente alfa de Cronbach (0,72). CONCLUSÕES: A versão brasileira do Índice de Capacidade para o Trabalho mostrou propriedades psicométricas satisfatórias quanto à validade de construto, de critério e de confiabilidade, representando uma opção adequada para avaliação da capacidade para o trabalho em abordagens individuais e inquéritos populacionais.

Desigualdade social, crescimento urbano e hanseníase em Manaus: abordagem espacial

OBJETIVO: Analisar a epidemiologia de hanseníase segundo a distribuição espacial e condições de vida da população. MÉTODOS: Estudo ecológico baseado na espacialização da hanseníase em Manaus (AM), entre 1998 e 2004. Os 4.104 casos obtidos do Sistema de Informações de Agravos de Notificação foram georreferenciados de acordo com a localização dos endereços em 1.536 setores censitários urbanos, por meio de quatro técnicas: correios (73,7% dos endereços encontrados); Programa de Cadastro de Logradouros (7,3%); Programa Saúde da Família (2,1%) e folhas de coleta do Instituto Brasileiro de Geografia e Estatística (1,5%). Para cálculo do coeficiente de detecção utilizou a população de 2001. Na análise espacial foi aplicado o método bayesiano empírico local para produzir uma estimativa do risco da hanseníase, suavizando o efeito da flutuação das taxas, quando calculadas para pequenas áreas. Para análise da associação entre espacialização e fatores de risco empregou-se a regressão logística, tendo como variáveis explicativas a ocorrência de casos em menores de 15 anos (indicador de gravidade) e o Índice de Carência Social construído a partir das variáveis do Censo 2000. RESULTADOS: O coeficiente de detecção apresentou-se hiperendêmico em 34,0% dos setores e muito alto em 26,7%. A medida de associação (odds ratio) referente às variáveis explicativas foi significativa. A combinação de baixa condição de vida e ocorrência em menores de 15 anos foi adotada para identificar as áreas prioritárias para intervenção. CONCLUSÕES: A análise espacial da hanseníase mostrou que a distribuição da doença é heterogênea, atingindo mais intensamente as regiões habitadas por grupos em situação de maior vulnerabilidade.

Análise espacial da dengue e o contexto socioeconômico no município do Rio de Janeiro, RJ

OBJETIVO: Analisar a epidemia de dengue em relação ao contexto socioeconômico segundo áreas geográficas. MÉTODOS: Foi realizado estudo ecológico no município do Rio de Janeiro (RJ), em áreas delimitadas como bairros, a partir de informações de casos de dengue notificados em residentes no município. Foi calculada a taxa de incidência média de dengue entre as semanas epidemiológicas: 48ª de 2001 a 20ª de 2002. A ocorrência de dengue foi correlacionada com variáveis socioeconômicas utilizando-se o coeficiente de correlação de Pearson. Utilizou-se o Índice de Moran global e local para avaliar a autocorrelação espacial da dengue e das variáveis correlacionadas significativamente com a doença. O modelo de regressão linear múltipla e o modelo espacial condicional auto-regressivo foram usados para analisar a relação entre dengue e contexto socioeconômico. RESULTADOS: Os bairros da zona oeste do município apresentaram elevadas taxas de incidência média de dengue. Apresentaram correlação significativa as variáveis: percentual de domicílios ligados à rede sanitária geral, domicílios com lavadora de roupas e densidade populacional por área urbana. O índice de autocorrelação espacial Moran revelou dependência espacial entre a dengue e variáveis selecionadas. Os modelos utilizados apontaram o percentual de domicílios ligados à rede sanitária geral como única variável associada significativamente à doença. Os resíduos de ambos os modelos revelaram autocorrelação espacial significativa, com índice de Moran positivo (p<0,001) para o de regressão e negativo (p=0,005) para o espacial condicional auto-regressivo. CONCLUSÕES: Problemas relacionados ao saneamento básico contribuem decisivamente para o aumento do risco da doença.

Fatores psicossociais e organizacionais na adesão às precauções-padrão

OBJETIVO: Analisar a influência de fatores psicossociais e organizacionais na adesão às precauções-padrão para prevenir a exposição a material biológico em hospital. MÉTODOS: Estudo transversal realizado com 270 profissionais médicos e de enfermagem de um hospital universitário em São Paulo, SP, em 2002. Após seleção por amostragem casual simples, os participantes responderam um questionário sobre variáveis psicossociais na forma de escala de Likert. Foram avaliadas a validade de constructo (análise fatorial) e a confiabilidade (coeficiente alfa de Cronbach). A associação dos fatores psicossociais com a adesão às precauções-padrão foi obtida por meio de análise de regressão logística múltipla, com eliminação retrógrada de variáveis não significativas. RESULTADOS: As escalas mostraram validade e confiabilidade satisfatórias (coeficiente alfa de Cronbach entre 0,67 e 0,82). Fatores individuais relativos ao trabalho e organizacionais explicaram 38,5% do índice global de adesão às precauções-padrão. Esse índice global apresentou associação significativa entre adesão e pertencer ao grupo profissional de médicos, ter recebido treinamento em precauções-padrão no hospital, perceber menos intensamente os obstáculos para seguir as precauções-padrão, perceber mais intensamente a carga de trabalho, o feedback das práticas de segurança e as ações gerenciais de apoio à segurança. CONCLUSÕES: Fatores individuais, relativos ao trabalho e organizacionais influenciam conjuntamente a adesão às precauções-padrão. Programas de prevenção da exposição ocupacional a material biológico devem considerar os obstáculos para seguir as precauções-padrão na prática clínica e enfatizar políticas organizacionais de apoio à segurança no trabalho.

Análise TG-ROC de testes de imunofluorescência no diagnóstico de leishmaniose visceral canina

OBJETIVO: Analisar a acurácia do diagnóstico de dois protocolos de imunofluorescência indireta para leishmaniose visceral canina. MÉTODOS: Cães provenientes de inquérito soroepidemiológico realizado em área endêmica nos municípios de Araçatuba e de Andradina, na região noroeste do estado de São Paulo, em 2003, e área não endêmica da região metropolitana de São Paulo, foram utilizados para avaliar comparativamente dois protocolos da reação de imunofluorescência indireta (RIFI) para leishmaniose: um utilizando antígeno heterólogo Leishmania major (RIFI-BM) e outro utilizando antígeno homólogo Leishmania chagasi (RIFI-CH). Para estimar acurácia utilizou-se a análise two-graph receiver operating characteristic (TG-ROC). A análise TG-ROC comparou as leituras da diluição 1:20 do antígeno homólogo (RIFI-CH), consideradas como teste referência, com as diluições da RIFI-BM (antígeno heterólogo). RESULTADOS: A diluição 1:20 do teste RIFI-CH apresentou o melhor coeficiente de contingência (0,755) e a maior força de associação entre as duas variáveis estudadas (qui-quadrado=124,3), sendo considerada a diluição-referência do teste nas comparações com as diferentes diluições do teste RIFI-BM. Os melhores resultados do RIFI-BM foram obtidos na diluição 1:40, com melhor coeficiente de contingência (0,680) e maior força de associação (qui-quadrado=80,8). Com a mudança do ponto de corte sugerido nesta análise para a diluição 1:40 da RIFI-BM, o valor do parâmetro especificidade aumentou de 57,5% para 97,7%, embora a diluição 1:80 tivesse apresentado a melhor estimativa para sensibilidade (80,2%) com o novo ponto de corte. CONCLUSÕES: A análise TG-ROC pode fornecer importantes informações sobre os testes de diagnósticos, além de apresentar sugestões sobre pontos de cortes que podem melhorar as estimativas de sensibilidade e especificidade do teste, e avaliá-los a luz do melhor custo-benefício.

Escala de Depressão Pós-natal de Edimburgo para triagem no sistema público de saúde

OBJETIVO: Avaliar a utilização da Escala de Depressão Pós-natal de Edimburgo como instrumento de triagem no sistema público de saúde. MÉTODOS: A Escala foi administrada entre o 40º e 90º dia do pós-parto, a 245 mulheres que tiveram parto em uma maternidade privada no município de Belo Horizonte (MG), entre 2005 e 2006. As participantes foram submetidas a uma entrevista psiquiátrica estruturada (Mini-Plus 5.0) utilizada como padrão-ouro para diagnóstico de depressão. Foram calculadas sensibilidade e especificidade da escala e utilizou-se a curva ROC para achar o melhor ponto de corte. Foi utilizado o teste t de Student para comparação das variáveis numéricas e o qui-quadrado para as variáveis categóricas. A confiabilidade foi aferida pelo coeficiente de consistência interna á de Cronbach. RESULTADOS: Foram diagnosticadas 66 mulheres com o quadro depressivo pós-parto (26,9% da amostra). Não houve diferença entre as mulheres com e sem depressão pós-parto em relação à idade, escolaridade, número de partos anteriores e estado civil. Utilizando-se o ponto de corte de 10, a sensibilidade da escala foi 86,4, a especificidade 91,1 e o valor preditivo positivo 0,78. CONCLUSÕES: As propriedades psicométricas da Escala a caracterizam como um bom instrumento de triagem da depressão pós-parto e seu uso disseminado no Sistema Único de Saúde poderia repercutir positivamente com aumento significativo na taxa de reconhecimento, diagnóstico, e tratamento da depressão pós parto.

Estimativas da cobertura de mamografia segundo inquéritos de saúde no Brasil

OBJETIVO: Inquéritos populacionais constituem ferramenta fundamental para monitorar a cobertura de mamografia e os fatores associados à sua realização. Em inquéritos baseados na população residente em domicílios com telefone as estimativas tendem a ser superestimadas. O estudo teve por objetivo estimar a cobertura de mamografia com base em pesquisas de base populacional. MÉTODOS: A partir das coberturas por mamografia em mulheres de 50 a 69 anos, com e sem telefone fixo, observadas na Pesquisa Nacional por Amostra de Domicílios (PNAD) 2003, calcularam-se as razões entre elas e o respectivo coeficiente de variação. A razão de cobertura foi multiplicada pela cobertura estimada pelo sistema de Vigilância de Fatores de Risco e Proteção para Doenças Crônicas por Inquérito Telefônico (VIGITEL), permitindo estimar a cobertura entre mulheres sem telefone em 2007. Essas estimativas foram aplicadas à população de mulheres, com e sem telefone, obtidas a partir da PNAD 2006, obtendo-se assim as estimativas finais para as capitais. RESULTADOS: Em 2007, para o conjunto das capitais, estimou-se a cobertura de mamografia em aproximadamente 70%, variando de 41,2% em Porto Velho (RO) a 82,2% em Florianópolis (SC). Em 17 municípios a cobertura foi maior que 60%; em oito, de 50%-60%; e em dois, a cobertura foi inferior a 50%. Em termos absolutos, a diferença entre as coberturas do VIGITEL e as estimadas foi de 6,5% para o conjunto dos municípios, variando de 3,4% em São Paulo (SP) a 24,2% em João Pessoa (PB). CONCLUSÕES: As diferenças nas magnitudes das estimativas da cobertura de mamografia por inquéritos populacionais são em grande parte reflexo dos desenhos dos estudos. No caso específico da mamografia, seria mais apropriado estimar sua cobertura combinando dados do VIGITEL com aqueles de outros inquéritos, que incluam informações sobre mulheres com e sem telefone fixo, especialmente em municípios de baixa cobertura de telefonia fixa.

Tendências temporais e fatores associados à duração do aleitamento materno em Pernambuco

OBJETIVO: Analisar as tendências temporais de duração do aleitamento materno exclusivo e os fatores de proteção. MÉTODOS: Foram analisadas as prevalências de amamentação total e o aleitamento materno exclusivo de lactentes aos quatro meses de idade. Os dados foram obtidos de inquéritos populacionais realizados em 1991, 1997 e 2006, com 935, 2.081 e 1.568 crianças, respectivamente. Os dados foram coletados por entrevista realizada com o responsável pela criança. As prevalências foram analisadas por regressão de Poisson em relação a: condições ambientais, comportamentais e socioeconômicas, variáveis maternas e fatores biológicos da criança. RESULTADOS: A duração mediana do aleitamento total elevou-se de 89 dias (1991) para 106 dias (1997) e, finalmente, 183 dias (2006). A mediana do aleitamento materno exclusivo manteve-se estacionária em torno de 30 dias entre 1997 e 2006. Na análise multivariada, das dez variáveis analisadas, apenas a escolaridade e idade maternas, situação do domicílio e a criança ser do sexo feminino se mantiveram no modelo explicativo final. CONCLUSÕES: Apesar do aleitamento total ter tido importante aumento de duração, o mesmo não ocorreu com a duração do aleitamento materno exclusivo.

Análise espacial e sazonal da leptospirose no município de São Paulo, SP, 1998 a 2006

OBJETIVO: Analisar a distribuição espacial e sazonal da leptospirose, identificando possíveis componentes ecológicos e sociais para a sua transmissão. MÉTODOS: Foram georreferenciados 2.490 casos em cada distrito do município de São Paulo, SP, registrados de 1998 a 2006. Os dados foram obtidos do Sistema de Informação de Agravos de Notificação. Foram realizados mapas temáticos com as variáveis taxa de incidência, letalidade, taxa de alfabetização, renda média mensal, número de moradores por domicilio, abastecimento de água e rede de esgoto. Para identificar o padrão espacial (disperso, em aglomerado ou randômico), foram analisadas pelo Índice de Moran global e local. Foi utilizado o coeficiente de correlação de Spearman para testar associações entre as variáveis com padrão espacial em aglomerados. RESULTADOS: O padrão espacial em aglomerados foi observado nas variáveis taxa de incidência de leptospirose, taxa de alfabetização, renda média mensal, número de moradores por domicílio, abastecimento de água e rede de esgoto. Foram notificados 773 casos no período seco e 1.717 no úmido. A incidência e a letalidade estão correlacionadas com as condições socioeconômicas da população, independentemente do período. CONCLUSÕES: A leptospirose está distribuída por todo o município de São Paulo e sua incidência aumenta no período das chuvas. No período seco, os locais de aparecimento dos casos coincidem com as áreas de piores condições de moradia e, durante o período úmido, também aumenta em outros distritos, provavelmente devido à proximidade de rios e córregos.

Risco de infecção tuberculosa em agentes comunitários de saúde

OBJETIVO: Estimar o risco de infecção tuberculosa em agentes comunitários de saúde envolvidos no controle da doença. MÉTODOS: Foi seguida uma coorte prospectiva, de abril de 2007 a maio de 2008, no município de Cachoeiro de Itapemirim, ES. A coorte foi composta por 61 agentes comunitários, divididos em não-expostos (n=37) e expostos (que acompanharam pacientes com tuberculose, n=24). Durante os 12 meses de seguimento, foi realizado teste tuberculínico, utilizando a tuberculina PPD RT23. Foi calculado o risco relativo e intervalo com 95% de confiança e foi avaliada a correlação entre a viragem tuberculínica e a história ocupacional dos agentes por meio do coeficiente de correlação de Pearson. RESULTADOS: A incidência da viragem foi de 41,7% no grupo dos expostos e 13,5% no grupo dos não expostos. O risco anual de infecção foi de 52,8% no grupo dos expostos e de 14,4% no grupo dos não expostos (p= 0,013). Observou-se associação entre viragem tuberculínica e exposição a paciente com tuberculose (RR= 3,08; IC 95%: 1,201;7,914). CONCLUSÕES: Os agentes que acompanharam pacientes com tuberculose em suas rotinas de serviço apresentaram risco de infecção maior que aqueles que não acompanharam pacientes com essa doença. A implementação de medidas administrativas de biossegurança de rotina, entre as quais a prova tuberculínica, devem ser priorizadas, considerando o alto risco de infecção tuberculosa entre os agentes comunitários de saúde.

Correção da magnitude da mortalidade por câncer do colo do útero no Brasil, 1996-2005

OBJETIVO: Desenvolver uma metodologia para correção da magnitude dos óbitos por câncer do colo do útero no Brasil. MÉTODOS: Os dados sobre os 9.607.177 óbitos foram obtidos do Sistema de Informação sobre Mortalidade, para o período de 1996 a 2005. Para a correção do sub-registro, foram utilizados os fatores de expansão gerados pelo Projeto Carga Global de Doença no Brasil - 1998. Para correção das categorias de diagnósticos desconhecidos, incompletos ou mal definidos de óbitos, foi aplicada redistribuição proporcional. Os dados ausentes de idade foram corrigidos por imputação. As correções foram aplicadas por Unidade Federativa e os resultados apresentados para o Brasil, região e áreas geográficas (capital, demais municípios das regiões metropolitanas e interior) por meio do percentual de variabilidade da magnitude das taxas, antes e após a correção dos óbitos. O comportamento das correções foi analisado por modelo de regressão linear multivariada com termos de interação entre região do País e área geográfica. RESULTADOS: As taxas corrigidas de mortalidade por câncer do colo do útero no Brasil mostraram um acréscimo de 103,4%, variando de 35% para as capitais da região Sul a 339% para o interior da região Nordeste. A redistribuição dos óbitos por câncer de útero sem especificação de localização anatômica promoveu os maiores acréscimos na magnitude das taxas. Os percentuais de correção, segundo ano de ocorrência do óbito, mostraram tendência estacionária no Brasil. CONCLUSÕES: Os resultados permitem concluir que a metodologia proposta foi adequada para corrigir a magnitude das taxas de mortalidade por câncer do colo do útero no País, mostrando que a mortalidade por esse câncer é ainda maior do que o observado nos informes oficiais.

Aumento na sobrevida de crianças de grupos de peso baixo ao nascer em Santa Catarina

OBJETIVO: Analisar fatores associados à sobrevida no primeiro ano de vida. MÉTODOS: Estudo de coorte histórica foi realizado com dados dos sistemas de informação de nascimento e mortalidade sobre 90.153 registros de nascidos vivos e 1.053 registros de óbitos de menores de um ano em hospitais de Florianópolis e São José, SC, entre 1999 e 2006. Foram estimadas curvas de sobrevida (Kaplan-Meier) para grupos de peso ao nascer, período (quadriênios) e tipo de maternidade. Foram calculadas razões de riscos proporcionais para óbitos utilizando regressão de Cox. RESULTADOS: A sobrevida (98,8%) não se modificou entre os grupos de peso, mas aumentou nos grupos de menos de 2.000 g (77,7% para 81,2%, p = 0,029) entre os quadriênios de 1999 a 2002 e 2003 a 2006. Houve aumento de menores de 2.000 g no segundo quadriênio estudado. O tipo de hospital foi associado significativamente à probabilidade de sobrevida. CONCLUSÕES: Há maior probabilidade de sobrevida entre nascidos em hospitais privados e no hospital de ensino para todos os grupos de peso e para o grupo de menos de 2000 g. A sobrevida dos grupos de peso abaixo de 2000 g aumentou no quadriênio mais recente. Entretanto, o coeficiente de mortalidade infantil não diminuiu nesse período, pois a prevalência dos nascidos em grupos de menor peso também aumentou.

Riscos à saúde em áreas próximas a aterros de resíduos sólidos urbanos

OBJETIVO: Avaliar a associação entre residência próxima a aterros de resíduos sólidos e a ocorrência de câncer e malformações congênitas nessas populações vizinhas. MÉTODOS: Foram selecionados e geocodificados óbitos ocorridos no período de 1998 a 2002 entre residentes do município de São Paulo, SP, segundo causas selecionadas. No período avaliado ocorreram 351 óbitos por câncer de fígado, 160 de bexiga e 224 por leucemia em adultos, 25 óbitos por leucemia entre crianças e 299 por malformação congênita nas áreas próximas aos aterros. Buffers com raios de 2 km em torno de 15 aterros delimitaram as áreas expostas. Razões de mortalidade padronizadas de cada desfecho foram analisadas em modelos espaciais bayesianos. RESULTADOS: De modo geral, os maiores valores das razões de mortalidade padronizadas localizam-se em áreas mais centrais do município e os aterros, nas áreas mais periféricas. As razões de mortalidade padronizadas não indicaram excesso de risco para os residentes nas áreas próximas aos aterros de resíduos sólidos no município de São Paulo. Para aterros em funcionamento encontrou-se risco aumentado para câncer de bexiga, fígado e para mortes por malformações congênitas, porém, sem significância estatística. CONCLUSÕES: Não se encontrou aumento no risco de câncer ou de malformações congênitas nas áreas vizinhas aos depósitos de resíduos urbanos do município de São Paulo. As fracas associações e a imprecisão das estimativas obtidas não permitem estabelecer relação causal.