Estoy construyendo una matriz (20x20
) para unos indicadores que se tienen que emplear más adelante en un Balance Scorecard.
Para efectuar las estimaciones correspondientes, se toma la diagonal
cuyos valores son unos rendimientos provenientes de históricos.
Con dicha diagonal
se obtienen promedios
para obtener una matriz triangular inferior: lower triangular matrix
Para que se visualice lo anterior aquí un pequeño esquema:
Dado que la matriz es relativamente pequeña (20x20
) es válido emplear un ataque de fuerza bruta
promedioMatrizTriangular <- function(diagonal){
len <- length(diagonal)
m <- matrix(0, nrow=len, ncol = len)
m <- diag(diagonal)
for (j in 1:len) { # columnas
for (i in 1:len){ # filas
if(i>j) m[i,j]<-mean(diagonal[j:i])
}
}
round(m, digits = 1)
}
Corriendo la función anterior con la diagonal mostrada en el esquema previo, se tiene:
promedioMatrizTriangular(c(9,11,8,5))
[,1] [,2] [,3] [,4]
[1,] 9.0 0.0 0.0 0
[2,] 10.0 11.0 0.0 0
[3,] 9.3 9.5 8.0 0
[4,] 8.2 8.0 6.5 5
Es decir funciona, pero... si se tuviera una matriz de
1000x1000
eso implicaría ejecutarif(i>j)
un millón de veces.
Entonces para la respuesta, cualquier enfoque que se emplee y que conlleve a una reducción de ese millón de operaciones, será aceptada como valida.
Nota: Algunas aplicaciones prácticas de la matriz triangular: