Tengo unas matrices
del tipo
A <- matrix(c(1, 0, 0
1, 0.1,0
1, 0.2,0
1, 0.9 0), ncol = 3, byrow = TRUE) #En mis datos reales es una matriz 121*3
B <- matrix(c(0.06273800 -0.08007714 -0.03178301
-0.08007714 1.29460000 0.20540957
-0.03178301 0.20540957 0.23629000), ncol = 3, byrow = TRUE)
Y debo operar de tal forma que la primera fila de A la pueda multiplicar por B y luego por la traspuesta de la primera fila de A. Luego la segunda fila de A, por B, por a segunda fila de A, y así sucesivamente fila a fila Ejemplo
[1, 0, 0] %*% B * [1
0
0] = 0.062738 #Este resultado lo quiero almacenar ordenadamente
[1, 0.1, 0] %*% B * [1
0.1
0] = 0.05966857 Y así sucesivamente y separadamente por cada fila de la matriz.
De esta forma, cada operación que esté realizando involucra a la primera fila de A, que tiene una dimensión de 1x3, a B que tiene 3x3, y a la traspuesta de A que tiene 3x1, de tal forma que el resultado es un escalar [1,1].
Y lo quiero escribir en una salida:
result <- data.frame(c(0.062738, 0.05966857, 0.08249114, 0.9672252))
Intenté hacerlo con el siguiente bucle, pero como pueden ver en mis comentarios, hay partes que no funcionan
resul <- data.frame()
for(i in seq_len(nrow(A))){ #Compruebo que recorre todas las filas de todas las columnas
print(combi[,i])
for(j in seq_len(ncol(A))){
print(combi[i,j]) #Compruebo que va leyendo la primera fila de las 3 columnas
#bec = as.vector([i,j]) #Intento crear un vector que me almacene la 1,0,0 y resuelva
"Operación". No funciona
#operacion = bec *%* P * t(bec)
result = rbind (result, operacion)
}
}
Cómo debo hacerlo?