Aplicação de regressão linear para correção de dados dietéticos

OBJETIVO: Estratégias metodológicas vêm sendo desenvolvidas para minimizar o efeito do erro de medida da dieta. O objetivo do estudo foi descrever a aplicação de uma estratégia para correção da informação dietética pelo erro de medida. MÉTODOS: Foram obtidos dados de consumo alimentar pela aplicação do Questionário de Freqüência Alimentar a 79 adolescentes do Município de São Paulo em 1999. Os dados dietéticos obtidos foram corrigidos por meio de regressão linear, após o ajuste pela energia usando-se o método dos resíduos. O método de referência utilizado foi o recordatório de 24 horas, aplicado em três momentos distintos. RESULTADOS: Os valores corrigidos aproximaram-se dos valores de referência. O fator de correção lambda foi de 0,89 para energia. Para os macronutrientes, os fatores foram de 0,41; 0,22 e 0,20, para carboidratos, lipídios e proteínas, respectivamente. CONCLUSÕES: As médias e desvios-padrão dos valores corrigidos denotam que houve uma correção do erro de medida. Apesar disso, debate-se o desempenho desses métodos, que são notoriamente imperfeitos quando seus pressupostos teóricos não são atendidos, o que é comum nos estudos da dieta que usam instrumentos de medida baseados no relato dos indivíduos.

Uso de análise exploratória de dados e de regressão robusta na avaliação do crescimento de espécies comerciais de terra firme da Amazônia

Esta pesquisa originou-se de uma base de dados de dois inventários a 100%, realizados em 1984 e 2000 sobre uma mesma área de 576 ha de floresta tropical primária, localizada na Floresta Nacional do Tapajós, Belterra, Pará. O objetivo da pesquisa foi utilizar a análise exploratória de dados e a regressão robusta para modelar o crescimento em diâmetro e área basal. Nos dois inventários as circunferências à altura do peito (CAP) foram medidas com fitamétrica, enquanto as alturas comerciais em 1984 foram medidas com o hipsômetro de Weiss e em 2000, estimadas com o uso de varas; o DAP mínimo em 1984 foi de 55 cm para todas as espécies e em 2000 foi adotado o DAP mínimo de 35 cm. As análises estatísticas foram precedidas de análise exploratória de dados (AED), em que foram utilizados o box plot (caixa-de-bigodes) na detecção de outliers (observações discrepantes) e o gráfico stem-and-leaf (tronco-e-folhas) para filtrar as observações extremas. Utilizou-se a regressão robusta para ajustar os modelos na presença de outliers. A pesquisa mostrou que, apesar das variações intra e interespecíficas, as seis espécies, em conjunto, apresentaram taxas de crescimento medianas observadas e esperadas em diâmetro e área basal de 4,56 e 4,28 mm/ano e 13,00 e 13,09 cm²/ha/ano, respectivamente. A pesquisa também mostrou que o uso da análise exploratória de dados e da regressão robusta viabilizou a análise e a determinação dos incrementos periódicos em diâmetro e área basal em bases consistentes. A metodologia empregada no caso específico dos dados disponíveis e tipos de variáveis mostrou ser mais eficiente com o uso de regressão robusta pelo método Least Trimmed Square do que pelo método dos mínimos quadrados ordinários.

Modelos de regressão lineares para estimativa de produtividade da soja no oeste do Paraná, utilizando dados espectrais

O trabalho teve o objetivo de avaliar modelos lineares de regressão entre resposta espectral e produtividade em soja, na escala regional. Para isso, foram monitorados 36 municípios do oeste do Paraná, utilizando cinco imagens do satélite Landsat 5/TM da safra de 2004/2005. Foram realizados os procedimentos de transformação radiométrica e correção atmosférica nas imagens, determinando valores físicos das refletâncias aparente e de superfície. Posteriormente, foram calculados os índices de vegetação NDVI e GVI, os quais, por meio de regressões lineares simples e múltiplas, compararam-se com as produtividades oficiais dos municípios, obtidas das estatísticas IBGE. Aplicou-se também uma análise de diagnóstico, para detectar pontos influentes e de colinearidade. Os resultados mostraram que a média dos valores de NDVI e GVI de todas as imagens foi mais bem relacionada com a produtividade do que para cada data separadamente. O uso de regressões múltiplas com os dois índices, em todas as datas, propiciou melhores resultados de relação com a produtividade.

