library(readr)
## Warning: package 'readr' was built under R version 4.4.3

2018

conc_2018 <- read_csv("concentradohogar.csv")
## Rows: 74647 Columns: 126
## ── Column specification ────────────────────────────────────────────────────────
## Delimiter: ","
## chr   (5): folioviv, ubica_geo, est_dis, upm, educa_jefe
## dbl (121): foliohog, tam_loc, est_socio, factor, clase_hog, sexo_jefe, edad_...
## 
## ℹ Use `spec()` to retrieve the full column specification for this data.
## ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.
# 2. Entidad federativa
conc_2018$entidad <- substr(conc_2018$folioviv, 1, 2)
inpc_base_2018   <- 100.000   # base de la ENIGH 2018
inpc_dic_2024    <- 135.38   # AVG de todo el año . 
factor_defl_2018      <- inpc_dic_2024 / inpc_base_2018
# Tabla de ingresos
ingresos_hogar_2018 <- data.frame(
  folioviv      = conc_2018$folioviv,
  foliohog      = conc_2018$foliohog,
  entidad       = conc_2018$entidad,
  factor        = conc_2018$factor,
  tot_integ     = conc_2018$tot_integ,
  ing_cor_2018  = (conc_2018$ing_cor/3),
  ing_cor_2024  = (conc_2018$ing_cor/3) * factor_defl_2018,
  ing_pc_2024   = ((conc_2018$ing_cor/3)*factor_defl_2018) / conc_2018$tot_integ)
ingresos_hogar_2018$pobre_extremo <- as.integer(ingresos_hogar_2018$ing_pc_2024 < 1514.91)
ingresos_hogar_2018$pobre_moderado <- as.integer(
  ingresos_hogar_2018$ing_pc_2024 >= 1514.91 & ingresos_hogar_2018$ing_pc_2024 < 3275.6
)
# Solo Estado de México (clave 15)
ingresos_edomex_2018 <- ingresos_hogar_2018[ingresos_hogar_2018$entidad == "15", ]

