[Mortality attributable to influenza in pre-vaccination and post-vaccination periods in Argentina: an ecological study (2002-2016)Mortalidade atribuível à gripe no período pré-vacinal e pós-vacinal na Argentina: estudo ecológico (2002-2016)].

OBJECTIVES: Compare mortality from severe acute respiratory infection (SARI) attributable to influenza between pre-vaccination (pre-V) and post-vaccination (post-V) periods, to determine the historical evolution and seasonality of time series between 2002 and 2016, and to estimate the risk of death in children between 6 and 23 months of age, using a statistical model.
METHODS: Time-series study using official mortality data from the official statistical database on deaths. ICD-10 codes between J09-18.9 and J22X were considered to represent SARI. Crude rates and age-adjusted rates (AAR) were calculated, and pre-V (2002-2009) and post-V (2010-2016) periods were compared using the chi-squared (χ2) test. The best statistical model was the one that compared deaths from SARI in children during 2002 with other years. The data were analyzed with R programming (p <0.05).
RESULTS: 4.6% of deaths (301,747) were from SARI, with a median age of 82 years. The percentage of deaths under age 2 declined in the post-V period (from 2.34% to 0.99%, p < 0.05). Marked seasonality was observed in winter. The AAR in persons over age 64 rose from 259.8 per 100,000 population (pre-V) to 328.6 (post-V) (p < 0.05). In children, the crude rate dropped significantly. Compared with the year 2002, there was a significantly lower estimated risk of dying from SARI during the three years post-V.
CONCLUSIONS: The reduction in mortality from influenza in Argentina was more pronounced in children, with an estimated 3.5 fewer child deaths from SARI per month.

La gripe es una enfermedad respiratoria aguda, autolimitada, causada por el virus influenza y posee una distribución mundial. El cuadro clínico varía desde un catarro de vías aéreas superiores hasta una neumonía que puede llevar a la muerte; se presenta en epidemias anuales de magnitud y gravedad variadas (1). Un reciente estudio realizado (2), estima que, en el año, se producen en todo el mundo entre 260 000 y 650 000 muertes por causas respiratorias; esto se calcula sin considerar las muertes causadas por la descompensación de enfermedades preexistentes, como las cardiovasculares. Las tasas más elevadas se presentan en la población que tiene más de 64 años de edad (2), los niños menores de dos años, embarazadas y personas con factores de riesgo (3-5).

Durante el invierno, en Argentina, las infecciones respiratorias agudas son uno de los motivos más frecuente de consulta e internación, y se encuentran entre las cinco primeras causas de muerte (6). Según las recomendaciones de la Organización Mundial de la Salud (OMS), Argentina inició la vacunación gratuita a los mayores de 64 años a partir del año 2003. En 2010 se realizó una Campaña Nacional, y desde 2011 fue incorporado al Calendario Nacional de Vacunación (7) para los niños de 6 a 23 meses, embarazadas, madres de niños menores de seis meses, personal de salud y personas entre 2 a 64 años con factores de riesgo. Según los datos publicados por el Programa Nacional de Inmunización del Ministerio de Salud Pública (PNI), durante la Campaña Nacional de Vacunación antigripal, se colocaron más de 8,7 millones de dosis de vacunas y se calcula una cobertura alrededor del 93% de la población estimada (8). De esta manera, la vacunación antigripal pasó a ser obligatoria y gratuita para toda la población, y debe ser reconocida por las obras sociales y servicios de salud prepagos.

Debido a que la decisión de incorporar una nueva vacuna al Calendario Nacional implica un elevado costo en la salud pública del país, nos propusimos evaluar el efecto de la vacunación antigripal en la población. Se plantearon los siguientes objetivos: comparar las tasas de mortalidad por infección respiratoria aguda grave atribuible a gripe (IRAG) entre los años previos y posteriores a la implementación de la vacunación; conocer la evolución histórica y la estacionalidad de las tasas de mortalidad por IRAG en Argentina entre los años 2002 y 2016; y estimar el riesgo de morir en niños entre 6 y 23 meses mediante un modelo estadístico para evaluar modificaciones de la mortalidad por IRAG.

