Nuestro estudio se llevó a cabo en el Parque Nacional Valle Nuevo e incluyó cuatro clases de hábitat: pastizal de tussock (Sabapa), bosque de pinos (Pinoc), bosque de hojas anchas (Boslat) y ecosistema agrícola (Ecoag). Tuvimos dos objetivos principales: describir comparativamente las comunidades de milpiés y determinar las relaciones entre la densidad/diversidad poblacional y las variables fisicoquímicas del suelo. La investigación fue transversal y no manipulativa, con un alcance descriptivo y correlacional; el muestreo siguió un diseño sistemático estratificado, con ocho transectos y 32 cuadrantes de 1 m, cubriendo 21.7 km. Encontramos un suelo franco arenoso con un pH extremadamente ácido. La mayor densidad poblacional de milpiés se registró en Sabapa, y la más baja en Ecoag. La mayor diversidad alfa se compartió entre Boslat (Margalef = 1.72) y Pinoc (Shannon = 2.53); Sabapa y Boslat mostraron la mayor similitud de Jaccard (0.56). La prueba de hipótesis nula utilizando el índice de Shannon ponderado reveló una diferencia estadísticamente significativa en la diversidad entre los pares Boslat-Sabapa y Pinoc-Sabapa. Dos de las especies registradas mostraron valores indicativos (IndVal) altamente significativos para dos clases de hábitat. Encontramos correlaciones significativas (< 0.05) entre varias variables fisicoquímicas del suelo y la densidad y diversidad de milpiés.