Epi Info Version 6 Guia Rápida
Epi Info Version 6 Guia Rápida
Epi Info Version 6 Guia Rápida
EPIINFO es la denominaci6n de un conjunto logico definido como "programa procesador de palabras, base de dates y estadfstico para Salud Publica". Tiene un grado de difusi6n elevado por ser de gran utilidad, estar avalado por la Organizaci6n Mundial de la Salud asf como por el Centre de Control de Enfermedades de Atlanta (EUA) y porque se obtiene en las oficinas respectivas sin costos para el interesado. Este tipo de instrumentos de trabajo va asociado a ideas de perfecci6n y de cierta infalibilidad en todo desempeno, lo que puede inducir a su utilizacidn sin prevenci6n alguna. El objetivo de este estudio fue verificar hasta que punto el programa mencionado se ajusta, en la parte de analisis estadistico de resultados, a ciertas normas generalmente aceptadas y, al mismo tiempo, establecer los aspectos en los cuales muchos usuarios pudieran requerir de cierto cuidado.
Material y Metodo
Se dispuso del programa EPIINFO version 6.0 de abril de 1994 y se revisaron, ofreciendo distintos tipos de problemas, parte de los subprogramas estadfsticos STATCALC Y EPITABLE, escogiendo en forma arbitraria dentro de ellos las rutinas probablemente rnas.empleadas por clfnicos, con la intenci6n de hacer al lector algunas advertencias que pudieran ser utiles al analizar con EPIINFO los datos de una investigacidn.
Resultados y Cornell tario STATCALC, en tabla de 2 x 2, indica ventaja o razon de productos cru/ados ("odds ratio, OR") 1 ' 2 y entrega los resultados de la prueba de Ji cuadrado con correcci6n de Yates3 y sin ella. Si hay una frecuencia esperada menor de 5, recomienda efectuar la prueba de Fisher-Irwin de probabilidades exactas4, pero cuando la tabla es 2 x 3 o mas grande, el usuario puede encontrar que en algunas tablas donde hay cifras esperadas menores de 5 el programa advierte sobre este hecho e indica que la prueba de Ji cuadrado no es va"lida, pero en otras tablas no
1. Universidad de Chile, Facultad de Medicina, Departamento de Pediatria, Campus Norte. Con fmanciamiento de Proyecto Fondecyt N 1950656.
341
342
Duffau G,
ocurrira" esto y se entregara" el resultado de la prueba. Asf por ejemplo, en la tabla siguiente:
0 2 9 8
14
1 1 7 9
0
se encuentran cinco celdas con valores esperados inferiores a 5 y el programa no calcula Ji cuadrado, expresando que "hay un valor esperado menor que 5 en una celda". Sin embargo, en esta otra:
0 2 9 8 15 1 1 7 9 0
hay 4 celdas con cifras esperadas inferiores a 5 y pese a ello se indica que Ji cuadrado es 12,94 con 4 grados de libertad y p = 0,011. Esta misma tabla analizada por EPITABLE senala, en cambio, "invalid input" porque ningun valor esperado puede ser inferior a 1. EPITABLE en tabla de 2 x 2 entrega los resultados sin correccio'n por continuidad, excepto cuando la cifra esperada en una celda es inferior a 5 (ver opcion "Help" (Ayuda)), lo que resulta sorprendente porque la correccion de Yates no persigue resolver este problema. Si hay una cifra esperada menor que la unidad, no efectua la prueba y senala "invalid input". Para tablas mayores muestra el porcentaje de valores esperados inferiores a 5 y calcula Ji cuadrado. Ante cifras esperadas inferiores a 1, como se explic6, no realiza la prueba. La correcci6n de Yates se considera importante en tabla de 2 x 2, al aproximar la distribucion de Ji cuadrado a la normal5. Una situaci6n similar se da con la distribuci6n binomial. Es cierto que la correccion hace la prueba mas conservadora y tambie'n lo es que no todos la aplican 6 ' 7 , pero lo claro es que EPIINFO 6.0 la incluye o no segiin el programa que este corriendo. Senalaremos que para la ejecucion de la prueba de Ji cuadrado es generalmente aceptado que no debe haber ningun valor esperado menor que 1 ni mas de 20% de celdas con valor esperado inferior a 5 (entonces, en tabla de 2 x 2 no
puede haber ninguna). La solucidn aqui es reagrupar categorias vecinas1"5. Las ventajas ("OR") son calculadas en STATCALC y se entregan sus h'mites de confianza segiin Cornfield1, aunque se pueden obtener valores exactos. Cuando en la tabla hay un cero, el programa entrega resultado de ventajas o no lo hace (declarando OR "indefmido") dependiendo de la posicion del cero. Asi es como, si este esta ubicado en la celda "a" o "d" muestra valores, cosa que no hace si el cero se encuentra en las celdas "b" o "c". Es mas, si se ponen dos ceros uno en "a" y el otro en "d", tambien da un resultado. Esto es raro si consideramos que para el calculo de ventajas en presencia de un cero se ban propuesto soluciones que no justificarian expresarlo como "indefinido"1'8. Por otra parte, el problema persiste independientemente de la Iocalizaci6n del cero, ya que "OR" resultaria un valor incalculable (a x d/ b x 0 o a x d/0 x c) o cero (a x 0/b x c o 0 x d/b x c) y por ultimo, con dos ceros se dan situaciones que no siempre tienen sentido. Cuando se desea efectuar un analisis estratificado de datos, STATCALC determina la totalidad de las etapas pero no incluye la muy importante prueba de homogeneidad de Mantel y Haenszel1. EPITABLE pide -para el ca~lculo de intervalo de confianza (95%) de un promedio el tamano de la poblaci6n de la que se extrajo la muestra, ademas del tamafio muestral, el promedio y la desviaci6n estandar. Esto es bastante inusual ya que raramente se dispone del tamano de la poblaci6n. De todos modos, al faltar el dato el programa supone una poblaci6n de 99.999.999 efectivos. La idea es que si la muestra tiene un tamano ms cercano al de la poblaci6n, el intervalo de confianza es mas estrecho que si ocurre lo opuesto. Para el analisis de concordancia en la evaluacio"n por dos observadores, el programa calcula el estadistico kappa corregido (para efecto del azar) y el error estdndar para la eventualidad que kappa = O9. Ademas senala el valor de "p" asociado. Estos calculos los efectua de la misma manera en tabla de 2 x 2 y otras mayores. De hecho, no toma en consideraci6n la importancia del grado de discordancia, que tiene s61o una expresi6n si la tabla es 2 x 2, puesto que dos de las cuatro celdas corresponden a acuerdo entre los observadores y las restantes dos a desacuer-
Volumen 66 Numero 6
Epiinfo6.0
343
do. Sin embargo, en una tabla mayor, existen grades de desacuerdo que el programa no contempla. De hecho, en una tabla como la que se muestra a continuacidn :
20 2 0 0
5 0 30 3 3 40 0 5 0 0 2 10
se obtiene un indice kappa practicamente igual que en la siguiente, pese a que la discordancia, siendo igual el numero de casos, es absolutamente diferente:
20 0 0
10
0 0 30 0 0 40 0 0
10 0 0 10
Para rcsolver esto se debe calcular el estadistico kappa ponderado9, que reflejara" mejor la situacion, ya que en el primer caso senalara un valor de 0,83 y en el segundo s61o 0,56. For otra parte, tampoco se entrega el error estandar de kappa encontrado, lo que hace engorroso para el usuario el comparar dos de estos estadisticos, debiendo efectuar el calculo por si mismo. Recordaremos que el error estandar que entrega solo permite decir cudn diferente de cero es el valor entregado de kappa. Al comparar una distribucidn determinada con otra te6rica por medio de Ji cuadrado en la opci6n "bondad del ajuste" ("Goodness of fit"), el numero de grades de libertad no aparece correctamente si hay cifras esperadas inferiores a 5 y -si se coloca un cero- simplemente el programa se devuelve al comienzo de EPIINFO 6.0. Si se desea comparar promedios no aparece otra posibilidad que el analisis de varianza que -"Help" explica- es, para dos muestras, una prueba equivalente a "t" de Student 10 y que para poder efectuarla en propiedad se requiere que las varianzas no sean estadfsticamente diferentes. Si lo son, empleara el procedimiento de Satterthwaite 11 que corrige los grados de libertad, haciendo mas exigente el mdtodo, esto es, requiriendo un valor de "t" mayor para declarar diferente los promedios. Sobre esto hay que aclarar que cuando se dan las circunstancias, aparece en la pantalla que las varianzas son di-
ferentes, se indica el "p" asociado y luego viene la comparaci6n de promedios. Queda senalado, ademas, que se emplear el metodo de Satterthwaite. Sin embargo, aunque las varianzas sean diferentes, cosa comprobable con la rutina "comparaci6n de varianzas", si los grupos en estudio tienen 30 o mas casos, no aparece ninguna de las indicaciones senaladas y solo se destaca la comparacidn de promedios. Esto no esti bien, ya que dos grupos pueden ser diferentes por sus varianzas, sus promedios o ambos7, lo que nos llevarfa, entonces, al punto de partida correcto, es decir, a la rutina ofrecida como comparacion de varianzas, prueba "F'10. Si la razon varianza mayor /varianza menor indica que estas son significativamente distintas, tambien lo son los grupos que se esta~n comparando. Si se insiste en contrastar los promedios, se requerira" una correccion de los grados en libertad, que pudiera ser la propuesta o bien aquella sugerida por Welch 12 . EPITABLE puede proporcionar el intervalo de confianza de una proporcidn, pero es necesario tener presente que cuando e"sta es muy pequena, cercana a cero, el tiempo que demora el programa en dar una respuesta se va dilatando progresivamente. En el estudio de eficacia de vacunas por medio de cohortes se calcula el riesgo relativo1 (RR), pero tambin se incluye OR, lo que no corresponde, puesto que OR representa una aproximacitfn al RR cuando la situaci6n estudiada tiene poca frecuencia y se trata de un diseno caso-control, precisamente porque una investigaci6n de cohortes no ha sido posible. La eficacia de una vacuna aparece defmida como (1-RR) x 100 o bien, para estudio casocontrol, (1 -OR) x 100. En cohortes, "Help" presenta : Enfermedad
Vacunados
(SI)
(no)
(si) a c a+c
(no) b d b+d
a+b c+d N
En donde, el riesgo de enfermar en los sujetos vacunados es a/(a+b) y en los no vacunados c/(c+d). El RR sera" la razon a/(a+b) cuyo valor b/(c-i-d) inferior a uno, se utiliza justamente para determinar eficacia de vacunas en la f6rmula
344
Duffau G.
senalada ms arriba. Ahora bien, en "Help" se propone calcular RR de la siguiente manera: a x (b + d) lo que es b x (a + c), claramente un error. Como sea, si se construye un ejemplo, el programa genera un resultado correcto pero la instrucci6n para hacerlo da un valor diferente y equivocado. Conclusions En lo principal, al analizar una tabla de 2 x 3 o mayor con Ji cuadrado por STATCALC se obtienen resultados del estadistico y de "p" asociado a este aunque el valor esperado sea inferior a uno. Por EPITABLE ello no es posible. EPITABLE, en tabla de 2 x 2, calcula Ji cuadrado con correcci6n de Yates solo si hay un valor esperado < 5. STATCALC determina "OR" con ceros, o no lo hace, dependiendo de la ubicacion de 6stos en la tabla. Incluso acepta y calcula con 2 ceros siempre que est6n en la diagonal principal. En analisis estratificado, no incluye la prueba de homogeneidad de Mantel y Haenszel. EPITABLE, en estudio de concordancia, no proporciona el error estndar del indice kappa calculado y, cuando corresponde, no determina kappa ponderado. En comparacion de promedios por andlisis de varianza (2 grupos) -si e"stas difieren- lo advierte, excepto si hay ms de 30 casos en cada grupo. En el estudio de eficacia de vacunas por medio de cohortes, la opci6n "ayuda" propone una soluci6n equivocada. Rusu men En un estudio critico de las rutinas estadisticas del programa EPIINFO versi6n 6.0, 1994 se advierten posibles problemas para el usuario en la prueba de Ji cuadrado, ya que en tablas de
2 x 3 o mayores -con STATCALC- se obtienen resultados estadisticos y de "p" asociado aunque el valor esperado sea menor que uno, lo que no es posible conseguir con EPITABLE. Tambie'n se detectan fallas en la correccidn de Yates, calculo de ventajas ("odds ratios") con ceros, analisis estratificado, comparaci6n de dos promedios, estudio de concordancia y en evaluaci6n de vacunas por medio de cohortes. (Palabras clave: biometria, estadistica, metodos epidemiologicos, programas computadorizados.) Referencias
1. Fleiss J: Satatistical methods for rates and proportions. 2nd. edition, New York; John Wiley & Sons, 1980: 56-82. 2. Matthews D, Farewell V: Estadfstica me"dica. Aplicacidn e interpretacion. Barcelona: Salvar Editores, S.A., 1988: 11-66. 3. Siege! S: Nonparametric statistics for the behavioral sciences. New York: McGraw-Hill Book Co., 1956: 95-110. 4. Saunders BD, Trapp R: Bioestadfstica me'dica. Mexico DF: Manual Moderno, S.A., 1993: 173-180. 5. ingelfinger J, Mosteller F, Thibodeau L, Ware J: Biostatistics in clinical medicine. 2nd edition. New York: Macmillan Publishing Co., 1987: 168-180. 6. Cotton T: Estadfstica en medicina. Barcelona: Sal vat Editores, S.A., 1979: 183-190. 7. Schwartz D: Mdtodos estadisticos para me'dicos y bidlogos, Barcelona: Herder 1985: 87-105. 8. Valenzuela C: Dos soluciones para la estimacidn de odds ratios con ceros. Rev Med Chil 1993: I21-: 1441-1444. 9. Kramer M, Feinstein A: Clinical bio statistics. LIV. The biostatistics of concordance. Clin Pharmacol Ther I9S1;29: 111-123. 10. Latnotte M: Estadfstica biolo"gica. Barcelona: TorayMassonS.A, 1965:92-100. 11. Steel R, Torrie J: Bioestadfstica: Principios y procedimientos. Madrid: McGraw-Hill, 1990: 102-105. 12. Lentner C: Introduction to statistics in: Scientific tables, 8th edition, Switzerland: Ciba-Geigy Ltd, 1982:209.