MATERIALES Y MÉTODOS

Se realizó un estudio ecológico de serie de tiempo (ST) entre 2002 y 2016. Para evaluar el primer objetivo, el efecto de la vacunación de la población en la mortalidad por IRAG atribuible a gripe. Se realizaron estudios de series históricas de las tasas y se compararon las tasas entre los períodos antes y después de la vacunación antigripal. Para estimar el riesgo de morir en niños, se trabajó con modelos estadísticos hasta conseguir los mejores indicadores.

El trabajo fue aprobado por el comité de Ética e Investigación del Instituto Nacional de Epidemiología Dr. Juan H. Jara en noviembre de 2014. La base de datos de mortalidad no poseía datos de identificación de las personas, las otras bases empleadas fueron realizada a partir de datos poblacionales y de acceso gratuito.

Los datos fueron anónimos, procedentes de las siguientes bases:

Informe Estadístico de Defunción (IED): base de datos suministrado por la Dirección de Estadísticas e Información de Salud de Argentina (DEIS). Esta base se elabora a partir del IED cargado por cada provincia y luego compilados por la DEIS, según las normativas de la Ley Nacional 26413, el IED y el certificado de defunción (CD) que son llenados por el médico o personal de salud idóneo quienes certifican la muerte. Esto es obligatorio de realizar en todas las personas, ya sea que haya fallecido en una institución de salud, en el domicilio o en la vía pública. Los IED y CD son entregados al Registro Civil de las personas de cada jurisdicción. El personal que trabaja en Estadísticas Vitales de cada provincia, recolecta todos los IED para codificarlos y digitalizarlos en la oficina central de cada provincia. En la base de datos no figura información referida a la identificación del fallecido (número de documento de identidad, nombre y apellido). En el IED, se registran tres tipos de causas de muerte: directa, intermedia y básica. Desde el año 1997, se siguen las normas de las Clasificación Internacional de Enfermedades, décima revisión (CIE-10) para la codificación y selección de causa de muerte. En el análisis, se tuvo en cuenta la causa de muerte básica, definida como la enfermedad que inicia la cadena de acontecimientos patológicos que conducen a la muerte. Esta causa es la que considera el Ministerio de Salud Pública de Argentina (MSAL) para elaborar estadísticas de mortalidad atribuida a alguna enfermedad. Con base en las normativas de Vigilancia Epidemiológica del Ministerio de Salud, se consideraron muertes atribuidas a gripe los códigos de causa de muerte básica entre J09 y J18.9 (gripe y neumonías) y J22X (infección aguda sin especificar de las vías respiratorias inferiores).

Número de vacunas aplicadas y porcentaje de población vacunada o cobertura de la vacuna, por año y provincia: estos datos fueron proporcionados por el PNI y de datos oficiales publicados (6, 7).

Número de personas según edad y provincia obtenido de las publicaciones del Censo Nacional de las Personas realizado por Instituto Nacional de Estadísticas y Censos de Argentina (INDEC) para los años 2001 y 2010. Se utilizaron las proyecciones por provincia y grupo de edad publicadas por INDEC (9, 10). La población mensual se estimó mediante una proyección lineal entre el año en estudio y el posterior, según provincia y grupo de edad.

Población en niños 6 a 23 meses: a partir de la población censada en 2001 y en 2010, se obtuvo el total de niños de uno y dos años. Para estimar el grupo entre 6 a 12 meses de edad, se calculó la mitad de la población de niños menores de un año. Luego se sumaron a los niños que tenían un año, y entre los puntos 2001 y 2010 se realizó una proyección lineal simple entre 2001 y 2016.

Tipo de cepa viral circulante y cepa vacunal colocada: se realizó una base de dato con la información disponible en línea del MSAL y del Grupo Colaborativo GROG (11). Se consideró dicotómica y afirmativa cuando había por lo menos una coincidencia.

Se calcularon las tasas de mortalidad específica por IRAG utilizando la siguiente fórmula:

Donde:

M = tasa de mortalidad por IRAG.