head(ingresos_edomex_2018)
##         folioviv foliohog entidad factor tot_integ ing_cor_2018 ing_cor_2024
## 34498 1500006701        1      15   2232         5    15043.263    20365.570
## 34499 1500006702        1      15   2232         4    15896.703    21520.957
## 34500 1500006703        1      15   2232         2     8213.110    11118.908
## 34501 1500006704        1      15   2232         2     4747.867     6427.662
## 34502 1500006705        1      15   2232         6    29265.563    39619.720
## 34503 1500006706        1      15   2232         5     9404.090    12731.257
##       ing_pc_2024 pobre_extremo pobre_moderado
## 34498    4073.114             0              0
## 34499    5380.239             0              0
## 34500    5559.454             0              0
## 34501    3213.831             0              1
## 34502    6603.287             0              0
## 34503    2546.251             0              1
vector_2018 <- paste0(
  "[",
  paste(ingresos_edomex_2018$ing_cor_2018, collapse = ","),
  "]"
)
vector_2018
## [1] "[15043.2633333333,15896.7033333333,8213.11,4747.86666666667,29265.5633333333,9404.09,11419.8866666667,9145.42,11299.1166666667,15838.3333333333,8834.42,14337.0133333333,10034.0033333333,12253.2333333333,12020.08,10651.25,13892.83,10701.96,14760.3733333333,15061.7333333333,14231.8466666667,16364.8266666667,9619.85,14920.9433333333,7048.51666666667,8095.10333333333,11874.74,7749.73666666667,40650.4166666667,74524.5866666667,6654.10666666667,17721.3033333333,11005.7133333333,39073.7633333333,10427.06,20812.5433333333,9592.51333333333,18652.6866666667,27364.6233333333,26316.2533333333,12138.8133333333,10980.71,150016.293333333,19512.6166666667,34502.3466666667,7785.05333333333,11500.1233333333,6332.39666666667,9430.22,60481.76,2367.79333333333,26295.0733333333,30111.4733333333,7685.24,18146.8033333333,16324.58,15152.45,5702.08666666667,11606.5533333333,12588.5233333333,12949.7333333333,30055.8966666667,64985.4333333333,22619.91,15903.5733333333,10754.1,25842.3833333333,10102.1633333333,20012.27,35714.7433333333,17736.63,23184.78,10759.5666666667,2705.8,58128.5366666667,6629.83333333333,12807.34,28721.41,9701.08333333333,13146.7266666667,20315.2133333333,9870.43666666667,20677.7066666667,12103.8733333333,37811.9533333333,24059.0033333333,51052.4566666667,18483.7833333333,22514.7533333333,5459.01333333333,42147.5366666667,13205.8066666667,22372.8066666667,5914.75333333333,30085.54,5333.60666666667,23470.09,4172.19666666667,8374.43666666667,21688.5233333333,21704.9133333333,12732.5766666667,9138.94666666667,17860.6533333333,14886.6466666667,14185.2366666667,23777.5066666667,39161.3233333333,20190.6,36349.17,13572.1533333333,6797.06,18443.31,10472.09,10621.5,21861.4666666667,11971.8533333333,23362.2766666667,42432.54,73158.6033333333,28638.1333333333,10701.34,17462.2633333333,11489.3966666667,7172.59333333333,14836.9533333333,27338.7033333333,27391.3,22499.9966666667,5831.84333333333,19780.4166666667,13304.5466666667,10597.7866666667,19431.8433333333,16659.0233333333,15031.05,34991.6066666667,30719.8666666667,18866.5366666667,20612.52,10271.7366666667,16008.5133333333,15722.28,14182.2633333333,52369.7133333333,10131.3666666667,10533.6533333333,4825.82,8543.47,6897.21,7988.94,3763.57,17765.3366666667,15942.32,14991.7966666667,14867.7633333333,8922.07666666667,17008.3733333333,7655.73666666667,13834.56,6405.66,13056.2766666667,5457.21666666667,13206.52,20897.7166666667,21552.1333333333,10036.7333333333,5287.23,5825.38,9659.01,12039.29,6320.45333333333,10299.6133333333,10200.2233333333,20844.2533333333,9204.91666666667,10289.9133333333,3997.76,27573.76,12040.98,9733.94666666667,6639.34,15151.2033333333,14764.1333333333,17862.5433333333,23098.48,3472.37,10939.0666666667,14364.61,9210.63666666667,18250.87,10897.61,19075.8833333333,22745.6066666667,18113.58,18072.63,8321.27333333333,8495.08666666667,11222.2766666667,28025.3133333333,3284.58333333333,19764.8333333333,16448.6366666667,21576.95,10745.0866666667,13058.7466666667,11167.4366666667,20292.6333333333,5965.42666666667,16331.6,15803.1133333333,7185.29333333333,16140.42,8599.43,3805.74666666667,8440.91333333333,6481.16,6147.25666666667,9084.14666666667,12140.1533333333,4671.71666666667,11726.5033333333,13571.5433333333,7388.39333333333,7030.13666666667,13967.21,4521.68333333333,1967.21,21827.8633333333,30482.9,10131.7766666667,10167.0433333333,15253.0166666667,6356.43,5018.35333333333,17103.78,3956.51333333333,18608.7,5955.06666666667,10618.3766666667,25622.24,6937.23333333333,24437.5766666667,13260.98,21957.2133333333,3398.03333333333,18632.8633333333,7043.47333333333,19495.2566666667,5852.73333333333,4289.61333333333,11315.7366666667,11316.2633333333,9375.52,8793.82333333333,5967.38666666667,14123.2433333333,6266.04,7619.98666666667,9566.75,12338.1466666667,18150.7933333333,14816.2866666667,4308.19333333333,6338.73,4131.14666666667,6319.66666666667,15429.4933333333,13799.7433333333,11045.08,9285.48,3016.82,8906.40333333333,17424.6133333333,12239.8966666667,8126.08333333333,15482.0233333333,11566.9833333333,10434.4233333333,4918.03,15819.6666666667,8081.96666666667,6314.75333333333,18108.1333333333,14356.9633333333,22718.2433333333,6849.17666666667,31659.6166666667,9027.87333333333,12120.92,8196.35,19307.8366666667,13898.5033333333,6767.35333333333,4448.87333333333,16483.5966666667,4234.21,6837.93666666667,15631.1366666667,19434.42,11663.08,5997.56,14789.3333333333,15557.37,13424.5866666667,9723.66333333333,7059.00666666667,44565.5633333333,23909.83,30573.7666666667,7654.72666666667,10266.6233333333,19137.59,12029.7433333333,10935.8166666667,7030.22666666667,21566.2033333333,10258.83,14810.19,4740.33,18175.46,9823.79,19837.0433333333,12049.39,7551.67,8315.21333333333,8606.88,14507.1166666667,19744.9366666667,16389.2766666667,10722.31,5611.49333333333,4400.64,4908.24,5696.52333333333,6824.65666666667,16422.94,19365.3533333333,7696.38333333333,8307.05666666667,8135.98666666667,10227.15,7287.68,3373.77,11556.4766666667,39357.37,12526.2366666667,9632.18,23244.05,10778.2,7243.07,11121.5666666667,8943.08666666667,15356.06,7051.63,17981.65,20162.4233333333,13086.3866666667,4842.90666666667,7805.04666666667,9245.72333333333,8783.30666666667,12815.21,14381.7266666667,19646.1266666667,10225.52,39021.7633333333,4914.22333333333,6311.47,9270.29,10942.62,6540.98,11363.1033333333,34728.3133333333,15466.2133333333,5045.93,11918.0266666667,13669.35,14424.6066666667,14760.28,10999.76,10435.5366666667,14779.3833333333,5546.56,4009.13,5096.68,9506.18333333333,23583.0866666667,10471.8633333333,7076.04666666667,17080.2333333333,14520.31,11429.4366666667,17529.7066666667,13327.8533333333,20385.8366666667,11642.26,14581.96,5663.93333333333,10147.5366666667,15940.2566666667,9196.06666666667,26120.6033333333,13185.48,12254.1333333333,5828.52666666667,6527.32333333333,10693.31,8793.82666666667,16147.53,13327.86,31999.9833333333,14648.3733333333,12827.8633333333,7770.48666666667,13327.3,50212.6466666667,20330.2366666667,6754.47,7980.52333333333,6645.18,5744.91333333333,8972.69333333333,6103.63333333333,7480.15333333333,11047.27,10176.5633333333,8090.16,9401.09666666667,14549.1633333333,14508.1933333333,16044.6966666667,7183.84,6327.86666666667,23903.1033333333,8074.74666666667,9574.03,8728.61666666667,4597.87333333333,4531.26666666667,5323.21666666667,14522.82,21191.1233333333,8107.38,17683.2066666667,13420.0733333333,6752.84666666667,13761.3533333333,14503.5233333333,23751.94,17647.08,8309.21,3130.43,5022.78333333333,14481.09,13197.13,15620.61,10432.04,5537.70333333333,4161.39333333333,8459.01333333333,9134.42333333333,8401.63333333333,5879.38333333333,15081.9566666667,14348.87,9513.05333333333,22551.01,23148.2566666667,12623.2966666667,8901.63666666667,18573.7666666667,14028.9766666667,21530.9866666667,7434.77666666667,20294.0833333333,8516.00666666667,12883.0033333333,10008.1933333333,12836.0566666667,8000.98333333333,27836.06,9940.31333333333,31426.4633333333,23641.3,27338.7033333333,15962.9966666667,12294.5966666667,32999.4566666667,8996.25333333333,21314.91,22617.35,10614.1233333333,14642.35,10097.8633333333,12582.2166666667,11499.31,27663.9233333333,7672.12666666667,16734.42,8196.72,8901.63666666667,32786.88,16309.8333333333,9631.14333333333,15803.2766666667,32704.91,40409.83,10691.9366666667,5220.99,8486.18666666667,14585.6266666667,17651.43,80349.3966666667,14839.2733333333,4890.56,2185.30333333333,9097.82,6756.03333333333,9030.23666666667,4121.76,6044.76333333333,29573.4533333333,16285.0633333333,6432.24333333333,4909.54666666667,4744.56,17611.7233333333,19660.8966666667,14084.23,11500.82,15848.7033333333,13286.8733333333,5516.38333333333,6901.86333333333,7424.50333333333,18598.5166666667,15854.6066666667,15655.5133333333,12005.25,12093.1266666667,12584.32,9427.81333333333,10385.4633333333,17191.4466666667,14694.0333333333,6462.97666666667,18148.7,14034.3966666667,9721.32666666667,9061.93,13256.58,9391.3,6772.61,28725.4666666667,10565.2,19049.7833333333,22983.6033333333,45417.0133333333,29016.38,22544.0933333333,62860.2933333333,10859.0233333333,24754.9366666667,52448.6566666667,47195.5866666667,12788.3866666667,48130.35,19408.8033333333,20533.21,13265.4,5247.63,32339.1033333333,31989.9133333333,19978.5,15642.35,8791.50666666667,11803.2766666667,12542.79,7919.62,15340.3966666667,32550.4933333333,10327.86,18658.4866666667,14327.8633333333,16877.13,11804.2133333333,3683.03,56518.79,18571.2533333333,8253.87,8843.60666666667,29295.2266666667,36307.4333333333,29347.8233333333,26677.8466666667,34133.9333333333,18953.37,28369.56,14276.3566666667,10567.3166666667,32266.3,28899.45,20117.2233333333,12293.9933333333,10840.7533333333,28979.6566666667,6855.97,22294.8733333333,8242.33,12381.5133333333,9865.75,18840.3133333333,11228.9533333333,8680.43,10923.91,21864.1333333333,6622.95,17540.9733333333,14726.6933333333,10099.87,12336.7733333333,11571.32,12370.26,8322.47333333333,26602.6366666667,73174.64,7986.49333333333,26474.05,12131.9833333333,4506.22666666667,7984.43666666667,21434.42,20279.08,25143.4933333333,8836.06333333333,7868.85,18647.5233333333,7177.46,4549.72,14449.29,18234.71,12610.85,7790.05333333333,8204.41666666667,3853.58666666667,23926.74,9823.10333333333,15597.96,14169,7712.84333333333,16773.4866666667,0,10549.43,13098.6966666667,21018.02,15910.41,24024.01,5266.30333333333,14184.7766666667,14827.4866666667,22091.07,11712.8266666667,11707.57,3505.95666666667,23745.89,8072.12333333333,9765.40666666667,26311.47,13450.81,15770.48,17677.0666666667,44785.3933333333,27780.2966666667,12446.5233333333,14722.82,5532.06333333333,16146.3333333333,11820.6466666667,21786.88,15106.5533333333,23327.8566666667,9991.79333333333,21704.9166666667,9829.5,18560.6533333333,19365.3533333333,11135.8666666667,12076.7866666667,22772.4366666667,9537.51666666667,7220.19333333333,4523.61333333333,16573.68,17522.96,11761.7433333333,16367.2033333333,8467.20666666667,8934.42333333333,7631.14333333333,5824.59,8606.55,36434.41,9836.05666666667,8852.45666666667,6713.11333333333,7139.33666666667,14190.5366666667,8196.72,19301.6333333333,8775.40666666667,48780.7833333333,25922.9833333333,21679.9766666667,7540.22666666667,14298.5166666667,13883.8866666667,6327.13333333333,9349.22333333333,6610.09666666667,32061.57,10218.09,20605,20590.16,24601.1366666667,8681.96666666667,14605.1033333333,28129.9533333333,27259.0033333333,12396.71,9180.39333333333,8288.91666666667,8617.31333333333,7611.49666666667,22549.9566666667,6195.65,12533.3066666667,13248.1766666667,8632.31,5521.30333333333,18490.01,13242.18,8360.65,35966.8433333333,6927.45666666667,17660.1733333333,24190.09,11183.28,5839.47333333333,14247.2533333333,34037.6366666667,10173.91,13057.5633333333,8726.28333333333,11331.5133333333,9691.29666666667,7441.49333333333,69926.3866666667,6826.78333333333,59826.22,52392.6933333333,9586.95333333333,30788.8333333333,40981.7633333333,13982.7566666667,13499.99,20424.5866666667,19262.29,34147.5366666667,16639.34,20204.9133333333,13048.0933333333,15918.0233333333,22492.21,15318.0266666667,6603.38666666667,13976.3866666667,13755.14,12992.78,10257.9266666667,22032.78,6794.39,12717.8533333333,9231.56666666667,6010.89,48999.8133333333,10291.4533333333,5054.86666666667,17630.5133333333,7695.64666666667,35639.89,9131.54666666667,17546.89,7326.06,61974.9166666667,51367.7733333333,27410.85,27143.6133333333,107136.286666667,9070.60333333333,36586.17,20816.3933333333,19944.25,19963.4566666667,15126.22,2664.13,17078.7966666667,20242.1666666667,20326.7733333333,139153.563333333,30474.6566666667,18058.1633333333,8874.78,10641.9933333333,10271.7333333333,29837.6266666667,13751.95,6066.26,59405.2366666667,29040.1366666667,32045.9266666667,28930.6533333333,6491.91666666667,5187.86333333333,6353.62,12038.6266666667,17634.9933333333,21055.0366666667,17560.0333333333,23960.4066666667,30087.2033333333,4531.87666666667,21504.2233333333,6834.79666666667,6929.34333333333,7774.87,4356.51666666667,6974.27666666667,4695.65,9988.02333333333,11778.2166666667,16237.28,8302.17,10930.2166666667,9509.11,14036.87,16537.4566666667,18753.7633333333,11608.8933333333,14262.2933333333,17146.36,25318.3466666667,17111.98,35282.88,25145.0666666667,5308.19,15232.5766666667,15008.72,21669.9666666667,7606.55666666667,16340.9766666667,7213.11333333333,12959.0033333333,7365.57,4677.97,12281.9633333333,13508.1933333333,14961.86,3097.55,2983.42333333333,27596.6766666667,12149.1666666667,15355.0033333333,24306.0933333333,10675.66,81513.3133333333,38482.86,20622.9433333333,19445.95,30235.3433333333,72557.3633333333,25736.9766666667,5508.19666666667,27459.0066666667,9754.09333333333,14434.42,30019.9433333333,11514.33,33966.2933333333,40698.6166666667,50018.2566666667,18633.2333333333,109191.786666667,13664.0833333333,81492.5466666667,9017.35,32116.9266666667,8355.10333333333,24515.7033333333,19623.1033333333,26556.2033333333,27728.1833333333,17217.39,26176.43,11773,57554.3366666667,26023.84,29125.1733333333,57338.6933333333,35129.0633333333,16308.01,8790.54,13278.6833333333,4286.88333333333,14848.6933333333,6706.49,15071.5566666667,30707.9133333333,9237.10333333333,19381.13,12283.0433333333,25188.89,8159.39666666667,8458.16,9432.95,9130.42666666667,7204.32333333333,11161.9366666667,8137.58,14394.4866666667,70704.76,9489.95666666667,19924.3933333333,10750.3466666667,31727.6766666667,21097.8233333333,10255.4333333333,12207.22,5457.22,10739.83,23915.0566666667,9063.05,10563.7366666667,7805.04666666667,15554.3466666667,6831.48333333333,33668.4666666667,22077.1966666667,7033.69,34411,5268.82666666667,48789.4366666667,9867.17,2689.42666666667,27064.63,5863.19333333333,9822.66333333333,20113.2733333333,37536.4166666667,21827.04,7897.61,15804.6133333333,25379.5933333333,10962.8233333333,17584.7,19017.1666666667,11799.6066666667,28551.1866666667,13510.19,9136.16666666667,9899.66666666667,11436.0566666667,12459.0133333333,11593.57,9385.24333333333,2459.01333333333,28748.9466666667,6479.80333333333,18636.0366666667,10924.96,14286.8466666667,31467.1966666667,16981.2733333333,5583.46,23098.3466666667,30332.5233333333,15516.1966666667,4372.98,13060.6633333333,19889.98,12335.28,19280.9866666667,11770.4833333333,42852.98,11213.11,12999.9933333333,14134.4166666667,26524.5766666667,16908.1233333333,32262.6133333333,18085.7933333333,12843.6133333333,49362.5366666667,6067.34,13782.42,32779.7,13686.7466666667,7736.04333333333,18318.7633333333,18465.42,7842.91,11941.7866666667,8593.38333333333,9997.14333333333,7902.33333333333,12622.9433333333,3167.1,7655.73333333333,9434.42333333333,18008.1933333333,36322.8733333333,18737.7,12008.19,7527.15,26950.8133333333,14754.0966666667,16190.1533333333,19282.63,13557.31,22813.1466666667,10826.8066666667,24294.1766666667,11535.91,14673.7533333333,7305.39666666667,15053.2133333333,11707.5666666667,5074.33333333333,21082.5666666667,11340.76,5826.85666666667,19161.0166666667,11467.7366666667,7086.66,10931.47,10076.8133333333,14855.9133333333,8630.44,9695.84333333333,9593.44,8797.2,10426.2233333333,5620.51,5105.70666666667,1880.99666666667,6855.03,8954.71666666667,8976.36333333333,4888.71333333333,23403.0333333333,11775.4166666667,11925.85,17343.0833333333,8625.52,11909.5533333333,18962.4733333333,9799.43333333333,5086.95333333333,13524.0033333333,11163.27,21961.2166666667,14869.56,11466.8533333333,9696.59,16876.5533333333,39467.2033333333,13221.69,14327.53,6434.64666666667,16093.44,40256.4766666667,15378.77,7262.47666666667,54239.8033333333,11522.4366666667,13069.6166666667,56116.09,4967.42333333333,13685.86,12728.4266666667,5343.16,28194.98,22447.4,10084.0533333333,4976.90666666667,13105.07,53637.1933333333,9024.58666666667,9661.37333333333,6604.31,4755.76,8411.46,28327.54,20051.2333333333,10278.6866666667,16219.9866666667,10796.85,10598.3066666667,16679.3366666667,15077.3066666667,7027.82333333333,17089.78,20129.49,25351.64,28424.82,18448.31,17079.1466666667,36024.5866666667,19106.5433333333,42352.45,16270.49,13197.1333333333,13481.3466666667,4292.50666666667,8624.99,20881.3933333333,8023.65333333333,16930.74,13887.3,11103.2533333333,12938.4666666667,8778.04666666667,14124.5766666667,22011.7,5148.70666666667,8003.25666666667,9918.59666666667,14777.0433333333,16179.1766666667,7209.83333333333,11892.9166666667,17760.6533333333,3888.52333333333,19855.73,6158.25666666667,38849.7266666667,5873.81666666667,6691.85,7749.95,8780.87333333333,11590.28,29488.5066666667,34185.47,6918.82666666667,5581.07333333333,39851.1533333333,23343.6166666667,19517.8733333333,32052.7266666667,47441.18,20398.5666666667,50613.0866666667,4839.21,12095.7,5945.99666666667,8942.98,13501.3233333333,13909.3866666667,12572.83,27258.43,5763.92666666667,3467.44333333333,15437.6133333333,18481.29,20017.5133333333,21035.9333333333,9769.71666666667,16672.85,18644.5266666667,23119.9666666667,18038.9166666667,41446.69,6510.31666666667,5176.62666666667,29771.05,13551.8566666667,6016.47666666667,7007.65666666667,9348.17,6092.69666666667,48655.32,16630.4333333333,16518.4466666667,4320.4,15993.44,9602.48,13428.7633333333,7401.63666666667,3681.05,8103.78,7258.84666666667,17547.1533333333,21995.08,12225.99,5703.27666666667,17857.9833333333,4147.54,21945.61,11555.73,12643.7566666667,9984.99666666667,49306.4466666667,23854.3,15680.04,7063.16333333333,21244.5566666667,32450.1133333333,16604.1333333333,38200.5,7223.49666666667,10679.8666666667,4978.60666666667,11178.91,3601.58333333333,31538.9066666667,27679.3,49809.7766666667,13336.0566666667,19475.3966666667,26967.3433333333,13634.4166666667,35557.36,12954.0933333333,9229.05,9582.93666666667,8811.47,10754.0933333333,4459.01333333333,6926.22666666667,12309.8533333333,11500.0066666667,12409.83,7049.17333333333,6276.22666666667,10814.92,7962.83,9352.62333333333,9391.3,12370.26,11399.5666666667,6225.09,14643.3766666667,3956.51333333333,19683.4366666667,17369.7266666667,13287.3833333333,12533.78,20654.8433333333,13466.3333333333,4132.88,20006.4733333333,10250.3233333333,18327.13,7885.87666666667,21528.04,16662.8966666667,6467.08,16662.6766666667,7868.85,8578.21666666667,4557.37333333333,5470.26,18805.22,24343.4333333333,7071.34666666667,16779.2866666667,9416.73,46931.9666666667,17264.2466666667,22981.4966666667,51215.8933333333,8103.78,9444.25333333333,7098.17333333333,6323.56,9389.79333333333,10926.1566666667,14516.7266666667,23216.2233333333,29835.6333333333,5623.15333333333,7060.23666666667,21518.4466666667,6397.46,8903.92333333333,39236.08,14694.3033333333,30437.7033333333,10596.3233333333,9752.96333333333,17022.7866666667,13212.7066666667,4948.52333333333,6631.13333333333,7081.86666666667,20886.9133333333,25236.7466666667,16232.1833333333,14482.0333333333,16131.7,12227.1566666667,13876.5233333333,15892.99,9848.58,8986.51666666667,22689.5266666667,10566.3866666667,17145.33,16300.02,17309.7433333333,9210.37666666667,19133.3,7219.78333333333,6486.75666666667,13244.7666666667,9082.04,6478.60666666667,5570.82333333333,11076.06,2527.69666666667,2837.74,6799.34,3853.97333333333,4828.13,5923.89,36959.2266666667,2168.08666666667,11974.09,34729.5033333333,6485.93666666667,7766.65,25621.62,7527.34666666667,16359.0333333333,8386.21666666667,4133.67,5023.34333333333,6897.26,11788.56,6903.82666666667,13307.76,24611.4033333333,10900.4266666667,22563.36,12686.5266666667,27311.01,14743.44,12012.6566666667,16213.6966666667,8624.99666666667,10341.3466666667,6468.96,14264.7966666667,11220.4333333333,22703.5666666667,2998.43666666667,12304.7533333333,4144.03333333333,5325.70666666667,12084.5433333333,8788.47333333333,9918.01333333333,4853.10666666667,6854.91333333333,8273.76333333333,8051.63666666667,5324.57666666667,7849,7780.32,3437.7,6901.63666666667,6416.68666666667,6078.68333333333,3893.44,15557.3666666667,32426.2166666667,8790.48666666667,13324.58,5330.55,4399.79666666667,13529.1733333333,10430.04,6578.70666666667,5573.70666666667,11350.7033333333,8373.44333333333,20393.5566666667,6858.02333333333,13557.2933333333,16610.5266666667,14577.0666666667,9925.87666666667,9784.69,6192.53,8416.40666666667,22066.6933333333,7268.02,2183.60333333333,9546.45666666667,12662.94,12851.39,19027.3166666667,12012.0533333333,9048.11333333333,6206.32333333333,16062.41,19355.0866666667,4640.37333333333,13028.21,41015.9366666667,13768.42,7929.69333333333,12984.3533333333,4385.08666666667,5143.75333333333,18265.5666666667,21475.4033333333,11868.8433333333,11951.73,24329.4966666667,19754.0933333333,16413.74,5451.43333333333,11511.9166666667,13455.6566666667,15443.9933333333,10245.44,11991.7966666667,10025.5633333333,6147.53666666667,11603.2733333333,8199.99333333333,5114.75333333333,7774.43333333333,14016.09,5818.50666666667,23836.4966666667,14511.3833333333,21897.11,11061.18,15627.86,8732.12666666667,11693.6533333333,13114.7466666667,4888.52,35821.2966666667,19977.0333333333,14721.2966666667,52524.5733333333,9892.9,17940.4066666667,7993.64,10890.1533333333,16113.6933333333,14041.17,11347.8233333333,10470.8033333333,10108.1666666667,14039.17,15625.34,10105.4866666667,11787.7233333333,5987.79,7291.20333333333,10590.1633333333,6491.46,26662.03,14710.4266666667,30000.0266666667,15630.1,11029.5033333333,7847.53666666667,10232.4466666667,28393.1733333333,7096.72,9459.01333333333,5644.44666666667,7272.49333333333,14966.8433333333,5558.36666666667,5257.35,28971.5,23486.0566666667,22093.2633333333,19894.8433333333,10649.48,23229.2033333333,33827.8333333333,6498.84333333333,9284.95,19574.1433333333,4065.75333333333,18508.1933333333,28046.2333333333,15881.3733333333,49934.3933333333,41201.17,6557.37333333333,3893.63,16892.9233333333,13823.97,8990.72666666667,7278.91,11421.4566666667,15009.9266666667,583.796666666667,4432.37666666667,12686.0333333333,33507.36,6442.63,21629.96,12906.9233333333,11626.2266666667,11753.6,10803.42,6231.15666666667,6969.15,16701.4266666667,3921.22,23177.68,10591.1533333333,3292.65,34430.2833333333,2257.34666666667,10021.4533333333,9663.45333333333,9514.88,23900.12,11674.7066666667,6325.69666666667,2749.74,10517.7333333333,10891.1633333333,7025.40666666667,8153.46666666667,4884.7,9909.64,20561.3433333333,9464.02666666667,7669.34666666667,16508.0966666667,9163.56333333333,6095.45666666667,5163.21666666667,11142.81,13180.22,17837.99,10367.7166666667,26510.8966666667,13909.93,10266.9966666667,27020.24,7597.60333333333,5641.28666666667,19487.5966666667,11336.8333333333,11739.1266666667,18547.28,2421.71666666667,3700.87,13430.0433333333,113579.56,185553.99,6358.69333333333,18595.3333333333,7899.62666666667,13883.8266666667,12585.65,9532.91333333333,76393.4333333333,13495.79,13821.8766666667,21963.53,11120.6133333333,10839.0033333333,13701.1366666667,25507.52,10365.2366666667,7016.35666666667,11783.2966666667,20047.9566666667,15319.9733333333,9193.39333333333,11845.8566666667,6470.64,8738.59333333333,16547.8533333333,7924.62333333333,9785.64666666667,8025.23666666667,16220.97,15322.47,10035.58,27896.94,19278.68,20944.26,19675.9,14925.1866666667,17082.7466666667,3204.36333333333,8478.78,4960.48,10401.96,9277.69666666667,9204.32333333333,10882.6466666667,15720.22,10747.3266666667,10548.8666666667,16310.65,9383.57333333333,4263.20333333333,6802.49666666667,13012.97,9464.4,6481.23333333333,10986.0733333333,23684.6666666667,10286.4533333333,7973.06333333333,15063.2633333333,20239.3366666667,41448.0366666667,8984.58666666667,12792.8833333333,6215.39666666667,8651.93,13224.86,4419,8029.61,5129.34,8398.51,4990.17333333333,8535.07666666667,9847.82,12401.78,5097.96666666667,9691.93333333333,16027.2066666667,16471.73,8734.01666666667,12997.98,16430.0633333333,2032.85333333333,20971.07,5663.39,5735.36,7721.03333333333,14624.01,16712.2533333333,16978.6033333333,12351.6633333333,2770.23666666667,12873.9033333333,6346.28666666667,6689.97,5142.33666666667,7073.76666666667,20258.72,5835.07666666667,10469.7466666667,5361.82333333333,5270.48666666667,2774.93333333333,8211.49333333333,18557.3666666667,5579.73333333333,7126.77666666667,3641.23666666667,12997.79,38918.9933333333,4587.94666666667,13723.52,18408.65,9766.26333333333,11815.3533333333,9878.4,16696.6933333333,9121.3,24665.57,10942.62,4445.89666666667,8215.91333333333,22540.9633333333,9824.81666666667,6933.45,12166.0266666667,8032.80666666667,4730.16666666667,2970.05333333333,2867.29,3006.41333333333,7910.62,15326.0833333333,10424.62,7434.03333333333,14244.7266666667,11532.9533333333,2898.16666666667,5920.76,11900.03,7790.15333333333,8480.31666666667,6114.75333333333,300,10613.8133333333,4678.60666666667,5309.54,14122.9666666667,25806.4833333333,5646.5,7924.95333333333,17476.4,1286.98333333333,11314.8233333333,5236.43,13740.1333333333,7901.63666666667,7447.92666666667,10042.3733333333,12399.4066666667,7820.75333333333,12450.81,8550.97333333333,14153.2933333333,7035.47666666667,14983.1666666667,13833.12,30347.9933333333,8729.99333333333,4435.03333333333,13909.5466666667,1489.83666666667,3749.92666666667,6992.33,7128.14,5087.99666666667,6409.25,5106.32333333333,3345.22666666667,11511.1166666667,5180.12666666667,15203.0233333333,17041.8966666667,10605.37,21672.1266666667,6288.28666666667,9314.75,17492.9466666667,6849.17666666667,11498.92,8319.42,10620.38,11425.2866666667,23660.3,7450.43666666667,9397.08,22593.3466666667,9885.25666666667,10113.3133333333,3682.41,8005.95333333333,6646.11333333333,8831.77,12413.5566666667,8648.73333333333,10817.31,23756.3266666667,5547.54,3454.1,9277.69333333333,9361.21333333333,7491.81666666667,5533.42333333333,5856.93333333333,9641.21,5753.38333333333,4469.51666666667,2222.06333333333,18111.5866666667,21191.4033333333,3525.87666666667,32695.28,15681.4166666667,11768.27,11965.06,4347.16666666667,19575.7033333333,11340.15,10811.56,5400.20666666667,4611.79,9135.26666666667,9986.06333333333,9891.28333333333,5398.98333333333,3860.21666666667,7027.83666666667,6095.08,22322.4233333333,8233.82333333333,10006.38,4219.31,6146.56666666667,15420.75,22734.5533333333,9025.94666666667,14030.41,8974.98333333333,14457.0333333333,17652.45,3509.83,20032.7766666667,22368.8433333333,7526.22,9180.59333333333,4918.22666666667,36808.2566666667,41919.5566666667,25011.9166666667,9981.10333333333,17790.1366666667,19592.13,7821.34666666667,9553.92333333333,4105.89333333333,8713.11333333333,12196.71,15100.79,20458.1833333333,8762.99,14831.6933333333,26550.3033333333,9759.93333333333,5786.17333333333,7529.07666666667,3625.71666666667,10717.7366666667,6461.45333333333,21535.2966666667,3690.23333333333,9109.18666666667,20903.1633333333,8325.65,4074.55333333333,5297.85333333333,6226.06666666667,3649.48,9486.56333333333,1790.48333333333,3301.25333333333,5777.28333333333,3129.45,3249.08666666667,9724.60666666667,4400.22,8615.59333333333,17942.5633333333,11574.7966666667,5131.13,6797.42666666667,7312.12666666667,2736.84,5501.42333333333,8263.96666666667,11160.4366666667,6499.01,10540.2766666667,14768.4833333333,6508.41,4406.55333333333,10285.7433333333,8099.4,17340.65,12084.0633333333,5914.21,6688.51666666667,7988.21,13458.3666666667,10019.42,5087.48333333333,7789.37333333333,32031.7133333333,4461.77,8573.09333333333,11880.2333333333,16950.24,5312.15,7640.14,8695.20666666667,7060.21666666667,34623.3233333333,9672.12,5526.04333333333,20708.4166666667,2693.42333333333,4316.38333333333,10938.35,13221.3,17042.7533333333,7823.92333333333,10557.3666666667,18700.17,14434.2,6591.66333333333,10049.2166666667,8630.52666666667,10090.2766666667,15579.6833333333,8354.34,2904.69333333333,5636.04,7768.98666666667,7365.97333333333,4230.38,2012.51333333333,3363.68333333333,9556.85666666667,6883.55666666667,5341.01666666667,2675.39,4799.25333333333,8075.95333333333,7285.67,4822.65666666667,3185.12666666667,3511.21333333333,1801.21333333333,1959.41,6349.74666666667,3898.05333333333,5434.42333333333,2996.71666666667,23409.83,3289.83,12093.4333333333,27596.7133333333,10886.29,39166.9033333333,10613.4666666667,8363.38333333333,4781.17666666667,7948.1,7470.23666666667,4489.47666666667,2748.63333333333,23745.9633333333,8785.40666666667,9959.3,10962.4,5663.38666666667,7795.08,42117.3766666667,5745.9,10188.5233333333,6326.08,8695.99666666667,16237.2233333333,8305.9,11647.0833333333,7761.77666666667,7453.98666666667,6105.62333333333,5487.26666666667,21536.7933333333,17044.12,4577.88333333333,2038.56333333333,10072.9233333333,3614.30333333333,2340.56333333333,16742.3466666667,8072.65666666667,6854.81,12258.23,5710.43666666667,8310.48,10150.1566666667,21537.0466666667,15976.1033333333,4393.53,8045.89,4247.06666666667,5030.30333333333,8216.42333333333,7691.18,5068.01666666667,6056.02333333333,15031.8666666667,19337.71,12287.9066666667,16405.0833333333,10508.2833333333,8401.63333333333,24081.96,9065.57,15619.5466666667,11819.6666666667,7265.57,9645.89,16762.2866666667,9836.06333333333,7532.78666666667,3123.18,6235.59,6024.16666666667,5718.09,12793.2333333333,6256.34,8484.80333333333,12467.6633333333,14814.2833333333,5281.02666666667,25484.74,1894.99666666667,7932.21666666667,14012.4566666667,3704.76333333333,19523.14,10333.8966666667,21180.3133333333,6247.02333333333,7116.03333333333,36983.61,14280.4933333333,12856.0466666667,9982.03666666667,7562.58,10077.87,6390.64333333333,4712.47333333333,7776.45333333333,12711.3633333333,9197.16333333333,9761.57,4113.75666666667,9350.38,4509.98666666667,37233.7333333333,7970.68,13233.6933333333,5508.19333333333,11791.32,10983.6033333333,5721.31,3642.84666666667,15092.6733333333,30392.77,7902.33333333333,13633.92,18466.75,7278.91333333333,2187.14666666667,8784.82666666667,5715.46,12079.1233333333,11702.9833333333,5760.16333333333,11126.2833333333,31794.19,13130.0666666667,16208.8766666667,15392.51,7557.37666666667,12372.1166666667,14872.1233333333,5540.98333333333,16081.96,4166.40333333333,21795.07,9478.60666666667,27609.3733333333,5324.4,17962.6533333333,10060.2966666667,5362.13,4453.71,5308.09,8521.30333333333,5169.53,11437.6966666667,16713.1066666667,12801.6366666667,1926.22333333333,5804.58666666667,24165.5966666667,2285.12666666667,10373.8633333333,4124.38333333333,6859.91333333333,7537.64666666667,7483.13,8809.44,5684.57,5880.82333333333,9429.04666666667,6459.61333333333,5775.81666666667,6806.31333333333,8581.96333333333,6065.57,9464.22333333333,8509.57333333333,9231.3,11080.62,24803.45,17280.23,10626.8766666667,10117.9533333333,11921.3966666667,9201.88333333333,6020.26666666667,9224.80333333333,6351.25333333333,10244.1133333333,12255.51,7989.75,3874.34333333333,3665.69333333333,3559.63666666667,6390.29333333333,10487.29,11850.8033333333,25594.2066666667,11123.0266666667,16779.7366666667,8734.49,6981.88666666667,7270.48,6323.70333333333,15761.97,4799.25666666667,26615.2166666667,5808.16,4848.86666666667,2102.52,10277.51,5651.28,17908.7866666667,11063.2166666667,5101.31,13147.0766666667,3584.39333333333,6127.86666666667,11395.77,6383.27,8039.99333333333,11050.4866666667,15500.6866666667,7663.56,8025.51666666667,5963.71333333333,11198.9533333333,10534.9433333333,7099.60666666667,4266.67666666667,4102.42333333333,4861.31666666667,9761.57,17372.2566666667,7151.73,6397.08666666667,21412.6866666667,8133.95333333333,3802.59333333333,9815.91333333333,10281.0233333333,13281.4,7954.56333333333,8188.52333333333,3860.65,10795.0766666667,12808.1933333333,15997.0466666667,30582.1433333333,16254.2,12555.0833333333,7009.97666666667,8167.51666666667,7629.06,13013.49,6669,11913.74,4674.61333333333,11913.74,6273.74666666667,7762.78333333333,13190.2133333333,13724.9833333333,13541.8933333333,25351.78,16053.04,12228.2566666667,22096.58,7451.60333333333,7336.95333333333,7411.79333333333,17655.9733333333,76088.6933333333,8587.65333333333,12478.4533333333,5074.33,5668.8,17317.7333333333,20906.8166666667,4951.78666666667,15982.37,2755.95666666667,5632.78333333333,25950.8033333333,4205.38,28654.3833333333,4990.80666666667,10888.0466666667,10487.8533333333,11174.5833333333,9766.03333333333,10788.0533333333,8605.52,7314.53,12235.7566666667,26439.1233333333,10862.1433333333,7730.17,20945.88,26258.6766666667,13100.3333333333,6212.47666666667,11812.7866666667,7842.25666666667,10701.5233333333,17356.8433333333,2924.26,5255.53666666667,4576.18,10492.6666666667,7443.24333333333,7948.89,5570.08333333333,19398.0233333333,9563.92,12212.9366666667,11055.5566666667,9446.97666666667,12911.0633333333,13457.49,4003.02666666667,19291.2066666667,21239.4333333333,6629.83333333333,10912.49,21798.2466666667,15688.64,7258.74333333333,3995.93,6446,8156.61666666667,7445.13666666667,5670.69666666667,7190.21333333333,15022.8966666667,17426.2166666667,8114.75333333333,12567.9066666667,6410.77666666667,7368.20666666667,8793.82666666667,16455.2433333333,20010.4733333333,14403.04,5951.32333333333,10010.3633333333,23946.01,12331.79,9170.66,3039.94333333333,7294.27,7663.69,22568.8466666667,11499.1033333333,9817.51333333333,15091,11914.9566666667,12861.2866666667,16628.6133333333,18340.8066666667,6558.91333333333,12166.72,5568.72,3609.04333333333,18113.6033333333,12040.9,5750.68,3077.05666666667,10439.1033333333,7610.82666666667,8960.33,18357.8233333333,4844.26,18794.9766666667,8581.83333333333,6811.47333333333,14809.1533333333,9214.37,6599.34666666667,24858.6733333333,5437.66,13533.8766666667,12710.8266666667,12722.3233333333,9243.34333333333,25799.7733333333,21670.6866666667,25284.6433333333,31930.74,15221.2133333333,11685.0066666667,26762.2866666667,26596.4633333333,21648.3466666667,12440.74,24647.5233333333,9959.01,19696.7133333333,28844.2533333333,9893.51666666667,5648.39666666667,7158.18,10185.8133333333,4139.61,5221.17,7981.49666666667,2817.58666666667,16492.03,3721.98333333333,2244.08333333333,3548.24,7820.82,4446.10333333333,2774.63,4960.99333333333,3628.63333333333,4377.04666666667,4480.65,11426.22,19385.9933333333,12890.48,3876.14666666667,26001.7466666667,18026.5,28642.5233333333,10697.43,13673.6733333333,15297.4766666667,16691.1766666667,22977.9366666667,6813.88,16257.8633333333,12293.44,8502.61666666667,14508.1933333333,13454.3733333333,17403.9666666667,13700.6466666667,19549.9566666667,8779.37666666667,13917.4366666667,6064.01,8294.17666666667,10799.68,7782.95,6968.79,2014.46333333333,12770.2833333333,3838.56,4230.22333333333,16047.5366666667,5213.11333333333,7081.96666666667,7049.17333333333,11749.1733333333,8046.47666666667,13719.6866666667,12122.9466666667,17960.32,13144.7166666667,10003.23,11275.8733333333,10776.16,19269.9133333333,13658.3633333333,3638.65,4110.06666666667,1480.64,6143.40666666667,3418.38666666667,2842.00333333333,2388.07,3332.55666666667,3089.28,2307.29,10121.2333333333,5194.50333333333,8236.46,6160.10333333333,6451.39333333333,1942.44,7382.94333333333,5532.15,5885.32666666667,8092.52,4080.98333333333,2678.15333333333,3958.35333333333,11866.5133333333,21735.07,14155.5833333333,19024.7566666667,6856.85,1567.84333333333,14367.97,13385.3066666667,11271.78,9890.21,27229.5,9652.75,19538.8433333333,7633.24,7910.52666666667,10534.42,16551.1633333333,27253.8033333333,10203.27,10050.01,13777.74,10497.5366666667,4692.49,11101.21,13031.3166666667,9547.87333333333,6869.71666666667,2190.34333333333,1362.51666666667,2847.43,2263.54333333333,2507.65666666667,3323.44666666667,16524.58,319896.293333333,4655.73333333333,15172.7833333333,8652.48,3241.04333333333,6174.04,16895.56,17777.2266666667,5609.2,25713.0033333333,4614.79666666667,24563.11,20649.0433333333,10547.7366666667,4476.31,49837.4666666667,18879.7833333333,5873.31666666667,12599.88,29723.5333333333,9445.54666666667,4699.32333333333,7345.65666666667,16515.2333333333,13444.1866666667,7049.86666666667,8963.59666666667,4839.99666666667,10188.5166666667,21245.8866666667,10745.0833333333,13157.3966666667,24588.0033333333,17825.7666666667,18580.2166666667,42663.7166666667,31082.8366666667,9185.13,6411.28666666667,5180.28,8464.58333333333,6652.17,15449.15,13437.93,17404.44,9639.33333333333,22704.5833333333,14882.2933333333,16171.56,20794.62,4902.87333333333,3353.06333333333,14028.0433333333,16652.83,15574.0633333333,13507.5566666667,11092.2633333333,7641.62666666667,4374.29,5810.32666666667,4366.83666666667,5789.02666666667,8791.79666666667,5874.09666666667,13644.26,8942.62,6254.09666666667,7274.46333333333,5379.38,8196.01666666667,10037.34,14125.16,9590.16333333333,25508.8833333333,3550.01666666667,13217.7333333333,3350.04,2976.14,14663.39,6594.79333333333,8612.84333333333,13932.73,7979.08,4631.14333333333,11987.9366666667,6390.15666666667,9826.22333333333,13182.7466666667,8334.33666666667,6393.47666666667,5210.44666666667,3807.02,4335.83333333333,18997.0833333333,4004.55333333333,28897.7466666667,12225.14,8407.86,45087.1333333333,11283.13,13477.69,4456.86333333333,10253.85,14344.26,10843.8066666667,6150.81666666667,25675.49,6972.13,5486.87666666667,9062.51666666667,8940.39,8228.85,5952.68666666667,8030.54666666667,8648.18,8683.92,14346.6366666667,7183.20666666667,6410.62,1357.61333333333,24012.61,21014,89472.2233333333,18358.6933333333,42102.0233333333,46553.1066666667,36067.3033333333,11800.14,38241.5766666667,59301.7966666667,91321.1,10112.0633333333,24382.8766666667,94887.6166666667,7970.92666666667,5538.34333333333,6777.00666666667,11218.8166666667,7669.55333333333,31452.13,7564.47333333333,7579.93666666667,13954.1433333333,16232.2633333333,23382.44,5517,16840.3766666667,7690.16333333333,9942.94333333333,15967.2066666667,7934.42333333333,21785.24,13795.0733333333,11732.4533333333,21324.9666666667,7191.22666666667,10654.2466666667,10547.5766666667,8381.35666666667,13984.5366666667,1818.39,8631.70666666667,12703.6733333333,3187.68,5654.97666666667,5017.73666666667,2076.38666666667,14508.79,8374.32666666667,5780.49,5830.63333333333,7157.34,7753.23333333333,3506.16,8934.03666666667,11108.88,9473.87,15226.7133333333,18130.8433333333,16614.65,4434.22,3481.5,6994.82,5658.46,13085.5666666667,14996.21,8687.31333333333,2175.19333333333,21728.4266666667,9457.36333333333,7307.20333333333,4364.27666666667,4016.38666666667,2336.38666666667,11755.72,7237.7,13824.4,13147.5366666667,4058.19333333333,4545.41333333333,3891.3,4465.74,8296.69333333333,3648.31,9688.46,4528.37333333333,8587.43,5727.86333333333,9468.07666666667,6939.09333333333,9877.79666666667,6799.75,10518.34,8725.72333333333,6347.74333333333,4800.81,2610.48333333333,8310.28,1663.87,10048.84,5914.76666666667,14278.4633333333,14562.0933333333,6865.08,12888.07,28103.95,7022.43333333333,28431.3433333333,21445.2566666667,12202.3633333333,14192.9233333333,10443.0866666667,8930.04333333333,14124.9766666667,11424.0433333333,14888.2633333333,7267.20333333333,14416.7166666667,6459.97,7207.87666666667,8630.78,9677.52,5628.67666666667,3014.09333333333,2106.97666666667,8127.03666666667,14612.11,5179.52,4882.88666666667,3543.82333333333,2241.58333333333,1959.01,5721.31,3547.54,3881.96,4881.84333333333,4791.58,23455.66,1922.66,5031.80666666667,8178.33,16281.72,11327.0433333333,10348.52,9239.82666666667,4875.52,15726.06,8045.79666666667,10628.69,18332.4933333333,5146.36,12254.4133333333,9948.5,23189.2433333333,9792.10666666667,46586.68,11464.2566666667,16204.9133333333,12485.94,12092.7166666667,15573.76,10418.63,16345.03,8359.75,15098.8633333333,8039.61333333333,11052.3633333333,14091.0633333333,8905.91666666667,7653.71333333333,11741.22,16949.7633333333,10998.57,24402.61,12527.4133333333,16355.3666666667,11929.1433333333,7916.15,7407.76333333333,12089.24,5609.21333333333,16194.3466666667,5344.26,5255.26666666667,10542.15,12963.2566666667,14858.2966666667,7047.29,8671.41,6581.96333333333,7360.65,5087.41,1346.65666666667,8686.8,3629.19,2515.56333333333,15085.1666666667,7237.76333333333,7963.12,8350.55,10405.97,8353.23333333333,6803.26666666667,6491.79666666667,5052.13666666667,2039.44,3660.23333333333,5055.34,2880.08,4000.34,7309.93666666667,6494.87666666667,23185.8133333333,4272.35333333333,3153.92666666667,16546.0766666667,5944.7,11625.2266666667,5847.67,5083.2,3252.46,15563.8466666667,4289.28333333333,9831.82,7647.53666666667,8619.00666666667,11455.7366666667,6820.72,11065.5733333333,21469.1633333333,7032.78666666667,18411.78,14094.0966666667,11717.2633333333,17143.94,12687.9,5235.96333333333,7314.61666666667,13723.7533333333,11432.7166666667,3123.68,31354.3,8064.65,15483.0066666667,8635.43666666667,21419.43,8380.23666666667,2718.16666666667,12544.2566666667,20053.1033333333,17311.4633333333,8688.52333333333,7549.17666666667,7401.63666666667,12086.8733333333,9199.99,6819.67,7868.85,10424.58,4362.56666666667,6257.85,35593.4266666667,26590.54,10729.31,14284.9966666667,10956.52,31418.7366666667,12091.4266666667,8504.86666666667,14642.6166666667,10889.92,22537.74,14377.0433333333,12375.52,1540.02,4270.26333333333,14150.33,10245.44,4045.26]"
vector_2018_2024precios <- paste0(
  "[",
  paste(ingresos_edomex_2018$ing_cor_2024, collapse = ","),
  "]"
)