Elaboração e validação de intervalos de referência longitudinais de peso fetal com uma amostra da população brasileira

OBJETIVOS: Elaborar modelos de predição de peso fetal e de percentis longitudinais de peso fetal estimado (PFE) com uma amostra da população brasileira. MÉTODOS: Estudo observacional prospectivo. Dois grupos de gestantes foram recrutados: Grupo EPF (estimativa de peso fetal): pacientes para elaboração (EPF-El) e validação (EPF-Val) de um modelo de predição de peso fetal; Grupo IRL (intervalos de referência longitudinais): gestantes para elaboração (IRL-El) e validação (IRL-Val) de intervalos de referência longitudinais de PFE. Regressão polinomial foi utilizada com os dados do subgrupo EPF-El para gerar o modelo de predição de peso fetal. O desempenho deste modelo foi comparado com os de outros disponíveis na literatura. Modelos lineares mistos foram usados para elaboração de intervalos longitudinais de PFE com os dados do subgrupo IRL-El. Os dados do subgrupo IRL-Val foram usados para validação destes intervalos. RESULTADOS: Quatrocentos e cinqüenta e oito pacientes compuseram o Grupo EPF (EPF-El: 367; EPF-Val: 91) e 315 o Grupo IRL (IRL-El: 265; IRL-Val: 50). A fórmula para cálculo do PFE foi: PFE=-8,277+2,146xDBPxCAxCF-2,449xCFxDBP². Os desempenhos de outras fórmulas para estimativa de peso fetal em nossa amostra foram significativamente piores do que os do modelo gerado neste estudo. Equações para predição de percentis condicionais de PFE foram derivadas das avaliações longitudinais do subgrupo IRL-El e validadas com os dados do subgrupo IRL-Val. CONCLUSÕES: descrevemos um método para adaptação de intervalos de referência longitudinais de PFE, sendo este obtido por meio de fórmulas geradas em uma amostra da população brasileira.

Razões de proporções e uso do método delta para intervalos de confiança em regressão logística

OBJETIVO: Apresentar uma aplicação da regressão logística para estimar razões de prevalência ou risco relativo, e o uso do método delta para a estimativa de intervalos de confiança. MÉTODOS: Utilizou-se o método delta, visto que esta técnica é adequada para estimar variâncias de funções não lineares de variáveis aleatórias, do tipo razão de proporções. O método baseia-se na expansão da função em série de Taylor, provendo uma aproximação para o cálculo das variâncias. Para estimar as razões de proporções e as respectivas variâncias, que permitem a estimação dos intervalos de confiança, desenvolveu-se programa utilizando-se o módulo matricial do SAS. A aplicação prática do modelo proposto é feita utilizando-se dados de um estudo transversal com uma amostra de 611 mulheres, para testar a hipótese de que a falta de compartilhamento do trabalho doméstico associa-se com altos escores de sintomas psicológicos medidos através de um questionário validado. As razões de prevalência calculadas através da regressão logística, e os intervalos de confiança estimados pelo método delta foram comparadas com os resultados obtidos pela técnica de Mantel-Haenszel. RESULTADOS: As estimativas pontuais das razões de prevalência brutas, cruas ou ajustadas, obtidas através da regressão logística, foram semelhantes às encontradas na análise tabular. Também os limites dos intervalos de confianças para as razões de prevalência através do método delta foram praticamente iguais em relação aos encontrados através da técnica de Mantel-Haenszel. CONCLUSÕES: Os resultados apóiam o uso do método delta para estimar intervalos de confiança para razões de proporção. Esse método deve ser visto como uma alternativa para situações em que o efeito em estudo (variável dependente) não é raro, situação onde a "odds ratio" se afasta das medidas epidemiológicas caracterizadas como razão de proporções.

