SciELO - Brazil - Razões de proporções e uso do método delta para intervalos de confiança em regressão logística Razões de proporções e uso do método delta para intervalos de confiança em regressão logística
Artigos Originais • Rev. Saúde Pública 31 (1) • Fev 1997 • https://doi.org/10.1590/S0034-89101997000100012 copiar
Ratio of proportions and the use of the delta method for confidence interval estimation in logistic regression
Logistic regression; Confidence intervals; Relative risk; Delta method
Nelson F. de Oliveira, Vilma S. Santana e Antonio Alberto Lopes
Instituto de Matemática da Universidade Federal da Bahia. Salvador, BA - Brasil (N.F.O.); Instituto de Saúde Coletiva da Universidade Federal da Bahia. Salvador, BA - Brasil (V.S.S.); Departamento de Medicina da Universidade Federal da Bahia. Salvador, BA - Brasil (A.A.L.)
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.
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.
Regressão logística. Intervalo de confiança. Risco relativo. Método delta.
Objective To present an application of logistic regression modelling to estimate ratios of proportions, such as prevalence ratio or relative risk, and the Delta Method to estimate confidence intervals.
Method The Delta Method was used because it is appropriate for the estimation of variance of non-linear functions of random variables. The method is based on Taylor's series expansion and provides a good approximation of variance estimates. A computer program, utilizing the matrix module of SAS, was developed to compute the variance estimates. A practical demonstration is presented with data from a cross-sectional study carried out on a sample of 611 women, to test the hypothesis that the lack of housework sharing is associated with high scores of psychological symptoms as measured by a validated questionnaire.
Results Crude and adjusted prevalence ratio estimated by logistic regression were similar to those estimated by tabular analysis. Also, ranges of the confidence intervals of the prevalence ratio according to the Delta Method were nearly equal to those obtained by the Mantel-Haenszel approach.
Conclusions The results give support to the use of the Delta Method for the estimation of confidence intervals for ratios of proportions. The method should be seen as an alternative for situations in which the need to control a large number of potential confounders limits the use of stratified analysis.
Logistic regression. Confidence intervals. Relative risk. Delta method.
O uso de modelagem matemática, especialmente a regressão logística, tem se tornado uma prática quase rotineira na análise epidemiológica. Aparentemente, isso se relaciona com a popularidade de pacotes estatísticos de utilização muitas vezes fácil e rápida, que permitem o manejo de múltiplas variáveis simultaneamente. Todavia, enquanto nos estudos de caso-controle a estimativa da "odds ratio" (OR) é obtida pela simples exponenciação do coeficiente de regressão da variável independente de interesse, nos estudos transversais e coortes com populações fixas, as medidas de associação de maior interesse epidemiológico são, respectivamente, a razão de prevalências (RP) e a razão de incidência cumulativa14 (risco relativo ou RR). Este é provavelmente o motivo principal para a crescente apresentação de "odds ratios" de prevalência (ORP) ou "odds ratio" de incidência (ORI) como medidas de associação em estudos que não os de caso-controle.
As ORs têm tido o seu significado epidemiológico per se questionado, apontando-se para o seu valor interpretativo apenas quando se aproxima das estimativas de razões de proporções ou de funções de densidade de incidências8. Vários estudos têm apontado as limitações das razões de produtos cruzados, as ORs, como estimadores de RPs. Por exemplo, aproximação matemática entre RPs e ORs ocorre apenas em certas situações específicas, como a raridade do evento predito (doença) entre expostos e não expostos, e estabilidade da prevalência da exposição11, 12. Exceção para os estudos de caso-controle, coorte por amostragem baseada na densidade de incidência, nos quais os controles são selecionados do conjunto de não-casos a cada vez que ocorre um caso6.Isto tem levado a que alguns autores4, 11 recomendem o uso de OR apenas nas situações onde esta medida pode ser tomada como uma aproximação não viezada de razões de densidade de incidência, embora outros3, 8 indiquem o uso de OR para estudos transversais. Discussão mais ampla desta questão vem sendo recentemente desenvolvida por Lee e Chia9, Axelsson e col.1, Stromberg15 e Lee10.
Enquanto é possível, com alguns cálculos, se estimar as razões de proporções com base nos coeficientes de regressão logística, mesmo quando o ajustamento para variáveis de confusão se faz necessário, a estimação de intervalos de confiança para estas medidas requer procedimento de difícil operacionalização. Neste trabalho, usa-se o modelo de regressão logística para estimar razões de proporções e desenvolve-se, baseado no método delta2, um procedimento para estimação de variâncias e, conseqüentemente, intervalos de confiança para as razões de proporções. Este procedimento envolve cálculo matricial e pode ser implementado mediante programas computacionais que disponham de recursos para manipulação de matrizes. Deve ser notado que este método não se aplica para os estudos de função de risco (hazard function) ou de densidade de incidência, quando a medida de associação não é uma razão de proporções ou de probabilidades.
O objetivo do presente artigo é apresentar uma aplicação da regressão logística para estimar razões de prevalências ou risco relativo, e o uso do método delta2 para a estimativa de intervalos de confiança.
APLICAÇÃO DA MODELAGEM LOGÍSTICA EM ANÁLISE EPIDEMIOLÓGICA
A regressão logística tem a sua aplicação mais imediata para a estimativa de ORs. Isso se deve ao fato de que com o antilogaritmo do coeficiente de regressão (eb) de uma variável independente obtém-se uma estimativa não viesada da OR para a associação entre essa variável e a dependente. Assim, para a associação entre uma variável dependente D e uma variável independente E, considerando-se todas as variáveis como dicotômicas e definidas como de desenho, "dummy variables", com a seguinte codificação: presença do fator=1 e ausência do fator=0, o modelo para o logito da probabilidade de ter o efeito D (D = 1) dada a exposição E (E = 1), na ausência de interação, com erro aleatório e, pode ser escrito como
Assim, a estimativa da OR não ajustada ou bruta, (OR(D|E)), poderá ser calculada da seguinte forma:
Quando covariáveis de efeito de confusão são identificadas, na ausência de interação, é recomendado o ajustamento da medida da associação principal para essas covariáveis. Para estimar medidas ajustadas da OR para o efeito de confusão de covariáveis, por exemplo, da covariável C, o ajustamento advirá da manutenção dessa covariável no modelo final. Isto é
e a estimativa da "odds ratio" ajustada será dada pela mesma expressão em (2), com sendo estimado através do modelo (3).
Estes são, portanto, os caminhos básicos para estimativa da medida de associação para estudos de caso-controle com resultados provenientes da regressão logística. Estas equações aplicam-se para estudos não pareados, ou seja, se prestam à situação de modelo de regressão logística não condicional, considerando-se o atendimento a todos os pressupostos desse modelo, observando-se que inferência sobre b0 não pode ser feita a menos que se conheça a fração amostral de casos7.
ESTIMATIVA DE RAZÕES DE PROPORÇÕES COM REGRESSÃO LOGÍSTICA
É bem conhecido que através de modelos de regressão logística (1) é possível determinar probabilidade de ocorrência da variável dependente D, de acordo com cada categoria da variável independente7. Em epidemiologia, a variável dependente é comumente representada por um evento (por exemplo, doença) e a variável independente por um fator de exposição. Em um estudo de corte transversal, as proporções da variável dependente entre expostos e não expostos permitem determinar a razão de prevalência (RP). Semelhantemente, o risco relativo (RR) poderá ser determinado num estudo de coorte fixa, situação onde a proporção de indivíduos que mudam o estado de não-doença para o de doença representa o risco ou incidência cumulativa.
Na Tabela 1, do tipo 2 x 2, apresentam-se os dados para uma situação onde as variáveis relativas à associação principal são dicotômicas. As colunas representam os dois níveis da variável dependente D(1 = doente, 0 = não doente) e as linhas, os dois níveis da variável independente principal, ou de exposição E(1 = exposto, 0 = não exposto). Tomando-se os coeficientes e fornecidos pelo modelo de regressão logística, determina-se a probabilidade (a prevalência ou o risco cumulativo, de acordo com a situação particular) de doença entre os expostos como
Entre os não expostos, a probabilidade de doença é dada por
A razão de probabilidades (a razão de prevalência ou o risco relativo) é o resultado da divisão da probabilidade do evento entre expostos pela probabilidade do evento entre não expostos, ou seja
ESTIMAÇÃO DA VARIÂNCIA DA RAZÃO DE PREVALÊNCIA E INTERVALO DE CONFIANÇA
Nesta seção, para simplicidade da notação, as variáveis independentes, sejam de confusão para a associação principal ou a de exposição, serão denotadas por X1, X2, etc. A variável X1(1 = exposto, 0 = não exposto) é dicotômica, enquanto que as variáveis de confusão poderão ser de qualquer natureza. Será considerado um modelo geral com r variáveis independentes. Neste contexto, o estimador de Pr(D = 1| X1, X2, ..., Xr), é dado por
O estimador da razão de prevalência denotado por , é dado por
ou simplificadamente . Assume-se que ln( ) tem distribuição aproximadamente normal, com média E[ln( )] e variância Var[ln( )]. Para estimar Var[ln( )] será empregado o método delta que se baseia no desenvolvimento, em série de Taylor, da função ln( ). De modo geral, esse método provê uma boa aproximação para médias, variâncias e covariâncias de funções não lineares de uma ou mais variáveis aleatórias. Essencialmente, o método consiste na expansão em série de Taylor da função de interesse, sendo usual tomar-se a aproximação de primeira ordem para a média e de segunda ordem para a variância. No caso de mais de uma variável, o método envolve o cálculo de derivadas parciais da função. Usando-se a notação g(X1, X2, ..., Xr) para função genérica de interesse, E(Xi) = µi para as médias, Var(Xi) para as variâncias das r variáveis e Cov (Xi, Xj) para a covariância entre duas variáveis, obtém-se as seguintes aproximações.
Do mesmo modo obtém-se a seguinte aproximação para a covariância entre duas funções não lineares g1e g2 de variáveis aleatórias
onde o símbolo representa a derivada parcial da função g com relação a variável Xi, em cujo resultado substituiu-se Xi po µi. Estes resultados são aplicados no presente trabalho (veja Apêndice 1), tanto no caso multivariado como no caso univariado. Após a estimação de Var[ln( )], pode-se determinar um intervalo de (1 - a) X 100% de confiança para a razão de prevalência, através de
onde za /2 é o ponto da distribuição normal padronizada acima do qual há uma probabilidade de a /2.
Para a estimação da razão de incidências cumulativas, ou risco relativo, e as suas respectivas variâncias, que permitirão a estimação dos intervalos de confiança, desenvolveu-se um programa utilizando o módulo matricial "Interactive Matrix Language" do "SAS"13 (IML/SAS), apresentado no Apêndice 2. Para sua aplicação, em uma situação concreta de análise, indica-se a realização das seguintes etapas:
1) definir o modelo final através do processo de modelagem logística comum;
2) estimar os parâmetros do modelo da regressão logística (coeficientes de regressão) e a matriz de variância-covariância do vetor de parâmetros;
3) estimar a razão de prevalências para a associação principal baseada no modelo final definido;
4) especificar os vetores de dados;
5) estimar Var[ln( )], de acordo com a expressão matricial (9) desenvolvida no Apêndice 1;
6) estimar o intervalo de confiança com as especificações desejadas, através de (8).
EXEMPLO EM UMA SITUAÇÃO DE ANÁLISE EPIDEMIOLÓGICA
Com dados de um estudo transversal, realizado em uma amostra aleatória de 470 domicílios de um bairro de Salvador, analisou-se a hipótese de que a falta de compartilhamento do trabalho doméstico era um fator de risco para número de sintomas psicológicos em mulheres. A população de estudo era composta de 611 mulheres com idades entre 14 e 70 anos, para as quais dispunha-se de informações sociodemográficas, número de mulheres que desenvolviam atividades domésticas na mesma família e respostas a um questionário sobre sintomas psicológicos, o Questionário de Morbidade Psiquiátrica de Adultos (QMPA). Para essa análise, foi definida como variável de exposição o Compartilhamento do Trabalho Doméstico (COMP), do tipo dicotômica (1 = não compartilhamento, 0 = compartilhamento) e como variável dependente, o escore do QMPA (SCOR) codificado como: 1 = escore alto (³ 7 respostas positivas), 0 = escore baixo (< 7 respostas positivas). Durante a análise, identificou-se que o estado civil (ECIV) e nível de educação formal (EDC), ambas categorizadas como dicotômicas (1 = casadas ou de menor escolarização e 0 = solteiras ou de maior escolarização) eram de confusão da associação principal. Assim, o modelo final pôde ser escrito da seguinte forma:
logito Pr(SCOR = 1) = b0 + b1(COMP) + b2(ECIV) + b3(EDC)
Na análise tabular, utiliza-se o Método de Mantel-Haenszel para estimação da razão de prevalências sendo que os intervalos de confiança são baseados no método de Greenland e Robins (1985)5, enquanto que com a análise logística emprega-se estimativas baseadas na máxima verossimilhança. Encontrou-se uma prevalência global de escores altos do QMPA de 24,60%, considerando-se, para simplicidade de análise, a amostragem como aleatória simples. Para esta magnitude de ocorrência do efeito, a OR não é um bom estimador da razão de prevalências, medida de escolha para a análise desse tipo de estudo, ensejando uma situação para o uso de parâmetros da regressão logística para a estimativa adequada.
De acordo com os dados da Tabela 2, evidencia-se uma associação positiva entre a falta de compartilhamento do trabalho doméstico e escores altos do QMPA, estatisticamente significante ao nível de 5%. Como se esperava, as estimativas pontuais para as associações evidenciam um padrão regular de superestimação das ORs em comparação às RPs. Além disso, a amplitude dos intervalos de confiança (IC) correspondentes, tanto para as medidas brutas, estrato-específicas, ou ajustadas para cada uma das variáveis de confusão, separadamente ou simultaneamente, são maiores para as ORs quando comparadas às RPs . Assim, pode-se dizer que para este efeito, não raro (prevalência de 24,6%), o uso de OR em vez de RP superestima a magnitude da medida de associação. Quanto à precisão, os resultados apontam para uma maior imprecisão das ORs quando comparadas às RPs. Portanto, conclusões poderão ser distintas, de acordo com a medida de associação utilizada.
Na Tabela 3 podem ser vistos os resultados baseados nos coeficientes da regressão logística. Consistentemente, as medidas para a associação e intervalos de confiança para OR e RP diferem, seguindo o mesmo padrão observado na análise tabular. Quando se comparam os resultados obtidos para as estimativas pontuais das ORs ou RPs - estas últimas calculadas com base nas equações derivadas do modelo logístico (Tabela 1, Eq. 1) - nota-se que os valores são essencialmente iguais àqueles calculados na análise tabular. Exceção para pequenas diferenças entre as medidas ajustadas para cada uma das variáveis de confusão, separadamente ou ajuste múltiplo.
Quanto às estimativas dos intervalos de confiança para as RPs, também mostradas na Tabela 3, evidencia-se que as amplitudes obtidas com o método delta2 são muito semelhantes, como era de se esperar, àquelas obtidas na análise estratificada tabular, onde o método delta também é usado para a estimativa das variâncias das ORs. Notam-se pequenas diferenças nas estimativas das RPs ajustadas para uma ou duas variáveis de confusão, quando comparadas àquelas que foram obtidas com as ORs. Portanto, pode-se sugerir o uso desse método, desde que os valores estimados, pontuais ou intervalares são essencialmente iguais aos que advêm da análise tabular. Com isto, uma mesma medida de associação poderá ser utilizada para comparação de resultados provenientes de modelagem logística e análise tabular, evitando-se desse modo o uso inadequado de OR em estudos transversais, que ocasionam superestimação das estimativas pontuais e amplificação da imprecisão evidente nos intervalos de confiança.
Nas Tabelas 4 e 5 apresentam-se os resultados obtidos com o múltiplo controle das variáveis de confusão para as ORs e RPs com a análise tabular e modelagem logística. Observou-se similaridade das estimativas obtidas através dos dois procedimentos, ambos baseados no método delta para o cálculo das variâncias. Os valores ajustados mostram pequenas diferenças uma vez que as estimativas pontuais são calculadas por diferentes métodos, Mantel- Haenszel e regressão logística.
Os fundamentos teóricos e resultados mostrados no presente trabalho apóiam o uso do método delta para estimar intervalos de confiança para razões de proporções com as matrizes de covariância geradas com a regressão logística, que produz valores aproximadamente iguais àqueles que resultam do método tabular. Razões de proporções podem, por sua vez, ser obtidas através do método convencional indicado nos livros-textos3, 7, que se baseia na estimativa das probabilidades esperadas, conforme é mostrado nas equações apresentadas na Tabela 1.
Atenção deve ser dada aos limites da exploração feita no presente artigo: é analisada uma situação com uma prevalência especificada, na ausência de interação e com apenas duas variáveis de confusão, ambos categoriais e dicotômicos. É necessário, a verificação da aplicação e adequação desse método com diferentes valores de prevalência do efeito, desde que estas são conhecidas como positivamente associadas com o desvio entre OR e razões de proporções. Situações onde interações estejam presentes poderão ser abordadas através de uma análise estrato-específica. É também possível apontar que aproximações dos valores ajustados para confusão apresentem, consistentemente, afastamentos porque os pesos conferidos para a obtenção de medidas sumarizadas ajustadas modifica-se de acordo com o método empregado3. Por fim, recomenda-se o uso cuidadoso de técnicas de modelagem em análise epidemiológica, que deve se limitar a situações onde seja necessário o ajuste simultâneo de variáveis de confusão, na condição de que a análise estratificada seja inviável. A análise estratificada é mais vantajosa pela transparência e clareza dos procedimentos e, freqüentemente, por ser de fácil interpretação, não requerendo a adoção dos muitos pressupostos como na modelagem matemática. Desse modo, os métodos apresentados no presente artigo tornam viável a estimação, com as vantagens inerentes à modelagem, da razão de prevalência e intervalos de confiança em estudos do tipo corte transversal ou coorte.
Ao Dr. Sander Greenland, da University of Califórnia (USA). Ao Dr. Bikas K. Sinha, do Indian Statistical Institute, Índia, ao Dr. Dana Loomis, University of North Carolina at Chapel Hill, USA, pelo estímulo à realização desse trabalho.
Correspondência para/Correspondence to: Vilma S. Santana - Instituto de Saúde Coletiva da Universidade Federal da Bahia. Rua Padre Feijó, 29 - 4º andar - 40110 -170 Salvador, BA - Brasil. E-mail: vilma@sunrnp.ufba.br
Recebido em 27.12.1995. Reapresentado em 9.7.1996. Aprovado em 2.9.1996.
1. AXELSON, O.; FREDERIKSON, M.; EKBERG, K. Use of the prevalence ratio v prevalence odds ratio as a measure of risk in cross-sectional studies. Occup. Environ. Med, 51:143-4, 1994.
2. BISHOP, Y.M.M.; FIENBERG, S.E.; HOLLAND, P.W. Discrete multivariate analysis: theory and practice Cambridge, MA, The MIT Press, 1975.
3. CHECKOWAY, H.; PEARCE, N.; CRAWFORD-BROWN, D. J. Research methods in occupational epidemio New York, Oxford University Press, 1989.
4. GREENLAND, S. Interpretation and choice of effect measures in epidemiologic analysis. Am. J. Epidemiol., 125:761-8, 1987.
5. GREENLAND, S. & ROBINS J.M. Estimation of a common effect parameter from sparse follow-up data. Biometrics, 41:55-68, 1985.
6. GREENLAND, S.& THOMAS, D.C. On the need for rare disease assumption in case control in case-control studies. Am. J. Epidemiol, 116:547-53, 1982.
7. HOSMER JR., D. W. & LEMESHOW, S. Applied logistic regression New York, John Wiley &Sons, 1989.
8. KLEINBAUM, D.G; KUPPER, L.L; MORGENSTERN, H. Epidemiologic research: principles and quantitative methods. New York, Van Nostrand Reinhold, 1982.
9. LEE, J. & CHIA, K.S. Estimation of prevalence rate ratios for cross-sectional data: an example in occupational epidemiology. Br. J. Ind. Med, 50:861-2, 1993.
10. LEE, J. Odds ratio or relative risk for cross-sectional data? Int. J. Epidemiol,23:201-3,1994.
11. MIETTINEN, O. Estimability and estimation in case-referent studies. Am. J. Epidemiol,103:226-35, 1976.
12. RODRIGUES, L. & KIRKWOOD, B. R. Case-control designs in the study of common diseases: updates on the demise of the rare disease assumption and the choice of sampling scheme for controls. Int. J. Epidemiol, 19:205-13, 1990.
13. SAS. SAS/ IML Software, Version 6 Cary,NC, SAS Institute Inc., 1990.
14. SAVITZ, D.A. Measurements, estimates and inferences in reporting epidemiologic study results. Am. J. Epidemiol, 135:223-24, 1992.