vector_2018_pc <- paste0(
  "[",
  paste(ingresos_edomex_2018$ing_pc_2024, collapse = ","),
  "]"
)
mean(ingresos_edomex_2018$ing_pc_2024)
## [1] 5765.114
names(conc_2018)
##   [1] "folioviv"   "foliohog"   "ubica_geo"  "tam_loc"    "est_socio" 
##   [6] "est_dis"    "upm"        "factor"     "clase_hog"  "sexo_jefe" 
##  [11] "edad_jefe"  "educa_jefe" "tot_integ"  "hombres"    "mujeres"   
##  [16] "mayores"    "menores"    "p12_64"     "p65mas"     "ocupados"  
##  [21] "percep_ing" "perc_ocupa" "ing_cor"    "ingtrab"    "trabajo"   
##  [26] "sueldos"    "horas_extr" "comisiones" "aguinaldo"  "indemtrab" 
##  [31] "otra_rem"   "remu_espec" "negocio"    "noagrop"    "industria" 
##  [36] "comercio"   "servicios"  "agrope"     "agricolas"  "pecuarios" 
##  [41] "reproducc"  "pesca"      "otros_trab" "rentas"     "utilidad"  
##  [46] "arrenda"    "transfer"   "jubilacion" "becas"      "donativos" 
##  [51] "remesas"    "bene_gob"   "transf_hog" "trans_inst" "estim_alqu"
##  [56] "otros_ing"  "gasto_mon"  "alimentos"  "ali_dentro" "cereales"  
##  [61] "carnes"     "pescado"    "leche"      "huevo"      "aceites"   
##  [66] "tuberculo"  "verduras"   "frutas"     "azucar"     "cafe"      
##  [71] "especias"   "otros_alim" "bebidas"    "ali_fuera"  "tabaco"    
##  [76] "vesti_calz" "vestido"    "calzado"    "vivienda"   "alquiler"  
##  [81] "pred_cons"  "agua"       "energia"    "limpieza"   "cuidados"  
##  [86] "utensilios" "enseres"    "salud"      "atenc_ambu" "hospital"  
##  [91] "medicinas"  "transporte" "publico"    "foraneo"    "adqui_vehi"
##  [96] "mantenim"   "refaccion"  "combus"     "comunica"   "educa_espa"
## [101] "educacion"  "esparci"    "paq_turist" "personales" "cuida_pers"
## [106] "acces_pers" "otros_gas"  "transf_gas" "percep_tot" "retiro_inv"
## [111] "prestamos"  "otras_perc" "ero_nm_viv" "ero_nm_hog" "erogac_tot"
## [116] "cuota_viv"  "mater_serv" "material"   "servicio"   "deposito"  
## [121] "prest_terc" "pago_tarje" "deudas"     "balance"    "otras_erog"
## [126] "smg"        "entidad"

número de personas en pobreza

cat(sum(ingresos_edomex_2018$tot_integ * ingresos_edomex_2018$factor * ingresos_edomex_2018$pobre_extremo))
## 733403
cat(sum(ingresos_edomex_2018$tot_integ * ingresos_edomex_2018$factor * ingresos_edomex_2018$pobre_moderado))
## 4862417

tasa

tasa_extrema_2018 <- (sum(ingresos_edomex_2018$tot_integ * ingresos_edomex_2018$factor * ingresos_edomex_2018$pobre_extremo))/(sum(ingresos_edomex_2018$tot_integ * ingresos_edomex_2018$factor))
tasa_moderada_2018 <- (sum(ingresos_edomex_2018$tot_integ * ingresos_edomex_2018$factor * ingresos_edomex_2018$pobre_moderado))/(sum(ingresos_edomex_2018$tot_integ * ingresos_edomex_2018$factor))

cat(tasa_extrema_2018)
## 0.04352249
cat(tasa_moderada_2018)
## 0.2885514

pgr

z_2018 <- 1514.91  # línea de pobreza extrema 
#z_2018 <- 3275.6  # línea de pobreza moderada 
# Solo se cuenta la brecha de quienes están por debajo de z; los demás aportan 0
brecha_2018 <- ifelse(ingresos_edomex_2018$ing_pc_2024 < z_2018,
                  (z_2018 - ingresos_edomex_2018$ing_pc_2024) / z_2018,
                  0)
peso_2018 <- ingresos_edomex_2018$tot_integ * ingresos_edomex_2018$factor
n_2018 <- sum(peso_2018)  # población total expandida
pgr_2018 <- sum(brecha_2018 * peso_2018) / n_2018
pgr_2018
## [1] 0.008716165

Gini

library(DescTools)
gini_2018 <- Gini(x = ingresos_edomex_2018$ing_pc_2024,
     weights = ingresos_edomex_2018$factor * ingresos_edomex_2018$tot_integ,
     unbiased = FALSE)
gini_2018
## [1] 0.3975647
writeLines(vector_2018_2024precios, "vector_2018_a_precios_2024.txt")
writeLines(vector_2018, "vector_2018.txt")
writeLines(vector_2018_2024precios, "vector_2018_a_precios_2024.txt")
writeLines(vector_2018, "vector_2018.txt")

atkinson

ingresos_edomex_2018_sin0 <- ingresos_edomex_2018[ingresos_edomex_2018$ing_pc_2024 > 0, ]
atkinson_2018 <- Atkinson(x = ingresos_edomex_2018_sin0$ing_pc_2024,
         n = ingresos_edomex_2018_sin0$factor * ingresos_edomex_2018_sin0$tot_integ,
         parameter = 1.5,
         na.rm = TRUE)
atkinson_2018
## [1] 0.3123059

2024

library(readr)

conc_2024 <- read_csv("concentradohogar24.csv")
## Rows: 91414 Columns: 126
## ── Column specification ────────────────────────────────────────────────────────
## Delimiter: ","
## chr   (5): folioviv, ubica_geo, est_dis, upm, educa_jefe
## dbl (121): foliohog, tam_loc, est_socio, factor, clase_hog, sexo_jefe, edad_...
## 
## ℹ Use `spec()` to retrieve the full column specification for this data.
## ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.
# Entidad federativa
conc_2024$entidad <- substr(conc_2024$folioviv, 1, 2)

# Base de ingresos (ya en pesos de 2024)
ingresos_hogar_2024 <- data.frame(
  folioviv     = conc_2024$folioviv,
  foliohog     = conc_2024$foliohog,
  entidad      = conc_2024$entidad,
  tot_integ    = conc_2024$tot_integ,
  factor       = conc_2024$factor,
  ing_cor_2024 = conc_2024$ing_cor/3,
  ing_pc_2024   = (conc_2024$ing_cor/3)  / conc_2024$tot_integ
)
  ingresos_hogar_2024$pobre_extremo <- as.integer(ingresos_hogar_2024$ing_pc_2024 < 2328.59)
ingresos_hogar_2024$pobre_moderado <- as.integer(
  ingresos_hogar_2024$ing_pc_2024 >= 2328.59 & ingresos_hogar_2024$ing_pc_2024 < 4528.71
)

