Skip to content

Instantly share code, notes, and snippets.

@timriffe
Created July 7, 2023 07:05
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save timriffe/4d21788219b1c6ec0dc0acec9cd7f2fa to your computer and use it in GitHub Desktop.
Save timriffe/4d21788219b1c6ec0dc0acec9cd7f2fa to your computer and use it in GitHub Desktop.
test data for Coale's r estimation function
srb <- 1.05
pfb <- 1/ (1 + srb)
fxf <- c(0.00000,0.00003,0.00038,0.00146,0.00436,0.01006,0.01963,0.03576,
0.06036,0.09112,0.12532,0.15873,0.18765,0.20418,0.20929,0.20586,
0.19741,0.18698,0.16816,0.14313,0.12684,0.11769,0.11308,0.10262,
0.08831,0.07622,0.06532,0.05493,0.04629,0.03623,0.02455,0.01515,
0.00881,0.00494,0.00305,0.00219,0.00174,0.00140,0.00124,0.00065,
0.00019,0.00005,0.00001,0.00000) * pfb
fxf <- c(rep(0,12),fxf,rep(0,55))
Lx <- c(97891, 97468, 97363, 97290, 97230, 97185, 97145, 97108, 97073,
97042, 97012, 96982, 96952, 96924, 96897, 96868, 96837, 96800,
96756, 96712, 96669, 96623, 96574, 96520, 96464, 96402, 96333,
96266, 96195, 96123, 96047, 95970, 95894, 95809, 95716, 95611,
95498, 95382, 95253, 95106, 94953, 94799, 94624, 94433, 94255,
94057, 93830, 93589, 93321, 93029, 92718, 92370, 91970, 91552,
91113, 90631, 90105, 89542, 88941, 88283, 87561, 86788, 85938,
84983, 83935, 82759, 81473, 80077, 78537, 76830, 74909, 72822,
70557, 68089, 65380, 62386, 59151, 55750, 52125, 48461, 44798,
41078, 37225, 33219, 29244, 25384, 21743, 18271, 15190, 12509,
10045, 7981, 6233, 4690, 3416, 2403, 1681, 1145, 758, 488, 305,
185, 109, 62, 34, 19, 10, 5, 2, 1, 1) / 1e5
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment