Hola David, Si entiendo correctamente, lo siguiente deberia hacer lo que describes:
set.seed(123) A <- matrix(rpois(20, 4), ncol = 4) B <- matrix(rpois(274*3, 4), ncol = 3) # el resultado es una matriz de orden 5x3 sapply(1:ncol(B), function(.col) sapply(A[, .col + 1], function(x) sum(x > B[, .col]))) # [,1] [,2] [,3] #[1,] 5 253 243 #[2,] 112 119 69 #[3,] 243 164 3 #[4,] 112 119 69 #[5,] 112 18 260 Este resultado indica que el elemento A[1, 2] es 5 veces mayor que todos los elementos de B[,1]; A[2, 2] es 112 veces mayor; A[3,2] es 243 veces mayor; A[4,2] es 112 veces mayor y finalmente que A[5,2] es 112 veces mayor que B[, 1]. Similarmente, A[5, 4] es 260 veces mayor que los elementos de B[,3]. Un por uno seria algo como R> sapply(A[,2], function(x) sum(x > B[,1])) [1] 5 112 243 112 112 R> sapply(A[,3], function(x) sum(x > B[,2])) [1] 253 119 164 119 18 R> sapply(A[,4], function(x) sum(x > B[,3])) [1] 243 69 3 69 260 Saludos cordiales, Jorge.- 2015-03-09 16:04 GMT+11:00 David Contreras <[email protected]>: > Buen d�a amigos, > > Antes que nada, reitero mis agradecimientos por sus respuestas tan precisas > y oportunas, sigo en mi trabajo de grado y tengo dificultad con un proceso: > > Tengo dos matrices A y B, A es de dimensi�n 5x4 y B es de dimensi�n 274x3. > Requiero de un proceso en el que tome el primer elemento de la Matriz A > columna 2 y lo compare con todos los elementos de la Matriz B columna 1 y > me diga en un vector de dimensi�n (1x5) cuantas veces fue el elemento de A > mayor que cada uno de los 274 elementos de B y asi hasta llegar a los 5 > elementos de A, es decir el resultado de este proceso sera un vector con 5 > elementos y en cada uno me dice cuentas veces fue mayor el elemento [i,j] > mayor que los 274 elementos de B. > > Agradezco su colaboraci�n con el tema que como siempre es muy valiosa. > > > Saludos, > > David Contreras. > > [[alternative HTML version deleted]] > > _______________________________________________ > R-help-es mailing list > [email protected] > https://stat.ethz.ch/mailman/listinfo/r-help-es > [[alternative HTML version deleted]]
_______________________________________________ R-help-es mailing list [email protected] https://stat.ethz.ch/mailman/listinfo/r-help-es