# Solo Estado de México (clave 15)
ingresos_edomex_2024 <- ingresos_hogar_2024[ingresos_hogar_2024$entidad == "15", ]
head(ingresos_edomex_2024)
##         folioviv foliohog entidad tot_integ factor ing_cor_2024 ing_pc_2024
## 42306 1505018801        1      15         4   4340    24678.677    6169.669
## 42307 1505018803        1      15         2   4340   112830.360   56415.180
## 42308 1505018804        1      15         2   4340    18260.883    9130.442
## 42309 1505019801        1      15         3   1286     6533.440    2177.813
## 42310 1505019802        1      15         1   1286     7073.770    7073.770
## 42311 1505019803        1      15         2   1286     7147.537    3573.768
##       pobre_extremo pobre_moderado
## 42306             0              0
## 42307             0              0
## 42308             0              0
## 42309             1              0
## 42310             0              0
## 42311             0              1
vector_2024 <- paste0(
  "[",
  paste(ingresos_edomex_2024$ing_cor_2024, collapse = ","),
  "]"
)
cat(vector_2024)
## [24678.6766666667,112830.36,18260.8833333333,6533.44,7073.77,7147.53666666667,14754.0933333333,25737.6966666667,24508.1866666667,14101.65,11676.81,21254.0933333333,50606.5466666667,19496.3866666667,23672.1233333333,17337.2966666667,12372.65,13032.6733333333,27044.1766666667,15875.6966666667,15086.2466666667,12957.22,10194.8966666667,12964.76,34284.3133333333,65172.1133333333,21365.1166666667,49230.0666666667,24330.2166666667,19803.2733333333,16805.0433333333,28674.6033333333,16575.48,33786.35,14317.4266666667,21631.95,10918.6766666667,9091.44666666667,32428.63,26023.7333333333,11139.5533333333,33426.6566666667,166435.37,58269.6266666667,93505.4266666667,19117.5633333333,65117.4566666667,6427.09,5384.66333333333,30393.0233333333,37676.7,11739.1266666667,36676.35,20115.3533333333,9256.13333333333,13928.3766666667,29085.63,23534,39865.6933333333,20430.92,47945.02,48672.5433333333,9599.18333333333,44712.69,22088.8933333333,85965.4533333333,50030.28,49921.8,58852.4533333333,105048.703333333,112811.473333333,217540.97,69098.3466666667,80758.7666666667,85316.95,34916.54,52599.8466666667,29901.7366666667,99076.25,32391.3,17090.1466666667,14210.06,19471.83,6724.9,34688.52,31031.1333333333,16503.82,16544.26,10792.4533333333,10852.4566666667,43524.5733333333,43348.3333333333,12287.6666666667,8008.14666666667,11868.13,85544.24,18256.6533333333,17338.27,25408.48,20942.26,36761.7333333333,26404.6166666667,29660.1933333333,57566.4266666667,103948.086666667,15001.78,12290.2133333333,35771.5133333333,54592.71,102414.576666667,64438.6766666667,55979.13,51019.84,72882,5790.32666666667,6289.85333333333,10685.2233333333,9456.74666666667,9804.39,23520.4533333333,12715.9266666667,13725.9333333333,21592.1566666667,36189.5633333333,9326.83333333333,5282.60666666667,9264.58666666667,11711.4533333333,9068.84666666667,4670.72,31508.1833333333,4852.45,26743.7,19195.39,15701.6033333333,18265.9,15742.9633333333,11359.3933333333,10243.1666666667,11597.2,4849.92333333333,16389.9933333333,15867.79,6998.50666666667,8788.04,16891.0833333333,34399.95,11130.3966666667,19833.11,19416.3833333333,10852.4566666667,11418.03,16473.1,36726.59,8265.48,46257.5366666667,12713.63,21570.0033333333,57460.7133333333,18296.62,32698.1,41404.6233333333,30272.5366666667,11393.4366666667,22638.1533333333,25082.5433333333,11868.8466666667,14144.26,49362.6233333333,23999.3933333333,21915.1233333333,52508.1833333333,50852.4366666667,37885.9433333333,15540.9833333333,43459.6933333333,13017.1766666667,16690.8,13809.7733333333,11674.9533333333,18758.5633333333,20744.9666666667,15430.2,8364.05,15981.2433333333,33932.9566666667,20011.2033333333,86012.18,131144.786666667,25335.63,31531.5466666667,49493.5733333333,9235.84666666667,30763.71,44016.38,21508.1966666667,16762.2933333333,33704.9033333333,41315.5733333333,26930.56,47010.73,35982.5233333333,17347.57,33655.7233333333,34090.1533333333,21903.96,15510.02,23205.57,24836.06,36815.21,33008.1333333333,3387.09666666667,30302.89,43251.22,17404.8866666667,27151.6233333333,30001.03,21125.6866666667,27562.23,19588.35,44673.6066666667,30255.8733333333,29274.1866666667,7989.13,9224.04666666667,17512.72,23698.7733333333,45732.3633333333,17170.6366666667,20289.4333333333,16459.0066666667,38236.9366666667,25024.54,13670.1966666667,28358.54,27359.7433333333,14201.0833333333,29242.6333333333,48615.42,46012.0666666667,23539.2666666667,9230.36,28306.4466666667,17849.4633333333,29550.26,11478.2566666667,36454.03,40056.09,26759.5166666667,10844.82,9956.04,9727.64,9289.51666666667,22292.62,7120.87333333333,36961.18,11035.4366666667,24376.35,31606.0966666667,32531.0033333333,29613.1133333333,28410.1266666667,50009.2933333333,26699.33,12825.8733333333,43171.5766666667,57537.1866666667,25616.44,49682.6333333333,19912.2333333333,15088.7266666667,19437.21,14343.1466666667,21903.0266666667,9063.73,15334.01,21813.3666666667,34080.4,59397.7866666667,15836.04,8207.38666666667,18454.6666666667,9677.67333333333,40800.41,17606.55,34426.2133333333,9155.73666666667,25169.7633333333,6442.61666666667,10533.36,27958.56,26001.16,22231.6333333333,26005.1333333333,20852.43,10385.24,12262.6533333333,24385.2366666667,23517.7033333333,42444.8966666667,15431.59,22647.5366666667,16046.1766666667,10434.7033333333,37439.3,14471.7633333333,54133.0066666667,21775.7633333333,16348.87,22426.6233333333,52142.35,3244.48,6908.56666666667,11999.9966666667,17048.66,41526.2433333333,24343.43,13316.2266666667,20934.4166666667,24290.0033333333,30040.9666666667,20795.0766666667,7096.40666666667,23068.1233333333,30533.7,38343.92,21547.24,19713.0033333333,13298.4,30946.6966666667,13593.0866666667,21776.6633333333,38089.6066666667,9599.67,13095.07,17210.4,57717.7633333333,17951.34,10038.1833333333,15885.54,12359.7433333333,8829.50333333333,38310.95,19666.8433333333,12206.1466666667,4227.67,37253.4566666667,13840.9433333333,27877.3133333333,5934.42,26306.19,19154.6566666667,43254.0933333333,18578.68,16008.1933333333,3992.50333333333,20696.7066666667,19967.2066666667,10972.8266666667,10090.16,58592.6933333333,12772.0633333333,28230.4166666667,14627.6133333333,23025.6766666667,26052.66,12829.4333333333,10506.5833333333,13973.5433333333,44641.47,9638.03666666667,17604.9166666667,14086.9533333333,19827.13,18126.2233333333,10173.91,15163.04,8315.21666666667,23369.0033333333,19113,10367.7166666667,23209.02,33012.1166666667,30106.5533333333,39762.2866666667,22745.8966666667,20991.7966666667,22721.3066666667,81934.42,19823.9766666667,23652.8633333333,7432.12666666667,10047.19,26865.3533333333,33207.2166666667,13221.3033333333,6816.39333333333,11580.8333333333,9836.06333333333,24222.65,18062.2933333333,20049.1666666667,22908.1966666667,16899.99,10819.6666666667,13721.3066666667,13479.0133333333,11291.75,8938.23333333333,15613.63,14800.7366666667,14744.91,14216.18,18610.0033333333,32440.4133333333,6732.66333333333,13505.4933333333,37110.7133333333,17853.9333333333,24911.9266666667,13940.73,15718.43,22685.86,10299.5033333333,13036.37,7808.75333333333,14976.8533333333,11913.5866666667,36544.34,20310.16,6418.01666666667,10234.42,30236.2233333333,16975.17,14934.4233333333,11889.5166666667,14316.2633333333,16308.86,11052.2933333333,7039.19666666667,17256.3066666667,22475.0033333333,14684.1466666667,16265.5033333333,25310.1633333333,13001.77,26957.9033333333,8865.56,33848.3366666667,10661.22,18610.0933333333,13846.6366666667,44016.39,13980.3233333333,21639.3366666667,3970.25333333333,15590.1633333333,12168.76,6116.63,9599.13666666667,19022.38,4963.93333333333,13378.9033333333,12203.2766666667,14169,34242.5366666667,27417.44,15978.2566666667,11934.78,16604.13,19512.6133333333,33110.96,17556.0933333333,11099.5733333333,19158.65,12701.6066666667,18071.5266666667,10159.7466666667,27354.65,31805.0433333333,37431.49,32740.1733333333,12844.26,13708.7466666667,22536.05,11311.4733333333,28453.0266666667,27839.3166666667,28670.5833333333,6628.33666666667,18137.2633333333,25888.3266666667,15147.26,25116.9966666667,37968.4433333333,14996.27,8869.18666666667,6001.05666666667,19214.45,13743.5233333333,59046.7266666667,27339.48,11828.3933333333,14788.0366666667,12704.6533333333,24158.4166666667,20546.6,70394.0133333333,18865.3366666667,15403.15,22770.4833333333,19754.0866666667,15208.8333333333,17364.5033333333,27525.2266666667,5491.82,19534.7466666667,12572.12,12407.08,20086.11,14510.6033333333,10868.4,12631.0366666667,22349.57,15516.9966666667,29263.6666666667,13185.4766666667,4891.30333333333,15310.3,31631.1266666667,20145.8933333333,23274.4766666667,31703.26,42760.2433333333,12945.6433333333,52136.3366666667,35001.7466666667,34864.39,14777.6833333333,12792.5966666667,16182.9366666667,11414.31,12137.94,40024.54,17328.8866666667,23492.96,28754.8066666667,42012.6133333333,22447.8833333333,12939.5533333333,20841.91,29316.7233333333,32964.37,36189.36,29688.5166666667,27008.1866666667,89016.3766666667,39532.7733333333,45918.03,19303.27,25973.8666666667,14836.9533333333,10108.6933333333,21541.1966666667,10067.93,30295.0733333333,40147.5166666667,28803.2666666667,62944.3933333333,7999.99666666667,28393.43,29959,13918.03,12336.0566666667,22065.5666666667,15797.5166666667,12868.8466666667,9147.54,42377.04,6098.36,16456.71,13153.92,23425.66,32981.0633333333,33214.77,18048.3633333333,43972.7366666667,13695.65,20511.9066666667,19095.0233333333,20033.3,29916.68,13871.0533333333,11476.1433333333,18745.2533333333,16706.8366666667,16040.5766666667,9995.91333333333,18995.92,24815.66,13144.5933333333,4723.52666666667,5909.21666666667,10082.92,9995.61,14463.94,88550.9966666667,26659.7066666667,7826.08666666667,23327.8666666667,15122.94,21672.13,26344.2533333333,16160.64,25680.3166666667,18939.3366666667,12612.91,13213.1033333333,44909.83,27950.8133333333,17282.84,9761.57,14824.4333333333,20917.59,45915.14,33901.1566666667,18527.07,17923.16,18125.2666666667,42828.92,12354.6066666667,23515.9866666667,17119.56,23798.4966666667,18422.78,28195.5766666667,46329.3866666667,21129.1733333333,19234.1833333333,106114.313333333,14663.9233333333,19614.2766666667,21496.94,12426.2233333333,33540.97,8852.45333333333,3937.7,30057.37,11860.11,29798.5633333333,11688.5166666667,22539.44,36322.39,24403.9233333333,36805.7433333333,20403.57,16304.3466666667,24403.9233333333,47264.8433333333,11948.45,93526.9833333333,33991.93,31199.15,13797.43,13389.5633333333,11936.62,7859.84333333333,10539.4333333333,4464.75333333333,30294.52,9593.17,22857.7066666667,43598.87,51458.6166666667,5085.32333333333,26037.61,42110.7,7809.15666666667,23227.62,47386.5166666667,26479.3066666667,18167.6566666667,16630.4266666667,13359.1266666667,54242.93,7281.34,11739.13,31201.78,57572.0266666667,29009.5733333333,18187.4866666667,76947.1766666667,17478.2566666667,21113.9,18941.77,21214.46,18317.4,15309.21,27341.9266666667,10660.9766666667,22682.91,28021.97,43080.46,37031.01,23901.09,122282.6,52699.8533333333,143895.33,175008.75,21098.87,11800.3933333333,29872.04,13321.0066666667,10337.4766666667,24442.6133333333,19805.0666666667,23680.7633333333,31473.5133333333,20646.1266666667,30920.3533333333,95221.0066666667,142699.39,88626.0033333333,55447.8733333333,48084,88206.52,28727.2066666667,116393.433333333,79262.2933333333,18233.4466666667,18758.2266666667,13814.2433333333,25122.99,11533.8266666667,18381.4033333333,39743.87,27153.6766666667,11030.4466666667,18501.1066666667,41744.1433333333,30613.1066666667,35887.1,23124.0666666667,20838,24814,20287.2433333333,8035.29,20028.0466666667,15217.6533333333,13417.4266666667,21617.22,18878.8533333333,16513.2566666667,21956.0033333333,11871.8,31048.9,7848.27666666667,47754.8633333333,35010.6933333333,207869.616666667,24122.9433333333,40039.3433333333,16819.6633333333,7868.85,21568.2933333333,21024.5766666667,12595.5233333333,27518.9533333333,27362.63,38611.2966666667,13061.85,17254.5166666667,15951.28,14748.2133333333,10991.7933333333,32983.6,13762.2866666667,7803.27333333333,69512,16655.28,10948.1466666667,9193.24,6995.08666666667,8655.99333333333,15833.68,28244.52,17090.85,30320.7066666667,15452.6133333333,8275.77333333333,10215.7533333333,30924.0533333333,13595.9166666667,30749.9933333333,17535.06,16577.8366666667,23485.0833333333,5338.5,13990.1733333333,17396.7933333333,11344.2566666667,13146.1133333333,17127.84,35698.0966666667,10565.2166666667,45549.3433333333,12501.7466666667,14389.7266666667,13555.74,9943.69333333333,12266.65,12339.1233333333,30696.3366666667,72132.5266666667,19164.6233333333,15377.0466666667,18155.7333333333,30885.24,27959.7,22229.5,32614.7433333333,13631.1433333333,32852.4533333333,124852.456666667,0,39960.2866666667,66019.7666666667,49020.9066666667,19401.6333333333,3983.60333333333,12836.0633333333,8902.43666666667,47983.5966666667,30516.38,22001.6333333333,47385.23,73999.9833333333,44981.2466666667,21545.89,36749.6766666667,10750.3466666667,19415.8733333333,33559.3066666667,32846.5133333333,15431.9233333333,15691.6133333333,67952.2966666667,63807.11,35506.64,53657.08,72713.43,17930.2,26005.74,12857.1633333333,26136.9033333333,30560.2533333333,9836.56333333333,16858.65,6659.23666666667,18338.8766666667,26338.8633333333,18797.88,21372.7766666667,32960.8666666667,10722.3966666667,15211.8766666667,56616.1233333333,22099.2266666667,6023.26333333333,44325.2666666667,13763.0466666667,27981.5933333333,9801.31,30180.3266666667,47579.5633333333,23042.6066666667,122700.453333333,14475.4033333333,50086.6366666667,49336.05,19891.3033333333,19401.5133333333,12010.8233333333,14073.88,26031.7166666667,16255.43,51077.3,53046.97,32311.5266666667,28955.98,44517.8433333333,82414.9133333333,33153.7466666667,10530.9966666667,21095.01,49096.94,36807.0533333333,47496.65,20351.5466666667,16222.2966666667,27649.96,29744.59,53272.8133333333,21216.83,27963.4333333333,14098.5233333333,14213.8933333333,20799.2833333333,5537.7,33877.0433333333,12163.93,28032.78,10836.0633333333,12425.69,18789.0866666667,23513.5933333333,6457.25333333333,21843.08,27179.95,17791.72,56440.39,11410.1933333333,18177.7466666667,18691.6766666667,30161.1033333333,89318.2933333333,21326.6066666667,14120.08,15822.0533333333,21264.6766666667,18867.5433333333,30132,32810.9266666667,33590.1533333333,28948.67,8577.04333333333,9049.17666666667,24834.9833333333,19426.22,11032.7833333333,12722.6333333333,5192.21666666667,36229.4966666667,21451.7766666667,14883.65,18501.69,13950.4666666667,32516.1233333333,80962.6433333333,36984.57,58490.5266666667,45938.8466666667,19024.73,53033.8033333333,16646.73,44642.34,22744.5566666667,33751.0433333333,32608.6933333333,73927.7166666667,55097.64,19616.6433333333,26056.4433333333,13940.21,31249.9533333333,15861.33,16484.2133333333,20253.7833333333,22972.4,59855.9566666667,18184.9133333333,49531.6,47426.6733333333,52885.2366666667,85245.8966666667,61803.2733333333,60500.15,27189.9366666667,35736.4866666667,24730.0033333333,18527.5166666667,6929.34666666667,17217.82,87923.5433333333,41046.3433333333,40406.3566666667,136462.77,40913.12,65753.8466666667,58354.83,48531.4633333333,60420.7533333333,15702.4766666667,22622.27,20652.8633333333,17485.4866666667,11543.47,13523.6733333333,14836.9533333333,18861.6666666667,28511.5666666667,14376.6266666667,23136.3833333333,12131.14,30107.3733333333,21885.2433333333,6999.99666666667,15702.9433333333,15446.8766666667,19016.3866666667,42455.7433333333,12644.02,50410.3833333333,68537.0066666667,13010.8666666667,10258.5966666667,26359.2633333333,16570.9933333333,27445.0433333333,15535.7066666667,13415.7433333333,17779.5833333333,11983.6933333333,10907.6066666667,12412.3666666667,12207.22,16794.52,27848.8733333333,24822.01,44137.9433333333,32586.5833333333,29918.03,8368.37666666667,9647.54,16581.5166666667,22028.2066666667,19523.14,18055.7433333333,21384.3066666667,16381.19,32393.4366666667,16205.9733333333,8934.42,66721.31,38814.84,19268.85,44526.05,32618.9366666667,36580.5033333333,17532.78,15558.4666666667,17956.48,30049.1733333333,29213.1,28819.66,29663.9233333333,47058.8066666667,23051.58,11446.6966666667,15586.9466666667,9751.04666666667,26818.0166666667,61606.55,22163.9233333333,17052.3266666667,23712.3866666667,36810.68,32135.86,29510.07,20904.3833333333,30828.8866666667,50780.12,14740.1766666667,21267.7,18725.5333333333,38507.1733333333,22790.2566666667,33381.8266666667,23338.3566666667,12456.7633333333,19943.0533333333,71501.1033333333,23004.1133333333,34557.05,62464.0466666667,12024.4533333333,42309.77,10846.3733333333,50960.4466666667,18330.03,15946.6966666667,20706.52,14300.4866666667,21495.9633333333,15576.0566666667,18647.96,15393.8833333333,13605.73,19862.2933333333,25000.92,35237.69,37092.7866666667,31035.4966666667,25723.5233333333,12740.1533333333,24416.4666666667,10011.39,17103.78,12237.6766666667,8157.31333333333,6725.27333333333,23980.85,16934.8066666667,12649.7633333333,6592.63,25036.7166666667,27822.0466666667,28501.04,11298.9066666667,49859.5866666667,44160.6466666667,24134.5033333333,5759.38333333333,22380.4166666667,34508.63,14120.5933333333,12310.3,34421.1433333333,58695.65,42744.8766666667,15510.4133333333,59081.9533333333,27811.4666666667,22704.9133333333,39093.8766666667,20947.3966666667,8951.25,16130.7833333333,16751.4,5141.13666666667,16719.7,19467.3266666667,39364.7433333333,26332.9366666667,81792.7933333333,8625.52,37537.27,24256.9566666667,79802.8233333333,15001.17,16824.6033333333,10122.0366666667,11106.7533333333,15945.92,11954.37,21890.18,18346.6233333333,27755.06,31286.5233333333,14391.8033333333,19195.3033333333,14496.8666666667,32791.08,14315.0366666667,12240.1166666667,9860.65,16577.5033333333,71264.97,12919.8933333333,8901.63333333333,15081.9633333333,23163.92,8868.85,11004.9166666667,14811.47,9868.84333333333,18461.3066666667,15371.8733333333,11503.9066666667,15147.0166666667,11403.45,32176.2466666667,29697.17,28793.4466666667,9505.89333333333,54548.6633333333,13372.9766666667,19915.9833333333,48112.3,23173.3533333333,32386.8833333333,44737.4133333333,8555.99666666667,14181.1166666667,21996.0066666667,57989.6666666667,17639.34,46013.57,22442.6066666667,21922.3066666667,15491.8,29508.1966666667,6247.80333333333,21297.7666666667,11547.7766666667,9881.58333333333,28924.53,38411.0633333333,6299.99666666667,38449.1766666667,1967.21,29527.8666666667,32245.89,19868.8466666667,16419.5233333333,19494.6933333333,22870.26,9544.63666666667,26771.7333333333,17062.3066666667,53769.72,11345.91,14508.2266666667,29508.5833333333,17177.5666666667,16877.0366666667,29459.0066666667,26066.2133333333,31049.1666666667,42852.4533333333,16808.9333333333,22688.5133333333,13196.7166666667,14868.8466666667,5901.63666666667,20424.44,14106.5533333333,12919.3366666667,13353.3866666667,16372.0466666667,13580.2533333333,21884.2,15174.61,16266.9966666667,25648.5133333333,11257.0933333333,10729.31,21376.6866666667,19688.5433333333,13707.7366666667,9189.45666666667,5354.15333333333,17881.4866666667,30607.15,24095.3933333333,26585.2766666667,12406.8566666667,30924.5966666667,24115.1733333333,18642.1766666667,24043.0633333333,10411.3733333333,19523.1333333333,5033.09333333333,45550.6566666667,20180.3233333333,44005.78,27303.2666666667,25745.8933333333,23011.2833333333,13334.18,22416.29,11009.5366666667,6816.26666666667,17294.17,16006.13,16048.0533333333,13051.3666666667,17713.58,25233.1033333333,25166.79,11464.24,21913.2466666667,25629.5333333333,10540.98,13991.7966666667,82985.4366666667,53593.9666666667,56119.0366666667,11843.14,44270.84,71340.35,21609.62,13994.2333333333,13573.3833333333,17832.2133333333,43857.8066666667,52555.1666666667,21483.5966666667,15351.59,41135.4033333333,38267.62,24920.7233333333,8687.74333333333,34119.6733333333,29803.7833333333,8612.52666666667,20739.1266666667,15280.5866666667,10739.8266666667,15234.5633333333,29454.78,22114.75,27934.42,15281.66,23991.8,19565.57,17885.86,20852.99,17417.4833333333,21061.0033333333,17384.2766666667,17344.26,9172.13,19131.1366666667,18262.2866666667,12240.9766666667,20295.6433333333,10886.44,37972.8166666667,16979.03,22744.56,12160.6333333333,21334.8733333333,38075.62,21971.19,22770.4866666667,29636.7533333333,32713.1133333333,25655.73,7245.9,10278.03,60409.8266666667,33724.2166666667,19162.0966666667,58911.5466666667,8920.04666666667,32734.18,40592.86,27698.7833333333,14583.0633333333,45694.86,28809.7733333333,10676.7133333333,16124.9933333333,74072.75,25350.6266666667,12700.59,13713.5166666667,7594.96666666667,9720.41666666667,14603.3366666667,22457.3866666667,13529.1466666667,20686.1466666667,12943.01,11218.84,16960.7466666667,11352.3933333333,16941.9966666667,40324.43,38852.76,31821.0033333333,44033.01,57040.7633333333,12450.1966666667,42907.8466666667,12717.3866666667,13674.6133333333,9961.95,37950.5033333333,15535.65,24175.1333333333,17803.2666666667,21688.5133333333,13335.22,4688.05,16016.39,15329.8533333333,25212.8266666667,22767.5766666667,8067.84333333333,39283.7166666667,31784.1833333333,21877.6866666667,11568.16,37732.5433333333,25353.55,18892.5233333333,12294.52,10973.2666666667,25907.3833333333,17598.6966666667,14169,20522.4333333333,31882.6466666667,15484.39,24701.0766666667,20748.32,14549.7833333333,10265.42,22057.1466666667,50313.8533333333,19207.5633333333,10250.6966666667,19708.1466666667,8255.27666666667,17741.98,15283.5533333333,17294.1766666667,20815.91,18713.11,13803.2766666667,17165.0966666667,13849.51,67484.57,30060.9733333333,11780.8166666667,85030.98,31900.3466666667,65544.1833333333,42778.9666666667,4528.21,10327.8633333333,14765.5633333333,10677.7833333333,21262.2833333333,24581.9366666667,10190.2133333333,23525.3933333333,21329.5433333333,22473.6933333333,8284.16,3028.77,18494.75,11762.2666666667,15432.06,28466.0633333333,18932.7566666667,13648.31,13376.4,20234.2133333333,16687.9033333333,14198.16,28183.84,15641.65,18908.2966666667,34841.0666666667,22561,29283.0433333333,13751.39,41967.2,55358.9033333333,15236.8,13893.4366666667,18745.5566666667,16098.2633333333,20005.4266666667,17882.1833333333,16891.2966666667,9684.06333333333,7837.19,13891.8733333333,7043.57666666667,14298.1133333333,9857.13666666667,17556.0933333333,23359.39,26370.9566666667,51523.3,19568.0466666667,29959.9733333333,17587.6533333333,9293.47666666667,17000.6933333333,7410.32333333333,14839.9966666667,11212.6533333333,49963.0466666667,17603.4033333333,84893.0333333333,16539.8833333333,15065.2066666667,18520.15,19463.1733333333,23772.7766666667,43409.0033333333,65153.3466666667,119625.453333333,16685.8433333333,12149.3033333333,23098.3566666667,26639.34,52918.03,17901.1,26065.57,16318.0233333333,17823.6866666667,43101.6266666667,16261.4666666667,26114.74,26918.03,36308.3733333333,13249.4033333333,14836.0633333333,15124.7033333333,18018.5933333333,85070.26,25225.4766666667,31229.5033333333,31475.4066666667,28159.92,15740.9233333333,35438.28,12685.83,21118.2833333333,28416.2333333333,32665.9766666667,12170.4,18309.6966666667,35549.2533333333,29249.9933333333,48664.4066666667,9782.60666666667,11640.8533333333,12300.53,14957.6033333333,13885.8133333333,3716.55333333333,20429.4833333333,11147.5366666667,12718.0266666667,33302.2866666667,14016.3833333333,77262.29,24901.6366666667,76819.6633333333,9836.06333333333,9662.11,13207.94,16000.9766666667,12250.5,11713.13,29206.55,31724.8133333333,21606.5466666667,10295.0766666667,14631.8333333333,33037.2633333333,33886.74,15326.0833333333,12771.88,16001.6366666667,20598.3566666667,16081.96,23114.74,22081.9533333333,11795.0766666667,50339.6666666667,33276.01,22129.1766666667,36099.3933333333,15652.1666666667,25230.9733333333,20477.5733333333,25397.9633333333,40991.22,3815.66666666667,6560.12333333333,30780.7533333333,17340.8666666667,7614.59333333333,16881.83,20457.21,9598.52333333333,13421.89,11347.8233333333,6316.76666666667,16088.7033333333,27518.05,18572.2766666667,42082.3,13377.0466666667,34750.8,40470.4866666667,25683.5966666667,9877.55666666667,20122.94,33607.0633333333,32311.5233333333,10597.8233333333,24399.9233333333,12115.1766666667,31131.3266666667,10966.9533333333,17159.1133333333,16380.8866666667,18684.8566666667,14246.1933333333,11902.6933333333,9913.19333333333,9272.43666666667,20625.2266666667,14033.3033333333,22349.6133333333,6468.45,13674.13,18278.88,35556.3066666667,28589.3166666667,17106.88,26163.93,45358.1266666667,18522.4466666667,23473.0633333333,18639.9366666667,31437.6933333333,20870.3366666667,17112.64,25606.5533333333,15353.4966666667,28715.08,17527.3033333333,17172.9133333333,3653.22,16403.22,13142.3466666667,20171.8466666667,20909.1866666667,19997.2633333333,5470.13,26676.06,54437.0633333333,32665.09,57192.48,30778.95,9337.65333333333,9486.69666666667,35914.53,8756.14666666667,41355.55,15090.5433333333,30976.33,10437.27,16967.2033333333,15982.7466666667,27375.8733333333,38084,76222.61,67672.12,102934.32,19754.0966666667,47829.4966666667,25508.72,36014.1266666667,18482.02,15117.5366666667,22384.6833333333,23700.13,24468.05,27532.5166666667,13201.7866666667,39762,16981.2166666667,52819.6566666667,19252.91,11126.22,32316.37,26696.7,10301.63,19569.63,10374.7466666667,19428.5566666667,62604.07,12276.22,11100.1266666667,8700.50333333333,27651.7366666667,20148.3133333333,8241.58,6859.02333333333,29155.1566666667,17729.98,47958.9633333333,45292.17,21820.46,17474.05,16079.66,15156.1366666667,15426.14,20398.7566666667,28406.4633333333,11893.6633333333,31031.4733333333,32326.7833333333,29878.7566666667,36385.7533333333,45107.3333333333,38890.23,25721.28,26299.9566666667,16689.96,22171.7966666667,22881.3033333333,48179.3333333333,33073.7633333333,28588.24,13335.9433333333,79570.5566666667,15053.8466666667,16121.3066666667,28763.46,29344.3833333333,50510.0233333333,8326.78,37859.54,35568.2366666667,9524.89,22024.68,31046.28,11234.2166666667,25571.8833333333,16240.02,9307.55666666667,18072.93,15620.61,11666.9133333333,17278.6733333333,25061.6266666667,22844.2533333333,31622.94,13729.5033333333,11175.0066666667,6885.24333333333,23704.9033333333,10278.68,20946.6966666667,48075.6733333333,14315.2066666667,10104.99,13039,9798.53666666667,19677.22,10989.99,13962.5933333333,8364.05333333333,26205.5433333333,11539.5166666667,12022.4566666667,18747.0666666667,16278.2366666667,7377.45333333333,23107.83,52844.4833333333,50028.8966666667,8739.13,16100.6666666667,18974.6466666667,12546.0766666667,24308.24,53558.13,48647.5633333333,39752.9133333333,8853.87333333333,41168.4766666667,16493.5766666667,44453.0133333333,45586.9533333333,36708.44,20469.84,17762.27,22778.7433333333,6847.82333333333,26594.49,11555.0466666667,14642.35,18907.6466666667,63807.6833333333,23955.8333333333,37981.7833333333,31230.2233333333,22483.6566666667,23770.49,17426.2233333333,13893.4366666667,17934.4233333333,11959.0066666667,12540.97,21132,8039.09,16115,19763.17,33922.23,39196.69,37080.32,22688.5066666667,49180.3266666667,45491.7966666667,31927.8566666667,71229.49,26810.3133333333,17479.42,21374.9966666667,20252.7666666667,60254.61,10856.5166666667,37273.8366666667,25707.7333333333,14504.8633333333,27449.2366666667,21211.7966666667,25641.72,25147.5366666667,39163.9233333333,33934.4166666667,15740.98,31721.3066666667,14426.22,28663.9266666667,12836.0533333333,21331.7933333333,9385.24,20819.6666666667,10365.57,13073.7666666667,18065.57,19809.78,15375.9933333333,59050.6533333333,43681.44,31467.39,12376.7733333333,19202.31,14778.0366666667,24363.36,21892.5233333333,18728.3,21414.44,51876.74,12485.97,29474.65,25800.31,19914.4366666667,44483.0333333333,39730,78761.8033333333,26550.8133333333,25818.77,84626.21,35754.0866666667,34965.5766666667,14088.8166666667,19819.3433333333,29686.4966666667,8581.43333333333,26306.82,119463.52,100399.46,114503.09,112308.926666667,126003.146666667,81702.7433333333,35110.1833333333,28714.0533333333,19941.9366666667,12391.8233333333,15033.7766666667,4992.39,22343.2633333333,31409.87,18109.5933333333,25553.4766666667,24339.33,13836.0566666667,20160.6233333333,16876.29,9049.17666666667,22778.6766666667,16843.5366666667,22580.31,11144.79,8922.53333333333,13696.7166666667,15795.0766666667,7401.82,19704.9033333333,17533.2333333333,20991.79,13353.78,18370.2633333333,19220.4833333333,39587.9933333333,11719.37,19515.46,22928.0733333333,8271.19666666667,10066.2966666667,17572.2366666667,23050.1366666667,15739.07,10926.4433333333,8009.39,16480.42,17786.2433333333,9000.18666666667,22818.15,30538.39,22764.97,7167.20666666667,12709.73,44781.7333333333,45690.57,21691.75,11432.76,10614.7333333333,12538.5633333333,11291.8066666667,7999.95333333333,12956.0166666667,14340.6566666667,25047.1966666667,21130.69,17292.0266666667,41027.0866666667,38987.33,15803.27,30842.6433333333,31973.77,31803.2766666667,331162.336666667,53541.37,22794.0333333333,36778.8866666667,32947.8766666667,57573.49,17487.9166666667,10136.98,53903.2,34854.48,16190.2133333333,39570.64,38596.1166666667,23135.8633333333,40743.5,13936.54,12304.94,51755.5633333333,9933.44,9807.68666666667,18886.08,19889.98,15906.5166666667,15126.2266666667,28093.4333333333,14968.4366666667,10885.2433333333,9049.18,17737.69,21229.5033333333,5631.14333333333,12295.08,22538.6766666667,72881.87,49229.9666666667,10575.3833333333,17606.55,22836.0533333333,32403.25,39531.9633333333,17493.2833333333,20909.8233333333,44715.52,34586.45,12071.3066666667,25778.6733333333,6235.58333333333,17182.6633333333,22818.0266666667,18172.12,9691.47,13401.6366666667,9360.65333333333,14917.31,14604.2066666667,26301.6333333333,14303.2766666667,12983.6,9262.37333333333,6995.78,20921.3033333333,22008.1933333333,17564.1233333333,53580.92,34859.4633333333,13216.7933333333,83573.7566666667,25385.2433333333,25868.8466666667,23122.84,30234.5633333333,20949.7666666667,38902.1633333333,27196.9133333333,19421.1233333333,18487.98,20920.6233333333,21383.93,11089.7133333333,5869.56333333333,10469.63,13674.6133333333,19230.7066666667,54359.2133333333,9459.33333333333,25689.9333333333,26035.4033333333,23377.02,12817.5133333333,24310.13,30057.3733333333,26803.76,6926.52333333333,16893.4366666667,8057.37333333333,34459.0066666667,14685.2433333333,12708.1966666667,3436.90666666667,46454.86,28484.4966666667,9989.05,11376.4133333333,35780.6733333333,16180.2166666667,19946.5233333333,21203.01,5324.25666666667,12727.3633333333,12446.1333333333,19519.19,57606.8833333333,14320.8133333333,39622.9433333333,78209.8266666667,14969.8233333333,3431.14666666667,10394.32,12715.73,14558.4333333333,11488.1866666667,3134.21333333333,4065.56333333333,7334.3,24943.3166666667,13363.9333333333,9612.92333333333,27077.82,12315.0933333333,14636.06,14286.8833333333,8959.01333333333,17442.6166666667,12879.9033333333,7968.09,11302.5866666667,7166.64666666667,20026.7666666667,28768.7333333333,20400.4733333333,19197.6333333333,13344.56,5140.59333333333,3706.87,26370.8233333333,7189.69333333333,16422.0766666667,10714.4666666667,15909.1766666667,17756.5433333333,3327.86,9762.29333333333,8098.35333333333,3641.64666666667,9245.08666666667,8772.78666666667,10739.8233333333,9725.85333333333,11908.0066666667,23813.14,11656.7633333333,22304.0866666667,8365.90333333333,28213.6233333333,8925.15,29122.74,4905,12749.1766666667,12375.1166666667,9853.85333333333,3978.54,16962.7633333333,7434.39666666667,12734.4233333333,13590.1566666667,24593.0966666667,12163.8966666667,14093.6366666667,11559.25,11028.0433333333,15957.2133333333,5907.95666666667,7326.96,9135.72333333333,14151.73,35856.2933333333,32826.4266666667,30695.5766666667,18823.2766666667,9888.29666666667,13126.5566666667,13601.24,4713.98,11002.8,10865.4266666667,4576.71333333333,11336.9933333333,11207.49,20242.3266666667,19031.3733333333,14387.06,5871.2,15389.5833333333,6679.10666666667,12546.0733333333,54668.7566666667,16539.4366666667,14884.1666666667,22963.1833333333,15773.33,16301.3333333333,23590.1566666667,17245.8933333333,22529.5,13931.8866666667,8926.22666666667,7278.68666666667,5950.81666666667,7222.47666666667,18819.6666666667,19942.6133333333,13606.5533333333,19496.9033333333,11643.2766666667,14076.39,5891.52333333333,13542.52,10240.53,6699.21,6165.87666666667,6230.29666666667,8990.33,32169.52,10525.25,6405.23333333333,22880.7766666667,21007.8266666667,14633.7666666667,10433.5033333333,13432.4866666667,4297.21,3958.84666666667,4605.93,7897.08,2690.21333333333,4425.19,18111.8633333333,7125.50666666667,10641.16,11283.4233333333,14548.41,29854.3133333333,8187.50666666667,3671.68,7862.37333333333,9816.10333333333,16963.4166666667,8288.29,9875.69333333333,7401.63666666667,10081.96,13786.8766666667,8436.06333333333,5934.42,15425.7433333333,14944.4266666667,15599.3433333333,15458.51,3527.47666666667,15510.14,10520.9433333333,5702.00333333333,29231.7933333333,45797.71,11592.91,9888.73666666667,18402.8633333333,29658.6466666667,9993.81666666667,13322.6866666667,8442.39333333333,16918.03,16968.91,9668.85,6372.25,13196.7133333333,6395.73333333333,34108.2133333333,24310.5066666667,14092.87,26643.1866666667,15888.5033333333,6869.07666666667,7771.63333333333,11612.5133333333,12126.57,12826.5066666667,12276.2333333333,8350.98,33159.6466666667,40626.71,8639.95333333333,8346.89,9112.66333333333,22713.1133333333,41153.7366666667,35270.6233333333,5532.95666666667,16286.22,10739.8266666667,7462.64666666667,23000.03,7799.78,33140.2066666667,27624.22,15168.86,15097.8233333333,18267.7,23762.2666666667,34565.21,17540.84,22059.78,39402.14,59386.7333333333,21783.1233333333,23722.82,21541.4533333333,36765.35,45065.55,14898.8233333333,9131.14333333333,39045.41,10145.3733333333,14820.7133333333,8511.38,20592.0166666667,12479.66,4388.57333333333,8668.22333333333,9103.60666666667,11186.9733333333,11489.8266666667,4779.56,8490.88,25399.7866666667,8996.71333333333,10418.03,20934.42,9540.97666666667,9306.55666666667,13090.1566666667,15120.97,8183.6,16909.83,11729.5,8022.47333333333,18591.9466666667,11362.9866666667,20250.71,10462.3066666667,17548.0033333333,18419.1666666667,22513.8133333333,11216.9433333333,24471.3333333333,7468.52666666667,28926.5033333333,10220.6533333333,23508.1933333333,12125.47,17915.3066666667,16908.6566666667,20652.86,40302.4933333333,15546.9766666667,15525.1466666667,11565.5633333333,6327.48,12250.3433333333,25543.35,16398.9933333333,28781.7466666667,20424.6033333333,22159.1733333333,16787.24,30688.0033333333,12082.4333333333,19695.1733333333,14725.2566666667,11616.84,14517.8666666667,18982.85,48743.7866666667,20445.6466666667,6172.59666666667,7009.82666666667,16272.1266666667,14731.4633333333,42139.3366666667,18316.6366666667,23243.8666666667,9473.23,27795.38,51649.7566666667,7451.60666666667,7287.91,39950.84,12284.63,19601.2366666667,27936.1733333333,5257.28333333333,4937.69666666667,28244.91,24577.0433333333,20590.16,7050.81,11636.0633333333,7208.66,15373.4333333333,26270.14,14503.3133333333,9887.33333333333,19614.2566666667,31620.9566666667,49956.0933333333,9580.93333333333,15436.53,17832.2166666667,9766.16666666667,31320.5033333333,39239.3333333333,18360.6533333333,28462.9833333333,29590.1633333333,48901.63,78363.2066666667,47123.0633333333,37205.7466666667,56787.6866666667,41741.6266666667,16625.0733333333,13448.0266666667,11501.6,9383.67333333333,16195.6233333333,8411.86,17136.9966666667,19042.5966666667,13766.3133333333,10301.4066666667,3531.12666666667,27500.6466666667,12490.78,16362.3433333333,7020.10666666667,13624.76,24307.06,14283.13,10548.0866666667,3908.30333333333,17233.86,29669.1666666667,8650.67,12434.4633333333,25080.45,14831.6933333333,11294.7,18239.8266666667,19669.2966666667,35886.37,27641.9566666667,6434.40666666667,13073.7633333333,12975.4033333333,16262.29,17639.3366666667,12336.0566666667,19821.0566666667,15854.91,19477.7966666667,14132.65,24082.26,15480.1033333333,9823.89333333333,8671.39,7971.57333333333,12643.6433333333,22646.4166666667,11703.8366666667,10310.13,4565.83666666667,7753.51666666667,9374.41333333333,5750.81666666667,12778.6833333333,8483.6,12151.5133333333,3459.01,16081.9533333333,21503.03,16897.6966666667,16130.8133333333,9033.72,10852.4533333333,23415.09,34842.81,20858.2233333333,21844.9833333333,17523.8833333333,18590.5133333333,29932.9333333333,6613.67666666667,3516.98666666667,212826.493333333,17969.9566666667,17997.23,10836.4666666667,23824.3433333333,21699.9166666667,10865.63,10762.9466666667,12247.7166666667,12693.19,31389.8133333333,17069.64,27387.54,17597.0666666667,18545.7933333333,7880.78,18173.24,30120.4266666667,13546.5766666667,26551.7233333333,16103.96,33575.2966666667,17145.1666666667,34309.7333333333,10144.2533333333,25577.2066666667,11819.6666666667,16930.7633333333,12353.16,53475.4033333333,41395.2066666667,24187.33,58762.2666666667,59972.4833333333,18571.2066666667,27211.54,54610.9466666667,19232.1133333333,12648.4766666667,14862.54,10619.71,33937.6533333333,15685.8266666667,7842.91,24708.2433333333,8280.58333333333,13117.3633333333,22510.4533333333,13967.81,17262.2866666667,4345.9,16621.7666666667,14884.79,21066.7333333333,12332.39,19604.81,11730.7133333333,8216.64666666667,4925.96,121209.506666667,14517.1833333333,12076.5933333333,19550.5066666667,43515.4666666667,12095.3266666667,11035.8666666667,6826.78,5190.7,14365.1766666667,13518.01,16311.0633333333,7712.63333333333,7470.48333333333,16241.79,8671,10289.21,23506.54,20141.5566666667,8566.72,9114.48333333333,9100.68666666667,15665.08,14844.84,10266.9966666667,7463.01666666667,11987.1133333333,11477.73,16171.6266666667,14872.3233333333,12256.5366666667,3179.21333333333,16609.7533333333,19250.6966666667,16632.95,12088.35,16333.2933333333,11927.47,10961.7566666667,8699.77666666667,24442.6166666667,11649.3933333333,5742.61333333333,7081.86666666667,42675.9933333333,20995.79,21330.0266666667,18416.9166666667,27033.6533333333,42258.7533333333,11951.6033333333,19803.7566666667,15815.6766666667,16713,5828.09333333333,11113.72,23522.4333333333,6995.34666666667,24183.45,25360.7866666667,20147.5333333333,11794.4233333333,8166.88,17514.4966666667,14869.3466666667,11677.2333333333,9847.56666666667,18364.2633333333,22835.6733333333,11118.3033333333,17996.3166666667,39668.99,11859.3233333333,5878.13333333333,14305.3033333333,18228.6666666667,12781.3033333333,7839.03333333333,47561.1666666667,12196.17,27932.7933333333,20022.25,39646.14,7738.09,7976.59666666667,16199.2566666667,21984.57,24356.24,6723.68333333333,15788.9,14579.1166666667,5103.23,7321.17333333333,1451.61,8212.29333333333,12360.9466666667,11555.0466666667,3881.48333333333,8721.31,22598.3566666667,9505.24,8327.86666666667,12111.4666666667,12784.49,6675.94,13286.1,16121.9333333333,7787.66,24411.4733333333,17554.0633333333,11548.1866666667,10829.86,26894.53,8256.23666666667,28337.4633333333,9819.88666666667,8450.35666666667,15604.2533333333,9129.28666666667,13740.3833333333,8538.74333333333,7703.30333333333,6197.55666666667,14902.6933333333,27433.1233333333,12291.85,5447.70333333333,15396.65,11416.5033333333,6144.37,13728.2533333333,21156.13,11109.9733333333,45696.33,31989.6466666667,6614.57,11269.95,5777.07,14360.6466666667,19880.32,12547.42,5357.36666666667,7386.38666666667,12600.2466666667,9544.20666666667,6286.76,21912.4566666667,5132.02,16109.7433333333,10470.2966666667,15360.7933333333,9527.18666666667,5273.95,3501.75666666667,29880.3266666667,15114.6466666667,5635.12666666667,13631.4766666667,10858.61,25010.7233333333,13442.4766666667,13776.4866666667,15190.3433333333,16442.62,11819.6666666667,8592.48,10022.0233333333,7434.42,19768.34,11733.0033333333,16776.9066666667,42332.08,21788.7066666667,29758.2766666667,8832.21666666667,7517.64666666667,13664.09,32285.0466666667,46357.1233333333,81441.72,20217.9333333333,19770.4866666667,12808.83,23967.39,7310.65666666667,10503.6733333333,10701.35,8654.5,6758.32666666667,13382.2633333333,22244.2433333333,11404.71,20165.37,23236.72,11299.4333333333,7661.78,7430.54666666667,25478.4533333333,12491.7966666667,15655.7333333333,14554.35,11719.99,10999.9933333333,6931.96666666667,10778.6833333333,8245.89666666667,8909.83,9786.88,8155.73333333333,13508.19,14363.91,48013.07,27076.7833333333,12587.71,17563.1866666667,21746.6533333333,12194.43,17207.4066666667,11090.1,10456.8633333333,22753.2166666667,20358.9566666667,12139.66,19834.42,8402.8,13384.07,18840.9733333333,19501.96,14022.9433333333,3950.81666666667,19167.8566666667,16163.4366666667,7987.72666666667,32341.0466666667,11630.58,7739.73666666667,11754.5333333333,16676.1533333333,11862.8633333333,11286.0933333333,11912.7433333333,13213.3466666667,22334.1733333333,10268.9833333333,7259.28333333333,39346.05,6466.27,3958.35333333333,6327.25666666667,6780.19666666667,4978.60333333333,15160.9133333333,5740.17333333333,6718.6,8490.88333333333,12230.1666666667,7205.06,9990.07,7484.57333333333,6702.34666666667,25990.17,11621.83,3902.52,14771.7266666667,18081.91,10321.37,7827.91666666667,50484.89,30675.5533333333,10115.4366666667,19605.8933333333,22498.3533333333,26587.6533333333,9907.48333333333,15848.2333333333,5603.09666666667,9843.08666666667,26987.2333333333,10132.6366666667,9663.74333333333,9405.26333333333,18821.5233333333,30987.45,22485.7133333333,9183.76666666667,23912.49,61001.9133333333,22491.7066666667,9608.19,29768.2333333333,14318.6333333333,32385.2433333333,13055.76,34807.1966666667,20885.0566666667,11860.3166666667,20662.5666666667,21927.8533333333,20181.9333333333,16394.5866666667,9742.68333333333,13234.42,11311.0066666667,16874.47,16106.1166666667,14395.39,27675.8966666667,6402.86333333333,12555.39,10434.78,35990.5266666667,38119.2566666667,16486.0433333333,6118.39666666667,4449.92333333333,7672.12666666667,10729.5033333333,14103.2533333333,14168.9966666667,28164.4366666667,19433.27,13325.2233333333,13322.89,22692.5533333333,11278.4566666667,16698.83,21290.8766666667,7061.19,26163.77,4152.77,26918.5233333333,28197.0466666667,26764.9066666667,10814.8866666667,11038.6366666667,14598.6366666667,26843.0633333333,16180.3166666667,13067.39,9045.89666666667,22911.83,1866.97333333333,8959.01333333333,11898.3466666667,11276.8466666667,20309.8133333333,14613.1933333333,12288.7766666667,14505.94,8744.91333333333,16463.7033333333,9521.47333333333,9876.22333333333,8273.14,15622.9466666667,8885.24333333333,5434.42333333333,35565.5633333333,8057.37333333333,26507.37,23163.9266666667,19008.1966666667,12622.9466666667,20269.3,19237.6966666667,21353.4266666667,24918.92,24943.8,3902.52,7136.04333333333,16419.52,19538.3866666667,21943.8,21815.7333333333,13066.9466666667,12220.6266666667,20609.2133333333,9840.44,12508.1966666667,14803.2766666667,12495.59,4754.09333333333,14319.3633333333,29114.8166666667,12844.7166666667,12746.66,11718.0866666667,9222.99,13752.8666666667,15908.8333333333,24659.01,19487.18,8701.13333333333,15706.79,49406.6366666667,25980.4166666667,12369.8833333333,25028.2166666667,9580.08,41608.0166666667,10199.09,8717.56333333333,8154.39666666667,10940.4533333333,39590.1533333333,16774.6966666667,3619.34,4244.02,14410.2933333333,20102.7,34603.1433333333,12616.7966666667,11159.3266666667,24966.3933333333,50992.5466666667,11168.4733333333,25772.43,49772.35,9440.26666666667,7160.57666666667,78180.3066666667,10750.4566666667,19607.8566666667,14593.66,40172.12,25967.21,42805.2066666667,15807.06,11378.4333333333,29631.83,6865.49,15782.3333333333,12894.3266666667,19211.0233333333,24056.42,14916.9866666667,7179.69,3418.65,24644.12,22859.63,39794.9966666667,9602.31333333333,8998.52,13388.2933333333,13834.5933333333,18424.9566666667,19497.5866666667,42080.97,20848.1133333333,19999.9866666667,22563.46,10098.3566666667,9948.23,18950.81,37286.3366666667,18152.9166666667,11281.7833333333,29940.6166666667,27351.97,51987.6333333333,35119.56,38978.2033333333,22244.68,23586.0633333333,6190.40666666667,16579.6633333333,10352.45,18019.94,13136.76,26309.85,15025.4,10249.3866666667,14351.7033333333,28211.7733333333,9544.14,20712.2966666667,15722.4933333333,55603.4933333333,13989.6466666667,11644.64,9782.60333333333,8432.19333333333,20308.6033333333,16218.6133333333,16514.7333333333,31947.2833333333,11280.5066666667,51811.46,11803.2766666667,8885.24333333333,14632.31,10276.2733333333,20610.6233333333,9418.03,15249.3866666667,8299.26,5738.72333333333,11889.2966666667,8550.04666666667,12731.27,39543.2033333333,7987.02,7023.19666666667,9483.25,11752.4,13485.8733333333,10889.3,7595.87,8871.22,4792.59333333333,23283.5366666667,3795.07666666667,4382.65666666667,25129.49,19738.1133333333,26310.0633333333,48791.6966666667,21763.9433333333,7743.48333333333,5838.54333333333,9813.85333333333,13437.51,16397.4033333333,23918.99,14384.2033333333,47658.11,18264.0766666667,19108.4033333333,15582.22,6347.53333333333,66163.9233333333,24747.9933333333,21304.9066666667,17017.56,19805.3633333333,21680.4866666667,19637.8733333333,30065.7366666667,16911.8,8225.40333333333,2266.82666666667,4089.68,12595.0633333333,4427.15,14768.7266666667,22439.9166666667,16828.1833333333,17214.84,6641.60333333333,15320.91,11388.0466666667,28303.14,6587.24666666667,20930.4033333333,12572.7166666667,12067.3133333333,13608.0833333333,18899.01,5381.96333333333,10402.6133333333,9573.00666666667,7259.33666666667,14052.16,4442.46333333333,23533.6666666667,5074.79333333333,20278.1733333333,7353.24,4776.49333333333,16356.87,4228.28,5640.16,10597.8,9405.23333333333,8334.08,16124.1333333333,8513.47333333333,22999.7466666667,22517.7833333333,26059.9433333333,12774.95,26402.5166666667,18364.0633333333,19872.1133333333,11426.7133333333,13631.4166666667,17144.2,15149.36,13130.78,4449.61,34115.5266666667,14699.4066666667,22526.8933333333,17006.48,9939.85,17821.69,19686.06,3442.62,8159.01333333333,5996.72,14811.1366666667,19319.6666666667,9085.81666666667,13794.7433333333,33475.4,21467.1933333333,7241.50666666667,19901.9233333333,17756.76,6144.32,21081.4366666667,19793.4766666667,14633.67,10995.3233333333,41284.3833333333,9315.56333333333,36797.5666666667,4821.96333333333,13822.9466666667,6401.63333333333,27893.44,11179.3633333333,20359.78,11355.7866666667,17500.34,17980.3966666667,13833.5933333333,14080.53,17274.7033333333,78796.9233333333,26411.1,10406.36,11050.86,7947.67666666667,14786.5533333333,15153.12,14119.1033333333,14085.47,16140.9866666667,7125.91,12688.01,9117.32666666667,36179.73,26590.1566666667,10471.6333333333,23803.27,3467.21,5890.06,3720.65,16531.14,35918.02,19803.74,95348.8666666667,26131.1233333333,14261.5633333333,37514.09,15053.4866666667,6078.32333333333,19285.5766666667,14271.8033333333,17679.6533333333,19752.52,11925.4133333333,7403.22,17579.4966666667,20698.76,8551.58333333333,25827.12,26918.0266666667,16309.8333333333,11475.4066666667,18327.2066666667,19672.1266666667,8738.38666666667,14605.5033333333,8055.64666666667,18814.0933333333,26302.3566666667,8745.82666666667,11440.6166666667,13745.0533333333,17775.92,29354.65,18513.3066666667,15571.17,14606.5533333333,10286.88,5770.49,12991.8,16442.62,15244.26,22725.2966666667,19132.59,13249.69,13386.3133333333,22409.0766666667,6227.63666666667,12623.6966666667,8605.15333333333,21098.8666666667,13647.4566666667,4672.56333333333,12538.69,6870.56,10868.1933333333,20503.11,23293.9266666667,22376.9233333333,10614.7466666667,14049.1733333333,5468.37666666667,11950.81,5734.18,4520.24,16418.0266666667,25138.5733333333,35693.3433333333,18035.2366666667,45114.2966666667,5462.52,17011.8433333333,6056.8,9173.75666666667,16416.9533333333,17934.0066666667,3060.27666666667,4856.32333333333,7540.75666666667,17992.9266666667,14537.47,24492.7866666667,16125.52,5859.04333333333,8470.62666666667,16935.1566666667,8931.51666666667,13576.7866666667,2570.15333333333,40772.61,11157.49,7571.58,10606.55,12015.2166666667,3450.9,18314.24,28690.4466666667,11088.6,16322.13,6981.7,7823.05333333333,19333.9233333333,6078.16333333333,14207.3933333333,6404.86333333333,8608.65666666667,11768.7866666667,14000.4033333333,12634.55,10777.74,10649.17,13137.76,23789.5466666667,9871.86333333333,5569.55666666667,6798.90666666667,13293.5566666667,19087.7233333333,7506.30333333333,9609.04333333333,12788.9133333333,13830.5233333333,19195.9833333333,21190.91,16185.79,34778.6366666667,11316.8433333333,10335.6233333333,16263.7333333333,28937.2666666667,3882.33,9852.79,5491.43666666667,21847.0633333333,13448.01,20559.09,5794.74666666667,28573.7666666667,9197.60333333333,24014.67,13262.29,7458.65,9203.53333333333,13655.3066666667,12957.7466666667,16927.12,10315.91,8859.04,7577.87333333333,19921.2666666667,34482.9533333333,8892.06,9032.60666666667,10647.5333333333,11336.06,20622.94,10311.47,11126.2233333333,9606.55333333333,11383.9233333333,21833.31,10500.93,17970.4866666667,12292.0266666667,13243.0833333333,7332.23,14877.7633333333,10880.7833333333,4831.87,12607.3633333333,12977,4615.06,7601.50666666667,12196.1733333333,6453.88,12991.9266666667,24729.48,9263.59666666667,11334.2533333333,8680.32666666667,39122.9433333333,20993.43,8204.76666666667,6409.67333333333,3077.00666666667,19443.4033333333,9549.55333333333,30522.2433333333,23752.1533333333,3128.63333333333,9587.6,18994.0966666667,18874.4433333333,9328.70333333333,11905.5766666667,3751.13666666667,12000.1,14621.2966666667,3582.19333333333,22684.7,5185.24,23413.0366666667,13653.4366666667,2842.73666666667,30237.5033333333,10544.4333333333,10221.8266666667,11249.9933333333,7405.94,6324.88,8153.05666666667,8837.28333333333,12424.84,4838.96666666667,6641.62333333333,10442.8566666667,18905.0366666667,26252.88,50489.8266666667,4444.78333333333,11493.4333333333,12729.5033333333,30271.6966666667,13851.1833333333,8043.97333333333,13522.94,6160.11,14362.2666666667,4432.42,11368.1633333333,11639.34,7653.64666666667,18543.3566666667,9581.69,3239.83,2962.69333333333,11958.0966666667,8540.36333333333,18095.4333333333,13587.9633333333,3336.44,5995.26,4513.42333333333,12012.2933333333,3527.51666666667,10345.9933333333,9468.08666666667,6754.01,9213.11,10081.96,10315.82,16551.4666666667,14007.8633333333,9557.37333333333,6501.52333333333,11979.4933333333,11867.95,17665.4466666667,34226.1033333333,7245.75,9052.25333333333,14736.9133333333,10729.31,6921.91,33721.6966666667,52090.01,18744.5766666667,13614.3333333333,7995.41,11935.54,14024.5833333333,34309.48,6860.92666666667,16576.7833333333,16060.2966666667,14437.4533333333,16653.6466666667,8272.81333333333,16686.9066666667,11184.2,33054.31,7870.26333333333,8939.54,13636.2133333333,15714.7533333333,11415.5633333333,9186.26333333333,11073.9666666667,29883.4066666667,9991.04333333333,6084.36333333333,11771.47,46008.91,15460.39,36778.92,18829.86,3973.93,22501.6266666667,2999.01,27178.6066666667,11633.4566666667,8906.78666666667,18781.4366666667,18898.7433333333,7517.42666666667,8865.35666666667,14152.0566666667,5962.12333333333,21457.5666666667,7065.57333333333,8921.30333333333,11027.86,19147.5366666667,9051.17666666667,5030.79666666667,7188.9,14190.1533333333,23034.9133333333,12881.9733333333,38524.46,20656.2066666667,25619.49,16577.8333333333,18961.0866666667,2973.63666666667,12353.43,14495.0866666667,15073.8433333333,26769.5166666667,14695.9933333333,13957.7,22940.94,9979.62,27146.8633333333,13280.74,14116.95,60702.9533333333,17519.7,11055.2133333333,4409.83333333333,18162.9866666667,6401.63666666667,40373.2066666667,16670.4833333333,4147.53666666667,32391.25,20821.17,12419.1766666667,10258.3066666667,17344.13,14662.3866666667,12011.5666666667,41283.65,9126.67666666667,13314.07,15289.8933333333,8579.00666666667,11373.7633333333,8704.91666666667,15727.4766666667,10620.98,4465.02,5197.88,10751.7966666667,14622.9466666667,37601.6233333333,23536.06,8641.96333333333,9491.79333333333,10041.9033333333,11539.5166666667,9114.58333333333,20349.21,11869.1866666667,11503.4766666667,9235.94,11911.9066666667,24379.0666666667,3289.00666666667,5206.17,7836.29666666667,36158.7633333333,8430.16,14695.7266666667,10989.6466666667,16842.9066666667,23385.0466666667,18509.3133333333,13038.3433333333,6929.98,12053.1966666667,20290.8333333333,13001.4433333333,20055.5066666667,31570.6733333333,37258.1066666667,47935.8233333333,6283.65666666667,28615.6433333333,14622.1233333333,37540.3033333333,17619.9866666667,14606.5433333333,23922.8433333333,12360.6533333333,7279.99666666667,9259.48333333333,27850.55,13821.35,27289.9433333333,7431.38333333333,41940.4566666667,1518.51]
writeLines(vector_2024, "vector_2024.txt")
mean(ingresos_edomex_2024$ing_pc_2024)
## [1] 7330.027