Metodologia caso-controle em epidemiologia de doenças cardiovasculares: II - Análise de dados

A larga difusão, nos últimos anos, de softwares estatísticos indicados para análises multivariadas, aliada à facilidade de seu manuseio, poderá conduzir os usuários à obtenção de estimativas enganosas caso não se conheçam os pressupostos teóricos necessários para o emprego destas técnicas. Com o propósito de se destacar alguns desses pressupostos, apresentam-se as etapas executadas para a análise dos dados de um estudo caso-controle, conduzido no Município de São Paulo, entre março de 1993 e fevereiro de 1994, cujo objetivo foi testar a associação do diabetes mellitus com a doença isquêmica do coração após o ajustamento para possíveis variáveis de confusão e/ou modificadoras de efeito. Para uma melhor visualização didática desses pressupostos metodológicos são destacadas quatro etapas fundamentais durante a fase de análise dos dados: a) construção do banco de dados; b) cálculo do poder estatístico; c) categorização e codificação das variáveis e d) a escolha do modelo de regressão logística multivariada.

Comparação de dados sobre homicídios entre dois sistemas de informação, Minas Gerais

INTRODUÇÃO: Existem muitas críticas à qualidade dos dados criminais existentes no Brasil, mas há pouco estudo sistemático para corroborar essas críticas. Nesse sentido, foi feito estudo com o objetivo de comparar o número de homicídios registrados entre dois sistemas públicos de informação. MÉTODOS: Foram analisados os óbitos registrados pelo Sistema de Informações de Mortalidade do Ministério da Saúde (SIM/MS) e pela Polícia Militar (PMMG), de 1979 a 1998. A comparação entre os dois sistemas foi feita pela razão entre os números absolutos, categorizando os municípios em cinco grupos de acordo com o tamanho da população em 1998, e pelas taxas de homicídio. Foram ajustados dois modelos de regressão linear das taxas no tempo, um para cada sistema. RESULTADOS: Nas cidades de menos de 100.00 habitantes, as informações de homicídios da PMMG são mais completas do que as do SIM/MS. Nas cidades de mais de 100.000 habitantes, o SIM/MS é capaz de recuperar mais eficientemente os óbitos, embora possam ser incluídos indevidamente óbitos de homicídios cometidos em outros municípios e óbitos por arma branca e de fogo de intenção indeterminada que não tenham sido devido a homicídio. O sub-registro da PMMG pode ser devido às mortes hospitalares tardias, que não são acompanhadas pela PMMG. CONCLUSÕES: Os sistemas do SIM/MS e PMMG representam limites superiores e inferiores do número real de homicídios ocorridos nas grandes cidades. Sugere-se que, quando disponíveis, as duas fontes sejam analisadas e comparadas; quando apenas uma delas estiver disponível, propõe-se que a análise deva considerar as diferenças.

Qualidade dos dados da vigilância epidemiológica da dengue em Belo Horizonte, MG

OBJETIVO: Avaliar a qualidade dos dados do sistema de vigilância epidemiológica na detecção de casos suspeitos de dengue internados em hospitais públicos e conveniados do Sistema Único de Saúde. MÉTODOS: O estudo foi realizado em Belo Horizonte, Estado de Minas Gerais, no período de 1996 a junho de 2002. Foram adotados os critérios de avaliação de qualidade dos dados de vigilância da dengue do Guidelines for Evaluating Public Health Surveillance Systems. Como padrão de referência foram utilizados os prontuários médicos revisados e validados dos pacientes internados e registrados na rede hospitalar do sistema público de saúde. Foram obtidos 266 (90%) do total de 294 prontuários selecionados, 230 (86,5%) preencheram o critério de caso suspeito de dengue. Para verificar associação entre o sub-registro e variáveis selecionadas, utilizou-se o odds ratio com intervalo de confiança de 95% em modelo de regressão logística. Para avaliar a sensibilidade do sistema de notificação, utilizou-se a proporção de casos internados no sistema hospitalar que estavam notificados; para o valor preditivo positivo, utilizou-se a proporção de casos confirmados por laboratório e registrados no sistema de notificação. RESULTADOS: Verificou-se sub-registro de 37% dos casos no período de 1997 a 2002, com cinco vezes mais chances de ocorrência nos três primeiros anos (OR=5,93; IC 95%: 2,50-14,04), oito vezes mais nas internações em hospitais conveniados que naqueles públicos (OR=8,42, IC 95%: 2,26-31,27). O sub-registro associou-se também aos casos clínicos internados sem manifestações hemorrágicas (OR=2,81; IC 95%: 1,28-6,15), e sem exames laboratoriais específicos para dengue no prontuário (OR=4,07; IC 95%: 1,00-16,52). A sensibilidade estimada do sistema de notificação foi de 63% e o valor preditivo positivo foi de 43%. CONCLUSÕES: Os casos de dengue registrados no sistema de notificação foram aqueles de evolução mais grave e não representaram a totalidade de casos internados no sistema público de saúde, superestimando a taxa de letalidade da doença. Os resultados indicam a necessidade de mudanças no modelo da vigilância e de implementação da capacitação dos profissionais de saúde, principalmente aqueles que trabalham em hospitais conveniados.

