library(dplyr)
## Warning: package 'dplyr' was built under R version 4.2.3
##
## Attaching package: 'dplyr'
## The following objects are masked from 'package:stats':
##
## filter, lag
## The following objects are masked from 'package:base':
##
## intersect, setdiff, setequal, union
library(readxl)
MIP_49 <- read_excel("~/Downloads/MIP_49.xlsm")
X <- MIP_49[,1]
DF <- MIP_49[,22]
x <- MIP_49[,c(2:21)]
tX <- t(X)
a <- sweep(x, 2, tX, "/")
Id = diag(1, nrow = 20)
aL <- solve(Id - a)
EDC = colSums(aL)
EDCk = EDC/mean(EDC)
EDCk
## [1] 1.0573726 1.0075159 0.9616262 1.2054900 1.0977869 0.9107823 0.9013177
## [8] 1.0035968 1.0725182 1.0580281 0.8156409 0.9305823 1.1251490 0.9336891
## [15] 0.8502004 1.0362611 1.0020936 1.0531228 0.9540608 1.0231653
EDR = rowSums(aL)
EDRk = EDR/mean(EDR)
t(EDRk)
## 11 - Agricultura, cría y explotación de animales, aprovechamiento forestal, pesca y caza
## [1,] 0.9389028
## 21 - Minería
## [1,] 0.8808883
## 22 - Generación, transmisión, distribución y comercialización de energía eléctrica, suministro de agua y de gas natural por ductos al consumidor final
## [1,] 0.9458046
## 23 - Construcción 31-33 - Industrias manufactureras
## [1,] 0.8686232 2.417271
## 43 - Comercio al por mayor 46 - Comercio al por menor
## [1,] 1.029677 1.028826
## 48-49 - Transportes, correos y almacenamiento
## [1,] 0.9885228
## 51 - Información en medios masivos 52 - Servicios financieros y de seguros
## [1,] 0.9259019 0.9949904
## 53 - Servicios inmobiliarios y de alquiler de bienes muebles e intangibles
## [1,] 1.081696
## 54 - Servicios profesionales, científicos y técnicos 55 - Corporativos
## [1,] 1.054527 0.9078511
## 56 - Servicios de apoyo a los negocios y manejo de residuos, y servicios de remediación
## [1,] 1.291764
## 61 - Servicios educativos 62 - Servicios de salud y de asistencia social
## [1,] 0.7554517 0.7837693
## 71 - Servicios de esparcimiento culturales y deportivos, y otros servicios recreativos
## [1,] 0.7487236
## 72 - Servicios de alojamiento temporal y de preparación de alimentos y bebidas
## [1,] 0.7798288
## 81 - Otros servicios excepto actividades gubernamentales
## [1,] 0.829237
## 93 - Actividades legislativas, gubernamentales, de impartición de justicia y de organismos internacionales y extraterritoriales
## [1,] 0.7477418
#Cambiar DF
DFn = DF
DFn[10,] = 1074748.784*1.1
Xn = as.matrix(aL)%*%as.matrix(DFn)
Xn
## DF - Demanda final
## 11 - Agricultura, cría y explotación de animales, aprovechamiento forestal, pesca y caza 1240730.6
## 21 - Minería 1461409.1
## 22 - Generación, transmisión, distribución y comercialización de energía eléctrica, suministro de agua y de gas natural por ductos al consumidor final 736572.0
## 23 - Construcción 2963444.4
## 31-33 - Industrias manufactureras 15977970.1
## 43 - Comercio al por mayor 2895139.2
## 46 - Comercio al por menor 2891109.3
## 48-49 - Transportes, correos y almacenamiento 2691735.5
## 51 - Información en medios masivos 667799.0
## 52 - Servicios financieros y de seguros 1489068.2
## 53 - Servicios inmobiliarios y de alquiler de bienes muebles e intangibles 2344134.0
## 54 - Servicios profesionales, científicos y técnicos 588738.1
## 55 - Corporativos 216124.6
## 56 - Servicios de apoyo a los negocios y manejo de residuos, y servicios de remediación 1043187.0
## 61 - Servicios educativos 1019413.9
## 62 - Servicios de salud y de asistencia social 855054.7
## 71 - Servicios de esparcimiento culturales y deportivos, y otros servicios recreativos 179152.1
## 72 - Servicios de alojamiento temporal y de preparación de alimentos y bebidas 930483.9
## 81 - Otros servicios excepto actividades gubernamentales 614141.3
## 93 - Actividades legislativas, gubernamentales, de impartición de justicia y de organismos internacionales y extraterritoriales 1305982.0
#aL coeficiente de Leontiff / interrelacion total
aL_df <- data.frame(aL)
if ("52 - Servicios financieros y de seguros" %in% rownames(aL_df)) {
row_values <- aL_df["52 - Servicios financieros y de seguros", ]
top_5_values <- head(sort(row_values, decreasing = TRUE), 5)
result <- data.frame(
Valor = top_5_values,
Fila = rep("52 - Servicios financieros y de seguros", length(top_5_values)) # Añadir el nombre de la fila
)
print(result)
} else {
print("La fila '52 - Servicios financieros y de seguros' no existe en el dataframe.")
}
## Warning in xtfrm.data.frame(x): cannot xtfrm data frames
## Warning in data.frame(Valor = top_5_values, Fila = rep("52 - Servicios
## financieros y de seguros", : row names were found from a short variable and
## have been discarded
## Valor.X10 Valor.X13 Valor.X6 Valor.X20 Valor.X7 Valor.X2 Valor.X9
## 1 1.091199 0.1185505 0.02026959 0.01985867 0.01208949 0.01092235 0.008229832
## 2 1.091199 0.1185505 0.02026959 0.01985867 0.01208949 0.01092235 0.008229832
## 3 1.091199 0.1185505 0.02026959 0.01985867 0.01208949 0.01092235 0.008229832
## 4 1.091199 0.1185505 0.02026959 0.01985867 0.01208949 0.01092235 0.008229832
## 5 1.091199 0.1185505 0.02026959 0.01985867 0.01208949 0.01092235 0.008229832
## 6 1.091199 0.1185505 0.02026959 0.01985867 0.01208949 0.01092235 0.008229832
## 7 1.091199 0.1185505 0.02026959 0.01985867 0.01208949 0.01092235 0.008229832
## 8 1.091199 0.1185505 0.02026959 0.01985867 0.01208949 0.01092235 0.008229832
## 9 1.091199 0.1185505 0.02026959 0.01985867 0.01208949 0.01092235 0.008229832
## 10 1.091199 0.1185505 0.02026959 0.01985867 0.01208949 0.01092235 0.008229832
## 11 1.091199 0.1185505 0.02026959 0.01985867 0.01208949 0.01092235 0.008229832
## 12 1.091199 0.1185505 0.02026959 0.01985867 0.01208949 0.01092235 0.008229832
## 13 1.091199 0.1185505 0.02026959 0.01985867 0.01208949 0.01092235 0.008229832
## 14 1.091199 0.1185505 0.02026959 0.01985867 0.01208949 0.01092235 0.008229832
## 15 1.091199 0.1185505 0.02026959 0.01985867 0.01208949 0.01092235 0.008229832
## 16 1.091199 0.1185505 0.02026959 0.01985867 0.01208949 0.01092235 0.008229832
## 17 1.091199 0.1185505 0.02026959 0.01985867 0.01208949 0.01092235 0.008229832
## 18 1.091199 0.1185505 0.02026959 0.01985867 0.01208949 0.01092235 0.008229832
## 19 1.091199 0.1185505 0.02026959 0.01985867 0.01208949 0.01092235 0.008229832
## 20 1.091199 0.1185505 0.02026959 0.01985867 0.01208949 0.01092235 0.008229832
## Valor.X12 Valor.X16 Valor.X4 Valor.X1 Valor.X8 Valor.X5
## 1 0.006268929 0.006191874 0.005381136 0.00478314 0.00469103 0.004417929
## 2 0.006268929 0.006191874 0.005381136 0.00478314 0.00469103 0.004417929
## 3 0.006268929 0.006191874 0.005381136 0.00478314 0.00469103 0.004417929
## 4 0.006268929 0.006191874 0.005381136 0.00478314 0.00469103 0.004417929
## 5 0.006268929 0.006191874 0.005381136 0.00478314 0.00469103 0.004417929
## 6 0.006268929 0.006191874 0.005381136 0.00478314 0.00469103 0.004417929
## 7 0.006268929 0.006191874 0.005381136 0.00478314 0.00469103 0.004417929
## 8 0.006268929 0.006191874 0.005381136 0.00478314 0.00469103 0.004417929
## 9 0.006268929 0.006191874 0.005381136 0.00478314 0.00469103 0.004417929
## 10 0.006268929 0.006191874 0.005381136 0.00478314 0.00469103 0.004417929
## 11 0.006268929 0.006191874 0.005381136 0.00478314 0.00469103 0.004417929
## 12 0.006268929 0.006191874 0.005381136 0.00478314 0.00469103 0.004417929
## 13 0.006268929 0.006191874 0.005381136 0.00478314 0.00469103 0.004417929
## 14 0.006268929 0.006191874 0.005381136 0.00478314 0.00469103 0.004417929
## 15 0.006268929 0.006191874 0.005381136 0.00478314 0.00469103 0.004417929
## 16 0.006268929 0.006191874 0.005381136 0.00478314 0.00469103 0.004417929
## 17 0.006268929 0.006191874 0.005381136 0.00478314 0.00469103 0.004417929
## 18 0.006268929 0.006191874 0.005381136 0.00478314 0.00469103 0.004417929
## 19 0.006268929 0.006191874 0.005381136 0.00478314 0.00469103 0.004417929
## 20 0.006268929 0.006191874 0.005381136 0.00478314 0.00469103 0.004417929
## Valor.X18 Valor.X14 Valor.X3 Valor.X17 Valor.X15 Valor.X19
## 1 0.004265286 0.003651411 0.003062974 0.002434707 0.002180567 0.00156035
## 2 0.004265286 0.003651411 0.003062974 0.002434707 0.002180567 0.00156035
## 3 0.004265286 0.003651411 0.003062974 0.002434707 0.002180567 0.00156035
## 4 0.004265286 0.003651411 0.003062974 0.002434707 0.002180567 0.00156035
## 5 0.004265286 0.003651411 0.003062974 0.002434707 0.002180567 0.00156035
## 6 0.004265286 0.003651411 0.003062974 0.002434707 0.002180567 0.00156035
## 7 0.004265286 0.003651411 0.003062974 0.002434707 0.002180567 0.00156035
## 8 0.004265286 0.003651411 0.003062974 0.002434707 0.002180567 0.00156035
## 9 0.004265286 0.003651411 0.003062974 0.002434707 0.002180567 0.00156035
## 10 0.004265286 0.003651411 0.003062974 0.002434707 0.002180567 0.00156035
## 11 0.004265286 0.003651411 0.003062974 0.002434707 0.002180567 0.00156035
## 12 0.004265286 0.003651411 0.003062974 0.002434707 0.002180567 0.00156035
## 13 0.004265286 0.003651411 0.003062974 0.002434707 0.002180567 0.00156035
## 14 0.004265286 0.003651411 0.003062974 0.002434707 0.002180567 0.00156035
## 15 0.004265286 0.003651411 0.003062974 0.002434707 0.002180567 0.00156035
## 16 0.004265286 0.003651411 0.003062974 0.002434707 0.002180567 0.00156035
## 17 0.004265286 0.003651411 0.003062974 0.002434707 0.002180567 0.00156035
## 18 0.004265286 0.003651411 0.003062974 0.002434707 0.002180567 0.00156035
## 19 0.004265286 0.003651411 0.003062974 0.002434707 0.002180567 0.00156035
## 20 0.004265286 0.003651411 0.003062974 0.002434707 0.002180567 0.00156035
## Valor.X11 Fila
## 1 0.001296497 52 - Servicios financieros y de seguros
## 2 0.001296497 52 - Servicios financieros y de seguros
## 3 0.001296497 52 - Servicios financieros y de seguros
## 4 0.001296497 52 - Servicios financieros y de seguros
## 5 0.001296497 52 - Servicios financieros y de seguros
## 6 0.001296497 52 - Servicios financieros y de seguros
## 7 0.001296497 52 - Servicios financieros y de seguros
## 8 0.001296497 52 - Servicios financieros y de seguros
## 9 0.001296497 52 - Servicios financieros y de seguros
## 10 0.001296497 52 - Servicios financieros y de seguros
## 11 0.001296497 52 - Servicios financieros y de seguros
## 12 0.001296497 52 - Servicios financieros y de seguros
## 13 0.001296497 52 - Servicios financieros y de seguros
## 14 0.001296497 52 - Servicios financieros y de seguros
## 15 0.001296497 52 - Servicios financieros y de seguros
## 16 0.001296497 52 - Servicios financieros y de seguros
## 17 0.001296497 52 - Servicios financieros y de seguros
## 18 0.001296497 52 - Servicios financieros y de seguros
## 19 0.001296497 52 - Servicios financieros y de seguros
## 20 0.001296497 52 - Servicios financieros y de seguros