número de personas en pobreza

cat(sum(ingresos_edomex_2024$tot_integ * ingresos_edomex_2024$factor * ingresos_edomex_2024$pobre_extremo))
## 1417752
cat(sum(ingresos_edomex_2024$tot_integ * ingresos_edomex_2024$factor * ingresos_edomex_2024$pobre_moderado))
## 5650189

tasa

tasa_extrema_2024 <- (sum(ingresos_edomex_2024$tot_integ * ingresos_edomex_2024$factor * ingresos_edomex_2024$pobre_extremo))/(sum(ingresos_edomex_2024$tot_integ * ingresos_edomex_2024$factor))
tasa_moderada_2024 <- (sum(ingresos_edomex_2024$tot_integ * ingresos_edomex_2024$factor * ingresos_edomex_2024$pobre_moderado))/(sum(ingresos_edomex_2024$tot_integ * ingresos_edomex_2024$factor))

cat(tasa_extrema_2024)
## 0.07997301
cat(tasa_moderada_2024)
## 0.3187177

pgr

#z_2024 <- 2328.59  # línea de pobreza extrema 
z_2024 <- 4528.71  # línea de pobreza moderada 
# Solo se cuenta la brecha de quienes están por debajo de z; los demás aportan 0
brecha_2024 <- ifelse(ingresos_edomex_2024$ing_pc_2024 < z_2024,
                  (z_2024 - ingresos_edomex_2024$ing_pc_2024) / z_2024,
                  0)