Lógica fuzzy e regressão logística na decisão para prática de cintilografia das paratiróides

OBJETIVO: Desenvolver e comparar dois modelos matemáticos, um deles baseado em regressão logística e o outro em teoria de conjuntos fuzzy, para definir a indicação para a realização do exame cintilográfico a partir de resultados dos exames laboratoriais. MÉTODOS: Foram identificados 194 pacientes que tiveram cálcio e paratormônio séricos medidos a partir da base de registros de cintilografia de paratiróides realizadas em laboratório de diagnóstico de São Paulo, no período de janeiro de 2000 a dezembro de 2004. O modelo de regressão logística foi desenvolvido utilizando-se o software SPSS e o modelo fuzzy, o Matlab. A performance dos modelos foi comparada utilizando-se curvas ROC. RESULTADOS: Os modelos apresentaram diferenças estatisticamente significantes (p=0,026) nos seus desempenhos. A área sob a curva ROC do modelo de regressão logística foi de 0,862 (IC 95%: 0,811-0,913) e do modelo de lógica fuzzy foi 0,887 (IC 95%: 0,840-0,933). Este último destacou-se como particularmente útil porque, ao contrário do modelo logístico, mostrou capacidade de utilizar informações de paratormônio em intervalo em que os valores de cálcio mostraram-se pouco discriminantes. CONCLUSÕES: O modelo matemático baseado em teoria de conjuntos fuzzy pareceu ser mais adequado do que o baseado em regressão logística como método para decisão da realização de cintilografia das paratiróides. Todavia, sendo resultado de um exercício metodológico, inferências sobre o comportamento do objeto podem ser impróprias, dada a não representatividade populacional dos dados.

Assistência à saúde e mortalidade de recém-nascidos de muito baixo peso

OBJETIVO: Analisar a mortalidade intra-hospitalar dos recém-nascidos de muito baixo peso, considerando a evolução clínica e os fatores associados à mortalidade. MÉTODOS: Estudo longitudinal que incluiu 360 recém-nascidos com peso entre 500 e 1.500g, em Londrina, Paraná, de 1/1/2002 a 30/6/2004. Os dados foram coletados por meio de entrevistas com as mães, análise dos prontuários e acompanhamento dos recém-nascidos. Para determinação de associação entre as variáveis utilizou-se o teste do qui-quadrado e análise de regressão logística com modelo hierarquizado, com nível de significância de 5%. RESULTADOS: A taxa de mortalidade foi de 32,5%. Na análise bivariada, as variáveis associadas ao óbito oram: não uso de corticosteróide antenatal, ausência de hipertensão arterial/pré-eclampsia, presença de trabalho de parto, parto normal, apresentação não cefálica, Apgar < 3 no primeiro e quinto minutos, Clinical Risk Index for Babies > 5, reanimação na sala de parto, sexo masculino, idade gestacional < 28 semanas, peso < 750g, síndrome do desconforto respiratório, pneumotórax, hemorragia intracraniana e ventilação mecânica. Após regressão logística, permaneceram como fatores de risco: baixa renda per capita, não uso de corticosteróide antenatal e não uso de pressão positiva contínua de vias aéreas. CONCLUSÕES: Mesmo com o uso de tecnologias, a mortalidade observada nos recém-nascidos de muito baixo peso foi alta quando comparada com os países desenvolvidos. A maior utilização do corticosteróide antenatal poderá diminuir a morbidade e mortalidade de recém-nascidos de muito baixo peso.

