Cuantitativo, observacional, retrospectivo, descriptivo comparativo y de análisis de series temporales

Las personas infectadas con VIH en Cali procedentes y residentes en Cali, entre los años 2012 y 2024.

Todas las personas confirmadas con VIH por laboratorio, reportadas al SIVIGILA de Cali, código INS No. 850, entre 2012-2019 y 2021-2024.

Todos los casos registrados en las bases de datos del SIVIGILA con código INS No. 850 de los años 2012 a 2019 y 2021 a 2024, que cumplieron los criterios de caso confirmado por laboratorio, según protocolos del INS, SIVIGILA y la guía de práctica clínica para el VIH de Colombia ²⁰.

1. Casos registrados sin número de identificación

2. Casos sin información en todas las variables del estudio.

3. Casos reportados en el año 2020.

Edad: tipo cuantitativa; nivel de medición razón; valores años. Año de reporte: tipo cualitativa; nivel de medición ordinal; valores numero del año calendario. Conteo de casos: tipo cuantitativa; nivel de medición razón; valores número de casos. Semana: tipo cualitativa; nivel de medición ordinal; valores número de la semana.

Las siguientes variables son de tipo cualitativo; nivel de medición, nominal; valores, dicotómicas: sexo heterosexual, otras orientaciones sexuales (agrupó a Homosexual y bisexual), uso de drogas ilícitas inyectables, materno-fetal y del lactante (agrupó a transmisión madre-hijo y lactancia), procedimientos percutáneos de baja complejidad (agrupó a piercing, tatuajes, centros estéticos y acupuntura), procedimientos de mediana y alta complejidad (nivel II y III, agrupó al accidentes de trabajo, las transfusiones sanguíneas, las hemodiálisis, los trasplantes de órganos), migrante, desmovilizados, desplazados, víctimas de la violencia, indígenas, discapacidad física, personas en centro psiquiátrico, privados de la libertad, gestantes, población infantil en centro de bienestar estatal (ICBF), madre comunitaria.

Se realizó una descripción univariada con mediana y rango intercuartílico para las variables cuantitativas y porcentajes para variables cualitativas.

