Данные "Inf" при преобразовании данных в таблицу жизни

Question

Я начинающий программист, так что наберитесь терпения.Я скачал некоторые данные о смертности из «Базы данных о человеческой смертности (HMD)».

Необходимые пакеты:

Packages <- c("StMoMo", "lifecontingencies", "fda", "demography")
invisible(lapply(Packages, library, character.only = TRUE))
remove(Packages)

Загрузка наборов данных:

# mortality
Germany_M <- hmd.mx("DEUTNP", "USERNAME", "PASSWORD", "Germany")
# life expectancy at birth
Germany_E <- hmd.e0("DEUTNP", "USERNAME", "PASSWORD")

Информация о наборе данных:

> dput(Germany_M$year[1:25])
1990:2014

> dput(Germany_M$age)
c(0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, ..., 108, 109, 110)

Ошибка, которую я объясню вам позже, возникает в данных 1991 года, поэтому я покажу вам только данные за 1990-1991-1992 годы (только для всего населения, а не мужчин и женщин)

# amount of population in 1990,1991,1992
> dput(Germany_M$pop$total[,1:3])
structure(c(896913.64, 896875.32, 898207.72, 879287.22, 850054.15, 
839359.64, 835409.03, 857772.58, 863924.38, 868836.24, 845586.95, 
812876.68, 812902.07, 798796.39, 782007.23, 783589.61, 801720.01, 
847608.69, 977237.46, 1062028.44, 1133724.83, 1235571.04, 1300948.78, 
1347307.2, 1379973.91, 1403661.77, 1421445.54, 1399885.26, 1365206.32, 
1336701.04, 1300938.62, 1258303.96, 1211491.66, 1190653.48, 1163796.6, 
1142818.59, 1123727.65, 1116924.22, 1114520.8, 1110869.09, 1109656.64, 
1060850.28, 974871.83, 945218.06, 778546.33, 894046.67, 1036240.34, 
1033737.34, 1131248.17, 1279443.67, 1309442.5, 1260110.5, 1183538.33, 
1137830.5, 1114928.5, 1068930.5, 939304.33, 844730.67, 856880.5, 
895016.5, 909947.33, 898300.83, 865206.67, 825870.17, 815721.83, 
782561.83, 745536.33, 755017, 784991.33, 788563.33, 677506.67, 
471913.83, 359179.17, 362213.17, 428866.67, 539703, 589021.67, 
570643.16, 528446.5, 486068.83, 459015.43, 420121.6, 373631.07, 
327781.69, 279963.02, 236851.03, 197042.15, 162373.22, 132051.91, 
101656.46, 76674.86, 57747.07, 42098.56, 30077.66, 20894.71, 
13971.05, 9212.02, 5833.84, 3748.51, 2420, 1247, 631.67, 332.67, 
188.33, 105.83, 53.5, 24.17, 12.17, 8.17, 12, 0.67, 878166.7, 
903652.22, 906643.96, 908339.03, 889577.67, 859791.16, 848243.07, 
844237.3, 866647.03, 873190.11, 877873.56, 854677.7, 821620.09, 
821373.15, 807857.66, 792207.99, 793725.15, 810755.79, 861172.96, 
993605.42, 1079892.35, 1150944.61, 1254152.18, 1318950.79, 1363418.17, 
1395209.56, 1418141.74, 1434127.06, 1412051.42, 1376656.39, 1348414.96, 
1312142.02, 1268380.13, 1221039.99, 1199684.17, 1172311.52, 1150503.25, 
1130011.9, 1122993.64, 1120285.07, 1115514.34, 1113711.39, 1063604.51, 
977238.3, 945861.27, 780180.11, 893915.17, 1034896.67, 1033592.33, 
1129190.83, 1277177.5, 1306247.17, 1257165, 1179783, 1133570.83, 
1109991.83, 1062374.67, 933873.67, 839141, 850556.67, 886601.67, 
900670.67, 888103.33, 854665.33, 814967.33, 803662.5, 769712.67, 
732528.17, 740167, 767443.83, 769106.33, 657860, 458371, 347102.5, 
348705.33, 411406.67, 514879.17, 559087.83, 538073.67, 494522.33, 
450771.23, 421943.82, 382135.14, 336218.29, 291370.9, 245419.75, 
204568.59, 167526.38, 135807.46, 108396.5, 81806.9, 60515.31, 
44586.85, 31746.87, 22125.63, 14957.19, 9717.87, 6271.02, 3890, 
2502.17, 1414.33, 748, 352.33, 178.17, 102.17, 59.5, 27.83, 13.5, 
7.67, 4.83, 1.5, 819166.01, 884885.18, 912927.66, 916670.95, 
918548.98, 899701.15, 868811.6, 856905.08, 853167.18, 875518.85, 
882093.24, 886644.95, 863148.37, 830115.3, 830202.15, 817822.28, 
802868.4, 804501.45, 825853.31, 879978.54, 1016710.9, 1102191.16, 
1175521.91, 1277677.14, 1341265.76, 1383512.77, 1413815.23, 1435263.4, 
1449903.7, 1426655.43, 1390845.92, 1361316.33, 1325004.75, 1279849.66, 
1231674.43, 1209391.04, 1181560.17, 1158992.2, 1137668.22, 1129746.64, 
1126258.25, 1120579.75, 1118148.52, 1066975.05, 979179.59, 947586.78, 
779828.85, 892880.17, 1033473.5, 1031123.67, 1126768, 1273645.33, 
1302395.33, 1252700.5, 1175046.17, 1128229.33, 1103574.83, 1055905.83, 
926832.5, 832059.17, 842035.83, 876989.5, 890203.33, 876978.17, 
842838.5, 802872, 790588.33, 756321.83, 718271.5, 724024.33, 
748737, 748023.33, 639067, 443302.33, 334242.83, 334350, 392509.83, 
489339.5, 528428.33, 504987, 459708.14, 415829.23, 385123.65, 
345240.98, 300198.95, 256778.57, 213202.91, 174806.42, 140850.72, 
112219.3, 87897.33, 64885.23, 47084.48, 33881.35, 23571.03, 15983.97, 
10565.02, 6748.53, 4228.51, 2574.33, 1475.67, 853.33, 443.33, 
204.67, 104, 60.83, 32.5, 12.5, 6.5, 3.83, 2.17), .Dim = c(111L, 
3L), .Dimnames = list(c("0", "1", "2", "3", "4", "5", "6", "7", 
"8", "9", "10", "11", "12", "13", "14", "15", "16", "17", "18", 
"19", "20", "21", "22", "23", "24", "25", "26", "27", "28", "29", 
"30", "31", "32", "33", "34", "35", "36", "37", "38", "39", "40", 
"41", "42", "43", "44", "45", "46", "47", "48", "49", "50", "51", 
"52", "53", "54", "55", "56", "57", "58", "59", "60", "61", "62", 
"63", "64", "65", "66", "67", "68", "69", "70", "71", "72", "73", 
"74", "75", "76", "77", "78", "79", "80", "81", "82", "83", "84", 
"85", "86", "87", "88", "89", "90", "91", "92", "93", "94", "95", 
"96", "97", "98", "99", "100", "101", "102", "103", "104", "105", 
"106", "107", "108", "109", "110+"), c("1990", "1991", "1992"
)))