F = número de fallecidos por IRAG según edad, grupo de edad, provincia, mes y año.

p = número de personas de igual edad, provincia, mes y año.

Se calcularon las siguientes tasas específicas de IRAG:

Tasa de mortalidad por periodo en la población global y en personas mayores de 64 años (gerontes). Se consideró como período anterior a la vacunación (AV) el comprendido entre los años 2002-2009 y como período posterior al inicio de la vacunación (DV) al que va desde 2010 hasta 2016. Se estratificó por provincias y se calculó la tasa de mortalidad bruta y la ajustada por edad (TAE) por método directo y comparación con la población mundial (12). Se estimaron los intervalos de confianza de 95% (IC95%) de las tasas ajustadas y se compararon las tasas AV y DV con la prueba de la χ2; se consideró P < 0,05 como el nivel de significancia.

Tasa de mortalidad por período de los niños entre 6 a 23 meses, calculada por cada 100 000. Se calcularon las tasas de mortalidad bruta de los períodos AV y DV de cada provincia y se las comparó con la prueba de la χ2, se consideró P < 0,05 como el nivel de significancia.

Tasa de mortalidad mensual en la población global y en gerontes por cada 100 000 habitantes y personas mayores de 64 años, respectivamente. Se estandarizó de manera indirecta por edad y mes de cada provincia y año de la serie de tiempo con la metodología propuesta por la OMS (12). Se utilizó como tasa estándar la estimada para el período AV, calculada a partir de la media aritmética del período 2002-2009 por mes y población en estudio (numerador) y, en el denominador, la población de mitad del período (año 2005), según mes y población específica. Se realizó una estandarización indirecta debido al escaso número de casos mensuales. Se utilizaron los datos para realizar el análisis de ST.

Tasa bruta de mortalidad mensual de niños de 6 a 23 meses por provincia y total del país por cada 100 000 niños de 6 a 23 meses.

Con las tasas mensuales brutas en niños y las TAE en la población global y gerontes, se realizaron los análisis de ST. Se descompuso y graficó la ST de cada uno de los grupos de edad (datos originales, data) y los componentes estacionales (seasonal), tendencia (trend) y residuos (remander) y se realizó la comparación relativa de sus componentes. Se utilizó el paquete Stats del programa estadístico R® (13).

Se realizaron gráficos de cajas del número de fallecidos por mes e histogramas de edad (por año de vida) entre los períodos AV y DV; en cada uno de ellos se calcularon las varianzas (Var) y las desviaciones estándares (DE). En las variables cuantitativas se sacaron sus medidas de tendencia central y se compararon con prueba de χ2; se consideró significativo un P < 0,05.

Se seleccionó el grupo de niños entre 6 a 23 meses para evaluar el efecto en esta población específicamente, ya que se tuvo acceso al porcentaje de la población vacunada en cada año. Debido a que no hay en el país un patrón oro para la comparación de las tasas de mortalidad por mes y grupo de edad, fue preciso realizar modelos estadísticos para evaluar su comportamiento. Para ello, se estimó el efecto de la vacunación con un modelo estadístico de tipo modelos aditivos generalizado (GAM, por sus siglas en inglés) (14), que son una extensión de los modelos lineales generalizados (GLM, por sus siglas en inglés). La ventaja es que permiten la inclusión de una función de suavizamiento no paramétrica para una o más variables independientes. Ello es útil cuando no hay una relación lineal entre la variable independiente y la variable dependiente (15). Si bien este método es muy utilizado para el modelaje de datos que relacionan patologías respiratorias contaminantes del aire (16), recientemente se han se han publicado estudios relacionados con influenza (15, 17).

Se consideró:

Donde:

Y = es la variable respuesta (niños fallecidos por IRAG).

X = son las variables predictoras: mes, año, dosis de vacunas en el año, porcentaje de población vacunada, período DV y coincidencia entre la cepa de la vacuna y la cepa viral en el año (variable dicotómica).

= es una constante.

ε = es el error.

= se estiman mediante el uso de suavizadores.