Se tabularon las distribuciones y luego se visualizaron mediante gráficas de nube de puntos, diagrama de cajas y bigotes. Se eliminaron los registros del año 2020 por el confinamiento pandémico de COVID-19, se creó la variable confinamiento y las categorías preconfinamiento (PreC, de 2012-2019) y el posconfinamiento (PostC 2021-2024). En el análisis comparativo se utilizó la prueba no paramétrica U de Mann-Whitney para comparar las medianas del preC y el postC. En el análisis de tendencia de la serie se utilizó el método de análisis de las pendientes de los segmentos generados por los puntos de inflexión, mediante regresión lineal ²¹. Para el análisis multivariado se usó la técnica estadística de regresión de Poisson ²², de conteo de casos y la variable respuesta fue el número de casos semanales reportados por año y como variable explicatoria, cada determinante. Se obtuvieron los estimados crudos del riesgo relativo (RR). Para la construcción de un modelo preliminar se seleccionaron las variables significativas del análisis comparativo (RR mayor que 1 y con intervalo de confianza del 95%, que no cruza el valor 1 y se incluyeron como variables explicatorias de la variable resultado que fue el conteo de casos. Se obtuvieron estimados de RR ajustados, después de controlar para confusión. Para el análisis se usó el software R versión 4.3.1 ²³.

Debido a que entre 2012 y 2014 se cambió la versión del instrumento de recolección de los datos (ficha epidemiológica v.2007.1 vs v2022.11.), se tuvieron que recodificar las categorías para resolver este problema. Se seleccionaron las variables y se inspeccionó la completitud de cada una de ellas, se depuraron las inconsistencias y se eliminaron los duplicados. Se cruzaron las bases de datos del SIVIGILA, los certificados de defunción y de las instituciones de salud, correspondientes al periodo 2012 a 2023, se adicionó la data del año 2024 y se eliminó la de año 2020.

El protocolo fue aprobado por los Comités de Ética y de Investigación de la Secretaría de Salud de Cali, con clasificación de riesgo mínimo, a la luz de la resolución número 8430 de octubre 4 de 1993, del Ministerio de Salud de Colombia y la “declaración de Helsinki de la Asociación Médica Mundial”, de octubre de 2024. Este trabajo hace parte del trabajo de grado de maestría patrocinado por el MSP de Colombia y el Ministerio de Sanidad de España y dirigido por la Universidad Rey Juan Carlos de España (URJC).

Cali es una ciudad de Colombia, América del sur, con 2,285,099 habitantes ²⁴. Se reportaron 16,954 casos nuevos de VIH entre 2012 y 2024, excluyendo los registros de 2020. Se destaca el aumento monotónico entre 2015 (1,019 casos) y 2019 (1,971 casos) y el aumento a 2,337 casos en 2022, año con el máximo número de casos reportados. Finalmente, se disminuyó el reporte de casos hasta 1,469 casos en el año 2024. El PreC (2012-2019) agregó 9,918 casos y el PostC (2021-2024) agregó 7,036 casos.

La mediana de edad fue 31 años (RIC: 25-43). La proporción de casos del sexo masculino fue del 78% en toda la serie. La proporción de mujeres disminuyó de 24 a 20% en PostC. La vía de transmisión heterosexual fue la más reportada (54%), seguida de otra orientación sexual (43%) y uso de drogas ilícitas inyectables (0.8%). Los determinantes, otra orientación sexual (8.0 %), uso de drogas ilícitas inyectables (1.0 %), vía percutánea nivel I (0.2 %) y percutáneo nivel II (0.1%) aumentaron su proporción en PostC. Según la condición de vulnerabilidad, ser migrante fue el determinante más reportado con 4.1% y el único que aumentó la proporción de reporte entre PreC y PostC (5.4%) (Tabla 1).

Tabla 1

Determinante		SEXO		p-value**
		F	M
	N: 16,954	N = 3,757*	N = 13,197*
año
2012	897 (5.3%)	250 (6.7%)	647 (4.9%)	<0.001
2013	720 (4.2%)	207 (5.5%)	513 (3.9%)
2014	697 (4.1%)	165 (4.4%)	532 (4.0%)
2015	1,019 (6.0%)	269 (7.2%)	750 (5.7%)
2016	1,362 (8.0%)	319 (8.5%)	1,043 (7.9%)
2017	1,406 (8.3%)	347 (9.2%)	1,059 (8.0%)
2018	1,846 (11%)	407 (11%)	1,439 (11%)
2019	1,971 (12%)	418 (11%)	1,553 (12%)
2021	1,700 (10%)	345 (9.2%)	1,355 (10%)
2022	2,337 (14%)	456 (12%)	1,881 (14%)
2023	1,530 (9.0%)	285 (7.6%)	1,245 (9.4%)
2024	1,469 (8.7%)	289 (7.7%)	1,180 (8.9%)
Edad (años)	31 (25, 43)	34 (26, 45)	31 (25, 42)	<0.001
Grupo_Poblacional
Discapacidad física/Mental	46 (0.3%)	13 (0.3%)	33 (0.3%)
Gestantes/Madre comunitaria/ población infantil ICBF	288 (1.7%)	263 (7.0%)	25 (0.2%)
Desmovilizado	2 (<0.1%)	1 (<0.1%)	1 (<0.1%)
Desplazado	23 (0.1%)	8 (0.2%)	15 (0.1%)
Indigena	233 (1.4%)	63 (1.7%)	170 (1.3%)
Victima de Violencia	13 (<0.1%)	3 (<0.1%)	10 (<0.1%)
Migrante	428 (2.5%)	74 (2.0%)	354 (2.7%)
Privado de Libertad	109 (0.6%)	6 (0.2%)	103 (0.8%)
Otros	15,459 (91%)	3,192 (85%)	12,267 (93%)
NA	353 (2.1%)	134 (3.6%)	219 (1.7%)
Modo de Transmisión
Heterosexual	9,097 (54%)	3,443 (92%)	5,654 (43%)	<0.001
Otra orientación sexual	7,282 (43%)	162 (4.3%)	7,120 (54%)
Usuarios drogas ilicitas	138 (0.8%)	19 (0.5%)	119 (0.9%)
materno-fetal y del lactante	60 (0.4%)	36 (1.0%)	24 (0.2%)
§ Percutanea nivel II y III	51 (0.3%)	17 (0.5%)	34 (0.3%)
‡ Percutanea I	36 (0.2%)	9 (0.2%)	27 (0.2%)
NA	290 (1.7%)	71 (1.9%)	219 (1.7%)

* Frequency (%)

La mediana de casos reportados semanalmente (CRS) aumentó a partir del año 2013 (13 casos) hasta el 2019 (35 casos). Aumentó del año 2021 (30 casos) a 2022 (40 casos), para finalmente disminuir hasta 29 casos en el 2024 (Figura 1).

Figura 1

El análisis de tendencia de la serie evidencia el aumento del riesgo anual desde 2014 hasta 2018, correspondiente al segundo segmento de la serie, con pendiente positiva y significativa estadísticamente. De 2018-2024 la pendiente es negativa, pero no es significativa estadísticamente en este 3er segmento. Figura 2.

Figura 2

El valor de la Mediana de CRS en el PreC (2012-2019) fue de 23 casos y del PostC (2021-2024) fue de 31 casos. (U de Mann-Whitney= 24886.5, p <0.0001 (Figura 3)

Figura 3

El análisis bivariado mediante la regresión de Poisson halló que los determinantes: ser varón y practicas heterosexuales inseguras, tuvieron un RR crudos mayores que 1 y con el intervalo de confianza estadísticamente significativo. Ser menor de 32 años, migrante, indígena, las vías de transmisión percutáneas como uso de drogas ilícitas inyectables, practicas cosméticas de baja complejidad, procedimientos médico quirúrgicos de mediana y alta complejidad y practicas no heterosexual inseguras tuvieron RR crudos menores que 1 y con intervalo de confianza estadísticamente significativo (Tabla 2). confinamiento no arrojo RR crudo estadísticamente significativo.

Tabla 2

Factor de Riesgo	IRR Crudo 1	95% CI	p-value
Sexo	1.89	1.82, 1.96	<0.001
No heterosexual	1.57	1.52, 1.62	<0.001
Heterosexual	0.80	0.77, 0.82	<0.001
Desmovilizado	0.33	0.22, 0.47	<0.001
Edad ( <32 /≥32)	0.96	0.93, 0.99	<0.001
Uso de drogas inyectadas	0.32	0.27, 0.38	<0.001
Procedimiento percutáneo de baja complejidad	0.30	0.21, 0.41	<0.001
Procedimiento percutáneo de mediana y alta complejidad	0.31	0.23, 0.40	<0.001
Materno-fetal y del lactante	0.31	0.24, 0.40	<0.001
Migrante	0.38	0.35, 0.41	<0.001
Indigena	0.32	0.28, 0.35	<0.001
Discapacidad	0.31	0.24, 0.39	<0.001
Desplazado	0.31	0.24, 0.40	<0.001
Privado de la libertad	0.34	0.29, 0.40	<0.001
Gestante	0.35	0.31, 0.39	<0.001
Población infantil ICBF	0.33	0.22, 0.47	<0.001
Institución psiquiatrica	0.32	0.24, 0.43	<0.001
Victima violencia	0.31	0.23, 0.41	<0.001
post confinamiento	0.99	0.96, 1.02	0.6

SIVIGILA JHRP

Los factores de riesgo significativos fueron incluidos como variables predictoras en un modelo de regresión de Poisson Preliminar con variable respuesta el conteo de casos por año (Tabla 3) Las variables del modelo preliminar con RR ajustados mayores que uno y significativos, se seleccionaron para el modelo definitivo que revelo el RR del modo de transmisión a través de prácticas heterosexuales inseguras fue 2.87 (IC 95%: 2.46-3.37); El RR asociado con prácticas no heterosexuales inseguras fue 3.34 (IC 95%: 2.86-3.92); El RR asociado con ser varón, 1.76 (IC 95%: 1.70-1.84) (Tabla 4).

Tabla 3

Characteristic	POISSON			QUASI-POISSON			BNM_NGT
	IRR	95% CI	p-value	IRR	95% CI	p-value	IRR	95% CI	p-value
Edad	1.00	0.97, 1.03	0.8	1.00	0.95, 1.04	0.8	0.97	0.93, 1.01	0.10
Sexo	1.76	1.70, 1.84	<0.001	1.76	1.67, 1.86	<0.001	1.74	1.65, 1.83	<0.001
Heterosexual	2.87	2.46, 3.37	<0.001	2.87	2.33, 3.59	<0.001	2.86	2.40, 3.42	<0.001
Otra orientación sexual	3.34	2.86, 3.92	<0.001	3.34	2.71, 4.18	<0.001	3.29	2.76, 3.95	<0.001
Desmovilizado	3.26	2.13, 4.77	<0.001	3.26	1.78, 5.45	<0.001	25.6	10.3, 67.0	<0.001
Madre comunitaria	115	6.14, 634	<0.001	115	1.04, 1,068	0.001	10.9	0.54, 75.7	0.035
Post Confinamiento	1.15	1.12, 1.19	<0.001	1.15	1.10, 1.20	<0.001	1.13	1.08, 1.18	<0.001
Uso de drogas ilícitas inyectables	0.97	0.77, 1.22	0.8	0.97	0.71, 1.33	0.9	1.03	0.79, 1.33	0.8
procedimientos percutáneos de baja complejidad	0.75	0.52, 1.07	0.12	0.75	0.44, 1.21	0.3	0.78	0.52, 1.16	0.2
procedimientos de mediana y alta complejidad	0.79	0.57, 1.08	0.15	0.79	0.50, 1.21	0.3	0.81	0.56, 1.14	0.2
materno-fetal y del lactante	0.94	0.69, 1.26	0.7	0.94	0.62, 1.41	0.8	0.91	0.65, 1.27	0.6
Migrante	0.31	0.29, 0.33	<0.001	0.31	0.28, 0.34	<0.001	0.32	0.29, 0.35	<0.001
Indigena	0.34	0.31, 0.39	<0.001	0.34	0.29, 0.40	<0.001	0.35	0.31, 0.40	<0.001
Discapacitado	0.36	0.28, 0.45	<0.001	0.36	0.25, 0.49	<0.001	0.35	0.26, 0.46	<0.001
Desplazado	0.51	0.39, 0.65	<0.001	0.51	0.35, 0.71	<0.001	0.48	0.34, 0.66	<0.001
Privado de la libertad	0.32	0.27, 0.37	<0.001	0.32	0.25, 0.39	<0.001	0.31	0.26, 0.38	<0.001
Gestante	0.52	0.46, 0.58	<0.001	0.52	0.44, 0.60	<0.001	0.52	0.46, 0.59	<0.001
Población infantil ICBF	1.64	1.10, 2.34	0.010	1.64	0.93, 2.65	0.062	2.60	1.33, 5.16	<0.001
Institución psiquiatrica	0.54	0.39, 0.71	<0.001	0.54	0.35, 0.79	0.003	0.51	0.35, 0.75	<0.001
Victima de violencia	0.39	0.29, 0.52	<0.001	0.39	0.25, 0.57	<0.001	0.37	0.25, 0.53	<0.001

SIVIGILA JHRP

Tabla 4

Determinante	POISSON_DF			QUASI-POISSON_DF			BNM_NGT_DF
	IRR	95% CI	p-value	IRR	95% CI	p-value	IRR	95% CI	p-value
Sexo	1.79	1.72, 1.86	<0.001	1.79	1.68, 1.90	<0.001	1.78	1.69, 1.87	<0.001
Heterosexual	2.92	2.66, 3.22	<0.001	2.92	2.53, 3.40	<0.001	2.90	2.59, 3.25	<0.001
Otra orientación sexual	3.47	3.15, 3.82	<0.001	3.47	3.00, 4.03	<0.001	3.38	3.02, 3.80	<0.001

SIVIGILA JHRP

El modelo final fue: (Principio de parsimonia): Y = 1.8*varón + 2.9*prácticas heterosexuales + 3.4*prácticas No heterosexuales

prácticas heterosexuales inseguras: RR = 2.90, IC 95%: 2.59-3.25

después de controlar para el efecto confusor de los otros determinantes explicatorios

El instrumento de recolección de la información varió su estructura desde 2015, no fue diseñado por los investigadores del estudio, el objetivo de su creación fue la vigilancia epidemiológica y no incluye hábitos de protección de barrera y farmacológicos, entre otros factores de riesgo. La recolección de los datos por parte de personal médico es acotada en el tiempo para una anamnesis amplia y empática. Por parte del paciente posibles sesgos debido a subregistros de información. En 290 registros (1.9%), no se halló información sobre la posible vía de transmisión y en 354 registros (2.4%) sobre vulnerabilidad social del caso. Además, en la fase de preprocesamiento se halló un Subregistro del 35% de casos, del esperado total de casos equivalente a la suma de 9,072 (subregistros) + 16,954 vigentes = 26,026 casos esperados. Los subregistros de la base de datos de mortalidad de VIH, fueron 1,431 registros (57.2%) y de la base de datos de IPS complementarias especializadas en la atención de casos VIH, 7,641 (51.8%).

Otro posible sesgo estaría relacionado con la mala clasificación de los MoT relacionados con la privacidad, confidencialidad, estigma y discriminación y la fase psicológica de negación. Un estudio en Ucrania ³⁰ entre 2014 y 2018 estimó una mala clasificación del MoT de personas VIH de hasta el 38% al momento del diagnóstico, que fue reduciéndose a medida que la cohorte avanzó en el programa de atención y posteriores visitas de seguimiento. Otro estudio en Alemania ³¹ en una comunidad de turcos migrantes HIV positivos identificó mediante técnicas genómicas y estadísticas que hasta un 59% de casos HIV positivos reportados como heterosexuales tenían genomas del VIH relacionados con genomas de virus de casos VIH positivos no heterosexuales.

Otro posible sesgo, estaría relacionado con una proporción de tamizajes predominante en determinados grupos de población de riesgo.

Se pueden relacionar las siguientes: contar con un instrumento de recolección de datos único nacional y oficial, talento humano médico para el diligenciamiento, digitación y transferencia electrónica en software único nacional, aseguran uniformidad, experticia y aseguramiento de calidad. Otra fortaleza del estudio es el uso de técnicas estadísticas de análisis multivariado que controlan posibles efectos confusores de otras variables. Se evaluaron supuestos como No equidispersión. Al ser detectada se corrieron modelos de regresión quasi-poisson y binomial negativa. Aunque el modelo de regresión binomial negativa arrojo el mejor ajuste según el criterio de información de Akaike, los estimados de asociación no variaron notablemente y siguiendo el principio de parsimonia; acogemos el modelo de Poisson inicial.

Los resultados de este estudio alertan sobre falsas percepciones de asociación, entre el aumento de la incidencia y algunos grupos vulnerables como por ejemplo los migrantes, aumentando erróneamente la discriminación y estigma. La mala clasificación del MoT, (38%-59%) se relacionan con largos periodos de transmisión de la infección de los contactos no tamizados, diagnósticos de infección en estadios avanzados de enfermedad de los contactos. La cuenta de alto costo (CAC) reporto que hasta el 40.5% de casos nuevos de infección por VIH fueron diagnosticados en estadio 2 y 34.7 en estadio 3 o SIDA, lo cual estaría relacionado con el fenómeno descrito. Aumento de la incidencia no detectada y los costos de tratamiento por infecciones oportunistas. Disminución de la calidad y su expectativa de vida. La magnitud del efecto sobre el conteo de CRS de los MoT sexual es FUERTE y pone de relieve la prioridad de las estrategias de control.

Determinación de factores administrativos y psico-sociales durante el confinamiento que contribuyeron al aumento de la incidencia (procesos de tamizaje, confirmación, acceso, seguimiento, adherencia). Cuantificar la mala clasificación de MoT mediante técnicas genómicas

Reducir el subregistro mediante el análisis periódico de reporte de casos de VIH por un EQUIPO de SIVIGILA, RUAF, y de clínicas de tratamiento. Entrenamiento en recolección de información.