peso_2024 <- ingresos_edomex_2024$tot_integ * ingresos_edomex_2024$factor
n_2024 <- sum(peso_2024)  # población total expandida
pgr_2024 <- sum(brecha_2024 * peso_2024) / n_2024
pgr_2024
## [1] 0.1255379

gini

library(DescTools)
gini_2024 <- Gini(x = ingresos_edomex_2024$ing_pc_2024,
     weights = ingresos_edomex_2024$factor * ingresos_edomex_2024$tot_integ,
     unbiased = FALSE)
gini_2024
## [1] 0.3875158

atkinson

ingresos_edomex_2024_sin0 <- ingresos_edomex_2024[ingresos_edomex_2024$ing_pc_2024 > 0, ]
atkinson_2024 <- Atkinson(x = ingresos_edomex_2024_sin0$ing_pc_2024,
         n = ingresos_edomex_2024_sin0$factor * ingresos_edomex_2024_sin0$tot_integ,
         parameter = 1.5,
         na.rm = TRUE)
atkinson_2024
## [1] 0.299208

2012

conc_2012 <- read_csv("concentradohogar2012.csv")
## Rows: 9002 Columns: 132
## ── Column specification ────────────────────────────────────────────────────────
## Delimiter: ","
## chr   (6): folioviv, ubica_geo, ageb, est_dis, upm, educa_jefe
## dbl (126): foliohog, tam_loc, est_socio, factor_hog, clase_hog, sexo_jefe, e...
## 
## ℹ Use `spec()` to retrieve the full column specification for this data.
## ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.
viv_2012  <- read_csv("viviendas2012.csv") 
## Warning: One or more parsing issues, call `problems()` on your data frame for details,
## e.g.:
##   dat <- vroom(...)
##   problems(dat)
## Rows: 8861 Columns: 58
## ── Column specification ────────────────────────────────────────────────────────
## Delimiter: ","
## chr  (9): folioviv, tipo_viv, num_dueno1, num_dueno2, tipo_finan, ubica_geo,...
## dbl (49): mat_pared, mat_techos, mat_pisos, antiguedad, antigua_ne, cocina, ...
## 
## ℹ Use `spec()` to retrieve the full column specification for this data.
## ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.
names(viv_2012)
##  [1] "folioviv"    "tipo_viv"    "mat_pared"   "mat_techos"  "mat_pisos"  
##  [6] "antiguedad"  "antigua_ne"  "cocina"      "cocina_dor"  "cuart_dorm" 
## [11] "num_cuarto"  "disp_agua"   "dotac_agua"  "excusado"    "uso_compar" 
## [16] "sanit_agua"  "drenaje"     "disp_elect"  "num_focos"   "combustible"
## [21] "estufa_chi"  "eli_basura"  "tenencia"    "num_dueno1"  "hog_dueno1" 
## [26] "num_dueno2"  "hog_dueno2"  "escrituras"  "tipo_finan"  "renta"      
## [31] "estim_pago"  "pago_viv"    "pago_mesp"   "tipo_adqui"  "viv_usada"  
## [36] "lavadero"    "fregadero"   "regadera"    "tinaco_azo"  "cisterna"   
## [41] "pileta"      "calentador"  "medidor_luz" "bomba_agua"  "tanque_gas" 
## [46] "aire_acond"  "calefacc"    "tot_resid"   "tot_hom"     "tot_muj"    
## [51] "tot_hog"     "ubica_geo"   "ageb"        "tam_loc"     "est_socio"  
## [56] "est_dis"     "upm"         "factor_viv"
names(conc_2012)
##   [1] "folioviv"   "foliohog"   "ubica_geo"  "ageb"       "tam_loc"   
##   [6] "est_socio"  "est_dis"    "upm"        "factor_hog" "clase_hog" 
##  [11] "sexo_jefe"  "edad_jefe"  "educa_jefe" "tot_integ"  "hombres"   
##  [16] "mujeres"    "mayores"    "menores"    "p12_64"     "p65mas"    
##  [21] "ocupados"   "percep_ing" "perc_ocupa" "ing_total"  "ing_cor"   
##  [26] "ing_mon"    "trabajo"    "sueldos"    "horas_extr" "comisiones"
##  [31] "otra_rem"   "negocio"    "noagrop"    "industria"  "comercio"  
##  [36] "servicios"  "agrope"     "agricolas"  "pecuarios"  "reproducc" 
##  [41] "pesca"      "otros_trab" "rentas"     "utilidad"   "arrenda"   
##  [46] "transfer"   "jubilacion" "becas"      "donativos"  "remesas"   
##  [51] "bene_gob"   "otros_ing"  "gasto_nom"  "autoconsum" "remu_espec"
##  [56] "transf_esp" "transf_hog" "trans_inst" "estim_alqu" "percep_tot"
##  [61] "percep_mon" "retiro_inv" "prestamos"  "otras_perc" "erogac_nom"
##  [66] "gasto_tot"  "gasto_cor"  "gasto_mon"  "alimentos"  "ali_dentro"
##  [71] "cereales"   "carnes"     "pescado"    "leche"      "huevo"     
##  [76] "aceites"    "tuberculo"  "verduras"   "frutas"     "azucar"    
##  [81] "cafe"       "especias"   "otros_alim" "bebidas"    "ali_fuera" 
##  [86] "tabaco"     "vesti_calz" "vestido"    "calzado"    "vivienda"  
##  [91] "alquiler"   "pred_cons"  "agua"       "energia"    "limpieza"  
##  [96] "cuidados"   "utensilios" "enseres"    "salud"      "atenc_ambu"
## [101] "hospital"   "medicinas"  "transporte" "publico"    "foraneo"   
## [106] "adqui_vehi" "mantenim"   "refaccion"  "combus"     "comunica"  
## [111] "educa_espa" "educacion"  "esparci"    "paq_turist" "personales"
## [116] "cuida_pers" "acces_pers" "otros_gas"  "transf_gas" "erogac_tot"
## [121] "erogac_mon" "cuota_viv"  "mater_serv" "material"   "servicio"  
## [126] "deposito"   "prest_terc" "pago_tarje" "deudas"     "balance"   
## [131] "otras_erog" "smg"
# 2. Entidad federativa
conc_2012$entidad <- substr(conc_2012$folioviv, 1, 2)
inpc_base_2012   <- 100.000   # base de la ENIGH 2012
inpc_dic_2024    <- 135.38   # AVG de todo el año . 
factor_defl_2012      <- inpc_dic_2024 / inpc_base_2012