Se utilizó una distribución binomial negativa debido a la diferencia entre la varianza y la media de los datos. Los grados de libertad se asignaron de manera automática, considerando el mejor valor posible. Se utilizaron suavizadores de tipo spline (s) para estimar funciones no paramétricas. El procedimiento del GAM contempló el valor de significancia (P < 0,05) y se seleccionó el modelo teniendo en cuenta el menor valor de máxima verosimilitud restringida (REML, por sus siglas en inglés) y el principio de parsimonia o navaja de Ockham. El modelo elegido arrojó un REML de 538,06 y una varianza explicada de 60,3%, relacionando en función del mes, el número de niños fallecidos, la población y compara cada año con los valores del año 2002 (intercepto):

Donde:

Caso: el número de niños 6 a 23 meses de edad, fallecidos por IRAG.

f (mes): al mes de la defunción.

factor (año): a cada año.

offset = log (pn): enlace de la función del logaritmo de la población de niños.

Los datos se analizaron con el programa R® y el paquete mgcv (18). Para el modelo seleccionado se calculó el riesgo relativo (RR) y el IC95% comparando el riesgo de morir en el año 2002 con los otros años y comparando los meses del año 2002 con los meses de los otros años. Se calculó la media aritmética de los períodos AV y DV de los fallecidos por mes en el modelo planteado.

RESULTADOS

De las 4 691 147 personas fallecidas entre los años 2002 y 2016, 15,2% fue por patología respiratoria (J00-J98.9) y el 6,4% (301 747) por IRAG; en especial por neumonía (J18), que representó 92,3% (31 219) de las IRAG denunciadas. El promedio anual de personas fallecidas por IRAG tuvo un incremento entre los períodos AV y DV, con 15 869 y 24 971 casos, respectivamente. La media de edad fue de 78 años, la mediana fue de 82 años y la mitad de los fallecidos tenía entre 72 y 88 años. En el histograma de edad se pudo observar una mayor concentración de datos en el período DV (edad del período AV: Var 372,5 y DE 19,3; edad período DV: Var 275,9 y DE 16,6). En comparación con el histograma del período AV, hubo una disminución del porcentaje de los casos en menores de dos años (P < 0,05) (figura 1), que pasó de representar el 2,34% del total de fallecidos por IRAG en el período AV (126 701) a ser el 0,99% de los fallecidos en el período DV (174 285).

FIGURA 1

Histograma comparativo entre el período antes de la vacunación y el período después de la vacunación de muertes por IRAG, Argentina

AV (período antes de la vacunación, 2002-2009, n = 126 701); DV (período después de la vacunación, 2010-2016, n = 174 285 casos).

Las TAE de la población total (tasas globales) y en personas mayores de 64 años (tasas en gerontes) aumentaron de manera significativa excepto en la Ciudad Autónoma de Buenos Aires (C.A.B.A.) (cuadro 1). La tasa bruta de mortalidad en niños entre 6 a 23 meses evidenció una disminución significativa en 14 de las 24 jurisdicciones (provincias y C.A.B.A.) del país (cuadro 1).

CUADRO 1

Tasas de mortalidad (e IC95%) por IRAG cada 100 000 habitantes, por provincia, grupo de edad, comparativo de los períodos anterior y posterior a la vacunación en Argentina

Provincia o distrito	Tasa de mortalidad global									Tasa de mortalidad en gerontes									Tasa de mortalidad en niños
	Período AV				Período DV				P	Período AV				Período DV				P	Período AV		Período DV		P
	% de la población	TB	TAE	IC95%	% de la población	TB	TAE	IC95%		% de la población	TB	TAE	IC95%	% de la población	TB	TAE	IC95%		% de la población	TB	% de la población	TB
