Relative sea level rate (RSL) estimation
BchronRSL(
BchronologyRun,
RSLmean,
RSLsd,
degree = 1,
iterations = 10000,
burn = 2000,
thin = 8
)Arguments
- BchronologyRun
Output from a run of
Bchronology- RSLmean
A vector of RSL mean estimates of the same length as the number of predictPositions given to the
Bchronologyfunction- RSLsd
A vector RSL standard deviations of the same length as the number of predictPositions given to the
Bchronologyfunction- degree
The degree of the polynomial regression: linear=1 (default), quadratic=2, etc. Supports up to degree 5, though this will depend on the data given
- iterations
The number of MCMC iterations to run
- burn
The number of starting iterations to discard
- thin
The step size of iterations to discard
Value
An object of class BchronRSLRun with elements
- BchronologyRun
The output from the run of
Bchronology- samples
The posterior samples of the regression parameters
- degree
The degree of the polynomial regression
- RSLmean
The RSL mean values given to the function
- RSLsd
The RSL standard deviations as given to the function
- const
The mean of the predicted age values. Used to standardise the design matrix and avoid computational issues
Details
This function fits an errors-in-variables regression model to relative sea level (RSL) data. An errors-in-variables regression model allows for uncertainty in the explanatory variable, here the age of sea level data point. The algorithm is more fully defined in the reference below
References
Andrew C. Parnell and W. Roland Gehrels (2013) 'Using chronological models in late holocene sea level reconstructions from salt marsh sediments' In: I. Shennan, B.P. Horton, and A.J. Long (eds). Handbook of Sea Level Research. Chichester: Wiley
See also
Examples
# \donttest{
# Load in data
data(TestChronData)
data(TestRSLData)
# Run through Bchronology
RSLrun <- with(TestChronData, Bchronology(
ages = ages,
ageSds = ageSds,
positions = position,
positionThicknesses = thickness,
ids = id,
calCurves = calCurves,
predictPositions = TestRSLData$Depth
))
#> ages not given as whole numbers - rounding occurred
#> ageSds not given as whole numbers - rounding occurred
#> ages not given as whole numbers - rounding occurred
#> ageSds not given as whole numbers - rounding occurred
#> ages not given as whole numbers - rounding occurred
#> ageSds not given as whole numbers - rounding occurred
#> Running Bchronology...
#>
|
| | 0%
|
| | 1%
|
|= | 1%
|
|= | 2%
|
|== | 3%
|
|== | 4%
|
|=== | 4%
|
|=== | 5%
|
|=== | 6%
|
|==== | 6%
|
|==== | 7%
|
|===== | 8%
|
|===== | 9%
|
|====== | 9%
|
|====== | 10%
|
|====== | 11%
|
|======= | 11%
|
|======= | 12%
|
|======== | 13%
|
|======== | 14%
|
|========= | 14%
|
|========= | 15%
|
|========= | 16%
|
|========== | 16%
|
|========== | 17%
|
|=========== | 18%
|
|=========== | 19%
|
|============ | 19%
|
|============ | 20%
|
|============ | 21%
|
|============= | 21%
|
|============= | 22%
|
|============== | 23%
|
|============== | 24%
|
|=============== | 24%
|
|=============== | 25%
|
|=============== | 26%
|
|================ | 26%
|
|================ | 27%
|
|================= | 28%
|
|================= | 29%
|
|================== | 29%
|
|================== | 30%
|
|================== | 31%
|
|=================== | 31%
|
|=================== | 32%
|
|==================== | 33%
|
|==================== | 34%
|
|===================== | 34%
|
|===================== | 35%
|
|===================== | 36%
|
|====================== | 36%
|
|====================== | 37%
|
|======================= | 38%
|
|======================= | 39%
|
|======================== | 39%
|
|======================== | 40%
|
|======================== | 41%
|
|========================= | 41%
|
|========================= | 42%
|
|========================== | 43%
|
|========================== | 44%
|
|=========================== | 44%
|
|=========================== | 45%
|
|=========================== | 46%
|
|============================ | 46%
|
|============================ | 47%
|
|============================= | 48%
|
|============================= | 49%
|
|============================== | 49%
|
|============================== | 50%
|
|============================== | 51%
|
|=============================== | 51%
|
|=============================== | 52%
|
|================================ | 53%
|
|================================ | 54%
|
|================================= | 54%
|
|================================= | 55%
|
|================================= | 56%
|
|================================== | 56%
|
|================================== | 57%
|
|=================================== | 58%
|
|=================================== | 59%
|
|==================================== | 59%
|
|==================================== | 60%
|
|==================================== | 61%
|
|===================================== | 61%
|
|===================================== | 62%
|
|====================================== | 63%
|
|====================================== | 64%
|
|======================================= | 64%
|
|======================================= | 65%
|
|======================================= | 66%
|
|======================================== | 66%
|
|======================================== | 67%
|
|========================================= | 68%
|
|========================================= | 69%
|
|========================================== | 69%
|
|========================================== | 70%
|
|========================================== | 71%
|
|=========================================== | 71%
|
|=========================================== | 72%
|
|============================================ | 73%
|
|============================================ | 74%
|
|============================================= | 74%
|
|============================================= | 75%
|
|============================================= | 76%
|
|============================================== | 76%
|
|============================================== | 77%
|
|=============================================== | 78%
|
|=============================================== | 79%
|
|================================================ | 79%
|
|================================================ | 80%
|
|================================================ | 81%
|
|================================================= | 81%
|
|================================================= | 82%
|
|================================================== | 83%
|
|================================================== | 84%
|
|=================================================== | 84%
|
|=================================================== | 85%
|
|=================================================== | 86%
|
|==================================================== | 86%
|
|==================================================== | 87%
|
|===================================================== | 88%
|
|===================================================== | 89%
|
|====================================================== | 89%
|
|====================================================== | 90%
|
|====================================================== | 91%
|
|======================================================= | 91%
|
|======================================================= | 92%
|
|======================================================== | 93%
|
|======================================================== | 94%
|
|========================================================= | 94%
|
|========================================================= | 95%
|
|========================================================= | 96%
|
|========================================================== | 96%
|
|========================================================== | 97%
|
|=========================================================== | 98%
|
|=========================================================== | 99%
|
|============================================================| 99%
|
|============================================================| 100%
#> Run completed!
# Now run through BchronRSL
RSLrun2 <- BchronRSL(RSLrun, RSLmean = TestRSLData$RSL, RSLsd = TestRSLData$Sigma, degree = 3)
#>
|
| | 0%
|
| | 1%
|
|= | 1%
|
|= | 2%
|
|== | 3%
|
|== | 4%
|
|=== | 4%
|
|=== | 5%
|
|=== | 6%
|
|==== | 6%
|
|==== | 7%
|
|===== | 8%
|
|===== | 9%
|
|====== | 9%
|
|====== | 10%
|
|====== | 11%
|
|======= | 11%
|
|======= | 12%
|
|======== | 13%
|
|======== | 14%
|
|========= | 14%
|
|========= | 15%
|
|========= | 16%
|
|========== | 16%
|
|========== | 17%
|
|=========== | 18%
|
|=========== | 19%
|
|============ | 19%
|
|============ | 20%
|
|============ | 21%
|
|============= | 21%
|
|============= | 22%
|
|============== | 23%
|
|============== | 24%
|
|=============== | 24%
|
|=============== | 25%
|
|=============== | 26%
|
|================ | 26%
|
|================ | 27%
|
|================= | 28%
|
|================= | 29%
|
|================== | 29%
|
|================== | 30%
|
|================== | 31%
|
|=================== | 31%
|
|=================== | 32%
|
|==================== | 33%
|
|==================== | 34%
|
|===================== | 34%
|
|===================== | 35%
|
|===================== | 36%
|
|====================== | 36%
|
|====================== | 37%
|
|======================= | 38%
|
|======================= | 39%
|
|======================== | 39%
|
|======================== | 40%
|
|======================== | 41%
|
|========================= | 41%
|
|========================= | 42%
|
|========================== | 43%
|
|========================== | 44%
|
|=========================== | 44%
|
|=========================== | 45%
|
|=========================== | 46%
|
|============================ | 46%
|
|============================ | 47%
|
|============================= | 48%
|
|============================= | 49%
|
|============================== | 49%
|
|============================== | 50%
|
|============================== | 51%
|
|=============================== | 51%
|
|=============================== | 52%
|
|================================ | 53%
|
|================================ | 54%
|
|================================= | 54%
|
|================================= | 55%
|
|================================= | 56%
|
|================================== | 56%
|
|================================== | 57%
|
|=================================== | 58%
|
|=================================== | 59%
|
|==================================== | 59%
|
|==================================== | 60%
|
|==================================== | 61%
|
|===================================== | 61%
|
|===================================== | 62%
|
|====================================== | 63%
|
|====================================== | 64%
|
|======================================= | 64%
|
|======================================= | 65%
|
|======================================= | 66%
|
|======================================== | 66%
|
|======================================== | 67%
|
|========================================= | 68%
|
|========================================= | 69%
|
|========================================== | 69%
|
|========================================== | 70%
|
|========================================== | 71%
|
|=========================================== | 71%
|
|=========================================== | 72%
|
|============================================ | 73%
|
|============================================ | 74%
|
|============================================= | 74%
|
|============================================= | 75%
|
|============================================= | 76%
|
|============================================== | 76%
|
|============================================== | 77%
|
|=============================================== | 78%
|
|=============================================== | 79%
|
|================================================ | 79%
|
|================================================ | 80%
|
|================================================ | 81%
|
|================================================= | 81%
|
|================================================= | 82%
|
|================================================== | 83%
|
|================================================== | 84%
|
|=================================================== | 84%
|
|=================================================== | 85%
|
|=================================================== | 86%
|
|==================================================== | 86%
|
|==================================================== | 87%
|
|===================================================== | 88%
|
|===================================================== | 89%
|
|====================================================== | 89%
|
|====================================================== | 90%
|
|====================================================== | 91%
|
|======================================================= | 91%
|
|======================================================= | 92%
|
|======================================================== | 93%
|
|======================================================== | 94%
|
|========================================================= | 94%
|
|========================================================= | 95%
|
|========================================================= | 96%
|
|========================================================== | 96%
|
|========================================================== | 97%
|
|=========================================================== | 98%
|
|=========================================================== | 99%
|
|============================================================| 99%
|
|============================================================| 100%
# Summarise it
summary(RSLrun2)
#> Posterior Medians with 95% credible intervals...
#> Power Lower Median Upper
#> mean -0.8312 -0.6579 -0.4859
#> linear -0.7778 -0.4595 -0.1477
#> quadratic -0.1478 0.3131 0.7867
#> cubic -0.8865 -0.3762 0.0877
#> Age RSL_2.5 RSL_50 RSL_97.5
#> 1 1956.000000 -1.8036736 -1.50260400 -1.25216279
#> 2 1936.111364 -1.7556459 -1.47602472 -1.24071976
#> 3 1916.222727 -1.7056543 -1.44771459 -1.22525394
#> 4 1896.334091 -1.6675277 -1.42108034 -1.20804675
#> 5 1876.445455 -1.6338059 -1.39494779 -1.18957857
#> 6 1856.556818 -1.5935126 -1.36991668 -1.16698062
#> 7 1836.668182 -1.5630950 -1.34511903 -1.14865261
#> 8 1816.779545 -1.5323828 -1.32162309 -1.12914324
#> 9 1796.890909 -1.4997452 -1.29920209 -1.10648669
#> 10 1777.002273 -1.4760987 -1.27572777 -1.08815014
#> 11 1757.113636 -1.4515341 -1.25273812 -1.07046175
#> 12 1737.225000 -1.4267056 -1.23048622 -1.05161674
#> 13 1717.336364 -1.4050828 -1.21129788 -1.03105197
#> 14 1697.447727 -1.3835620 -1.19292074 -1.01166284
#> 15 1677.559091 -1.3630812 -1.17535380 -0.99004677
#> 16 1657.670455 -1.3429682 -1.15702154 -0.96746418
#> 17 1637.781818 -1.3266718 -1.13759524 -0.95370455
#> 18 1617.893182 -1.3083973 -1.12045984 -0.94075945
#> 19 1598.004545 -1.2939574 -1.10161156 -0.92186618
#> 20 1578.115909 -1.2816099 -1.08472900 -0.90549469
#> 21 1558.227273 -1.2621487 -1.06823865 -0.89120386
#> 22 1538.338636 -1.2465085 -1.05412497 -0.88007037
#> 23 1518.450000 -1.2285063 -1.04003376 -0.86252146
#> 24 1498.561364 -1.2137494 -1.02645458 -0.85042173
#> 25 1478.672727 -1.1990293 -1.01278885 -0.84116163
#> 26 1458.784091 -1.1806451 -0.99929592 -0.82660587
#> 27 1438.895455 -1.1615166 -0.98587658 -0.81287782
#> 28 1419.006818 -1.1504033 -0.97358169 -0.80107132
#> 29 1399.118182 -1.1363702 -0.96144469 -0.78960248
#> 30 1379.229545 -1.1205081 -0.94933942 -0.78238990
#> 31 1359.340909 -1.1096544 -0.93819866 -0.77387987
#> 32 1339.452273 -1.0961770 -0.92708978 -0.76773203
#> 33 1319.563636 -1.0825131 -0.91673806 -0.75852289
#> 34 1299.675000 -1.0684370 -0.90658000 -0.74986437
#> 35 1279.786364 -1.0560002 -0.89681221 -0.74474386
#> 36 1259.897727 -1.0494471 -0.88692092 -0.74067880
#> 37 1240.009091 -1.0351117 -0.87766397 -0.73512341
#> 38 1220.120455 -1.0219713 -0.86931236 -0.72805392
#> 39 1200.231818 -1.0134990 -0.85888528 -0.72045728
#> 40 1180.343182 -1.0051541 -0.85090317 -0.71242913
#> 41 1160.454545 -0.9947714 -0.84215023 -0.70666050
#> 42 1140.565909 -0.9820526 -0.83385301 -0.69864792
#> 43 1120.677273 -0.9740254 -0.82639858 -0.69530807
#> 44 1100.788636 -0.9660214 -0.81786559 -0.68990870
#> 45 1080.900000 -0.9559408 -0.80980465 -0.68451756
#> 46 1061.011364 -0.9481109 -0.80228551 -0.67635652
#> 47 1041.122727 -0.9363558 -0.79573582 -0.66807457
#> 48 1021.234091 -0.9289819 -0.78901344 -0.65894688
#> 49 1001.345455 -0.9207791 -0.78104275 -0.65407145
#> 50 981.456818 -0.9146988 -0.77418375 -0.64510964
#> 51 961.568182 -0.9063048 -0.76786874 -0.63677160
#> 52 941.679545 -0.9003199 -0.76016323 -0.62892899
#> 53 921.790909 -0.8919792 -0.75258455 -0.61744758
#> 54 901.902273 -0.8869427 -0.74436152 -0.60630123
#> 55 882.013636 -0.8824253 -0.73624354 -0.59939275
#> 56 862.125000 -0.8754553 -0.72816195 -0.58840009
#> 57 842.236364 -0.8715165 -0.72003105 -0.57952856
#> 58 822.347727 -0.8595411 -0.71215078 -0.56820259
#> 59 802.459091 -0.8564357 -0.70611361 -0.55836018
#> 60 782.570455 -0.8538964 -0.69868942 -0.54594874
#> 61 762.681818 -0.8477299 -0.68933014 -0.52988639
#> 62 742.793182 -0.8449020 -0.68124612 -0.51679045
#> 63 722.904545 -0.8380421 -0.67273418 -0.50410703
#> 64 703.015909 -0.8337394 -0.66383095 -0.49273423
#> 65 683.127273 -0.8293988 -0.65401049 -0.48069285
#> 66 663.238636 -0.8232672 -0.64430150 -0.46907434
#> 67 643.350000 -0.8145400 -0.63388235 -0.45723126
#> 68 623.461364 -0.8052817 -0.62412110 -0.44720250
#> 69 603.572727 -0.7980852 -0.61398857 -0.43633836
#> 70 583.684091 -0.7890258 -0.60330362 -0.42407310
#> 71 563.795455 -0.7797042 -0.59263293 -0.41418780
#> 72 543.906818 -0.7706194 -0.58071990 -0.40349936
#> 73 524.018182 -0.7602079 -0.56872751 -0.39317196
#> 74 504.129545 -0.7488681 -0.55776812 -0.38089143
#> 75 484.240909 -0.7344616 -0.54563082 -0.36746641
#> 76 464.352273 -0.7184003 -0.53282718 -0.35433712
#> 77 444.463636 -0.7029787 -0.51795577 -0.34223253
#> 78 424.575000 -0.6869022 -0.50350040 -0.33135589
#> 79 404.686364 -0.6710036 -0.48905182 -0.32037449
#> 80 384.797727 -0.6563457 -0.47473714 -0.30928245
#> 81 364.909091 -0.6394329 -0.45908873 -0.29714978
#> 82 345.020455 -0.6158937 -0.44362105 -0.28289712
#> 83 325.131818 -0.5983914 -0.42739239 -0.27138473
#> 84 305.243182 -0.5799973 -0.41026519 -0.26058179
#> 85 285.354545 -0.5579820 -0.39337648 -0.24537076
#> 86 265.465909 -0.5327128 -0.37480191 -0.23073305
#> 87 245.577273 -0.5072470 -0.35643195 -0.21656536
#> 88 225.688636 -0.4812166 -0.33606213 -0.20118656
#> 89 205.800000 -0.4532858 -0.31655435 -0.18663565
#> 90 185.911364 -0.4260392 -0.29566908 -0.17146789
#> 91 166.022727 -0.3988338 -0.27384805 -0.15144894
#> 92 146.134091 -0.3712441 -0.25196163 -0.13625069
#> 93 126.245455 -0.3469415 -0.23025846 -0.11766079
#> 94 106.356818 -0.3215834 -0.20744046 -0.09693170
#> 95 86.468182 -0.2928731 -0.18180436 -0.07733801
#> 96 66.579545 -0.2669212 -0.15648784 -0.05139012
#> 97 46.690909 -0.2449200 -0.13194475 -0.01965832
#> 98 26.802273 -0.2250010 -0.10637602 0.01479978
#> 99 6.913636 -0.2090760 -0.07813643 0.05439399
#> 100 -12.975000 -0.1927903 -0.04812656 0.09682347
#> Age RSL_rate_2.5 RSL_rate_50 RSL_rate_97.5
#> 1 1956.000000 0.402750888 1.4380685 2.7324577
#> 2 1936.111364 0.405657378 1.3953101 2.6332725
#> 3 1916.222727 0.408723131 1.3568259 2.5259171
#> 4 1896.334091 0.413287236 1.3186393 2.4332780
#> 5 1876.445455 0.415961295 1.2794244 2.3424502
#> 6 1856.556818 0.427716011 1.2377065 2.2535147
#> 7 1836.668182 0.430091916 1.2005259 2.1663904
#> 8 1816.779545 0.436888486 1.1646959 2.0800483
#> 9 1796.890909 0.443871234 1.1280053 1.9944440
#> 10 1777.002273 0.454746610 1.0916530 1.9108723
#> 11 1757.113636 0.470996440 1.0581159 1.8290891
#> 12 1737.225000 0.475797706 1.0220744 1.7514715
#> 13 1717.336364 0.470652127 0.9892602 1.6776700
#> 14 1697.447727 0.471427217 0.9568134 1.6084295
#> 15 1677.559091 0.465226190 0.9238636 1.5360176
#> 16 1657.670455 0.458480541 0.8953063 1.4654445
#> 17 1637.781818 0.452277985 0.8674943 1.3943458
#> 18 1617.893182 0.454685185 0.8389779 1.3258874
#> 19 1598.004545 0.442766787 0.8121560 1.2615872
#> 20 1578.115909 0.435090009 0.7857876 1.2077775
#> 21 1558.227273 0.427796631 0.7615013 1.1569912
#> 22 1538.338636 0.411055063 0.7348155 1.1157552
#> 23 1518.450000 0.394429063 0.7101014 1.0748590
#> 24 1498.561364 0.374482210 0.6905867 1.0412744
#> 25 1478.672727 0.362265602 0.6690981 0.9948238
#> 26 1458.784091 0.350082731 0.6429821 0.9692869
#> 27 1438.895455 0.332084451 0.6255232 0.9465575
#> 28 1419.006818 0.314101514 0.6059742 0.9258154
#> 29 1399.118182 0.292158898 0.5897367 0.9033872
#> 30 1379.229545 0.265505188 0.5699963 0.8940844
#> 31 1359.340909 0.236330951 0.5527811 0.8769820
#> 32 1339.452273 0.219700243 0.5323060 0.8637156
#> 33 1319.563636 0.192473202 0.5171497 0.8511711
#> 34 1299.675000 0.168793863 0.5041730 0.8385973
#> 35 1279.786364 0.148789576 0.4900413 0.8247023
#> 36 1259.897727 0.129482426 0.4764323 0.8144459
#> 37 1240.009091 0.112072910 0.4614285 0.8163265
#> 38 1220.120455 0.087005000 0.4482778 0.8107683
#> 39 1200.231818 0.067052743 0.4381886 0.8023251
#> 40 1180.343182 0.047440742 0.4289908 0.7974021
#> 41 1160.454545 0.031391817 0.4193156 0.7937826
#> 42 1140.565909 0.016506956 0.4132341 0.7956122
#> 43 1120.677273 0.002716240 0.4043745 0.7929188
#> 44 1100.788636 -0.012684687 0.3972747 0.7885326
#> 45 1080.900000 -0.022707595 0.3918242 0.7862187
#> 46 1061.011364 -0.032174023 0.3864872 0.7801441
#> 47 1041.122727 -0.032565088 0.3819764 0.7773182
#> 48 1021.234091 -0.040364843 0.3815050 0.7773213
#> 49 1001.345455 -0.041839493 0.3793486 0.7735483
#> 50 981.456818 -0.042427641 0.3786798 0.7719160
#> 51 961.568182 -0.040992980 0.3787205 0.7757073
#> 52 941.679545 -0.040911026 0.3790360 0.7778294
#> 53 921.790909 -0.032706322 0.3798708 0.7773120
#> 54 901.902273 -0.027305251 0.3831243 0.7745862
#> 55 882.013636 -0.020368946 0.3869368 0.7663545
#> 56 862.125000 -0.005088542 0.3891886 0.7632407
#> 57 842.236364 0.009611830 0.3933050 0.7658508
#> 58 822.347727 0.017587527 0.3987294 0.7602823
#> 59 802.459091 0.031470022 0.4057973 0.7672717
#> 60 782.570455 0.048167784 0.4120263 0.7645398
#> 61 762.681818 0.068541901 0.4195324 0.7660238
#> 62 742.793182 0.088336184 0.4312021 0.7657355
#> 63 722.904545 0.110326741 0.4425296 0.7762669
#> 64 703.015909 0.134542408 0.4538858 0.7801340
#> 65 683.127273 0.161018005 0.4657147 0.7759823
#> 66 663.238636 0.189637991 0.4820390 0.7717290
#> 67 643.350000 0.221856392 0.4950589 0.7762717
#> 68 623.461364 0.247142450 0.5111254 0.7774120
#> 69 603.572727 0.272221470 0.5235115 0.7829993
#> 70 583.684091 0.297114109 0.5412542 0.7846777
#> 71 563.795455 0.334076805 0.5590309 0.7913973
#> 72 543.906818 0.355543729 0.5777675 0.8064734
#> 73 524.018182 0.383282104 0.5965315 0.8193558
#> 74 504.129545 0.405890871 0.6156350 0.8358226
#> 75 484.240909 0.432864279 0.6368730 0.8599072
#> 76 464.352273 0.457575709 0.6590579 0.8887440
#> 77 444.463636 0.469597563 0.6818418 0.9180242
#> 78 424.575000 0.478451156 0.7059256 0.9546628
#> 79 404.686364 0.490964167 0.7309619 0.9912194
#> 80 384.797727 0.495888370 0.7599122 1.0315777
#> 81 364.909091 0.490612725 0.7857276 1.0746569
#> 82 345.020455 0.508673719 0.8094808 1.1305968
#> 83 325.131818 0.490093272 0.8421728 1.1957835
#> 84 305.243182 0.492977388 0.8737566 1.2545488
#> 85 285.354545 0.492278935 0.9027287 1.3154709
#> 86 265.465909 0.483337617 0.9317276 1.3788827
#> 87 245.577273 0.473627238 0.9641029 1.4443447
#> 88 225.688636 0.466993435 0.9964051 1.5130257
#> 89 205.800000 0.468886127 1.0289292 1.5923249
#> 90 185.911364 0.460826511 1.0629890 1.6560581
#> 91 166.022727 0.451000687 1.0988001 1.7345373
#> 92 146.134091 0.445221504 1.1366794 1.8151148
#> 93 126.245455 0.441855554 1.1767943 1.9089876
#> 94 106.356818 0.433415053 1.2149198 1.9934682
#> 95 86.468182 0.413568836 1.2543069 2.0829127
#> 96 66.579545 0.411844273 1.2934464 2.1716075
#> 97 46.690909 0.397136704 1.3355517 2.2578115
#> 98 26.802273 0.380127053 1.3777284 2.3453865
#> 99 6.913636 0.368511555 1.4195652 2.4388328
#> 100 -12.975000 0.356836335 1.4634826 2.5324367
#> Age RSL_accel_2.5 RSL_accel_50 RSL_accel_97.5
#> 1 1956.000000 -5.2908248 -2.21235824 0.5420679
#> 2 1936.111364 -5.1850037 -2.17002426 0.5283138
#> 3 1916.222727 -5.0791825 -2.12502808 0.5145597
#> 4 1896.334091 -4.9733614 -2.08131404 0.5008055
#> 5 1876.445455 -4.8675403 -2.03415189 0.4870514
#> 6 1856.556818 -4.7642831 -1.98634784 0.4732972
#> 7 1836.668182 -4.6642825 -1.93910545 0.4595431
#> 8 1816.779545 -4.5643941 -1.89098958 0.4457890
#> 9 1796.890909 -4.4589286 -1.84366030 0.4341583
#> 10 1777.002273 -4.3497198 -1.79591535 0.4313297
#> 11 1757.113636 -4.2404774 -1.75104903 0.4285010
#> 12 1737.225000 -4.1313220 -1.70659783 0.4256724
#> 13 1717.336364 -4.0283841 -1.66437385 0.4228438
#> 14 1697.447727 -3.9250672 -1.62195588 0.4200152
#> 15 1677.559091 -3.8225884 -1.57033299 0.4171865
#> 16 1657.670455 -3.7139163 -1.52520738 0.4129682
#> 17 1637.781818 -3.6005714 -1.48206187 0.3947856
#> 18 1617.893182 -3.4930691 -1.43387574 0.3899574
#> 19 1598.004545 -3.3953708 -1.38608304 0.3804309
#> 20 1578.115909 -3.3055420 -1.34393997 0.3665383
#> 21 1558.227273 -3.1982984 -1.29724111 0.3425442
#> 22 1538.338636 -3.0974995 -1.25085649 0.3256466
#> 23 1518.450000 -3.0041612 -1.20461914 0.3130308
#> 24 1498.561364 -2.9029069 -1.16228729 0.2951909
#> 25 1478.672727 -2.8006985 -1.11809525 0.2773348
#> 26 1458.784091 -2.7018082 -1.07364717 0.2662974
#> 27 1438.895455 -2.6083010 -1.02679826 0.2589019
#> 28 1419.006818 -2.5080155 -0.98174785 0.2638182
#> 29 1399.118182 -2.3943995 -0.94024961 0.2577285
#> 30 1379.229545 -2.2961473 -0.89852243 0.2513295
#> 31 1359.340909 -2.1875275 -0.85361119 0.2499753
#> 32 1339.452273 -2.0742498 -0.80613099 0.2586458
#> 33 1319.563636 -1.9799444 -0.76429355 0.2530652
#> 34 1299.675000 -1.8809061 -0.72532498 0.2425230
#> 35 1279.786364 -1.7787080 -0.67235959 0.2417021
#> 36 1259.897727 -1.6747378 -0.63214354 0.2367564
#> 37 1240.009091 -1.5741997 -0.59338640 0.2525702
#> 38 1220.120455 -1.4814193 -0.55321120 0.2682161
#> 39 1200.231818 -1.3860141 -0.50675461 0.2547567
#> 40 1180.343182 -1.2929269 -0.46376905 0.2336956
#> 41 1160.454545 -1.2023331 -0.41934514 0.2307267
#> 42 1140.565909 -1.1107487 -0.37941691 0.2506890
#> 43 1120.677273 -1.0223477 -0.32631258 0.2536438
#> 44 1100.788636 -0.9301196 -0.29005577 0.2520736
#> 45 1080.900000 -0.8524455 -0.24886057 0.2765812
#> 46 1061.011364 -0.7669886 -0.20434092 0.2928207
#> 47 1041.122727 -0.6765624 -0.16315640 0.3222235
#> 48 1021.234091 -0.6228709 -0.11457857 0.3579591
#> 49 1001.345455 -0.5511985 -0.07595005 0.3932244
#> 50 981.456818 -0.4815112 -0.02703327 0.4316773
#> 51 961.568182 -0.4372059 0.01760296 0.4744578
#> 52 941.679545 -0.4158017 0.05707949 0.5247133
#> 53 921.790909 -0.3791185 0.09691811 0.5766326
#> 54 901.902273 -0.3465637 0.13857484 0.6405361
#> 55 882.013636 -0.3155000 0.18625210 0.7064286
#> 56 862.125000 -0.3177890 0.22754472 0.7774043
#> 57 842.236364 -0.3179365 0.27564723 0.8487656
#> 58 822.347727 -0.3115434 0.31305204 0.9395310
#> 59 802.459091 -0.3072050 0.36367229 1.0193068
#> 60 782.570455 -0.3059813 0.41606948 1.1142727
#> 61 762.681818 -0.3067069 0.46569441 1.2220345
#> 62 742.793182 -0.3026963 0.50996230 1.3154285
#> 63 722.904545 -0.3126772 0.55522662 1.4123742
#> 64 703.015909 -0.3022837 0.60110569 1.5203230
#> 65 683.127273 -0.2920453 0.64396280 1.6092574
#> 66 663.238636 -0.3080132 0.68639063 1.6876292
#> 67 643.350000 -0.3190743 0.73125618 1.7997242
#> 68 623.461364 -0.3341736 0.77211495 1.9138215
#> 69 603.572727 -0.3377805 0.81408895 2.0073454
#> 70 583.684091 -0.3604525 0.86111493 2.0963499
#> 71 563.795455 -0.3817217 0.90473349 2.1765200
#> 72 543.906818 -0.3724148 0.95039793 2.2739184
#> 73 524.018182 -0.3724986 0.99572231 2.3654738
#> 74 504.129545 -0.3974993 1.04181381 2.4609381
#> 75 484.240909 -0.4222890 1.08310311 2.5615740
#> 76 464.352273 -0.4470053 1.12739620 2.6622100
#> 77 444.463636 -0.4717215 1.16814840 2.7631509
#> 78 424.575000 -0.4964378 1.20823535 2.8700544
#> 79 404.686364 -0.5100574 1.25580574 2.9769578
#> 80 384.797727 -0.5199826 1.30237868 3.0838296
#> 81 364.909091 -0.5285832 1.34177647 3.1902945
#> 82 345.020455 -0.5371837 1.38799819 3.2973649
#> 83 325.131818 -0.5457492 1.43207684 3.4045716
#> 84 305.243182 -0.5542857 1.47684589 3.5116241
#> 85 285.354545 -0.5694933 1.52175838 3.6188933
#> 86 265.465909 -0.5802934 1.56806239 3.7260811
#> 87 245.577273 -0.5891016 1.61238713 3.8311833
#> 88 225.688636 -0.5979098 1.65890142 3.9264144
#> 89 205.800000 -0.6066443 1.70580401 4.0372636
#> 90 185.911364 -0.6162757 1.75199468 4.1517738
#> 91 166.022727 -0.6344693 1.79944736 4.2615034
#> 92 146.134091 -0.6418580 1.84622968 4.3586287
#> 93 126.245455 -0.6490832 1.89098764 4.4521822
#> 94 106.356818 -0.6574193 1.93622045 4.5457779
#> 95 86.468182 -0.6723823 1.98312037 4.6515661
#> 96 66.579545 -0.6873453 2.03180151 4.7573543
#> 97 46.690909 -0.7023083 2.07930978 4.8631425
#> 98 26.802273 -0.7172713 2.12582272 4.9689307
#> 99 6.913636 -0.7286941 2.16751526 5.0747189
#> 100 -12.975000 -0.7393445 2.21189872 5.1805072
# Plot it
plot(RSLrun2)
# }