Regressão logística ordinal em estudos epidemiológicos

Os modelos de regressão logística ordinal vêm sendo aplicados com sucesso na análise de estudos epidemiológicos. Entretanto, a verificação da adequação de cada modelo tem recebido atenção limitada. O artigo apresenta uma breve análise dos principais modelos de regressão logística ordinal e as estratégias para ajuste s, as técnicas de verificação de qualidade do ajuste, bem como os comandos para execução nos softwares R e Stata. A metodologia é ilustrada com aplicação dos dados do Second National Health and Nutrition Examination Survey (NHANES II), o conhecido levantamento de saúde e nutrição.

Modelo preditivo de óbito a partir de dados do Sistema de Informações Hospitalares

OBJETIVO: Desenvolver um modelo preditivo de óbito hospitalar com base nos dados do Sistema de Informações Hospitalares do Sistema Único de Saúde. MÉTODOS: Estudo transversal com dados de 453.515 autorizações de internação de 332 hospitais do Rio Grande do Sul no ano de 2005. A partir da razão entre óbitos observados e óbitos esperados elaborou-se um ranking ajustado dos hospitais que foi comparado ao ranking bruto da taxa de mortalidade. Utilizou-se regressão logística para desenvolvimento do modelo preditivo de probabilidade para óbito hospitalar segundo sexo, idade, diagnóstico e uso de unidade de terapia intensiva. Foram obtidos os intervalos com 95% de confiança para 206 hospitais com mais de 365 internações por ano. RESULTADOS: Obteve-se um índice de risco para mortalidade hospitalar. A ordenação dos hospitais utilizando apenas a taxa de mortalidade bruta diferiu da ordenação quando se utiliza o ranking ajustado pelo modelo preditivo de probabilidade. Dos 206 hospitais analisados, 40 hospitais apresentaram mortalidade observada significativamente superior à esperada e 58 hospitais com mortalidade significativamente inferior à esperada. Uso de unidade de terapia intensiva apresentou maior peso para a composição do índice de risco, seguida pela idade e diagnóstico. Quando os hospitais atendem pacientes com perfis muito diferentes, o ajuste de risco não resulta numa indicação definitiva sobre qual prestador é o melhor. Os hospitais de grande porte apresentaram, no conjunto, maior número de óbitos do que seria esperado de acordo com as características das internações. CONCLUSÕES: O índice de risco de óbito hospitalar mostrou-se preditor adequado para o cálculo dos óbitos esperados, podendo ser aplicado na avaliação do desempenho hospitalar. Recomenda-se que, ao comparar hospitais, seja utilizado o ajuste pelo modelo preditivo de probabilidade de risco, estratificando-se pelo porte do hospital.

Mudanças em dez anos das desigualdades sociais em saúde dos idosos brasileiros (1998-2008)