Argentinaa	100b	41	25,3	(25,2-25,5)	100c	60	31,6	(31,4-21,8)	< 0,01	100,00	338	259,8	(258,2-261,4)	100,00	468	328,6	(326,9-330,4)	< 0,01	100,0	14,2	100,00	10,0	< 0,01
Buenos Aires	37,00	39	23,6	(23,3-23,8)	38,30	60	32,0	(31,8-32,3)	< 0,01	40,20	300	233,1	(230,6-235,5)	39,40	450	330,0	(327,1-332,8)	< 0,01	37,1	13	39,50	8,1	< 0,01
Córdoba	8,40	32	17,9	(17,4-18,3)	8,40	59	29,8	(29,2-30,3)	< 0,01	13,00	241	185,3	(180,8-189,7)	13,50	437	322,0	(316,1-327,8)	< 0,01	7,9	6,5	7,70	4,7	< 0,01
Santa Fe	8,20	39	21,1	(20,7-21,6)	8,00	59	27,6	(27,1-28,2)	< 0,01	9,50	288	215,8	(211,1-220,6)	0,70	428	298,7	(293,1-304,3)	< 0,01	7,3	13,1	7,10	7,2	< 0,01
C.A.B.A.	7,80	108	34,9	(34,4-35,4)	7,60	107	26,7	(26,2-27,1)	< 0,01	9,30	584	381,2	(376,0-386,5)	1,70	515	280,9	(276,0-285,8)	< 0,01	5,0	10,1	5,10	2,2	< 0,01
Mendoza	4,40	41	25,4	(24,7-26,1)	4,10	43	22,4	(21,7-23,0)	0,12	4,20	356	277,7	(269,6-285,8)	0,90	349	250,8	(243,2-258,4)	0,24	4,5	11	4,70	6,8	0,3
Tucumán	3,70	28	22,7	(21,9-23,5)	3,70	66	46,7	(45,6-47,8)	< 0,01	2,80	311	260,5	(250,5-270,4)	1,80	676	529,0	(515,4-542,5)	< 0,01	4,1	8	3,70	4,5	< 0,01
Entre Ríos	3,20	32	20,4	(19,7-21,2)	3,10	50	28,0	(27,1-28,8)	< 0,01	3,10	264	206,5	(198,2-214,7)	9,10	392	294,7	(285,0-304,3)	< 0,01	3,2	5,9	2,80	5,4	< 0,01
Salta	3,00	23	23,7	(22,7-24,7)	3,10	48	42,7	(41,4-43,9)	< 0,01	1,90	249	217,7	(206,4-228,9)	3,10	515	422,1	(407,5-436,7)	< 0,01	3,9	29,6	3,70	29,4	< 0,01
Misiones	2,70	23	25,1	(23,9-26,2)	2,80	41	38,7	(37,4-40,1)	< 0,01	1,60	251	226,7	(213,9-239,5)	0,90	460	388,7	(373,1-404,3)	< 0,01	3,6	23,3	3,30	24,4	< 0,01
Chaco	2,70	32	33,6	(32,2-34,9)	2,70	50	44,8	(43,4-46,3)	< 0,01	1,70	347	322,2	(307,4-337,1)	1,20	543	459,1	(442,4-475,7)	< 0,01	3,1	11,7	2,70	7,1	0,12
Corrientes	2,50	29	25,6	(24,6-26,7)	2,40	44	32,3	(31,1-33,5)	< 0,01	1,80	313	259,3	(247,1-271,6)	0,90	418	329,6	(316,0-343,2)	< 0,01	2,9	41,6	2,50	23,4	< 0,01
Santiago del Estero	2,20	31	28,4	(27,2-29,7)	2,20	42	34,3	(33,0-35,7)	< 0,01	1,50	319	274,3	(260,4-288,3)	0,60	423	348,0	(332,8-363,2)	< 0,01	2,8	25,6	2,60	9,2	< 0,01
San Juan	1,70	64	48,3	(46,7-50,0)	1,70	81	51,6	(49,9-53,3)	0,14	1,40	661	551,4	(531,2-571,6)	4,10	776	612,7	(591,9-633,5)	0,07	2,0	20,6	1,90	11,1	< 0,01
Jujuy	1,70	18	17,8	(16,6-18,9)	1,60	44	36,2	(34,6-37,7)	< 0,01	1,10	199	177,0	(163,6-190,4)	1,00	466	385,8	(367,6-404,0)	< 0,01	1,9	9,7	1,80	6,3	0,13
Río Negro	1,50	19	14,8	(13,8-15,8)	1,50	25	16,2	(15,2-17,1)	0,39	1,20	189	158,8	(147,0-170,6)	1,70	233	180,7	(169,3-192,1)	0,23	1,6	7,9	1,60	7,8	0,2
Neuquén	1,40	22	22,5	(21,0-23,9)	1,40	31	24,9	(23,6-26,3)	0,21	0,80	318	267,9	(249,2-286,5)	0,90	374	299,4	(281,9-316,9)	0,19	1,5	10,4	1,60	7,1	0,26
Formosa	1,30	20	21,9	(20,3-23,4)	1,40	62	57,6	(55,3-59,9)	< 0,01	0,80	194	177,4	(161,4-193,3)	0,40	632	541,7	(516,2-567,2)	< 0,01	1,6	37,6	1,40	48,0	< 0,01
Chubut	1,20	33	30,0	(28,3-31,7)	1,20	71	51,3	(49,3-53,3)	< 0,01	0,80	371	327,3	(306,4-348,3)	1,40	786	627,6	(601,8-653,4)	< 0,01	1,3	9,9	1,40	5,0	0,09
San Luis	1,10	55	44,1	(42,0-46,2)	1,10	76	51,6	(49,5-53,6)	< 0,01	0,80	558	471,9	(447,5-496,3)	2,00	721	579,5	(554,6-604,3)	< 0,01	1,1	12,2	1,10	11,4	0,59
Catamarca	0,95	24	20,3	(18,7-21,8)	1,00	28	21,1	(19,6-22,6)	0,66	0,70	263	212,0	(194,0-229,9)	8,90	291	225,5	(207,8-243,1)	0,51	1,0	6	0,90	3,1	0,03
La Rioja	0,84	39	38,2	(35,8-40,6)	0,83	59	49,7	(47,0-52,3)	< 0,01	0,50	466	398,2	(370,1-426,2)	1,60	616	508,4	(477,6-539,2)	< 0,01	0,9	14,5	0,80	14,8	0,15
La Pampa	0,83	47	26,8	(25,2-28,4)	0,82	75	38,1	(36,2-40,1)	< 0,01	0,90	393	301,2	(282,7-319,7)	1,30	562	412,2	(391,3-433,1)	< 0,01	0,8	8,2	0,80	16,6	0,14
Santa Cruz	0,56	28	32,6	(29,7-35,5)	0,66	38	38,8	(36,0-41,6)	< 0,01	0,30	347	338,6	(302,7-374,5)	0,10	468	413,3	(378,4-448,2)	< 0,01	0,7	7,2	0,90	4,8	0,11
Tierra del Fuego	0,30	15,8	29,3	(24,5-34,2)	0,35	26	35,3	(30,9-39,6)	< 0,01	0,10	331	349,5	(283,5-415,4)	2,90	405	393,0	(336,4-449,6)	0,11	0,9	3,9	0,30	0,0	0,05

