El cambio climático genera diferentes respuestas en el patrón de distribución de la flora y fauna. Según las predicciones del Panel Intergubernamental del Cambio Climático (IPCC) por medio de modelos globales, las zonas altas son las que sufrirán efectos mayores, debido a los cambios drásticos en variables ambientales como la temperatura y precipitación. Por ello se realizaron modelos de nicho ecológico con el algoritmo MaxEnt (basado en máxima entropía) para las especies: Comarostaphylis arbutoides (Lindl.), Ilex pallida (Standl.), Myrrhidendron donnellsmithii (J.M. Coult. & Rose), Quercus costaricensis (Liebm.) y Vaccinium consanguineum (Klotzsch). A partir de las proyecciones al 2050 y 2070 de Worldclim (capas con información del clima global), se generaron los modelos en los escenarios de Representive Concentration Pathways de: RCP 2,6 (mitigación), RCP 4,5 (estabilización) y RCP 8,5 (pesimista). Las variables de mayor aporte fueron la temperatura máxima del periodo más caliente y la temperatura media en el trimestre más lluvioso; se utilizaron el test de Jackknife y el Área bajo la curva (AUC) para la validación y análisis de los modelos. La mayoría de las modelaciones muestran una tendencia a la reducción del área a mediano y largo plazo en comparación con nicho ecológico actual. Sin embargo, la especie Ilex pallida tiene proyecciones favorables en la colonización de nuevas áreas para los RCP 2,6 y RCP 4,5 en lugares de mayor elevación que las actuales.