Analisar mudanças nos diferenciais por renda nas condições de saúde e no uso de serviços de saúde por idosos brasileiros. MÉTODOS: Foram analisadas amostras representativas da população brasileira com 60 anos ou mais em 1998 e 2008 (n = 27.872 e 41.198, respectivamente), oriundas da Pesquisa Nacional por Amostra de Domicílios. As variáveis consideradas foram renda mensal domiciliar per capita, autoavaliação da saúde, capacidade funcional, consultas médicas e hospitalizações nos 12 meses precedentes e uso exclusivo do Sistema Único de Saúde. A análise dos dados foi baseada em estimativas de prevalência e em razões de prevalência obtidas por meio da regressão de Poisson. RESULTADOS: Em 1998 e 2008, as prevalências ajustadas por idade e sexo da autoavaliação da saúde como ruim, do comprometimento da mobilidade e da incapacidade para realizar atividades da vida diária apresentaram fortes gradientes com o quintil da renda domiciliar per capita, com pior performance entre aqueles com renda mais baixa. As razões de prevalência ajustadas por idade e sexo entre o quintil inferior (mais pobres) e o superior (mais ricos) de renda permaneceram estáveis para pior autoavaliação da saúde (RP = 3,12 [IC95% 2,79;3,51] em 1998 e 2,98 [IC95% 2,69;3,29] em 2008), comprometimento da mobilidade (RP = 1,54 [IC95% 1,44;1,65 e 1,69[IC95% 1,60;1,78], respectivamente) e incapacidade para realizar atividades da vida diária (RP = 1,79 [IC95% 1,52;2,11] e 2,02 [IC95% 1,78;2,29], respectivamente). Observou-se redução das disparidades por renda na realização de três ou mais consultas médicas e no uso exclusivo do Sistema Único de Saúde. Não foram observadas desigualdades entre os extremos de renda na ocorrência de hospitalizações no mesmo período. CONCLUSÕES: Apesar da redução das desigualdades por renda de indicadores do uso de serviços de saúde, a magnitude das disparidades nas condições de saúde não diminuiu. São necessários estudos longitudinais para um melhor entendimento da persistência dessas desigualdades entre idosos brasileiros.

Influencia de fatores psicossociais na cessacao do tabagismo: evidencias longitudinais no Estudo Pro-Saude

OBJETIVO Analisar a incidência de cessação de tabagismo e sua associação com fatores psicossociais. MÉTODOS Foram analisados dados de três fases consecutivas do Estudo Pró-Saúde, estudo longitudinal entre funcionários técnico-administrativos de uma universidade no Rio de Janeiro, Brasil. Os critérios de inclusão foram ter participado das Fases 1 e 3 e ser fumante na linha de base (Fase 1 – 1999). Foram excluídos indivíduos que, no seguimento (Fase 3 – 2007), tinham parado de fumar há menos de um ano. A população final de estudo foi composta de 661 funcionários (78% dos elegíveis). Os riscos relativos (RR) da cessação do tabagismo foram avaliados pelo modelo de regressão de Poisson com variância robusta. RESULTADOS A incidência acumulada de cessação do fumo em oito anos de seguimento foi de 27,7%. Entre os fatores psicossociais avaliados, ausência de experiência de violência física apresentou-se associada à maior cessação de tabagismo na análise multivariada (RR = 1,67; IC95% 1,09;2,55). CONCLUSÕES A incidência de cessação de tabagismo foi alta, e o fato de não terem sido encontradas associações com a maioria dos fatores avaliados sugere que grande parte do efeito encontrado seja decorrente do impacto relativamente homogêneo das políticas públicas implementadas no Brasil nas últimas décadas. A associação encontrada entre ausência de exposição à violência e maior incidência de cessação de tabagismo aponta para a importância desse fator nas políticas de controle do tabagismo.

Estudos silviculturais de uma população natural de Copaifera multijuga Hayne - Leguminosae, na Amazônia Central. IV. Interpretação de dados fenológicos em relação a elementos climáticos

O autor apresenta informações fenológicas mensais de 82 matrizes de Copaiba (Copaifera multljuga Hayne), obtidas na Reserva Florestal Adolfo Ducke no período de fevereiro de 1979 a dezembro de 1985. São mostradas as percentagens de ocorrência das fases de floração, frutificação e mudança foliar e uma análise de regressão múltipla das fenofases com sete variáveis do clima (temperatura máxima, temperatura minima e temperatura média, em C°; umidade relativa em %; precipitação total em milímetros; insolação total em horas e evaporação líquida total em milímetros. A plena floração ocorreu na estação de chuvas e a frutificação (frutos maduros) apresentou picos máximos no fim da estação chuvosa e inicio da estação seca, comportamento que. difere da maioria das espécies arbóreas tropicais. A espécie portou-se como perenifólia e a fase de folhas novas antecedeu a da floração.