Skip to main content
se añadieron 1035 caracteres en el cuerpo
Origen Enlace
Patricio Moracho
  • 61.1k
  • 12
  • 42
  • 72

Si entendí tu algoritmoLo que buscas es aplicar por fila un calculo como este:

(fila_x %*% B) %*% t(fila_x)

Ciertamente el resultado en este caso sería un matriz de 1 x 1, pero tal vez te encontraste con un problema por eso en tu formula indicas la multiplicación elemento a elemento fila_de_A %*% B * transpuesta(fila_de_A), y este es fácilun problema causado por un comportamiento por defecto de aplicarR, muchas veces cuestionado, que es la promoción, cuando corresponde, de un objeto más complejo a uno más simple, en este caso, cuando seleccionas una iteración porsola fila, por ej: A[1, ] obtienes un vector simple y no una matriz, por lo que la multiplicación matricial dejaría de funcionar como tal. Esto se puede corregir mediante el parámetro drop:

apply(AA[1, ]
[1] 1 0 0

A[1, FUN, =drop=FALSE]
 function(fila) {fila %*% B *[,1] t(fila)})
[,2] [,3]
[1,]    1    0 [,1]   0
 
dim(A[1, ])   [,2]       [,3]   # Un vector simple no tiene [,4]dimensiones
[1NULL

dim(A[1,] 0.062738, 0.054730286drop=FALSE]) 0.04672257# una matriz si
[1] 1 3

Entendiendo que siendo tal vez este el problema que moviliza tu pregunta, podría sugerirte algunos cambios en tú código:

result <-0 data.009331426frame()
[2for(i in 1:dim(A)[1]) {
  operacion = (A[i,] 0.000000, 0.004938286drop 0.03576857= FALSE] 0.976556574%*% B) %*% t(A[i, , drop=FALSE])
[3  result = rbind (result,] operacion)
}
result

          V1
1 0.00000006273800
2 0.00000000005966857
3 0.00000000 08249114
4 0.00000000096722515

Básicamente aplicamos la formulaSimplemente vamos iterando por fila, sin embargo, entiendo que por el ejemplo que estás dando, pareciera más bien que y aplicamos la formula es: Suma(fila_de_A %*% B * transpuesta(fila_de_A)),pero esta vez asegurándonos que básicamente podrías implementar desde la matriz generadasiempre trabajemos con matrices.

También puedes simplificar todo en una sola línea y evitas el potencial problema de usar applyrbind dentro de un ciclo que tiene una penalidad importante en performance:

colSumssapply(apply1:dim(A)[1], 1,
       FUN = function(filai) {fila(A[i, , drop = FALSE] %*% B) *%*% t(filaA[i, , drop=FALSE])}))

[1] 0.06273800 0.05966857 0.08249114 0.96722515

Si entendí tu algoritmo: fila_de_A %*% B * transpuesta(fila_de_A) es fácil de aplicar una iteración por fila:

apply(A, 1, FUN = function(fila) {fila %*% B * t(fila)})

         [,1]        [,2]       [,3]         [,4]
[1,] 0.062738 0.054730286 0.04672257 -0.009331426
[2,] 0.000000 0.004938286 0.03576857  0.976556574
[3,] 0.000000 0.000000000 0.00000000  0.000000000

Básicamente aplicamos la formula por fila, sin embargo, entiendo que por el ejemplo que estás dando, pareciera más bien que la formula es: Suma(fila_de_A %*% B * transpuesta(fila_de_A)), que básicamente podrías implementar desde la matriz generada con apply:

colSums(apply(A, 1, FUN = function(fila) {fila %*% B * t(fila)}))

[1] 0.06273800 0.05966857 0.08249114 0.96722515

Lo que buscas es aplicar por fila un calculo como este:

(fila_x %*% B) %*% t(fila_x)

Ciertamente el resultado en este caso sería un matriz de 1 x 1, pero tal vez te encontraste con un problema por eso en tu formula indicas la multiplicación elemento a elemento *, y este es un problema causado por un comportamiento por defecto de R, muchas veces cuestionado, que es la promoción, cuando corresponde, de un objeto más complejo a uno más simple, en este caso, cuando seleccionas una sola fila, por ej: A[1, ] obtienes un vector simple y no una matriz, por lo que la multiplicación matricial dejaría de funcionar como tal. Esto se puede corregir mediante el parámetro drop:

A[1, ]
[1] 1 0 0

A[1, , drop=FALSE]
     [,1] [,2] [,3]
[1,]    1    0    0
 
dim(A[1, ])             # Un vector simple no tiene dimensiones
NULL

dim(A[1, , drop=FALSE]) # una matriz si
[1] 1 3

Entendiendo que siendo tal vez este el problema que moviliza tu pregunta, podría sugerirte algunos cambios en tú código:

result <- data.frame()
for(i in 1:dim(A)[1]) {
  operacion = (A[i, , drop = FALSE] %*% B) %*% t(A[i, , drop=FALSE])
  result = rbind (result, operacion)
}
result

          V1
1 0.06273800
2 0.05966857
3 0.08249114
4 0.96722515

Simplemente vamos iterando por fila y aplicamos la formula pero esta vez asegurándonos que siempre trabajemos con matrices.

También puedes simplificar todo en una sola línea y evitas el potencial problema de usar rbind dentro de un ciclo que tiene una penalidad importante en performance:

sapply(1:dim(A)[1], 
       FUN = function(i){(A[i, , drop = FALSE] %*% B) %*% t(A[i, , drop=FALSE])})

[1] 0.06273800 0.05966857 0.08249114 0.96722515
Origen Enlace
Patricio Moracho
  • 61.1k
  • 12
  • 42
  • 72

Si entendí tu algoritmo: fila_de_A %*% B * transpuesta(fila_de_A) es fácil de aplicar una iteración por fila:

apply(A, 1, FUN = function(fila) {fila %*% B * t(fila)})

         [,1]        [,2]       [,3]         [,4]
[1,] 0.062738 0.054730286 0.04672257 -0.009331426
[2,] 0.000000 0.004938286 0.03576857  0.976556574
[3,] 0.000000 0.000000000 0.00000000  0.000000000

Básicamente aplicamos la formula por fila, sin embargo, entiendo que por el ejemplo que estás dando, pareciera más bien que la formula es: Suma(fila_de_A %*% B * transpuesta(fila_de_A)), que básicamente podrías implementar desde la matriz generada con apply:

colSums(apply(A, 1, FUN = function(fila) {fila %*% B * t(fila)}))

[1] 0.06273800 0.05966857 0.08249114 0.96722515