Hola,
Tengo una matriz numérica (A) con n columnas y k filas. Quiero construir una rutina que me permita ordenar cada columna de la matriz por separado manteniendo la información sobre el índice original [usando la función sort.int(…,index.return=T)] y que me devuelva una matriz de matrices (B) compuesta por n nuevas matrices, cada una con un numero de k filas y 2 columnas en que la primera columna contenga los valores ordenados y la segunda los índices originales.
Mi primer intento fue:
> B<-matrix(data=NA,0, 0)
> for (i in n) {B[,i]<-sort.int(A[,i], method=c("quick"), index.return=TRUE)}
Pero R me devolvió el siguiente error:
“Error en inherits(x, "factor") : subíndice fuera de los límites”
La cuestión es que no puedo generar la matriz B manualmente porque tengo que aplicar este procedimiento a un gran numero de matrices iniciales (A’s) con dimensiones variables.
¿Tenéis alguna idea de cómo hacerlo?
Muchas gracias!