Далее идут данные о смертности:

# mortality rate in 1990,1991,1992
> dput(Germany_M$rate$total[,1:3])
structure(c(0.007119, 0.000647, 0.00041, 0.000373, 0.000258, 
0.000251, 0.000218, 0.000245, 0.000225, 0.000165, 0.000175, 0.000187, 
0.000165, 0.000222, 0.00029, 0.000323, 0.000422, 0.000501, 0.000755, 
0.000782, 0.000766, 8e-04, 0.000763, 0.00077, 7e-04, 0.000782, 
0.000742, 0.000758, 0.000784, 0.000908, 0.00096, 0.000936, 0.001049, 
0.001146, 0.001228, 0.00133, 0.001402, 0.001577, 0.001624, 0.001782, 
0.001914, 0.001958, 0.002243, 0.002402, 0.002562, 0.003018, 0.003358, 
0.003487, 0.003914, 0.004087, 0.004717, 0.004996, 0.005402, 0.00617, 
0.006752, 0.007807, 0.008125, 0.009093, 0.009913, 0.011079, 0.012176, 
0.013241, 0.01393, 0.014956, 0.016064, 0.017474, 0.019039, 0.020692, 
0.02329, 0.025322, 0.031231, 0.028207, 0.03351, 0.037694, 0.03925, 
0.047793, 0.051304, 0.058191, 0.065135, 0.072105, 0.081741, 0.091005, 
0.101488, 0.113726, 0.126628, 0.141132, 0.157617, 0.172622, 0.191835, 
0.20896, 0.233755, 0.253642, 0.275071, 0.301754, 0.322428, 0.35452, 
0.380592, 0.393055, 0.395891, 0.428926, 0.55814, 0.573087, 0.643287, 
0.562832, 0.48189, 0.542056, 0.62069, 0.328767, 0.489796, 2.083333, 
6, 0.006503, 0.000614, 0.000386, 0.000334, 0.000273, 0.000238, 
0.000196, 0.000174, 0.000204, 0.000165, 0.000152, 0.000183, 0.000189, 
0.000173, 0.000213, 0.000294, 0.000454, 0.000577, 0.000756, 0.000828, 
0.000774, 0.000772, 0.000755, 0.000754, 0.000755, 0.000776, 0.000799, 
0.000826, 0.000846, 0.000923, 0.000951, 0.00103, 0.001122, 0.001152, 
0.001254, 0.001331, 0.001471, 0.001621, 0.001709, 0.001742, 0.001993, 
0.002116, 0.002215, 0.002409, 0.002662, 0.002966, 0.003305, 0.00356, 
0.003888, 0.004232, 0.00457, 0.00519, 0.005419, 0.006043, 0.006642, 
0.00711, 0.008186, 0.008741, 0.009777, 0.010714, 0.012058, 0.013017, 
0.014265, 0.015324, 0.016133, 0.017241, 0.018997, 0.02025, 0.02216, 
0.024794, 0.027147, 0.033103, 0.030128, 0.03639, 0.040891, 0.042846, 
0.05234, 0.055007, 0.06299, 0.070424, 0.07798, 0.087417, 0.099381, 
0.110708, 0.123358, 0.137141, 0.152633, 0.168511, 0.186389, 0.204628, 
0.221839, 0.245905, 0.265056, 0.29058, 0.31633, 0.341642, 0.362837, 
0.387338, 0.403856, 0.470392, 0.547961, 0.514706, 0.596026, 0.623012, 
0.626427, 0.605042, 0.57485, 0.592593, 0.652174, 0.206897, 0, 
0.006094, 0.000542, 0.000345, 0.000274, 0.000238, 0.000215, 0.000198, 
0.00017, 0.000163, 0.00014, 0.000164, 0.000162, 0.000162, 0.000164, 
0.000201, 0.000287, 0.00037, 0.000481, 0.000729, 0.000848, 0.000755, 
0.000756, 0.000776, 0.000735, 0.000703, 0.00075, 0.000778, 0.000749, 
0.000821, 0.00088, 0.00093, 0.000968, 0.001007, 0.001126, 0.001254, 
0.00134, 0.001389, 0.001491, 0.001646, 0.001793, 0.001907, 0.002147, 
0.002198, 0.002486, 0.002592, 0.002808, 0.00312, 0.003517, 0.003869, 
0.004032, 0.004632, 0.004681, 0.005419, 0.005692, 0.006297, 0.00699, 
0.007661, 0.008712, 0.009417, 0.010527, 0.011497, 0.012631, 0.013793, 
0.014925, 0.016234, 0.016933, 0.018151, 0.01982, 0.021382, 0.023815, 
0.026136, 0.029014, 0.035225, 0.032335, 0.038708, 0.044316, 0.046203, 
0.055532, 0.059261, 0.065984, 0.074739, 0.08383, 0.09276, 0.104582, 
0.116116, 0.130622, 0.143239, 0.161167, 0.17838, 0.192685, 0.2178, 
0.23634, 0.257813, 0.281541, 0.303635, 0.328079, 0.347846, 0.363042, 
0.395411, 0.44128, 0.499435, 0.485156, 0.509774, 0.478827, 0.375, 
0.575342, 0.615385, 0.88, 0.923077, 1.565217, 1.384615), .Dim = c(111L, 
3L), .Dimnames = ...))