# 4. Vector de dos columnas: ingreso nominal (2018) e ingreso deflactado (2024)
ingresos_hogar_2012 <- data.frame(
  folioviv       = conc_2012$folioviv,
  foliohog       = conc_2012$foliohog,
  entidad        = conc_2012$entidad,
  tot_integ      = conc_2012$tot_integ,
  ing_cor_2012   = conc_2012$ing_cor/3,          # nominal, como lo reporta la encuesta
  ing_cor_2024   = (conc_2012$ing_cor/3) * factor_defl_2012,  # deflactado a precios dic-2024
  ing_pcp = ((conc_2012$ing_cor/3) * factor_defl_2012) / conc_2012$tot_integ
)
ingresos_hogar_2012$pobre_extremo <- as.integer(ingresos_hogar_2012$ing_pcp < 1108.57)
ingresos_hogar_2012$pobre_moderado <- as.integer(
  ingresos_hogar_2012$ing_pcp >= 1108.57 & ingresos_hogar_2012$ing_pcp < 2568.83
)
ingresos_hogar_2012 <- merge(
  ingresos_hogar_2012,
  viv_2012[, c("folioviv", "factor_viv")],
  by = "folioviv",
  all.x = TRUE
)

# 5. Solo Estado de México (clave 15)
ingresos_edomex_2012 <- ingresos_hogar_2012[ingresos_hogar_2012$entidad == "15", ]
head(ingresos_edomex_2012)
##      folioviv foliohog entidad tot_integ ing_cor_2012 ing_cor_2024   ing_pcp
## 3959   151001        0      15         5     3473.023     4701.779  940.3558
## 3960   151003        0      15         2     8023.777    10862.589 5431.2944
## 3961   151004        0      15         8     2019.370     2733.823  341.7279
## 3962   151005        0      15         6     5534.053     7492.001 1248.6669
## 3963   151006        0      15         7     5111.610     6920.098  988.5854
## 3964   151007        0      15         8     8253.963    11174.216 1396.7769
##      pobre_extremo pobre_moderado factor_viv
## 3959             1              0       8470
## 3960             0              0       8470
## 3961             1              0       8470
## 3962             0              1       8470
## 3963             1              0       8470
## 3964             0              1       8470
mean(ingresos_edomex_2012$ing_pcp)
## [1] 5163.922
vector_2012 <- paste0(
  "[",
  paste(ingresos_edomex_2012$ing_cor_2012, collapse = ","),
  "]"
)
vector_2012
## [1] "[3473.02333333333,8023.77666666667,2019.37,5534.05333333333,5111.61,8253.96333333333,2665.38333333333,2550.03333333333,2094.91333333333,9146.78666666667,3524.32333333333,10326.48,14230.8333333333,17965.6133333333,3713.18,11081.7966666667,7864.08666666667,10908.42,18883.9166666667,12363.0566666667,8612.35,6459.54,13388.4333333333,8845.18666666667,35007.17,17495.7233333333,24857.96,3007.99,7181.81,11159.58,6065.21666666667,19648.2933333333,8624.99666666667,19332.7433333333,26713.3433333333,157517.34,5747.8,1864.47666666667,3417.83666666667,29458.4433333333,11984.3433333333,11290.6766666667,2131.89666666667,13613.8566666667,4007.34,5704.33333333333,5377.48666666667,1880.77666666667,20777.6266666667,1787.38,3331.07,2908.22666666667,4230.66333333333,3613.92666666667,18565.9166666667,7722.30333333333,8019.36666666667,4738.44666666667,4711.95,6936.18333333333,14807.3366666667,5273.14,1956.51666666667,2679.69333333333,9945.64666666667,6837.30333333333,15231.73,6748.4,5869.56333333333,16846.33,5282.60333333333,5076.96,5416.29666666667,7846.07,9214.14333333333,32338.75,38264.8733333333,57806.7966666667,22599.9266666667,24075.0166666667,13832.1766666667,10748.2433333333,7162.43666666667,7557.76,5253.67666666667,5478.26,6488.97666666667,2479.51666666667,16816.13,10586.9133333333,41309.5766666667,16474.2866666667,16907.3633333333,9557.04333333333,11739.13,13533.1266666667,7163.39,7721.69666666667,6867.71333333333,6531.04,13080.8366666667,6652.17333333333,11440.4633333333,8740.17666666667,12322.8366666667,12670.62,9420.22666666667,13701.7833333333,26348.71,48042.8333333333,16038.4166666667,24886.7366666667,2013.57666666667,6499.39666666667,16985.65,8862.95666666667,6368.69,15192.7333333333,3173.76666666667,5158.07666666667,3596.4,2655.56666666667,3714.98,2071.3,2223.69,16455.49,32262.2933333333,8644.53666666667,6921.88333333333,4509.82333333333,7950.89,4169.66333333333,14573.7633333333,31731.57,19321.3066666667,6828.63666666667,6066.42333333333,3626.23333333333,967.21,5290.39666666667,6153.03,18548.7933333333,5317.78666666667,4342.61666666667,5458.54,9021.30666666667,7868.84666666667,14711.3566666667,11945.6933333333,8207.95333333333,6885.24333333333,27491.7933333333,33459.01,46381.1433333333,21295.07,8826.22,7194.37333333333,6475.08,6401.63333333333,11196.7166666667,2481.92333333333,1999.99666666667,4303.27,2847.76666666667,6422.94,1841.08333333333,6615.89,6481.49333333333,1549.17666666667,131097.366666667,6118.07,5428.67666666667,6310.03,6403.74,27142.8466666667,4442.62,18157.85,7404.51,3286.88333333333,8893.43333333333,1355.73333333333,4663.76333333333,5173.16333333333,17713.1033333333,5114.75,5016.39,8990.89,9483.6,10975.4066666667,13114.7433333333,10819.67,12810.4666666667,9072.23,5614.76,18066.4033333333,4264.95333333333,13847.5333333333,19032.78,21131.5966666667,7815.92333333333,8856.19666666667,2572.36,6579.61,8278.68666666667,19922.2833333333,5846.59333333333,17742.9266666667,11903.33,4652.65333333333,7740.76333333333,7035.12,4951.43333333333,8274.32666666667,8705.84333333333,13620.2466666667,8886.69333333333,17651.62,5311.47333333333,5213.11,4478.68,9849.17333333333,42224.8,4934.41666666667,6891.97333333333,9940.17,6104.66666666667,18050.8066666667,16016.3833333333,12512.87,8956.76666666667,3417.79333333333,8416.39,15682.7833333333,8488.52333333333,1983.60333333333,3950.81666666667,3554.09666666667,4298.36,3850.81666666667,25459.01,16929.9533333333,11352.4566666667,8295.08,23770.49,11080.66,8402.52,6702.32,2243.09666666667,3993.1,11416.4933333333,2922.03,8630.77666666667,8262.76333333333,8652.86666666667,3092.56333333333,9186.85333333333,17196.8733333333,9243.49,3515.29666666667,10859.62,8625.52,9359.21666666667,7321.17333333333,18717.9066666667,7712.47666666667,12954.5866666667,12115.7,22941.7866666667,21888.5233333333,31792.47,7611.49666666667,5359.39333333333,11376.2266666667,4298.24333333333,16835.65,8489.03666666667,12609.7766666667,6248.24333333333,38870.0833333333,12208.74,2702.07666666667,9086.35,30847.8166666667,13122.42,5745.44,35217.3733333333,20284.0433333333,13147.3533333333,8137.83,11136.3033333333,7626.88666666667,6275.17,2187.37,5817.54333333333,6234.3,4345.99666666667,3860.1,16917.4366666667,14896.3066666667,7096.91,8449.33333333333,8390.13333333333,89115.36,44002.79,6440.21666666667,6031.62666666667,7840.39,9032.34,5663.39,3786.02666666667,4670.45666666667,4372.61333333333,21357.8566666667,8164.88666666667,15145.8866666667,16196.5133333333,6004.56,5282.60333333333,3418.65,8481.24,9185.12666666667,12589.0533333333,4452.30666666667,11580.4166666667,10255.95,6054.69333333333,8596.06666666667,13553.7433333333,7612.59666666667,10569.59,8804.34666666667,7218.09,13631.48,9603.25666666667,10553.4866666667,12539.09,4391.65,5469.84333333333,3418.65,2934.78,10266.4766666667,7000.34666666667,4096.07,3706.87,12544.98,9886.95,4435.30333333333,18253.2366666667,4378.27333333333,11175.7866666667,4934.07333333333,2104.78666666667,12370.2633333333,22888.29,11000.9533333333,9847.82333333333,9350.01333333333,8577.68666666667,4152.54,12980.0866666667,7915.11666666667,3042.38666666667,7949.52,5913.74,4695.65,10471.0666666667,8804.34333333333,4495.55333333333,5469.84333333333,6553.05333333333,4880.78333333333,4157.60666666667]"
vector_2012_2024precios <- paste0(
  "[",
  paste(ingresos_edomex_2012$ing_cor_2024, collapse = ","),
  "]"
)
cat(vector_2012_2024precios)
## [4701.77898866667,10862.5888513333,2733.823106,7492.00140266667,6920.097618,11174.2155606667,3608.39595666667,3452.23512666667,2836.09367066667,12382.9197893333,4771.22892866667,13979.988624,19265.7021666667,24321.8473306667,5026.903084,15002.5363273333,10646.4005293333,14767.818996,25565.0463833333,16737.1061153333,11659.39943,8744.925252,18125.2610466667,11974.6137093333,47392.706746,23685.7102486667,33652.706248,4072.216862,9722.734378,15107.839404,8211.09032333333,26599.8595146667,11676.5204873333,26172.6679246667,36164.5242046667,213246.974892,7781.37164,2524.12851133333,4627.06727933333,39880.8405846667,16224.4040046667,15285.3180713333,2886.16170733333,18430.4391553333,5425.136892,7722.52646666667,7280.04144933333,2546.19545133333,28128.7509813333,2419.755044,4509.602566,3937.15726133333,5727.47202066667,4892.53392133333,25134.5379833333,10454.4542526667,10856.6185933333,6414.90909733333,6379.03791,9390.20499666667,20046.1723793333,7138.776932,2648.73226333333,3627.76883466667,13464.4164573333,9256.34125266667,20620.716074,9135.98392,7946.21484066666,22806.561554,7151.58839266667,6873.188448,7332.58242733333,10622.009566,12474.1072446667,43780.19975,51802.9855186667,78258.8413273333,30595.7807213333,32592.7575633333,18726.0007713333,14550.9718246667,9696.50675933333,10231.695488,7112.42747133333,7416.468388,8784.77661133333,3356.76966333333,22765.676794,14332.5632706667,55924.9048913333,22302.8892893333,22889.1884806667,12938.3252646667,15892.434194,18321.1468813333,9697.797382,10453.6329473333,9297.51031066667,8841.721952,17708.8366793333,9005.71225866667,15488.0992606667,11832.4511713333,16682.6562793333,17153.485356,12753.1028613333,18549.4742766667,35670.883598,65040.3877666667,21712.8084833333,33691.6640993333,2725.98009133333,8798.88320733333,22995.17297,11998.6707353333,8621.932522,20567.9223866667,4296.64531333333,6983.00419133333,4868.80632,3595.10615333333,5029.339924,2804.12594,3010.431522,22277.442362,43676.6927146667,11702.9737393333,9370.84565666667,6105.39882866667,10763.914882,5644.89022066667,19729.9608006667,42958.199466,26157.1849653333,9244.60831933333,8212.72390866667,4909.19468666667,1309.408898,7162.13900733333,8329.972014,25111.3564146667,7199.21958933333,5879.03444333333,7389.771452,12213.0449653333,10652.8446173333,19916.2346553333,16172.0796346667,11111.9272226667,9321.24242466667,37218.3898146667,45296.807738,62790.7918446667,28829.265766,11948.936636,9739.74261866667,8765.963304,8666.53120666667,15158.1150233333,3360.02780866667,2707.59548733333,5825.766926,3855.30651333333,8695.376172,2492.45861666667,8956.591882,8774.64567466666,2097.27537133333,177479.614993333,8282.643166,7349.34247133333,8542.518614,8669.383212,36745.9858173333,6014.418956,24582.09733,10024.225638,4449.78265666667,12039.9300466667,1835.39178666667,6313.80280066667,7003.42852066667,23979.9992926667,6924.34855,6791.188782,12171.866882,12838.89768,14858.5055453333,17754.7395246667,14647.669246,17342.8097733333,12281.984974,7601.262088,24458.2968326667,5773.89382266667,18746.7906266667,25766.577564,28607.9555673333,10581.1970086667,11989.5190473333,3482.460968,8907.476018,11207.6860093333,26970.7871766667,7915.11805466667,24020.3741213333,16114.728154,6298.76208266667,10479.4454006667,9524.145456,6703.25044666667,11201.7834413333,11785.9707046667,18439.0899373333,12030.8054346667,23896.763156,7190.67259866667,7057.508318,6063.236984,13333.8108586667,57163.93424,6680.21328333333,9330.35349866667,13457.002146,8264.49773333333,24437.1820653333,21682.9797566667,16939.923406,12125.6707133333,4627.00861466667,11394.108782,21231.3520766667,11491.7628886667,2685.40219266667,5348.61560333333,4811.53606733333,5819.119768,5213.23560333333,34466.407738,22919.7708226667,15368.9558353333,11229.879304,32180.489362,15000.997508,11375.331576,9073.600816,3036.70426733333,5405.85878,15455.6486746667,3955.844214,11684.3454513333,11186.1290006667,11714.2508933333,4186.71224066667,12437.1620426667,23281.1271186667,12513.836762,4759.00862733333,14701.753556,11677.228976,12670.5075233333,9911.40445866667,25340.3020453333,10441.1509113333,17537.9194293333,16402.23466,31058.5907893333,29632.6828886667,43040.645886,10304.4441873333,7255.54669466667,15401.1356613333,5818.96182466667,22792.10297,11492.4578393333,17071.1156513333,8458.87182466666,52622.3188166667,16528.192212,3658.07139133333,12301.10063,41761.7742033333,17765.132196,7778.176672,47677.2800186667,27460.5378646667,17798.8869426667,11016.994254,15076.3274526667,10325.2791693333,8495.325146,2961.261506,7875.79016466667,8439.99534,5883.61028733333,5225.80338,22902.8257593333,20166.6199653333,9607.796758,11438.7074666667,11358.5625066667,120644.374368,59570.977102,8718.76532333333,8165.61618133333,10614.319982,12227.981892,7667.097382,5125.52290133333,6322.86423533333,5919.64393066667,28914.2663553333,11053.6235693333,20504.5013693333,21926.8397506667,8128.973328,7151.58839266667,4628.16837,11481.902712,12434.8244813333,17043.0604026667,6027.53276533333,15677.5680833333,13884.50511,8196.84383466667,11637.3550533333,18349.0577246667,10305.9333673333,14309.110942,11919.3245173333,9771.850242,18454.297624,13000.8888753333,14287.3102493333,16975.420042,5945.41577,7405.07390466667,4628.16837,3973.105164,13898.7561113333,9477.06931733333,5545.259566,5018.360606,16983.393924,13384.95291,6004.51365266667,24711.2317993333,5927.30643866667,15129.7799893333,6679.74847866667,2849.46018933333,16746.8625006667,30986.167002,14893.0906226667,13331.9832286667,12658.0480506667,11612.4722093333,5621.708652,17572.4413293333,10715.4849433333,4118.78306933333,10762.060176,8006.021212,6356.97097,14175.7300533333,11919.3200046667,6086.08010266667,7405.07390466667,8871.52360266667,6607.60447666667,5628.56790533333]
writeLines(vector_2012_2024precios, "vector_2012_a_precios_2024.txt")
writeLines(vector_2012, "vector_2012.txt")