Población en base a las proyecciones realizadas por el Instituto nacional de Estadísticas y Censos (INDEC) de los censos poblacionales de 2001 y 2010.

38 801 324 habitantes.

41 795 033 habitantes.

TB, tasa bruta; TAE, tasa ajustada por edad; AV, antes de la vacunación (2002-2009); DV, después de la vacunación (2010-2016); C.A.B.A., Ciudad Autónoma de Buenos Aires; IC95%, intervalo de confianza de 95%.

La descomposición de la ST de las tasas mensuales mostró, en los tres grupos estudiados, un marcado comportamiento estacional, con picos en el mes de julio (figura 2). La tendencia de la mortalidad en los niños parece descender a partir del año 2 004; mientras que en los adultos mayores y población global parece ascender en los últimos años. Luego de eliminado el efecto de la tendencia y estacionalidad, se evidenció un pico de la mortalidad en 2009 en los niños y en 2016 en los adultos mayores (figura 2).

FIGURA 2

Gráfico de caja de los casos de IRAG según el mes de defunción y ST de la tasa de mortalidad (por cada 100 000) por IRAG, según grupo en Argentina (2002-2016)

Para estimar el riesgo de morir en niños entre 6 y 23 meses se realizaron varios modelos estadísticos GAM. Si bien se incorporaron diferentes variables como el número y el porcentaje de vacunación y la coincidencia de la cepa viral y la de la vacuna del mismo año, el mejor modelo fue el más simple. Este comparó el riesgo de morir en el año 2002 con los otros años en el grupo de niños entre 6 a 23 meses, y fue significativamente menor en un año AV, el 2008, y en cuatro años del período DV (2010, 2012, 2014 y 2015) (cuadro 2). En base a los valores obtenidos del modelo, se estimó que el promedio mensual de niños fallecidos por IRAG disminuyó alrededor de 3,5 entre el período AV y DV, pasó de 11,89 a 8,19 defunciones mensuales, respectivamente (figura 3A). Se analizó el parámetro estacional del modelo planteado que comparó los meses de 2002 con los otros años; se observó que el RR de morir por IRAG en los niños de 6 a 23 meses fue de 0,46 (IC95%: 0,36-0,57) durante abril. A fines de junio el riesgo llegó a casi el doble (RR = 2,7; IC95%: 2,3-3,3) (figura 3B). El RR de morir por IRAG en este grupo de edad se mostró aumentado entre finales de mayo y julio, en comparación con el riesgo durante el año 2 002.

