ISSN2953-6367 Enero2026 http://revistainvestigo.com Vol.7No,18,PP.27-52 https://doi.org/10.56519/4ecdck16 RevistaCientíficaMultidisciplinariaInvestiGoRiobamba–EcuadorCel:+593979119620revisinvestigo@gmail.com 27 ANÁLISISESPACIALDELRIESGOSANITARIOPORCONCENTRACIÓNDEGRANJASAVICOLASSINBIOSEGURIDADENLAPROVINCIADETUNGURAHUA:UNAPROPUESTADEGESTIÓNPREVENTIVAPARALASOSTENIBILIDADEMPRESARIALSPATIALANALYSISOFTHEHEALTHRISKDUETOTHECONCENTRATIONOFPOULTRYFARMSWITHOUTBIOSECURITYINTHEPROVINCEOFTUNGURAHUA:APREVENTIVEMANAGEMENTPROPOSALFORBUSINESSSUSTAINABILITYSayakDamiánRamosCachiguango1,JulioMauricioVizuete-Muñoz2,MaríaDoloresGuamánGuevara3{sramos2562@uta.edu.ec1,jm.vizuete@uta.edu.ec2,md.guaman@uta.edu.ec3}Fechaderecepción:19/11/2025/Fechadeaceptación:11/12/2025/Fechadepublicación:06/01/2025RESUMEN:EntreLaausenciadeprácticasadecuadasdebioseguridadenlasgranjasavícolasdelaprovinciadeTungurahua,Ecuador,poneenriesgolasostenibilidadempresarial,lacompetitividadproductivaylasaludpúblicadeestesector.Anteestaproblemática,elpresenteestudiotienecomoobjetivoanalizarelriesgosanitarioasociadoalaconcentracióndegranjasavícolassinbioseguridad.Paraello,seaplicóunenfoqueespacialcuantitativoentresfases:enlaprimerafase,seefectuóunanálisisdeconglomeradosjerárquicosparaagruparparroquiasconsimilitudesproductivasysanitarias,enlasegundafase,secalculóelÍndicedeRiesgoSanitarioParroquial(IRSP)paramedirlaexposiciónasistemasdeficientesdedesinfección,controldeplagas,controldeanimalesybioseguridadyfinalmente,serealizóunaevaluacióndelaautocorrelaciónespacialmedianteelíndicedeMoranGlobalyLocal(LISA)con999permutacionesyunamatrizdepesosespacialesdecincovecinosmáscercanos.Losresultadosmostraronunaautocorrelaciónespacialpositivasignificativa(p<=0.05)identificandoclústeresdealtoriesgo.Estoshallazgosevidencianladesigualdadterritorialenlagestióndelriesgoysubrayanlanecesidaddeestrategiaspreventivasquefortalezcanlasostenibilidadproductiva,laseguridadalimentariaylacompetitividadempresarialdelsector 1 CarreraAdministracióndeEmpresas,FacultaddeCienciasAdministrativas,UniversidadTécnicadeAmbato-Ecuador(UTA).https://orcid.org/0009-0006-0330-9005 2 DocenteInvestigadorFacultaddeCienciasAdministrativas,UniversidadTécnicadeAmbato,Ecuador.https://orcid.org/0000-0003-3910-9292 3 DocenteInvestigadorFacultaddeCienciasAdministrativas,UniversidadTécnicadeAmbato,Ecuador.https://orcid.org/0000-0003-4771-6412
ANÁLISISESPACIALDELRIESGOSANITARIOPORCONCENTRACIÓNDEGRANJASAVICOLASSINBIOSEGURIDADENLAPROVINCIADETUNGURAHUA:UNAPROPUESTADEGESTIÓNPREVENTIVAPARALASOSTENIBILIDADEMPRESARIAL 28 avícolaecuatoriano.Enconclusión,elestudiodemuestralautilidaddelanálisisespacialparapriorizarintervencionesenlaszonascríticasidentificadas,loquepermitiríaoptimizarrecursosydiseñaraccionescorrectivasespecíficasparamitigarelriesgosanitarioyprotegerlaaviculturalocal.Palabrasclave:Análisisespacial,avicultura,bioseguridad,riesgosanitario,sostenibilidadempresarialABSTRACT:ThelackofproperbiosecuritypracticesinpoultryfarmsintheprovinceofTungurahua,Ecuador,jeopardizesbusinesssustainability,productivecompetitiveness,andpublichealthinthissector.Facedwiththisproblem,thepresentstudyaimstoanalyzethehealthriskassociatedwiththeconcentrationofpoultryfarmswithoutbiosecurity.Tothisend,aquantitativespatialapproachwasappliedinthreephases:thefirstphase,ahierarchicalclusteranalysiswasconductedtogroupparisheswithsimilarproductionandhealthcharacteristics,thesecondphase,theParishHealthRiskIndex(PHRI)wascalculatedtomeasureexposuretoinadequatedisinfectionsystems,pestcontrol,animalcontrol,andbiosecurity,andfinallyaspatialautocorrelationassessmentwasconductedusingtheGlobalandLocalMoran’sI(LISA)indexwith999permutationsandaspatialweightsmatrixofthefivenearestneighbors.Theresultsshowedasignificantpositivespatialautocorrelation(p<=0.05),identifyinghigh-riskclusters.Thesefindingshighlightterritorialinequalityinriskmanagementandunderscoretheneedforpreventivestrategiesthatstrengthenproductivesustainability,foodsecurity,andthebusinesscompetitivenessoftheEcuadorianpoultrysector.Inconclusion,thestudydemonstratestheusefulnessofspatialanalysisforprioritizinginterventionsinidentifiedcriticalareas,whichwouldoptimizeresourcesandenablethedesignofspecificcorrectiveactionstomitigatehealthrisksandprotectlocalpoultryfarming.Keywords:Biosecurity,businesssustainability,healthrisk,poultryfarming,spatialanalysisINTRODUCCIÓNLaalimentacióneslabasefundamentalparaqueelserhumanopuedaadquirirenergíaynutrientesnecesariosparaelmantenimientodelavidaylasalud.Elconsumobalanceadodealimentosricosenvitaminas,minerales,proteínasyentreotros,ayudaaqueelserhumanoseapropensoaenfermedadescrónicascomoporejemploladiabetes,laobesidad,laanemia,elcánceryotros(1).Enesteescenario,elconsumoadecuadodeproteínasesdesumaimportanciayaquebrindaalcuerpomacronutrientesquesonesencialesenlaformacióndetejidos,laregulacióndefuncionesfisiológicasyelfortalecimientodelsistemainmunológico(2).Alelegirunafuentedeproteínamásadecuadaparaelfortalecimientoydesarrollodelcuerpohumano,losexpertosmencionanquelasproteínasdeorigenanimalbrindanmayoresbeneficiosadiferenciadelasproteínasdeorigenvegetalyaquecuentanconunaciertacantidaddeaminoácidos(3).
ANÁLISISESPACIALDELRIESGOSANITARIOPORCONCENTRACIÓNDEGRANJASAVICOLASSINBIOSEGURIDADENLAPROVINCIADETUNGURAHUA:UNAPROPUESTADEGESTIÓNPREVENTIVAPARALASOSTENIBILIDADEMPRESARIAL 29 Enlasúltimasdécadaselsectordelaaviculturasehaconsolidadocomounodelosprincipalessectoresqueofertalaproteínadeorigenanimal,suprincipalproductoeslacarnedepolloelcualescomercializadaanivelmundial,actualmenteesteproductorepresentael40%delconsumomundialdecarne,posicionándose,asícomolaproteínapreferidaporlosconsumidoresgraciasasuvalornutricional,asequibilidadyversatilidadculinaria(4).Asimismo,laproducciónmundialdehuevosalcanzólos91millonesdetoneladasenelaño2023,loqueequivaleamásde1,7billonesdeunidades(huevos)producidas(5).EnAméricaLatina,elsectoravícolaesconsideradocomounodelossectoresconmayordesarrolloeconómicoysocialaportandosignificativamentealaindustriaalimentariayademásagenerandomillonesdeempleosdirectoeindirectos(6).EnelEcuador,elconsumopercápitadecarnedepolloasciendea34kgporpersonaalaño,mientrasquelaproducciónnacionaldehuevossuperalos3.800millonesdeunidadesanuales,cifraqueevidenciaelrolestratégicodelsectorenlaseguridadalimentariayenlacompetitividadagroindustrial(7).Porotraparte,elcrecimientoconstantedeestesectorysugranaportealaindustriaalimentariaestáligadaintrínsecamenteconelriesgosanitario,comoejemplocitaremoslaamenazadeenfermedadescomolainfluenciaaviar,lasalmonella,lapsitacosisentreotros.Sehanregistradoalrededorde365brotesdevirusen24paíseseuropeosocasionandoelsacrificioolamuertedemásde7,9millonesdeavesdomésticasysilvestresaestolesumamoslamuestrede4sereshumanoporelcontagiodeestasenfermedadestodoestoenlosmesesdemarzoajuniodel2025(8).Enelcontinenteamericanodurantelosañosde2022hastaprincipiosdel2025sehanregistradomásde4.700brotesdevirusenavesyalrededorde74infeccioneshumanasasociadasalvirusH5N1(9).Perodesdeelaño2003hastaelaño2025laOrganizaciónMundialdelaSaludyaharegistradomásde986casosdelvirusH5N1enlossereshumanos,deestetotalel48%yahanfallecido(10).DeacuerdoconAgrocalidad(11),existenalrededorde1900granjasavícolasquesededicanalacrianzaycomercializacióndeaves(gallinas,pavos,patos,entreotros).Sinembargo,notodaslasgranjasavícolascumplenconlosrequisitossanitariosbásicos,alrededorde1000granjasavícolascuentanconcertificaciónzoosanitariaydeestetotalapenas120cumplenconlasdistanciasmínimasdeaislamientoquefueronreguladasbajolaresolución0286emitidaporlaAgenciadeRegulaciónyControlFitoyZoosanitario,esunasituaciónpreocupanteyaqueel50%delasgranjasavícolasincumplenestándaresbásicosdebioseguridadycercadel90%nocumpleconlosaislamientosrecomendados.Esteescenariorepresentaunaltoriesgozoosanitariodondelafaltadeprevención,controlybuenasprácticasdejanexpuestoalsectorabrotesdeenfermedadesdegranimpacto.Anivelnacional,durantelosañosde2022y2023sehanregistrado27brotesdeinfluenciaaviarafectandoasíamásde1,25millonesdeavesyobligandoalsacrificiodecientosdemilesdeellos(12).Estasmedidastomadasporlospropietariosdelasgranjasrepresentanunagranpérdidatantoeconómicacomoemocionalyaqueeltiempoinvertidoenlacrianzadeestosanimalessefueronalabasura.
ANÁLISISESPACIALDELRIESGOSANITARIOPORCONCENTRACIÓNDEGRANJASAVICOLASSINBIOSEGURIDADENLAPROVINCIADETUNGURAHUA:UNAPROPUESTADEGESTIÓNPREVENTIVAPARALASOSTENIBILIDADEMPRESARIAL 30 LaregióncentraldelEcuadorparticularmentelasprovinciasdeCotopaxi,ChimborazoyTungurahuasonpuntosneurálgicosparaeldesarrolloeconómicoalconcentrarel85%delaproducciónavícola(13).PorloquelaprovinciadeTungurahuaseconsolidacomounadelascuencasavícolasmássignificativasenlasierraecuatoriana.Paraelaño2025,seproyectaquelapoblaciónavícolaenlaprovinciaalcanceaproximadamentelos5.8millonesdeavesloquerepresentacercadel15%delinventarionacionalavícolaestimadoen38millones(14).Enlaprovincia,laamenazadelainfluenciaaviaraltamentepatógena(IAAP)representaunriesgosanitariosignificativoparaelsectoravícola,debidoasualtadensidaddeexplotaciónavícolayalarecienteexperienciadebrotes.Enenerode2023,elMinisteriodeAgriculturayGanadería(MAG)reportólaaparicióndedosnuevosbrotesenlasgranjasdelaprovincia,loqueelevóasieteeltotaldebrotesconfirmadosenesemomentoanivelnacional.Estasituacióndemuestralavulnerabilidadespecíficadelaprovinciafrentealapropagacióndelviruscuandolasmedidasdeprevenciónnosonadecuadas(15).Unasuposición,silaprovinciaseveafectadasignificativamenteporestevirus,lasgranjasavícolasseveríanmuyperjudicadas,porejemplo,unagranjaquecontieneaproximadamente50.000gallinasponedorasgeneraríapérdidasdirectasqueascenderíanlos$385.000.Elobjetivodeestainvestigaciónfueanalizarlospatronesterritorialesderiesgosanitarioenelsectordelaaviculturadelaprovincia.Paraello,seplanteóunametodologíaqueestádivididaentresetapas:Enlaprimeraetapa,seclasificóalasparroquiassegúnsussimilitudesproductivasysanitariasmedianteelanálisisdeclúster.Enlasegundaetapa,secalculóelÍndicedeRiesgoSanitarioParroquial(IRSP)conlafinalidaddeidentificarnivelesderiesgoterritorial.Yporúltimo,seevaluólaautocorrelaciónespacialdelriesgomedianteelÍndicedeMoranGlobalyLocal(LISA).MATERIALESYMÉTODOSLainvestigaciónsedesarrollóbajounenfoquecuantitativo,descriptivoyespacialconundiseñonoexperimentalytransversal.Porello,esteestudiosecentróenanalizarlospatronesterritorialesdeconcentracióndelsectoravícoladelaprovinciadeTungurahua,Ecuador.Seempleóunenfoquecuantitativobasadoentresetapas.Enlaprimeraetapa,seaplicaunanálisisdeagrupamientonoespacial(clústerjerárquico)medianteelcálculodelmétododedistanciaeuclidianacuadradaentrelasparroquiasqueseencuentranenlaprovinciadeTungurahua,estemétodonospermitirádeterminarlasimilituddedosomásvariablesenunespacion-dimensional(16).Paraello,seutilizaronvariablesquereflejantantolacapacidadproductivacomolascondicionesdebioseguridaddecadaunadelasgranjasavícolas.Paradeterminaraquellasimilitudentrelasparroquiassehaconsideradoutilizarlassiguientesvariables:1)NúmerodeGalpones;2)CapacidaddeAvesInstalada;3)CapacidaddeAvesOcupada;4)SistemadeDesinfecciónparaPersonasqueIngresanalasGranjas;5)ControldeInsectos;6)ControldeRoedores;7)ControldeIngresodePerrosy/oGatosalosGalponesyÁreasdeCrianza;y8)ControldeBioseguridad.Pararealizardichoanálisisseaplicólasiguientefórmula:
ANÁLISISESPACIALDELRIESGOSANITARIOPORCONCENTRACIÓNDEGRANJASAVICOLASSINBIOSEGURIDADENLAPROVINCIADETUNGURAHUA:UNAPROPUESTADEGESTIÓNPREVENTIVAPARALASOSTENIBILIDADEMPRESARIAL 31 Posteriormenteseelaboróundendrogramaelcualserepresentacomoundiagramasemejanteaunárbolinvertidoquiennosfacilitalainterpretacióndelasrelacionesqueexistenentrelasvariables,cuandolaunióndedosramasseaenelnivelmásbajoestosignificaráunamayorsemejanzaentrelasvariablesoelementosconectados,perosilauniónocurreenelnivelmásaltosignificaráquelasvariablessonmuydistantesomenossimilares(17).Enlasegundaetapa,secalculóelÍndicedeRiesgoSanitarioParroquial(IRSP)conelobjetivodedimensionarelpesorelativodelasgranjasavícolasenriesgo(Sincondicionesadecuadasdebioseguridad)dentrodecadaparroquia,encomparaciónconelpromedioprovincial.Esteíndiceseelaboróapartirdelaadaptacióndelafórmuladelcocientedelocalización(LQ),demaneraqueparacadaparroquiadelaprovinciadeTungurahuasecalculólaproporcióndegranjasavícolasenriesgorespectoasunúmerototaldegranjas,tantoanivelparroquialcomoprovincial.Paraprofundizarelanálisis,seutilizóelsoftwareestadísticoSPSS,atravésdelcualserealizaroncontrastesdeproporcionesconelfindecomprobarsilasvariablesobservadasentrecadaparroquiayelpromedioprovincialpresentandiferenciasestadísticamentesignificativas.Graciasaesteprocedimiento,elíndiceelaboradonosolofacilitaunacomparaciónordenadaentrelasparroquias,sinoqueademáspermitereconocerconmayorclaridadcuálesdeellasconcentrannivelesmáselevadosderiesgosanitario. Donde:IRSP:ÍndicedeRiesgoSanitarioParroquial.GRp:Númerodegranjasconriesgosanitarioenlaparroquiap.GTp:Totaldegranjasenlaparroquiap.GRt:Númerototaldegranjasconriesgosanitarioenlaprovincia.GTt:Totaldegranjasenlaprovincia.Finalmente,seefectuóunAnálisisExploratoriodeDatosEspacialesconelobjetivodecomprenderladistribuciónterritorialdelriesgosanitarioasociadoalaproducciónavícolaenlaprovinciadeTungurahua.Paraelloseutilizóinformacióngeorreferenciadadelasgranjas,loquepermitióasignarunvalorderiesgoacadadivisiónadministrativa.Conestosinsumosfueposible
ANÁLISISESPACIALDELRIESGOSANITARIOPORCONCENTRACIÓNDEGRANJASAVICOLASSINBIOSEGURIDADENLAPROVINCIADETUNGURAHUA:UNAPROPUESTADEGESTIÓNPREVENTIVAPARALASOSTENIBILIDADEMPRESARIAL 32 construirunamatrizdepesosespacialesfundamentalesenlacontigüidadentreparroquias,empleandoparaelloelsoftwareGeoDa(18).Posteriormente,seestimóelíndicedeMoranGlobal,herramientaquepermiteevaluarelgradodeautocorrelaciónespacialpresenteenunaovariasvariablesdentrodeláreadeestudio(19).ElprocedimientocontinuóconelcálculodelestadísticodeMoranLocal(LISA),técnicaorientadaaidentificarpatronesdedependenciaespacialconunniveldedetallemásfino(20).Graciasaesteanálisis,selograronreconocersectoresespecíficosdondeelriesgosanitariotiendeaconcentrarse(alto-alto)yzonasconmenorincidencia(bajo-bajo).LainformaciónseobtuvodelCatastroAvícoladelInstitutoNacionaldeEstadísticaycensos(INEC)atravésdelsistemaESPAC(EncuestadeSuperficieyProducciónAgropecuariaContinua)correspondientealúltimoperiododisponible.RESULTADOSLosresultadosobtenidosenelpresenteestudiopermitenevidenciarlamaneraenqueelriesgosanitarioasociadoalaproducciónavícolasedistribuyedeformadiferenciadaenelterritoriodelaprovinciadeTungurahua.Enestasecciónsedetallanlosresultadosdeformasecuencial,siguiendolastresetapasmetodológicasdescritaspreviamente.AgrupamientonoespacialElanálisisdeconglomeradosmedianteladistanciaeuclidianacuadradapermitióidentificaryclasificarentresgruposoclústereslasparroquiasdelaprovinciadeTungurahuaquepresentanentresípatronesdesimilitudenrelaciónconelriesgosanitarioavícolaylacapacidadproductiva.Acontinuación,semuestraenlaFigura1,eldendrogramaconstruidoapartirdelaestandarizacióndelasvariablesdeestudioyelusodelmétododeenlaceentregrupos. G1 G2 G3
ANÁLISISESPACIALDELRIESGOSANITARIOPORCONCENTRACIÓNDEGRANJASAVICOLASSINBIOSEGURIDADENLAPROVINCIADETUNGURAHUA:UNAPROPUESTADEGESTIÓNPREVENTIVAPARALASOSTENIBILIDADEMPRESARIAL 33 Figura1.DendrogramadeAgrupamientodeParroquiasdelaProvinciadeTungurahua.Nota:Eldendrogramaresultantemuestralaagrupacióndelasparroquiasenfuncióndesusimilitud.CreadoporelautorutilizandoIBMSPSSStatistics.ElGrupo1estáconformadoporlasparroquiasquecuentanconunaaltahomogeneidadenlasbuenasprácticassanitarias,aquíseencuentranagrupadaslasparroquiasdeAmbato,Atahualpa,AugustoN.Martínez,Baños,Bolívar,Cevallos,ConstantinoFernández,Cotalo,Cunchibamba,ElTriunfo,GarcíaMoreno,Guambalo,HuachiGrande,JuanBenignoVela,Mocha,Montalvo,Patate,Pelileo,Pinguili,Quero,RioNegro,Rumipamba,SanBartolomédePinllo,SanMiguelito,SantaRosa,Sucre,Tisaleo,UlbayUnamunchoquienespresentanunnivelderiesgobajo,esdecir,estegrupocuentaconporcentajesreducidosdeincumplimientoporejemplo:enelcontroldedesinfecciónpresentóunamediade16.38%estoquieredecirqueel83.62%delasgranjasavícolasqueseencuentranubicadasenestasparroquiascuentanconuncontroldedesinfecciónrigurosoparalaspersonasqueingresanalosdiferentesgalpones,elcontroldeinsectosconel16,41%,elcontrolderoedoresconel18,62%yelcontroldebioseguridadconel27.79%perohayunaexcepciónenelcasodelcontroldeperrosygatosyaquepresentóunamediade77.45%estoquieredecirqueapenasel22.55%delasgranjasavícolascuentanconuncontroldelingresodelosperrosygatosalosgalpones,estopuedepresentarsecomoriesgoyaqueestosanimalitospuedenserconductoresdevariasenfermedadesyasuvezserdepredadoresquepuedanperjudicaralaproducciónycrecimientodelasgranjasavícolas.ElGrupo2estacaracterizadoporserelgrupoconmayorriesgosanitario,aquíseconcentranlassiguientesparroquiasBenítez,Izamba,MarcosEspinel,Picaihua,Píllaro,PresidenteUrbinayTotoras,todaslasvariables(Controldedesinfección81.86%,Controldeinsectos89.71%,Controlderoedores91.29%,Controldeperrosygatos70,71%yControldebioseguridad92,14%)superanel70%deincumplimientoloquesignificaqueapenasun30%omenosdelasmedidasestánsiendoimplementadascorrectamenteenlasgranjasavícolas.Estepatróndeincumplimientouniformeyelevadodenotaunaltoriesgosanitarioyproductivo,yaquetantolacapacidadinstaladacomolosmecanismosdecontrolpresentanfallassistemáticas.EnlaIlustración2sepuedeobservarquelasparroquiaspertenecientesaestegrupoconformanunbloquecompactoydiferenciadodelresto,evidenciandoquecompartenunaproblemáticaestructuralcomúnquedemandaintervencionesintegralesyurgentes.ElGrupo3conformadaporlasparroquiasLosAndesySanAndréspresentanunescenariomuyfavorableentérminosdebioseguridadycapacidadproductivayaquelamayoríadelasvariablesevaluadasmuestranunamediadel0%estoqueriendodecirquelasgranjasavícolasqueseencuentranenestasparroquiascumplenal100%conelcontroldedesinfección,elcontroldeinsectos,elcontrolderoedoresyconelcontroldebioseguridad,sedebepriorizaraccionesparacontrolarelingresodeperrosygatosalosgalponesyaqueeslavariableconstituyelaprincipaldebilidaddelgrupo.ÍndicedeRiesgoSanitarioParroquial(IRSP)
ANÁLISISESPACIALDELRIESGOSANITARIOPORCONCENTRACIÓNDEGRANJASAVICOLASSINBIOSEGURIDADENLAPROVINCIADETUNGURAHUA:UNAPROPUESTADEGESTIÓNPREVENTIVAPARALASOSTENIBILIDADEMPRESARIAL 34 ElanálisisdelÍndicedeRiesgoSanitarioParroquial(IRSP)identificólosnivelesdevulnerabilidadsanitariadelsectoravícolaenlaprovinciadeTungurahua.Laaltaconcentracióndeproducciónenestaregión,querepresentael15,9%deltotalnacional(11),seveexpuestaporriesgossanitariosyeconómicosasociadosalapropagacióndeenfermedadesinfecciosasenzonasdealtadensidadproductiva.Losresultadosobtenidossecalcularonparaslascincovariablesquesondesinfección,controldeinsectos,controlderoedores,controldeperrosy/ogatosybioseguridad.Estosvaloresobtenidosrevelancontrastessignificativosentrelasparroquias,tantoenlamagnituddelriesgocomoenladistribuciónespacialdelasdeficienciasdebioseguridad.DeacuerdoconlosvaloresdelIRSP,seidentificarontrescategoríasprincipalesderiesgo:Riesgoalto(IRSP>1)dondelaproporcióndegranjassinmedidasdebioseguridadsuperanelpromedioprovincial.Riesgomedio(0.70≤IRSP≤1)conunniveldecumplimientosanitariomoderado.Riesgobajo(IRSP<0.70)caracterizadoporlaadopcióndeprácticasmásconsolidadasdebioseguridadycontrolsanitario.LaTabla1presentalosresultadosdesagregadosporparroquiarevelandounaclaraheterogeneidadenlaadopcióndemedidaspreventivas.Tabla1.ÍndicedeRiesgoSanitarioParroquial. Parroquia IRSPDes. Cat.delRies. IRSPIns. Cat.delRies. IRSPRoe. Cat.delRies. IRSPP&G Cat.delRies. IRSPBio. Cat.delRies. Ambatillo 0.00 Bajo 0.00 Bajo 0.00 Bajo 0.00 Bajo 0.00 Bajo Ambato 0.65 Bajo 0.62 Bajo 0.59 Bajo 1.01 Alto 1.35 Alto Atahualpa 1.72 Alto 1.64 Alto 1.58 Alto 0.90 Medio 2.00 Alto AugustoN.Martínez 0.94 Medio 1.12 Alto 1.29 Alto 1.23 Alto 1.47 Alto Baños 1.72 Alto 2.46 Alto 2.37 Alto 1.13 Alto 1.20 Alto BaquerizoMoreno 0.00 Bajo 0.00 Bajo 0.00 Bajo 0.00 Bajo 0.00 Bajo Benítez 3.44 Alto 4.92 Alto 4.75 Alto 0.90 Medio 2.40 Alto Bolívar 0.57 Bajo 1.09 Alto 1.05 Alto 1.35 Alto 1.20 Alto
ANÁLISISESPACIALDELRIESGOSANITARIOPORCONCENTRACIÓNDEGRANJASAVICOLASSINBIOSEGURIDADENLAPROVINCIADETUNGURAHUA:UNAPROPUESTADEGESTIÓNPREVENTIVAPARALASOSTENIBILIDADEMPRESARIAL 35 Cevallos 1.29 Alto 0.62 Bajo 0.59 Bajo 1.01 Alto 1.35 Alto Chiquicha 0.00 Bajo 0.00 Bajo 0.00 Bajo 0.00 Bajo 0.00 Bajo ConstantinoFernández 1.55 Alto 0.98 Medio 0.95 Medio 1.08 Alto 1.08 Alto Cotalo 0.99 Medio 0.84 Medio 1.11 Alto 0.86 Medio 1.07 Alto Cuchibamba 0.00 Bajo 1.23 Alto 1.19 Alto 1.01 Alto 0.90 Medio ElRosario 0.00 Bajo 0.00 Bajo 0.00 Bajo 0.00 Bajo 0.00 Bajo ElTriunfo 0.00 Bajo 0.00 Bajo 0.00 Bajo 1.35 Alto 0.00 Bajo EmilioMaríaTerán 0.00 Bajo 0.00 Bajo 0.00 Bajo 0.00 Bajo 0.00 Bajo GarcíaMoreno 0.00 Bajo 4.92 Alto 4.75 Alto 1.35 Alto 0.00 Bajo Guambalo 0.76 Medio 1.08 Alto 0.81 Medio 0.99 Medio 1.05 Alto HuachiGrande 0.00 Bajo 0.00 Bajo 0.00 Bajo 1.35 Alto 0.00 Bajo Izamba 5.17 Alto 4.92 Alto 4.75 Alto 0.00 Bajo 3.59 Alto JuanBenignoVela 1.48 Alto 1.41 Alto 1.36 Alto 0.58 Bajo 2.57 Alto Lligua 0.00 Bajo 0.00 Bajo 0.00 Bajo 0.00 Bajo 0.00 Bajo LosAndes 0.00 Bajo 0.00 Bajo 0.00 Bajo 1.35 Alto 0.00 Bajo MarcosEspinel 5.17 Alto 4.92 Alto 4.75 Alto 1.35 Alto 3.59 Alto Mocha 1.72 Alto 0.55 Bajo 1.05 Alto 1.05 Alto 0.80 Medio Montalvo 0.00 Bajo 0.98 Medio 0.95 Medio 0.27 Bajo 1.44 Alto Pasa 0.00 Bajo 0.00 Bajo 0.00 Bajo 0.00 Bajo 0.00 Bajo Patate 0.94 Medio 0.00 Bajo 0.00 Bajo 0.61 Bajo 0.33 Bajo Pelileo 1.05 Alto 1.00 Alto 0.87 Medio 1.10 Alto 0.66 Bajo Picaihua 5.17 Alto 4.92 Alto 4.75 Alto 1.35 Alto 3.59 Alto Pilaguin 0.00 Bajo 0.00 Bajo 0.00 Bajo 0.00 Bajo 0.00 Bajo Píllaro 2.87 Alto 3.83 Alto 4.22 Alto 1.05 Alto 2.80 Alto Pinguili 0.00 Bajo 0.00 Bajo 1.19 Alto 1.35 Alto 0.00 Bajo
ANÁLISISESPACIALDELRIESGOSANITARIOPORCONCENTRACIÓNDEGRANJASAVICOLASSINBIOSEGURIDADENLAPROVINCIADETUNGURAHUA:UNAPROPUESTADEGESTIÓNPREVENTIVAPARALASOSTENIBILIDADEMPRESARIAL 36 PresidenteUribina 5.17 Alto 2.46 Alto 2.37 Alto 0.68 Bajo 3.59 Alto Quero 2.58 Alto 0.00 Bajo 1.19 Alto 1.01 Alto 2.70 Alto Quinchicoto 0.00 Bajo 0.00 Bajo 0.00 Bajo 0.00 Bajo 0.00 Bajo Quisapincha 0.00 Bajo 0.00 Bajo 0.00 Bajo 0.00 Bajo 0.00 Bajo RioNegro 1.29 Alto 1.23 Alto 1.19 Alto 1.35 Alto 0.90 Medio RioVerde 0.00 Bajo 0.00 Bajo 0.00 Bajo 0.00 Bajo 0.00 Bajo Rumipamba 0.00 Bajo 0.00 Bajo 0.00 Bajo 1.35 Alto 0.00 Bajo Salasaca 0.00 Bajo 0.00 Bajo 0.00 Bajo 0.00 Bajo 0.00 Bajo SanAndrés 0.00 Bajo 0.00 Bajo 0.00 Bajo 1.35 Alto 0.00 Bajo SanBartolomédePinllo 2.58 Alto 0.00 Bajo 0.00 Bajo 0.00 Bajo 3.59 Alto SanFernando 0.00 Bajo 0.00 Bajo 0.00 Bajo 0.00 Bajo 0.00 Bajo SanJosédePoalo 0.00 Bajo 0.00 Bajo 0.00 Bajo 0.00 Bajo 0.00 Bajo SanMiguelito 0.00 Bajo 0.00 Bajo 0.00 Bajo 1.35 Alto 0.00 Bajo SantaRosa 1.03 Alto 0.00 Bajo 0.00 Bajo 0.27 Bajo 2.16 Alto Sucre 0.00 Bajo 0.00 Bajo 0.00 Bajo 1.35 Alto 0.00 Bajo Tisaleo 1.72 Alto 1.64 Alto 1.58 Alto 1.35 Alto 1.20 Alto Totoras 2.58 Alto 4.92 Alto 4.75 Alto 1.35 Alto 3.59 Alto Ulba 0.00 Bajo 0.00 Bajo 0.00 Bajo 1.35 Alto 0.00 Bajo Unamuncho 0.00 Bajo 0.00 Bajo 0.00 Bajo 1.35 Alto 0.00 Bajo YanayacuMochapata 0.00 Bajo 0.00 Bajo 0.00 Bajo 0.00 Bajo 0.00 Bajo Nota:LosdatosfueronanalizadosporelautorutilizandoelsoftwareIBMSPSSStatistics.Losresultadosobtenidosreflejanlosiguiente,deacuerdoconlavariablededesinfecciónel35.85%delasparroquiascuentanconunnivelderiesgoaltoalnocontarconunsistemaeficazdedesinfecciónelcualpuedeincrementarlosriesgosdebrotessanitarios,porotraparte,aproximadamenteel33.96%delasparroquiasnohanimplementadoprogramasorecursosparamitigarocontrolarlasplagasdeinsectosocasionandoasíunentornollenodeenfermedadesyunriesgoalaproductividad.ParroquiascomoIzamba,MarcosEspinel,Picaihua,Benítez,Totoras,
ANÁLISISESPACIALDELRIESGOSANITARIOPORCONCENTRACIÓNDEGRANJASAVICOLASSINBIOSEGURIDADENLAPROVINCIADETUNGURAHUA:UNAPROPUESTADEGESTIÓNPREVENTIVAPARALASOSTENIBILIDADEMPRESARIAL 37 GarcíaMoreno,entreotrosrepresentanel37.74%deparroquiasquetienenunnivelderiesgoaltoenlavariabledecontrolderoedoresevidenciandoasíunavulnerabilidadmuyimportantealasaludpública.Asímismo,el50.94%delasparroquiasnotienenunrigurosocontroldelingresodelosperrosygatosalosgalponesloquecomprometealaintegridaddelasmedidassanitariasyfacilitaladiseminacióndeagentesinfecciosos.Finalmente,alrededordel41.51%delasparroquiascuentanconunnivelderiesgoaltoenelcontroldelabioseguridadloqueconfirmaunafaltadeprotocolossanitariosestandarizadosquelessirvaparacontrolar,supervisarymitigarlosbrotesepidémicosenlaproducciónavícola.LaTabla2presentaelIRSPGlobalobtenidocomopromediodelascincovariables.Esteíndiceintegradopermitióclasificaralasparroquiassegúnsuniveldevulnerabilidadsanitaria.Tabla2.ÍndicedeRiesgoSanitarioGlobal. Parroquia IRSPGlobal CategoríadelRiesgo Ambatillo 1.0 Bajo Ambato 1.8 Medio Atahualpa 2.8 Alto AugustoN.Martínez 2.8 Alto Baños 3.0 Alto BaquerizoMoreno 1.0 Bajo Benítez 2.8 Alto Bolívar 2.6 Alto Cevallos 2.2 Medio Chiquicha 1.0 Bajo ConstantinoFernández 2.6 Alto Cotalo 2.4 Alto Cuchibamba 2.4 Alto ElRosario 1.0 Bajo ElTriunfo 1.4 Medio EmilioMaríaTerán 1.0 Bajo GarcíaMoreno 2.2 Medio
ANÁLISISESPACIALDELRIESGOSANITARIOPORCONCENTRACIÓNDEGRANJASAVICOLASSINBIOSEGURIDADENLAPROVINCIADETUNGURAHUA:UNAPROPUESTADEGESTIÓNPREVENTIVAPARALASOSTENIBILIDADEMPRESARIAL 38 Guambalo 2.4 Alto HuachiGrande 1.4 Medio Izamba 2.6 Alto JuanBenignoVela 2.6 Alto Lligua 1.0 Bajo LosAndes 1.4 Medio MarcosEspinel 3.0 Alto Mocha 2.4 Alto Montalvo 1.8 Medio Pasa 1.0 Bajo Patate 1.2 Bajo Pelileo 2.4 Alto Picaihua 3.0 Alto Pilaguin 1.0 Bajo Píllaro 3.0 Alto Pinguili 1.8 Medio PresidenteUribina 2.6 Alto Quero 2.6 Alto Quinchicoto 1.0 Bajo Quisapincha 1.0 Bajo RioNegro 2.8 Alto RioVerde 1.0 Bajo Rumipamba 1.4 Medio Salasaca 1.0 Bajo SanAndrés 1.4 Medio SanBartolomédePinllo 1.8 Medio
ANÁLISISESPACIALDELRIESGOSANITARIOPORCONCENTRACIÓNDEGRANJASAVICOLASSINBIOSEGURIDADENLAPROVINCIADETUNGURAHUA:UNAPROPUESTADEGESTIÓNPREVENTIVAPARALASOSTENIBILIDADEMPRESARIAL 39 SanFernando 1.0 Bajo SanJosédePoalo 1.0 Bajo SanMiguelito 1.4 Medio SantaRosa 1.8 Medio Sucre 1.4 Medio Tisaleo 3.0 Alto Totoras 3.0 Alto Ulba 1.4 Medio Unamuncho 1.4 Medio YanayacuMochapata 1.0 Bajo Nota:Adaptadoderesultadosdesglosadosporvariable(Tabla2).Elaboradoporelautor.Losresultadosvisualizadosenlatabla2muestranquelasparroquiasdeBaños,MarcosEspinel,Picaihua,Píllaro,Tisaleo,entreotros,concentranel39.62%queregistranuníndicederiesgoalto.Estoseasociaconlaaltaconcentracióndegalpones,bajoniveldetecnificaciónylaescasaimplementacióndeprotocolosdebioseguridad,porotraparte,parroquiascomoAmbatillo,Pasa,Salasaca,ElRosario,entreotros,concentranel30.19%deparroquiasqueregistraníndicesbajosderiesgosanitarioloqueconfirmalaexistenciadegranjastecnificadasconestándaresdebioseguridadmásrobustas.Y,porúltimo,alrededordel30.19%deparroquiasmantienenunnivelderiesgomedio,sibienhanimplementadoprácticasdebioseguridadaúnpresentandeficienciasenaspectoscomoladesinfeccióndeinstalaciones,controldeinsectos,controlderoedores,entreotrosquegaranticenunaprotecciónsanitariaintegral.AnálisisespacialdelriesgosanitarioElanálisisdelaautocorrelaciónespacialparacadavariablederiesgosanitariopermitióidentificarlaexistenciadepatronesdeconcentraciónespacialdelosriesgosasociadosalasprácticasdebioseguridadenlasgranjasavícolasdelaprovinciadeTungurahua.ElíndicedeMoranGlobalpresentóvalorespositivos,estoimplicaquelasparroquiasconaltosvaloresderiesgotiendenaestarrodeadasdeparroquiasconcondicionessimilares,loquerefuerzalapresenciadezonascríticasquedemandanunagestiónpreventivamásfocalizada(21).Tabla3.ÍndicedeMoranGlobal. Variable ÍndicedeMoranGlobal(I) Valorp Interpretación
ANÁLISISESPACIALDELRIESGOSANITARIOPORCONCENTRACIÓNDEGRANJASAVICOLASSINBIOSEGURIDADENLAPROVINCIADETUNGURAHUA:UNAPROPUESTADEGESTIÓNPREVENTIVAPARALASOSTENIBILIDADEMPRESARIAL 40 Desinfección 0.069 0.002 Altaautocorrelaciónespacial ControldeInsectos 0.122 0.001 Altaautocorrelaciónespacial ControldeRoedores 0.118 0.001 Altaautocorrelaciónespacial ControldePerrosyGatos 0.154 0.001 Altaautocorrelaciónespacial ControldeBioseguridad 0.207 0.001 Altaautocorrelaciónespacial Nota.CálculorealizadoconelsoftwareGeoDa(18).Elaboradoporelautor.Losresultadosevidencianunpatróndeagrupamientoespacialpositivoencadaunadelasvariablesanalizadas.Estoindicaquelosnivelesderiesgosanitariopresentesenlasparroquiasnosiguenunadistribuciónaleatoria,porlocontrario,seidentificanáreasgeográficasdondelasdeficienciassanitariastiendenaconcentrarsedemaneraconsistente.Lavariablededesinfección,elvalorde0.069indicaquelasprácticasdelimpiezaydesinfecciónaúnnoseaplicandemaneraconstanteenlamayoríadelasgranjasloquereflejaunaejecuciónirregulardeestasmedidassanitarias.Estafaltadeestandarizaciónrepercuteenlaeficienciaoperativa,elevandoloscostosdemantenimientoygenerandoposiblespérdidaseconómicasderivadasdebrotessanitarios.Lasvariablesdeinsectosyroedores,losvaloresde0.112y0.118respectivamentereflejanqueestosriesgosestánfuertementeinterrelacionadosterritorialmente,loqueevidenciadeficienciasestructuralesenlagestióndelentornofísicodelasgranjasavícolas.Dichosfactorespuedenfacilitarlapropagacióndepatógenosycomprometerlasostenibilidadempresarialdelsectoravícola(22).Lavariabledecontroldeperrosygatospresentauníndicede0.154locualsugierequelapresenciadeanimalesdomésticosnocontroladostambiénsigueunpatróndeconcentraciónendeterminadasparroquias.Esteaspectopuedegenerarcontaminacióncruzadaydañosenlacadenaproductiva,afectandolatrazabilidadylacalidaddelproductofinal.Y,porúltimo,lavariabledebioseguridadmantieneuníndicede0.207loqueindicaquelasparroquiasmuestransimilitudeselevadasentresíencuantoalacarenciademedidasdecontrolinterno.Estopuedeincrementarelriesgodedifusióndeenfermedadesinfecciosasavícolas,afectandotantolaproductividadcomolarentabilidaddelosproductoreslocales.Elanálisislocaldeautocorrelaciónespacial(LISA)permitióidentificarlasparroquiasdelaprovinciadeTungurahuadondeelriesgosanitariopresentapatronesdeconcentraciónestadísticamentesignificativos,complementandolosresultadosdelíndicedeMoranGlobal.Atravésdeesteanálisis,sedeterminaronzonasderiesgoqueevidencianlaexistenciadeconglomeradosespaciales(High-HighyLow-Low)yáreasatípicas(Hig-LowyLow-High)lascualesaportanunavisióndetalladadeladistribucióndelriesgoavícola(23).
ANÁLISISESPACIALDELRIESGOSANITARIOPORCONCENTRACIÓNDEGRANJASAVICOLASSINBIOSEGURIDADENLAPROVINCIADETUNGURAHUA:UNAPROPUESTADEGESTIÓNPREVENTIVAPARALASOSTENIBILIDADEMPRESARIAL 41 Enesteanálisis,elcálculoserealizóutilizandoelsoftwareGeoDautilizandounamatrizbasadaenloscincovecinosmáscercanos(K=5)yuntotalde999permutacionesaleatorias,conelfindeasegurarunamayorsolidezestadística.Entodosloscasos,losvaloresdep≤0.05confirmaronunasignificanciaespacialelevada,loqueimplicaquelospatronesdeagrupamientoidentificadosnorespondenalazar,sinoqueaunaestructuraespacialreal.Figura2.ResumendeResultadosLISAparalavariabledeDesinfecciónNota:Análisisdeautocorrelaciónespaciallocalparalavariablededesinfección.CálculosrealizadosconelsoftwareGeoDa(18).Elaboradoporelautor.Figura3.MapadeSignificanciadelÍndicedeMoranLocalparalavariabledeDesinfección.Nota.MapadeclústeresestadísticamentesignificativosdelanálisisLISAparalavariablededesinfección.CálculosrealizadosconelsoftwareGeoDa(18).Elaboradoporelautor.
ANÁLISISESPACIALDELRIESGOSANITARIOPORCONCENTRACIÓNDEGRANJASAVICOLASSINBIOSEGURIDADENLAPROVINCIADETUNGURAHUA:UNAPROPUESTADEGESTIÓNPREVENTIVAPARALASOSTENIBILIDADEMPRESARIAL 42 Sehanidentificadotresnivelesdeprioridad(1:Intervencióninmediata;2:Acciónpreventivaymonitoreo;3:Vigilanciayestudioamedianoylargoplazo)enelriesgoAlto-Alto.LasparroquiasdePíllaroyCevallosalcontarconunasignificanciadep=0.001sonconsideradascomoprioridad1alserlosepicentrosenlamalagestiónycontroldedesinfecciónenlasgranjasavícolas,suintervencióndebeserinmediatayaquesinunprotocolodedesinfecciónpuedegenerarunambientepropicioparalaproliferaciónydiseminacióndeagentespatógenos,deigualmanerapuedeocasionarlapropagacióndeestasenfermedadesalasdemászonasaledañas.LasparroquiasdeIzambayPresidenteUrbinasonconsideradascomoprioridad2alcontarconunasignificanciadep=0.01,esdecir,estaszonasactúancomotransmisoresentrelosepicentrosprincipalesylasáreasdemenorriesgo,esdecir,dadoasuposiciónintermediaenelpaisajeepidemiológicolosconvierteenpotencialesamplificadoresdeladispersióndeenfermedadesalasdemászonassinriesgo,yesporelloqueserequieretomaraccionespreventivasparaconteneryevitarsupropagaciónalasdemászonas.LasparroquiasdeMarcosEspinel,BenítezyTotorassonconsideradascomoprioridad3debidoasusignificanciadep=0.05,enestaszonascuentanconunatendenciaincipientedelriesgoporloqueserecomiendarealizarintervencionestempranasmedianteestrategiaspreventivasantesdequeseconsolidencomofocosderiesgoestablecidos.Porotrolado,seidentificarondosclústeresdebajaincidencia(bajo-bajo)ubicadasenlasparroquiasdeGuambaló,Bolívar,PelileoySanFernando.Estosterritoriossedistinguenpormantenerprácticasconstantesdelimpiezaydesinfección,loquecontribuyeapreservarcondicionessanitariasestablesyasostenerdemaneracontinualaeficienciadesusprocesosproductivos.Figura4.ResumendeResultadosLISAparalavariabledeInsectos.
ANÁLISISESPACIALDELRIESGOSANITARIOPORCONCENTRACIÓNDEGRANJASAVICOLASSINBIOSEGURIDADENLAPROVINCIADETUNGURAHUA:UNAPROPUESTADEGESTIÓNPREVENTIVAPARALASOSTENIBILIDADEMPRESARIAL 43 Nota:Análisisdeautocorrelaciónespaciallocalparalavariabledeinsectos.CálculosrealizadosconelsoftwareGeoDa(18).Elaboradoporelautor.Figura5.MapadeSignificanciadelÍndicedeMoranLocalparalavariabledeInsectosNota:MapadeclústeresestadísticamentesignificativosdelanálisisLISAparalavariabledeinsectos.CálculosrealizadosconelsoftwareGeoDa(18).Elaboradoporelautor.LosresultadosmostraronlapresenciadeconglomeradosdeltipoAlto-Alto,enlosqueseagrupanparroquiasconniveleselevadosderiesgosanitariojuntoaotrasquepresentancondicionesparecidas,Estecomportamientoconfirmalaexistenciadeáreasespecíficasdondelavulnerabilidadsanitariaesparticularmentemarcadaanivelterritorial.Entreloshallazgosmásrelevantes,seidentificaronparroquiasconunniveldesignificanciaestadísticadep=0.001(MarcosEspinel,Píllaro,presidenteUrbinaeIzamba),representándosecomozonascríticascuyaconcentraciónderiesgonorespondealazar,sinoqueaunpatrónterritorialconsistenteyaltamenterobusto.Estosterritoriosconstituyenverdaderospuntoscalientesoepicentrosquerequierenatenciónprioritariaporpartedelasautoridadessanitariasylosdueñosdelasgranjasavícolas.Estepatrónsugierequelascondicionesambientales,elmanejoinadecuadoderesiduosorgánicosylalimitadaaplicacióndeproductosdecontrolcontribuyenacrearentornosfavorablesparalaproliferacióndeinsectosvectoresdeenfermedadesaviares.Deigualmanera,existenparroquiasconunasignificanciadep=0.01
ANÁLISISESPACIALDELRIESGOSANITARIOPORCONCENTRACIÓNDEGRANJASAVICOLASSINBIOSEGURIDADENLAPROVINCIADETUNGURAHUA:UNAPROPUESTADEGESTIÓNPREVENTIVAPARALASOSTENIBILIDADEMPRESARIAL 44 (Totoras)loquerefuerzalaevidenciadeconcentracióndelproblemaysubrayalanecesidaddeimplementarmedidaspreventivascoordinadasenelámbitointerparroquialalserconsideradasestaszonascomoconectoresotransmisoresdelasenfermedadesdesdeelepicentroalaszonasconbajosnivelesderiesgosanitario.Porotraparte,seidentificaron3zonas(Noroeste,SuryCentrodelaprovincia)quecuentanconclústeresderiesgoBajo-Bajoconfirmandoasíquesetratadezonasdondelaimplementacióndeprogramasregularesdecontrol,mantenimientosanitarioymayortecnificacióndelosgalponeshapermitidoreducirlosnivelesderiesgosanitario.Figura6.ResumendeResultadosLISAparalavariabledeRoedoresNota:Análisisdeautocorrelaciónespaciallocalparalavariablederoedores.CálculosrealizadosconelsoftwareGeoDa(18).Elaboradoporelautor.
ANÁLISISESPACIALDELRIESGOSANITARIOPORCONCENTRACIÓNDEGRANJASAVICOLASSINBIOSEGURIDADENLAPROVINCIADETUNGURAHUA:UNAPROPUESTADEGESTIÓNPREVENTIVAPARALASOSTENIBILIDADEMPRESARIAL 45 Figura7.MapadeSignificanciadelÍndicedeMoranLocalparalavariabledeRoedores.Nota:MapadeclústeresestadísticamentesignificativosdelanálisisLISAparalavariablederoedores.CálculosrealizadosconelsoftwareGeoDa(18).Elaboradoporelautor.Deacuerdoconlosresultadosobtenidos,lasparroquiasidentificadasconunasignificanciaestadísticadep=0.001(MarcosEspinel,Píllaro,presidenteUrbinaeIzamba)representanlosfocosdemayorvulnerabilidadsanitariaconrespectoaestavariable,lascualesreflejandebilidadesestructuralesenlagestiónsanitariaenlasgranjas.Estetipodefocosderiesgopuedentenerefectosdirectosenlapropagacióndepatógenoszoonóticos,especialmenteSalmonellayEscherichiacoli,entreotros,afectandotantolaproductividadavícolacomolasaludpública.Porotraparte,existenclústeresqueestánclasificadoscomoriesgoBajo-Altoconunasignificanciadep=0.001siendoestaslasparroquiasdeUnamuncho,SanAndrésyPicaihua,lasgranjasqueloconformanhanimplementadomedidasefectivasparacontrolaralosroedoresapesardeencontrasecercadezonasdealtaconcentracióndelriesgosanitario.También,existenzonasconunasignificanciadep=0.05clasificadacomoclústeresderiesgoAlto-Bajo(Cotalo,QueroySalasaca)loshallazgosmencionanqueestaszonassonconsiderascomounaalertatempranaparatomaraccionescorrectivasantesdequeelriesgollegueamayoresyseconsolidecomounriesgoAlto-Alto,porelloesimportanteincrementarlafrecuenciaderevisióndelosgalpones,mejorarsuinfraestructurayreforzarlosprotocolosdecontrolenestasgranjasespecíficas.Finalmente,seidentificaron82granjasavícolasdistribuidasen12parroquiasdelaprovinciadeTungurahuaquefueronclasificadasdentrodelacategoríaBajo-Bajo,loquereflejalapresenciadeterritoriosconprácticaspreventivasmásconsolidadas,unainfraestructuramásadecuadayuncontrolsanitariomásconsistente.
ANÁLISISESPACIALDELRIESGOSANITARIOPORCONCENTRACIÓNDEGRANJASAVICOLASSINBIOSEGURIDADENLAPROVINCIADETUNGURAHUA:UNAPROPUESTADEGESTIÓNPREVENTIVAPARALASOSTENIBILIDADEMPRESARIAL 46 Figura8.ResumendeResultadosLISAparalavariabledePerrosyGatos.Nota:Análisisdeautocorrelaciónespaciallocalparalavariabledeperrosygatos.CálculosrealizadosconelsoftwareGeoDa(18).Elaboradoporelautor.Figura9.MapadeSignificanciadelÍndicedeMoranLocalparalavariabledePerrosyGatos.Nota:MapadeclústeresestadísticamentesignificativosdelanálisisLISAparalavariabledeperrosygatos.CálculosrealizadosconelsoftwareGeoDa(18).Elaboradoporelautor.Losresultadosreflejanqueexistenclústeresconunasignificanciadep=0.001clasificadascomoriesgoAlto-Alto(Quero,Guambalo,Bolívar,Cunchibamba,UnamunchoyAugustoN.Martínez),indicandoasíunagranagrupaciónterritorialaltamentesignificativadegranjasavícolasquepresentandeficienciascríticasenelcontroldeanimalesdomésticos(perrosygatos)loscualesactúancomovectoresdeagentespatógenosdiseminadoresdeparásitoocontaminantesbiológicosafectandolabioseguridadylaeficienciadelossistemasproductivosavícolas.Deigualmanera,laparroquiadeTisaleomantieneunriegoAlto-Bajoyaquecuentaconunasignificanciadep=0.001,esdecir,queesteresultadoseñalalapresenciadegranjasavícolascongravesdeficienciasenelcontroldelasmascotas(perrosygatos)queseencuentrangeográficamenteinmersasenunazonadegranjasconbuenasprácticassanitarias,esdecir,queestasgranjasrepresentanaloseslabonesmásdébilesenelcontroldeestavariableyquecomprometesignificativamentelosesfuerzossanitariosdelasgranjasvecinasmejorgestionadas.
ANÁLISISESPACIALDELRIESGOSANITARIOPORCONCENTRACIÓNDEGRANJASAVICOLASSINBIOSEGURIDADENLAPROVINCIADETUNGURAHUA:UNAPROPUESTADEGESTIÓNPREVENTIVAPARALASOSTENIBILIDADEMPRESARIAL 47 Porotrolado,lasparroquiasdeCotaloySanBartolomédePinlloseclasificancomounriesgoBajo-Altoconunasignificanciadelp=0.05,esteresultadoindicalapresenciadegranjasconadecuadocontroldeperrosygatosquienesseencuentransituadasenzonasconpredominiodegranjascondeficienciasenelmanejoycontroldeestavariable.Finalmente,seidentificaron18granjasavícolasdistribuidasen8parroquiasdelaprovinciaquefueronclasificadasdentrodelacategoríaderiesgoBajo-Bajo.Estasunidadesproductivasseubicanenterritoriosdondepredominanprotocolosdeingresomásrigurosos,sistemasdecerramientoadecuadosycondicionesdemanejomejordesarrolladas,factoresquecontribuyenamantenerunamayorestabilidadsanitariayadisminuirlaprobabilidaddepérdidasenlaproducción.Figura10.ResumendeResultadosLISAparalavariabledeBioseguridad.Nota:Análisisdeautocorrelaciónespaciallocalparalosindicadoresdebioseguridad.CálculosrealizadosconelsoftwareGeoDa(18).Elaboradoporelautor.
ANÁLISISESPACIALDELRIESGOSANITARIOPORCONCENTRACIÓNDEGRANJASAVICOLASSINBIOSEGURIDADENLAPROVINCIADETUNGURAHUA:UNAPROPUESTADEGESTIÓNPREVENTIVAPARALASOSTENIBILIDADEMPRESARIAL 48 Figura11.MapadeSignificanciadelÍndicedeMoranLocalparalavariabledeBioseguridad.Nota:MapadeclústeresestadísticamentesignificativosdelanálisisLISAparalavariabledebioseguridad.CálculosrealizadosconelsoftwareGeoDa(18).Elaboradoporelautor.Sehanidentificadoclústeres(Píllaro,presidenteUrbina,Izamba,Atahualpa,AugustoN.Martínez,Ambato,SantaRosayJuanBenignoVela)conunriesgoAlto-Altoyconunasignificanciadep=0.001,estosresultadosevidencianlaexistenciadeunaagrupaciónterritorialaltamentesignificativadegranjasavícolascondeficienciascríticasygeneralizadasenbioseguridad,inmersasenunentornosimilarencaracterísticas.Laidentificacióndeestosclústeresdemandalaimplementacióninmediatadeunplandeemergenciasanitariaespecíficoparaestosterritorios,conasignaciónprioritariaderecursostécnicosyfinancieros,dadoquerepresentanlaprincipalamenazaparalasostenibilidaddelaaviculturaprovincialynacional.Porotraparte,sehaidentificadoclústeres(Quero,BañosyBolívar)conunriesgoAlto-Bajoyconunasignificanciadep=0.001,estosresultadosseñalanlapresenciadegranjasindividualescondeficienciascríticasenbioseguridadqueseencuentrangeográficamenteaisladasdentrodeterritoriosconadecuadosestándaresdebioseguridad.Asímismo,lasparroquiasdeConstantinoFernándezyUnamunchosonzonasqueseclasificancomounriesgoBajo-Altoyconunasignificanciadep=0.001,esteresultadoindicalapresenciadegranjasconexcelentesestándaresdebioseguridadqueseencuentraninmersasenterritoriosconpredominiodegranjasdealtoriesgosanitario.Finalmente,sehaidentificado171granjasavícolasubicadasen15parroquiasdelaprovinciaquefueronclasificadascomounriesgoBajo-Bajo,loqueindicaterritoriosconadecuadagestiónsanitaria,infraestructurafortalecidaycumplimientodenormasdebioseguridadloquefavorecelaeficienciaproductivaylaestabilidadeconómicadelsector.DISCUSIÓNLosresultadosobtenidosenesteestudiomuestranqueelriesgosanitarioenlaaviculturadelaprovinciadeTungurahuapresentapatronesterritorialesbiendefinidos,caracterizadosporlaconcentracióndeprácticasdeficientesdebioseguridad,especialmenteenvariablescomoelcontroldeinsectos,presenciaderoedoresyprocesosdedesinfección.Estoshallazgoscoincidenconestudiosinternacionalesquedemuestranquelabioseguridadnoactúacomounfactoraislado,sinocomounfenómenomulticausalinfluenciadoporelterritorio,elniveldetecnificaciónylasinteraccionesentreunidadesproductivas(24).
ANÁLISISESPACIALDELRIESGOSANITARIOPORCONCENTRACIÓNDEGRANJASAVICOLASSINBIOSEGURIDADENLAPROVINCIADETUNGURAHUA:UNAPROPUESTADEGESTIÓNPREVENTIVAPARALASOSTENIBILIDADEMPRESARIAL 49 ElanálisisdeclústeraplicadoenestainvestigaciónpermitióreconocerquelaszonasdemayorriesgosanitarioenTungurahuanoaparecendeformaaislada,sinoquerespondenaunadinámicaterritorialquetambiénhasidoobservadaendistintoscontextosinternacionales.DeacuerdoconKachigundaetal.(25),cuandolasfallasenbioseguridadserepitenenunterritorio,estastiendenaagruparseydarformaaáreascríticasquecompartencondicionessimilares,uncomportamientomuyparecidoalqueseevidencioenlasparroquiasconprácticasdeficientesdentrodelaprovincia.AlgosimilarocurrióenVietnam,despuésdeunlargoprocesosdevacunaciónagranescala,lapresenciadeclústeresderiesgosiguióestandovinculadaaterritoriosdondelaconcentraciónproductivaylascondicionesambientalesfacilitabanlapersistenciadelproblema(26).Alhaberaplicadouníndicederiesgosanitarioenlainvestigaciónpermitiósintetizardemaneraponderadalavulnerabilidadsanitariadecadaparroquia.HerramientascomolasdesarrolladasporVougatetal.(27)yNielsenetal.(28),muestrancómoesposibleasignarpuntajesaprácticascotidianas(desdelalimpiezayelcontroldeaccesohastalagestiónderesiduos)paraobtenerunalecturaobjetivadelriesgo.Otrostrabajos,comoeldeWielicketal.(29),demuestranqueestospuntajespuedenllevarseunpasomásallá,convirtiéndoseeníndicespredictivoscapacesdeanticiparlapresenciadeenfermedadesensistemaspecuarios.SeidentificóunaaltaautocorrelaciónespacialenlaconcentracióndegranjasavícolassinbioseguridadenlaprovinciadeTungurahua.Estoshallazgosmuestrancoherenciaconinvestigacionesprevias,deacuerdoconWilchesetal.(30),quienesreportaronunI=0.303entreelPM10ycasosdeCOVID-19enBogotá,detectandoasíclústeresAlto-Altoenzonassuroccidentalesdelaciudad,estosugierequelaconcentraciónlocaldeunriesgosueleasociarseconaumentoslocalesenelriesgosanitario.Asímismo,estudiosqueserealizaronenAragón–Españaevidenciaronqueelíndiceglobalvariaconeltiempo,losanálisisdeLISArevelanlaformacióndeconglomeradoslocalesquepuedenserefímerosopersistentesporelloesimportanterealizaranálisistemporalesdespuésdesuintervención(31).EsimportantedestacarquesielíndicedemoranhubiesetenidounaausenciadeautocorrelaciónesonohubieraimpedidolaidentificacióndeclústereslocalesdetectablesporLISA(32).Portanto,serecomiendapriorizarintervencionesfocalizadasenloscantonesconriesgoAlto-Altoidentificadosycomplementarelanálisisconindicadoresdemovilidadyrutasdecomercializaciónavícola.CONCLUSIONESEnconclusión,elanálisisdeclústerpermitióorganizaralasparroquiasdeTungurahuaengruposconcaracterísticasproductivasysanitariassemejantes.Estasegmentaciónpusoenevidenciaqueexistenterritoriosquecompartennivelessimilaresdecapacidadinstalada,cantidaddegalponesyprácticasdebioseguridad,locualfacilitódistinguirzonasconmayoresdesafíosyáreasconmenorvulnerabilidad.Dichoshallazgosconstituyenunrecursofundamentalparaorientaraccionesdeintervención,asistenciatécnicayplanificacióndiferenciadasegúnlasparticularidadesdecadaagrupaciónterritorial.
ANÁLISISESPACIALDELRIESGOSANITARIOPORCONCENTRACIÓNDEGRANJASAVICOLASSINBIOSEGURIDADENLAPROVINCIADETUNGURAHUA:UNAPROPUESTADEGESTIÓNPREVENTIVAPARALASOSTENIBILIDADEMPRESARIAL 50 LaelaboraciónyposterioraplicacióndelíndicedeRiesgoSanitarioParroquial(IRSP)permitiómediryordenarlosnivelesderiesgosanitariopresentesencadaunadelasparroquiasdelaprovinciadeTungurahua,loquehizoposibledistinguirzonasderiesgoalto,medioybajo.Losvaloresobtenidosevidenciarondebilidadesparticularmentemarcadasenelmanejoderoedores,lapresenciadeinsectosyelingresodeanimalesexternos,ademásderevelarinconsistenciasenlaimplementacióndelosprotocolosdedesinfección.Finalmente.Losresultadosobtenidosapartirdelosanálisisdeautocorrelaciónespacialconfirmaronqueleriesgosanitarionosedistribuyedemaneraaleatoria,sinoquetiendeaconcentrarseenzonasespecíficasdelaprovincia.TantoelíndicedeMoranGlobalcomoelestadísticoLocal(LISA)evidenciaronquelasparroquiasconmayoresnivelesdevulnerabilidadpresentanunaclaratendenciaaagruparse,generandoposiblesefectosdepropagaciónhaciaáreasvecinas.Estasituaciónponedemanifiestolaimportanciadegestionarelriesgodesdeunenfoqueterritorialycolectivo,superandolavisiónaisladadecadaunidadproductivaavícola.AGRADECIMIENTOSLosautoresagradecenalaDireccióndeInvestigaciónyDesarrollo(DIDE).EstetrabajodeTitulaciónespartedelosresultadosdelProyectodeInvestigaciónTitulado“FortalecimientocomercialparaeldesarrollosostenibledelasorganizacionesasociativasdeTungurahua,Ecuador”,conResoluciónNro.UTA-CONIN-2023-0330RPFCA29.Además,alosinvestigadoresdelgrupoDesarrolloFinancieroEmpresarial(DFE)ydelgrupoDesarrolloTerritorialEmpresaeinnovación(DeTEI)delaFacultaddeCienciasAdministrativasdelaUniversidadTécnicadeAmbato,quienespermitieronmejoraresteProyectodeInvestigación.REFERENCIASBIBLIOGRÁFICAS1.UNICEF.ReviewofnationalFood-BasedDietaryGuidelinesandassociatedguidanceforinfants,children,adolescents,andpregnantandlactatingwomen.NewYork.2.WHO.WorldHealthOrganization.[Online];2020.Disponibleen:https://www.who.int/news-room/fact-sheets/detail/healthy-diet.3.SpringmmannM,ClarkM,RaynerM,ScarboroughP,WebbP.Loscostosglobalesyregionalesdepatronesalimentariossaludablesysostenibles:unestudiodemodelado.LancetPlanetHealth.2021;5(11):p.797-807.4.OECD-FAO.AgriculturalOutlook2023-2032.Paris.5.FAO.WorldFoodandAgriculture-StatisticalYearbook.Roma:FoodandAgricultureOrganizationoftheUnitedNations;2024.6.CumbajinE,MenesesE,PilaM,QuevedoA,ToapantaA.DesarrollodelaaviculturaenEcuador,américayelmundoLatacunga:UniversidadTécnicadeCotopaxi;2022.7.MinisteriodeAgriculturayGanadería.Elhuevo,unsúperalimento,quegeneramásde300milempleosdirectos.Quito.
ANÁLISISESPACIALDELRIESGOSANITARIOPORCONCENTRACIÓNDEGRANJASAVICOLASSINBIOSEGURIDADENLAPROVINCIADETUNGURAHUA:UNAPROPUESTADEGESTIÓNPREVENTIVAPARALASOSTENIBILIDADEMPRESARIAL 51 8.EuropeanFoodSafetyAuthority,EuropeanCentreforDiseasePreventionandControl,EuropeanUnionReferenceLaboratoryforAvianInfl,AlexakisL,BuczkowskiH,DucatezM,etal.AvianinfluenzaoverviewMarch-June2025.EFSAJ.2025;25(7):p.9520-9585.9.PAHO-WHO.PanAmericanHealthOrganization.[Online].:PanAmericanHealthOrganization;2025..Disponibleen:https://www.paho.org/sites/default/files/2025-03/2025-mar-4-phe-epidupdate-avianinfluenza-eng-final.pdf.10.WHO.WorldHealthOrganization.[Online].;2025..Disponibleen:https://www.who.int/emergencies/disease-outbreak-news/item/2025-DON575.11.Agrocalidad.Resolución0286.Quito.12.SerranoC.ControlespermanentesrefuerzanlabioseguridadengranjasavícolasdelEcuador.Vistazo.28Octubre2024.13.Primicias.EscasezdemaízllevaallímiteapequeñasavícolasdeTungurahua,mientrassubeelpreciodelacubetadehuevos.Primicias.04Octubre2024.14.Avinews.Avinews.[Online].;2025.Acceso10deMarzode2025.Disponibleen:https://avinews.com/crecimiento-y-desafios-del-sector-avicola-de-postura-en-ecuador-perspectivas-para-2025/.15.ElUniverso.EngranjasdeTungurahuasedetectandosnuevosbrotesdeinfluenzaaviar.ElUniverso..16.ÁlavaI,GordonF.ModelodeaprendizajeautomatizadoparaelcálculodelaprimavehicularysegmentacióndelconsumidorGuayaquil:UniversidadCatólicadeSantiagodeGuayaquil;2025.17.FermínR.RevisióndemétodosdeclusteringparadatosfuncionalesCoruña:UniversidaddeCoruña;2024.18.AnselinL,SI,&KY.GeoDa:AnIntroductiontoSpatialDataAnalysis.GeographicalAnalysis.2006;38(1):p.5–22.19.PérezD,RodriguezH.AnálisisespacialdelpagodelimpuestopredialenlaciudaddeCartagena.Entramado.2022;18(1):p.1-18.20.ChacónC,RamírezM,RangelH.DinámicasespacialesdelíndicedepobrezamultidimensionalenColombia,2005-2018.Revistaciudades,estadosypolítica.2024;11(2):p.15-41.21.FernándezP,HerreraM.RegresionesSUREspaciales.AnálisisEspacio-temporaldelEmpleoSectorialenArgentina.RedNacionaldeInvestigadoresenEconomía.2023;(279):p.1-25.22.MoraM,GeldesC,MoralesC.Innovaciónenagronegociosencadenasdevalorenterritoriosconaltopotencialdecrecimiento.LoscasosdeEcuador,ParaguayyNicaraguaRoma:EconomíadeldesarrolloagrícoladelaFAO;2024.23.EscuderoM,OjedaS,MoyaL,EnriqueI.Lamallaestadísticacomounidaddeanálisisespacial.Razóndemortalidad,pobaciónyvivienda.EURE.2024;50(150):p.1-24.24.SánchezR,PérezI.Pertinenciadelconocimientoycumplimientodelabioseguridadparaelprofesionaldelasalud.HumanidadesMédicas.2021;21(1):p.239-258.25.KachigundaB,MengersenK,PereraD,CouplandG,VanderMerweJ,McKirdyS.Useofmixed-typedataclusteringalgorithmforcharacterizingtemporalandspatialdistributionofbiosecurityborderdetectionsofterrestrialnon-indigenousspecies.PLoSONE.2022;17(8).
ANÁLISISESPACIALDELRIESGOSANITARIOPORCONCENTRACIÓNDEGRANJASAVICOLASSINBIOSEGURIDADENLAPROVINCIADETUNGURAHUA:UNAPROPUESTADEGESTIÓNPREVENTIVAPARALASOSTENIBILIDADEMPRESARIAL 52 26.HenningJ,PfeifferD,VuleT.RiskfactorsandcharacteristicsofH5N1HighlyPathogenicAvianInfluenza(HPAI)post-vaccinationoutbreaks.VetRes.2009;40(3).27.VougatR,LeiteM,TilliG,LaconiA,MahmoodQ,ProdanovJ,etal.BiosecurityimplementationinpoultryfarmsacrossEuropeandneighboringcountries:asystematicreview.FrontiersinVeterinaryScience.2025;12.28.NielsenS,KemuntoN,MuloiD,BojesenA,KnightT,AyebareD,etal.Flexibleassessmentofbiosecurityinsmallandmediumscalepoultryfarmsinlowandmiddleincomecountries.ActaVeterinariaScandinavica.2025;69(9).29.WielickC,LudwigL,RibbensS,ThiryÉ,FaesC,SaegermanC.BiosecurityRiskFactorsandPredictiveIndexforHepatitisEVirusSerologicalStatusinBelgianPigFarms:ConventionalandFree-RangeSystems.Viruses.2025;17(432).30.WilchesJ,DíazH,CastilloM.AsociaciónespacialeimpactodevariablessocioambientalesenloscasosdeCOVID-19enBogotá(Colombia).InformaciónTecnológica.2022;33(2):p.67-76.31.EscolanoS,SalvadorJ.RedularidadglobalyvariabilidadlocaldelospatronesespaciotemporalesdelaCOVID-19enAragón(España).BoletíndelaAsociacióndeGeógrafosEspañoles.2022;1(93):p.1-31.32.RubioL,IpanaquéE,RojasM.AutocorrelacióngeoespacialdepartamentaldelCOVID-19asociadoconenfermedadesdecomorbilidad:Perú2020-2021.RevistadeInvestigaciónEstadística.2024;6(7):p.76-87.