número de personas en pobreza

cat(sum(ingresos_edomex_2012$tot_integ * ingresos_edomex_2012$factor_viv * ingresos_edomex_2012$pobre_extremo))
## 1223083
cat(sum(ingresos_edomex_2012$tot_integ * ingresos_edomex_2012$factor_viv * ingresos_edomex_2012$pobre_moderado))
## 5310334

tasa

tasa_extrema_2012 <- (sum(ingresos_edomex_2012$tot_integ * ingresos_edomex_2012$factor_viv * ingresos_edomex_2012$pobre_extremo))/(sum(ingresos_edomex_2012$tot_integ * ingresos_edomex_2012$factor_viv))
tasa_moderada_2012 <- (sum(ingresos_edomex_2012$tot_integ * ingresos_edomex_2012$factor_viv * ingresos_edomex_2012$pobre_moderado))/(sum(ingresos_edomex_2012$tot_integ * ingresos_edomex_2012$factor_viv))

cat(tasa_extrema_2012)
## 0.07557985
cat(tasa_moderada_2012)
## 0.3281497

pgr

z_2012 <- 1108.57  # línea de pobreza extrema 
#z_2012 <- 2568.83  # línea de pobreza moderada 
# Solo se cuenta la brecha de quienes están por debajo de z; los demás aportan 0
brecha_2012 <- ifelse(ingresos_edomex_2012$ing_pcp < z_2012,
                  (z_2012 - ingresos_edomex_2012$ing_pcp) / z_2012,
                  0)
peso_2012 <- ingresos_edomex_2012$tot_integ * ingresos_edomex_2012$factor_viv
n_2012 <- sum(peso_2012)  # población total expandida
pgr_2012 <- sum(brecha_2012 * peso_2012) / n_2012
pgr_2012
## [1] 0.02152469

gini

library(DescTools)
gini_2012 <- Gini(x = ingresos_edomex_2012$ing_pcp,
     weights = ingresos_edomex_2012$factor_viv * ingresos_edomex_2012$tot_integ,
     unbiased = FALSE)
gini_2012
## [1] 0.4758925

atkinson

atkinson_2012 <- Atkinson(x = ingresos_edomex_2012$ing_pcp,
         n = ingresos_edomex_2012$factor_viv * ingresos_edomex_2012$tot_integ,
         parameter = 1.5, # aversa a la desigualdad
         na.rm = TRUE)
atkinson_2012
## [1] 0.4234745

2006

library(foreign) 
## Warning: package 'foreign' was built under R version 4.4.3
conc_2006 <- read.dbf("concen.dbf", as.is = TRUE)
viv_2006  <- read.dbf("hogares.dbf", as.is = TRUE)
names(viv_2006)
##   [1] "FOLIO"      "ESTRATO"    "CONAPO"     "UBICA_GEO"  "MUROS01"   
##   [6] "TECHOS02"   "TECHOS03"   "PISOS04"    "ANTIGUA"    "CUA_COC"   
##  [11] "COC_DUER"   "CUA_DOR"    "NUM_CUA"    "COMBUS10"   "COMBUS11"  
##  [16] "TENENCIA12" "TENENCIA13" "TENENCIA14" "AGUA15"     "AGUA16"    
##  [21] "BANO17"     "BANO18"     "BANO19"     "BANO20"     "DRENAJE21" 
##  [26] "LUZ22"      "NO_FOCOS"   "BASURA24"   "BASURA25"   "EQV26_1"   
##  [31] "EQV26_2"    "EQV26_3"    "EQV26_4"    "EQV26_5"    "EQV26_6"   
##  [36] "EQV26_7"    "EQV26_8"    "EQV26_9"    "EQV26_10"   "EQV26_11"  
##  [41] "EQV26_12"   "EQV26_13"   "RESIDENTES" "HOGARES2"   "HOGARES3"  
##  [46] "HUESPEDES4" "HUESPEDES5" "HUESPEDES6" "DOMESTICO7" "DOMESTICO8"
##  [51] "DOMESTICO9" "SERV01_1"   "SERV01_2"   "SERV01_3"   "SERV01_4"  
##  [56] "SERV02_1"   "SERV02_2"   "SERV02_3"   "SERV02_4"   "VEHI04_1"  
##  [61] "VEHI04_2"   "VEHI04_3"   "VEHI04_4"   "VEHI04_5"   "VEHI04_6"  
##  [66] "VEHI04_7"   "VEHI04_8"   "VEHI04_9"   "VEHI05_1"   "VEHI05_2"  
##  [71] "VEHI05_3"   "VEHI05_4"   "VEHI05_5"   "VEHI05_6"   "VEHI05_7"  
##  [76] "VEHI05_8"   "VEHI05_9"   "EQH07_01"   "EQH07_02"   "EQH07_03"  
##  [81] "EQH07_04"   "EQH07_05"   "EQH07_06"   "EQH07_07"   "EQH07_08"  
##  [86] "EQH07_09"   "EQH07_10"   "EQH07_11"   "EQH07_12"   "EQH07_13"  
##  [91] "EQH07_14"   "EQH07_15"   "EQH07_16"   "EQH07_17"   "EQH07_18"  
##  [96] "EQH07_19"   "EQH07_20"   "EQH07_21"   "EQH07_22"   "EQH07_23"  
## [101] "EQH07_24"   "EQH07_25"   "EQH07_26"   "EQH07_27"   "EQH07_28"  
## [106] "EQH07_29"   "EQH07_30"   "EQH07_31"   "EQH07_32"   "EQH07_33"  
## [111] "EQH07_34"   "EQH08_01"   "EQH08_02"   "EQH08_03"   "EQH08_04"  
## [116] "EQH08_05"   "EQH08_06"   "EQH08_07"   "EQH08_08"   "EQH08_09"  
## [121] "EQH08_10"   "EQH08_11"   "EQH08_12"   "EQH08_13"   "EQH08_14"  
## [126] "EQH08_15"   "EQH08_16"   "EQH08_17"   "EQH08_18"   "EQH08_19"  
## [131] "EQH08_20"   "EQH08_21"   "EQH08_22"   "EQH08_23"   "EQH08_24"  
## [136] "EQH08_25"   "EQH08_26"   "EQH08_27"   "EQH08_28"   "EQH08_29"  
## [141] "EQH08_30"   "EQH08_31"   "EQH08_32"   "EQH08_33"   "EQH08_34"  
## [146] "TIENE_TBM"  "TIENE_TCM"  "TIENE_TBT"  "TIENE_TCT"  "TIENE_TBZ" 
## [151] "TIENE_TCZ"  "TIENE_TBSE" "TIENE_TCSE" "PAGO_TCRE"  "CLAVE_TB01"
## [156] "PAGO_TCOM"  "CLAVE_TC01" "FACTOR"
# 2. Entidad federativa
conc_2006$entidad <- substr(conc_2006$ubica_geo, 1, 2)
table(conc_2006$entidad)
## 
##   01   02   03   04   05   06   07   08   09   10   11   12   13   14   15   16 
##  378  559  312  477  596  412  713  850 1135  556 1889  762  541  657 1161  592 
##   17   18   19   20   21   22   23   24   25   26   27   28   29   30   31   32 
##  515  335  487  784  618  418  427  521  506  682  458  635  359 1745  431  364
inpc_base_2006   <- 100.000
inpc_dic_2024    <- 135.38
factor_defl_2006      <- inpc_dic_2024 / inpc_base_2006
ingresos_hogar_2006 <- data.frame(
  folio        = conc_2006$folio,
  entidad      = conc_2006$entidad,
  tot_integ = conc_2006$tot_resi,
  ing_cor_2006 = conc_2006$ingcor/3,
  ing_cor_2024 = (conc_2006$ingcor/3) * factor_defl_2006,
  ing_pcp      = ((conc_2006$ingcor/3) * factor_defl_2006) / conc_2006$tot_resi
)
ingresos_hogar_2006$pobre_extremo <- as.integer(ingresos_hogar_2006$ing_pcp < 786.58)
ingresos_hogar_2006$pobre_moderado <- as.integer(
  ingresos_hogar_2006$ing_pcp >= 786.58 & ingresos_hogar_2006$ing_pcp < 2013.21
)
ingresos_hogar_2006 <- merge(
  ingresos_hogar_2006,
  viv_2006[, c("FOLIO", "FACTOR")],
  by.x = "folio",
  by.y = "FOLIO",
  all.x = TRUE
)

# Solo Estado de México (clave 15)
ingresos_edomex_2006 <- ingresos_hogar_2006[ingresos_hogar_2006$entidad == "15", ]
head(ingresos_edomex_2006)
##            folio entidad tot_integ ing_cor_2006 ing_cor_2024   ing_pcp
## 9838 20061510010      15         1     8307.180    11246.260 11246.260
## 9839 20061510020      15         4     4507.617     6102.411  1525.603
## 9840 20061510030      15         4    10032.437    13581.913  3395.478
## 9841 20061510040      15         1     1965.597     2661.025  2661.025
## 9842 20061510050      15         2     5750.867     7785.523  3892.762
## 9843 20061510060      15         2     2456.583     3325.723  1662.861
##      pobre_extremo pobre_moderado FACTOR
## 9838             0              0    164
## 9839             0              1    164
## 9840             0              0    164
## 9841             0              0    164
## 9842             0              0    164
## 9843             0              1    164
mean(ingresos_edomex_2006$ing_pcp)
## [1] 4069.186

número de personas en pobreza

cat(sum(ingresos_edomex_2006$tot_integ * ingresos_edomex_2006$FACTOR * ingresos_edomex_2006$pobre_extremo))
## 365860
cat(sum(ingresos_edomex_2006$tot_integ * ingresos_edomex_2006$FACTOR * ingresos_edomex_2006$pobre_moderado))
## 5254731

tasa

tasa_extrema_2006 <- (sum(ingresos_edomex_2006$tot_integ * ingresos_edomex_2006$FACTOR * ingresos_edomex_2006$pobre_extremo))/(sum(ingresos_edomex_2006$tot_integ * ingresos_edomex_2006$FACTOR))
tasa_moderada_2006 <- (sum(ingresos_edomex_2006$tot_integ * ingresos_edomex_2006$FACTOR * ingresos_edomex_2006$pobre_moderado))/(sum(ingresos_edomex_2006$tot_integ * ingresos_edomex_2006$FACTOR))

cat(tasa_extrema_2006)
## 0.02510116
cat(tasa_moderada_2006)
## 0.3605199

pgr

z_2006 <- 786.58  # línea de pobreza extrema 
#z_2006 <- 2013.21  # línea de pobreza moderada 
# Solo se cuenta la brecha de quienes están por debajo de z; los demás aportan 0
brecha_2006 <- ifelse(ingresos_edomex_2006$ing_pcp < z_2006,
                  (z_2006 - ingresos_edomex_2006$ing_pcp) / z_2006,
                  0)
peso_2006 <- ingresos_edomex_2006$tot_integ * ingresos_edomex_2006$FACTOR
n_2006 <- sum(peso_2006)  # población total expandida
pgr_2006 <- sum(brecha_2006 * peso_2006) / n_2006
pgr_2006
## [1] 0.005680538

gini

library(DescTools)
gini_2006 <- Gini(x = ingresos_edomex_2006$ing_pcp,
     weights = ingresos_edomex_2006$FACTOR * ingresos_edomex_2006$tot_integ,
     unbiased = FALSE)
gini_2006
## [1] 0.4458356

Atkinson

atkinson_2006 <- Atkinson(x = ingresos_edomex_2006$ing_pcp,
         n = ingresos_edomex_2006$FACTOR * ingresos_edomex_2006$tot_integ,
         parameter = 1.5, # aversa a la desigualdad
         na.rm = TRUE)
atkinson_2006
## [1] 0.3698087
summary(ingresos_edomex_2006$tot_integ)
##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
##     1.0     3.0     4.0     4.2     5.0    15.0
names(conc_2006)
##   [1] "folio"      "hog"        "estrato"    "ubica_geo"  "tenencia"  
##   [6] "clase_hog"  "edad"       "ed_formal"  "tam_hog"    "hombres"   
##  [11] "mujeres"    "tot_resi"   "tot_hom"    "tot_muj"    "mayores"   
##  [16] "menores"    "p0a11"      "p12_64"     "p65mas"     "n_ocup"    
##  [21] "pering"     "perocu"     "ingtot"     "ingcor"     "ingmon"    
##  [26] "trabajo"    "sueldos"    "otra_rem"   "negocio"    "no_agrop"  
##  [31] "agrope"     "coopera"    "socie"      "rentas"     "transfer"  
##  [36] "jubila"     "beca_don"   "donativo"   "otros"      "gasnom"    
##  [41] "auto"       "pago"       "rega"       "esti"       "pertot"    
##  [46] "permon"     "retiro"     "prestamo"   "otras_pr"   "eronom"    
##  [51] "gastot"     "gascor"     "gasmon"     "alimentos"  "ali_dent"  
##  [56] "cereales"   "carnes"     "pescado"    "leche"      "huevo"     
##  [61] "aceites"    "tuberculo"  "verduras"   "frutas"     "azucar"    
##  [66] "cafe"       "especias"   "otro_ali"   "bebidas"    "fuera_hog" 
##  [71] "tabaco"     "vestido"    "vestido_c"  "ves_3ymas"  "ves_3men"  
##  [76] "calzado"    "vivienda"   "alquiler"   "agua"       "energia"   
##  [81] "limpieza"   "cuidados"   "enseres"    "salud"      "aten_pri"  
##  [86] "hospital"   "medica"     "transporte" "trans_ad"   "publico"   
##  [91] "foraneo"    "vehiculo"   "comunica"   "educacion"  "educa"     
##  [96] "esparci"    "personal"   "cuidado"    "accesorio"  "transfe"   
## [101] "erotot"     "eromon"     "cuota"      "material"   "servicio"  
## [106] "deposito"   "terceros"   "pago_tar"   "deudas"     "monedas"   
## [111] "casas"      "balance"    "otra_ero"   "smg"        "est_dis"   
## [116] "upm"        "entidad"

library(dplyr)
## Warning: package 'dplyr' was built under R version 4.4.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(tidyr)
## Warning: package 'tidyr' was built under R version 4.4.3
library(ggplot2)
## Warning: package 'ggplot2' was built under R version 4.4.3
# ---------------------------------------------------
# (1) Tasas de crecimiento por decil
# ---------------------------------------------------
nivel_por_decil <- function(x, peso) {
  ord <- order(x)
  x <- x[ord]; peso <- peso[ord]
  peso_acum <- cumsum(peso) / sum(peso)
  decil <- ceiling(peso_acum * 10)
  decil[decil > 10] <- 10
  sapply(1:10, function(d) weighted.mean(x[decil == d], peso[decil == d]))
}

decil_2006 <- nivel_por_decil(ingresos_edomex_2006$ing_pcp, peso_2006)
decil_2012 <- nivel_por_decil(ingresos_edomex_2012$ing_pcp, peso_2012)
decil_2018 <- nivel_por_decil(ingresos_edomex_2018$ing_pc_2024, peso_2018)
decil_2024 <- nivel_por_decil(ingresos_edomex_2024$ing_pc_2024, peso_2024)

data <- data.frame(
  Decil = 1:10,
  C_2006 = decil_2006,
  C_2012 = decil_2012,
  C_2018 = decil_2018,
  C_2024 = decil_2024
)

data <- data %>%
  mutate(Growth06_12 = (C_2012 - C_2006) / C_2006 * 100,
         Growth12_18 = (C_2018 - C_2012) / C_2012 * 100,
         Growth18_24 = (C_2024 - C_2018) / C_2018 * 100)

# ---------------------------------------------------
# (2) Tasas de crecimiento de la pobreza 
# ---------------------------------------------------
GrowthPov <- c(
  (tasa_extrema_2012 - tasa_extrema_2006) / tasa_extrema_2006 * 100,
  (tasa_extrema_2018 - tasa_extrema_2012) / tasa_extrema_2012 * 100,
  (tasa_extrema_2024 - tasa_extrema_2018) / tasa_extrema_2018 * 100
)
names(GrowthPov) <- c("GPov06_12", "GPov12_18", "GPov18_24")

# ---------------------------------------------------
# (3) Elasticidad pobreza-crecimiento
# ---------------------------------------------------
media_ingreso <- c(mean(ingresos_edomex_2006$ing_pcp),
                    mean(ingresos_edomex_2012$ing_pcp),
                    mean(ingresos_edomex_2018$ing_pc_2024),
                    mean(ingresos_edomex_2024$ing_pc_2024))

GrowthMeans <- c(
  (media_ingreso[2] - media_ingreso[1]) / media_ingreso[1] * 100,
  (media_ingreso[3] - media_ingreso[2]) / media_ingreso[2] * 100,
  (media_ingreso[4] - media_ingreso[3]) / media_ingreso[3] * 100
)

Elasticity_Pov_Growth <- GrowthPov / GrowthMeans
names(Elasticity_Pov_Growth) <- names(GrowthPov)

# ---------------------------------------------------
# (4) Elasticidad de la pobreza respecto a la desigualdad 
# ---------------------------------------------------
gini_vec <- c(gini_2006, gini_2012, gini_2018, gini_2024)

GrowthGini <- c(
  (gini_vec[2] - gini_vec[1]) / gini_vec[1] * 100,
  (gini_vec[3] - gini_vec[2]) / gini_vec[2] * 100,
  (gini_vec[4] - gini_vec[3]) / gini_vec[3] * 100
)

Elasticity_Pov_Gini <- GrowthPov / GrowthGini
names(Elasticity_Pov_Gini) <- names(GrowthPov)

# ---------------------------------------------------
# (5) Comparar crecimiento del quintil más bajo vs promedio (RAZÓN, no diferencia)
# Quintil 1 = deciles 1 y 2
# ---------------------------------------------------
quintil1 <- data %>%
  filter(Decil %in% c(1, 2)) %>%
  summarise(across(starts_with("Growth"), mean))

promedio <- data %>%
  summarise(across(starts_with("Growth"), mean))

Comparacion <- quintil1 / promedio

# ---------------------------------------------------
# (6) Curva de Incidencia del Crecimiento (CIC) — los 3 sexenios juntos
# ---------------------------------------------------
data_long <- data %>%
  select(Decil, starts_with("Growth")) %>%
  pivot_longer(-Decil, names_to = "Periodo", values_to = "Crecimiento")

ggplot(data_long, aes(x = Decil, y = Crecimiento, color = Periodo)) +
  geom_line(linewidth = 1.2) +
  geom_point() +
  scale_x_continuous(breaks = 1:10) +
  labs(title = "Curva de Incidencia del Crecimiento (CIC) — Estado de México",
       x = "Decil de ingreso",
       y = "Tasa de crecimiento (%)") +
  theme_minimal()

# ---------------------------------------------------
# Ver resultados
# ---------------------------------------------------
GrowthPov
## GPov06_12 GPov12_18 GPov18_24 
## 201.10108 -42.41523  83.75101
Elasticity_Pov_Growth
## GPov06_12 GPov12_18 GPov18_24 
##  7.475021 -3.643246  3.085374
Elasticity_Pov_Gini
##  GPov06_12  GPov12_18  GPov18_24 
##  29.829435   2.577002 -33.134241
Comparacion
##   Growth06_12 Growth12_18 Growth18_24
## 1   0.3307527    1.539833    1.080967