CUADRO 2

Número de casos y riesgo de fallecer por IRAG en niños de 6 a 23 meses comparando el riesgo en el año 2002 y cobertura de vacunación (n = 1 814)

Año	Número de casos	Población (hab.)a	RR	IC95%	Población vacunada
					Primera dosis		Segunda dosis
					%	Número de dosisb	%	Número de dosisb
2002	141	971 730	(Intercepto)		NA		NA
2003	208	977 113	1,3c	0,9-1,9	NA		NA
2004	132	982 496	0,9c	0,6-1,3	NA		NA
2005	120	987 878	0,8c	0,6-1,2	NA		NA
2006	126	993 261	0,9c	0,6-1,3	NA		NA
2007	141	998 644	0,9c	0,6-1,3	NA		NA
2008	98	1 004 027	0,6d	0,4-0,9	NA		NA
2009	159	1 009 409	0,9c	0,7-1,4	NA		NA
2010	98	1 014 792	0,7c	0,5-1,0	86,7	949 075	57,7	629 179
2011	107	1 020 175	0,7c	0,5-1,0	73,0	792 877	36,4	395 569
2012	78	1 025 557	0,5e	0,4-0,8	75,7	811 843	58,2	624 164
2013	109	1 030 940	0,7c	0,5-1,0	78,1	848 663	46,0	511 567
2014	89	1 036 323	0,6e	0,4-0,9	73,8	801 219	50,4	547 174
2015	79	1 041 706	0,5e	0,3-0,8	81,3	882 644	60,3	654 655
2016	129	1 047 088	0,8c	0,6-1,2	78,6	853 331	59,9	650 312

Estimaciones con base en el Censo 2001 y 2010.

Datos publicados por el Ministerio de Salud de la Nación Argentina.

P > 0,05.

P < 0,05.

P < 0,01.

RR, riesgo relativo; IC95%, intervalo de confianza de 95%; NA, no aplica.

FIGURA 3

Modelo estadístico planteado

DISCUSIÓN