Кажется, данные отсутствуют (однако есть некоторые аномальные данные по смертности).Тем не менее, когда я вычисляю команду lifetable() для создания жизненной таблицы из данных HMD (жизненной таблицы для всего населения), возникают некоторые проблемы:

# total population lifetable
 GER_t <- lifetable(Germany_M, series=names(Germany_M$rate)[3], years=Germany_M$year[1:25], ages=Germany_M$age,
                                          max.age=min(110, max(Germany_M$age)), type=c("period"))

Теперь созданная жизненная таблица имеет некоторые данные "inf", и я нене понимаю, почему

# Tx (total years lived by survivors at age x until age 110) in year 1991
> dput(GER_t$Tx[,2])
structure(c(Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, 
Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, 
Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, 
Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, 
Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, 
Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, 
Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, 
Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, 
Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf), .Names = c("0", 
"1", "2", "3", ...))

Также данные об ожидаемой продолжительности жизни представлены таким же образом:

# life expectancy for individuals at age x in 1991
> dput(GER_t$ex[,2])
structure(c(Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, 
Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, 
Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, 
Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, 
Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, 
Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, 
Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, 
Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, 
Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf, Inf), .Names = c("0", 
"1", "2", "3", ...))

Одна из моих целей - создать график (с X = годами, а Y = ожидаемой продолжительностью жизни).e_x) которая представляет кривую ожидаемой продолжительности жизни для каждого возраста на протяжении всех лет, включенных в набор данных.Для этого я должен правильно составить описанную выше жизненную таблицу.График должен выглядеть так, как показано на рисунке: График ожидаемой продолжительности жизни

Итак, мои вопросы: 1) Почему это происходит?2) Как я могу решить эту проблему?

Спасибо за помощь