estoy realizando un Análisis de Componentes Principales (ACP) en R para resumir los patrones de correlación mensual encontrados en 30 regiones muestreadas. Para ello, estoy utilizando como input una matriz de 10 (meses) x 30 (regiones) -descargar aquí-.
En un principio estaba trabajando con la función psych::principal
, pero al introducir la matriz obtengo el siguiente error:
pc <- principal(dat,nfactors = 2, residuals = FALSE, rotate="varimax", n.obs=NA, covar=TRUE,scores=TRUE, missing=FALSE, impute="median", oblique.scores=TRUE, method="regression")
Warning messages: 1: In cor.smooth(r) : Matrix was not positive definite, smoothing was done 2: In principal(dat) : The matrix is not positive semi-definite, scores found from Structure loadings
La función se ejecuta, el problema es que pc$weights es igual que pc$loadings. Cuando la matriz es cuadrada, o el número de columnas es menor al de las filas se obtienen los weights sin problema, no obstante, no es el caso de mis datos.
La función principal
, usando los datos de ejemplo de la función, me permite extraer directamente los loadings
, scores
y weights
. Me es importante extraer los weights
de las componentes principales pues necesito expresar los scores
en la misma unidad de correlación que la variable de entrada.
He probado la función stats::prcomp
con mis datos y no me alerta de ningún error, pero no devuelve a simple vista los weights tal cual se definen en principal
.
Agradezco de antemano cualquier ayuda facilitada. Gracias.
principal
una matriz de correlaciones? Eso parece viendo los datos (entre 1 y -1), sin embargo no es rectangular. En la llamada estás pasando el argumentocovar=FALSE
, que se usa para dar a la función los datos en crudo y que esta haga primero la matriz de correlación/covarianza y después el PCA. Segúnhelp(principal)
"covar: If false, find the correlation matrix from the raw data or convert to a correlation matrix if given a square matrix as input." Puede ser por ahí el problema.