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