Sobre la pregunta
Al ver la estructura del resultado final, se está hablando de una matriz triangular inferior.
Aquí se detalla cómo se implementó con más detalle en R lo anterior
Como obtener promedios de una matriz triangular sin emplear condicionales
La única diferencia sería calcular la expresión: ( x / x_previo - 1 ) * 100
, proporcionada en la pregunta, en lugar del promedio.
Datos
df <- read.table(text="
Codigo PXN
455-1 564
455-1 4568
455-1 213
455-1 564
321-2 678
321-2 1558
321-2 2314
321-2 12353
369-8 1524
369-8 978
369-8 234
369-8 2347
", header = TRUE)
Si bien los códigos se repiten igual cantidad de veces (4 por cada tipo) se está considerando el caso en que no lo sean.
Adicional a lo anterior hay que capturar los elementos PXN (estos serán las diagonales de las matrices triangulares inferiores) para cada código a través de un bucle for
codigos <- unique(as.vector(df$Codigo))
print(codigos)
# "455-1" "321-2" "369-8"
Matriz Triangular
Refactorizando el código del post que se mencionó más arriba
ratiosTriangularInferior <- function(diagonal, codigo){
len <- length(diagonal)
m <- matrix(0, nrow=len, ncol = len)
porcentaje <- c()
label <- c()
for (j in 1:(len-1)) { # columnas
for (i in (j+1):len){ # filas
print(paste0('(',diagonal[j],'/',diagonal[i],'-1)*100'))
m[i,j] <- diagonal[j]/diagonal[i] - 1
label <- c(label,paste0(diagonal[j],'-',diagonal[i]))
porcentaje<-c(porcentaje,m[i,j]*100)
}
}
print(round(m, digits = 2))
porcentaje <- round(porcentaje, digits = 2)
mm <- matrix(c(rep(codigo,times=length(label)),label,porcentaje),ncol=3)
mm
}
Los prints
que contienen se pueden omitir pero permiten visualizar el paso a paso:
# [1] "455-1"
# [1] "(564/4568-1)*100"
# [1] "(564/213-1)*100"
# [1] "(564/564-1)*100"
# [1] "(4568/213-1)*100"
# [1] "(4568/564-1)*100"
# [1] "(213/564-1)*100"
# [,1] [,2] [,3] [,4]
# [1,] 0.00 0.00 0.00 0
# [2,] -0.88 0.00 0.00 0
# [3,] 1.65 20.45 0.00 0
# [4,] 0.00 7.10 -0.62 0
# [1] "321-2"
# [1] "(678/1558-1)*100"
# [1] "(678/2314-1)*100"
# [1] "(678/12353-1)*100"
# [1] "(1558/2314-1)*100"
# [1] "(1558/12353-1)*100"
# [1] "(2314/12353-1)*100"
# [,1] [,2] [,3] [,4]
# [1,] 0.00 0.00 0.00 0
# [2,] -0.56 0.00 0.00 0
# [3,] -0.71 -0.33 0.00 0
# [4,] -0.95 -0.87 -0.81 0
# [1] "369-8"
# [1] "(1524/978-1)*100"
# [1] "(1524/234-1)*100"
# [1] "(1524/2347-1)*100"
# [1] "(978/234-1)*100"
# [1] "(978/2347-1)*100"
# [1] "(234/2347-1)*100"
# [,1] [,2] [,3] [,4]
# [1,] 0.00 0.00 0.0 0
# [2,] 0.56 0.00 0.0 0
# [3,] 5.51 3.18 0.0 0
# [4,] -0.35 -0.58 -0.9 0
Pero cómo el OP no indica la estructura que desea asumo que debe ser un dataframe
por tanto la respuesta final se dará en este formato
Bucle For
Se llama la función ratiosTriangularInferior
(3 veces en este caso) para obtener las combinaciones de los porcentaje por cada tipo de código
names <- c("Codigo","Ratio-PXN","Porcentaje")
ans <- matrix(names, ncol=3) # aca se almacenan los resultados parciales
for (i in 1:length(codigos)){
diag <- df[df$Codigo==codigos[i],]$PXN
print(codigos[i])
ans <- rbind(ans, ratiosTriangularInferior(diag, codigos[i]))
cat("\n")
}
Se le da la estructura al dataframe
ans <- as.data.frame(ans)
ans <- ans[-1, ]
colnames(ans) <- names
rownames(ans) <- NULL
print(ans)
# Codigo Ratio-PXN Porcentaje
# 1 455-1 564-4568 -87.65
# 2 455-1 564-213 164.79
# 3 455-1 564-564 0
# 4 455-1 4568-213 2044.6
# 5 455-1 4568-564 709.93
# 6 455-1 213-564 -62.23
# 7 321-2 678-1558 -56.48
# 8 321-2 678-2314 -70.7
# 9 321-2 678-12353 -94.51
# 10 321-2 1558-2314 -32.67
# 11 321-2 1558-12353 -87.39
# 12 321-2 2314-12353 -81.27
# 13 369-8 1524-978 55.83
# 14 369-8 1524-234 551.28
# 15 369-8 1524-2347 -35.07
# 16 369-8 978-234 317.95
# 17 369-8 978-2347 -58.33
# 18 369-8 234-2347 -90.03