“Inf” data when transforming data into life table

I'm a novice programmer so have a little patience with me. I downloaded some mortality data from "Human Mortality Database (HMD)".

Packages required:

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


Downloading datasets:

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


Dataset infos:

> 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)


The error I'll explain you next happens within 1991 data so I'll show you only data for year 1990-1991-1992 (only for total population, not male and female)

# 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"
)))


Next there are mortality rate data:

# 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 = ...))


It seems no data is missing (some anomalous data in death rates however). However, when I compute lifetable() command to produce a lifetable from HMD data (lifetable for total population) some problems occur:

lifetable()

# 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"))


Now the lifetable produced has some "inf" data, and I don't understand why

# 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", ...))


Also life expectancy data presents in the same way:

# 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", ...))


One of my objectives is to generate a plot (with X=years and Y=life expectancy e_x) that represents a life expectancy curve for each age throughout all the years included in the dataset. In order to do this I should produce the lifetable described above correctly. The plot should look like shown in the image:
Life expectancy plot

So my questions are:
1) Why this happens?
2) How can i fix this problem?

Thanks to all who'll try to answer.

dput()

0

Thanks for contributing an answer to Stack Overflow!

But avoid …

To learn more, see our tips on writing great answers.

Required, but never shown

Required, but never shown

By clicking "Post Your Answer", you acknowledge that you have read our updated terms of service, privacy policy and cookie policy, and that your continued use of the website is subject to these policies.