En la ST estudiada, la causa de muerte atribuible a IRAG más frecuente fue la neumonía por microorganismo no especificado (J18 de CIE-10), al igual que en otros estudios realizados en el país (19). Las tasas de mortalidad por IRAG mostraron un ascenso entre el período AV y DV en la población total de Argentina y en los mayores de 64 años. Esto coincide con la tendencia en ascenso observada en la descomposición de la serie histórica de los 15 años estudiados. Las tasas observadas en los mayores de 64 años fueron similares a las encontradas en otros estudios en el país (20) y en Brasil (21), pero mayores a las reportadas en Chile y otros países del hemisferio sur (22). Una revisión sistemática del Grupo Colaborativo Cochrane (23) estimó que los adultos mayores vacunados podrían tener un menor riesgo de enfermarse de gripe en la misma estación comparados con aquellos no vacunados (RR: 0,42; IC95%: 0,27-0,66). No obstante, el estudio no pudo estimar el efecto de la vacunación en la mortalidad debido a la poca información en los estudios seleccionados. En Argentina, no hay una cifra oficial del porcentaje de población de mayores de 64 años vacunados, debido a que se realiza en el sector público, privado y de seguridad social. Recién en el año 2 013 se realizó la primera encuesta poblacional sobre la cobertura de vacunación en el país (24), en la que se reportó que se había vacunado 51,6% (IC95%: 49,6-53,6) del total de personas encuestadas y 55,6% (IC95%: 52,8-58,3) de las personas mayores de 64 años sin comorbilidad. Si bien la vacunación antigripal en mayores de 64 años se realiza de forma obligatoria desde el año 2003, este estudio pretendió mostrar el posible efecto de la disminución de los susceptibles con la inmunidad de grupo.

Con respecto a las tasas de mortalidad bruta en niños entre 6 a 23 meses, se observó una disminución en la mayoría de las provincias entre el periodo AV y DV, similar a lo observado en una revisión sistemática del grupo Cochrane (25).

Los modelos estadísticos semiparamétricos, como el GAM, pueden reemplazar a los modelos armónicos estacionales convencionales para estimar las enfermedades respiratorias y la gripe. Un estudio realizado en Australia (15) modeló con GAM para evaluar la mortalidad atribuible a influenza en una serie histórica, como el presente trabajo.

En el análisis de la ST y en el modelo GAM realizado, el componente estacional fue bien marcado, con un claro incremento de la mortalidad entre junio y setiembre, con pico en julio, similar al comportamiento observado en otros estudios (26-29).

Este estudio tiene limitaciones, sobre todo por ser un estudio ecológico y la poca información pública disponible en el país. En primer lugar, se utilizó un grupo de enfermedades que se le puede atribuir a gripe, pero no se confirmó la etiología por virus influenza de las IRAG, que pueden haber sido causadas por otros virus respiratorios. En el país, la vigilancia de virus respiratorios se realiza en algunas provincias, en especial en la población pediátrica. Hay estudios que muestran que el comportamiento es diferente en niños y en mayores de 64 años (30). En segundo lugar, la fuente de datos fue el IED. Si bien se trata de datos oficiales y denunciados por profesionales, puede tener un margen de error en el registro, con sobreestimación o subestimación de las muertes por gripe. Otras medidas sanitarias no pudieron ser evaluadas, como el lavado de manos, la mejora en la atención de los pacientes, el inicio temprano de la medicación antigripal, la educación sanitaria, el avance tecnológico y el terapéutico durante los 15 años estudiados, la incorporación de la vacuna antineumocócica en el año 2012, las temperaturas en las diferentes regiones del país y la fecha de inicio de la vacunación con las nuevas cepas cada año.

Conclusiones

Las tasas de mortalidad por IRAG se mantuvieron similares en los adultos mayores de 64 años y en la población general entre los períodos AV y DV. En cambio, en los niños, se observó un claro descenso. Se estima que los niños entre 6 a 23 meses tuvieron menor riesgo de morir por IRAG en los tres años posteriores a la vacunación, en comparación con el año 2002.

Si bien la vacuna antigripal está disponible durante todo el año, la marcada estacionalidad que presenta la enfermedad hace necesario planificar la distribución de las partidas de la nueva vacuna desde principios del año, así podría estar disponible unos meses antes del ascenso de la epidemia anual.

Este trabajo pretende delinear un panorama general del efecto de la vacunación antigripal en la población luego de la implementación de una medida poblacional de salud pública. Esto debe ir acompañado de la educación de la población y profesionales de la salud, del avance tecnológico en la atención y tratamiento y un buen sistema de vigilancia epidemiológica que genere datos para la toma de decisiones políticas.