La forma más fácil de consultarlo es con summary(objeto_dudi)
, donde objeto_dudi
es el nombre del modelo que ajustaste. Y es consultarlo porque extrañamente el summary()
no produce un output propiamente dicho, así que no podemos asignarle un nombre y conservarlo como una estructura de datos.
De todos modos es muy fácil extraer la información directamente desde el objeto dudi.hillsmith
. La varianza explicada por un eje es el eigenvalor de ese eje como proporción de toda la varianza, es decir, de la suma de los eigenvalores. Entonces la obtendrías con objeto_dudi$eig/sum(objeto_dudi$eig)
. Creo que en el ejemplo queda más claro.
library(ade4)
data("dunedata") #datos de ejemplo
objeto_dudi <- dudi.hillsmith(dunedata$envir, scann = FALSE, nf=2)
summary(objeto_dudi)
Muestra esto en pantalla, pero no hay forma de pasarlo a una estructura de datos. Usa la terminología "inercia proyectada". Es normal, hay pocos acuerdos de lenguaje en este campo y cada autor/a usa su propia nomenclatura. Greenacre dice que
... podemos interpretar este porcentaje exactamente igual que, en regresión, explicamos el "porcentaje de varianza explicada". (p. 68)
Class: mix dudi
Call: dudi.hillsmith(df = dunedata$envir, scannf = FALSE, nf = 2)
Total inertia: 8
Eigenvalues:
Ax1 Ax2 Ax3 Ax4 Ax5
2.5421 1.8578 1.2306 0.9899 0.6927
Projected inertia (%):
Ax1 Ax2 Ax3 Ax4 Ax5
31.776 23.222 15.383 12.374 8.658
Cumulative projected inertia (%):
Ax1 Ax1:2 Ax1:3 Ax1:4 Ax1:5
31.78 55.00 70.38 82.76 91.41
(Only 5 dimensions (out of 8) are shown)
Extraer los datos del modelo
varianza_explicada <- objeto_dudi$eig/sum(objeto_dudi$eig)
# Le agrego nombres para facilitar la lectura
names(varianza_explicada) <- paste0("Eje", 1:length(objeto_dudi$eig))
varianza_explicada
Eje1 Eje2 Eje3 Eje4 Eje5 Eje6 Eje7 Eje8
0.31776060 0.23222407 0.15382594 0.12374192 0.08658363 0.05134938 0.02735787 0.00715659
Esto genera un vector con la varianza explicada por cada eje que se puede utilizar para incluirlo en tablas, gráficos, etc. Los resultados no son exactamente iguales al summary()
porque son proporciones en lugar de porcentajes y no están redondeados.