![[Experimental]](figures/lifecycle-experimental.svg)
A function that returns a dose() method that computes the dose reaching a
specific target value of a given independent variable, based on the model
specific parameters.
Usage
doseFunction(model, ...)
# S4 method for class 'GeneralModel'
doseFunction(model, ...)
# S4 method for class 'ModelPseudo'
doseFunction(model, ...)
# S4 method for class 'LogisticLogNormalOrdinal'
doseFunction(model, grade, ...)Value
A dose() method that computes doses.
Functions
doseFunction(GeneralModel): Return a dose function for aGeneralModel.doseFunction(ModelPseudo): Return a dose function for aModelPseudo.doseFunction(LogisticLogNormalOrdinal): Return a grade-specific dose function for aLogisticLogNormalOrdinalmodel.
Examples
my_model <- LogisticLogNormal(
mean = c(-0.85, 1),
cov = matrix(c(1, -0.5, -0.5, 1), nrow = 2),
ref_dose = 50
)
dose_fun <- doseFunction(my_model, alpha0 = 2, alpha1 = 3)
dose_fun(0.6)
#> [1] 29.38579
data_ordinal <- .DefaultDataOrdinal()
model <- .DefaultLogisticLogNormalOrdinal()
options <- .DefaultMcmcOptions()
suppressWarnings({
samples <- mcmc(data_ordinal, model, options)
})
doseFunction(
model,
alpha1 = samples@data$alpha2,
beta = samples@data$beta,
grade = 1L
)(x = 0.75)
#> [1] 87.21638 145883.27473 61.51115 67.42361 103.45997
#> [6] 76.45095 81.32959 162.93052 68.93592 98.68414
#> [11] 72.53608 6308.21941 61.38387 407.71674 154.63225
#> [16] 77.81467 72.73968 144.46059 72.47730 108.95528
#> [21] 76.74244 124.76481 72.58418 72.43440 103.62518
#> [26] 107.15173 95.91364 106.69285 113.44893 121.37328
#> [31] 84.55551 74.32140 118.42807 239.18557 79.35191
#> [36] 486.03796 69.29754 76.07944 83.78211 70.26298
#> [41] 74.07065 70.73024 69.36088 97.22546 83.48185
#> [46] 66.59772 69.49560 65.98788 62.53784 67.89687
#> [51] 75.75824 99.97549 108.95176 80.97078 81.12301
#> [56] 76.81664 266.40535 191.00187 67.80633 64.40575
#> [61] 75.54535 74.19220 188.74474 63.75622 70.97663
#> [66] 118.04091 247.65846 116.61838 92.17013 76.11969
#> [71] 104.44091 71.32964 109.87446 62.69247 114.66445
#> [76] 112.52920 87.86024 175.91918 742.36029 122.16484
#> [81] 65.85714 521.48547 1192.23883 73.29130 268.07394
#> [86] 68.92603 101.53562 75.37157 81.19490 85.28073
#> [91] 79.56663 84.80991 311.85311 75.02409 64.78874
#> [96] 68.99022 77.57387 393.68292 1699.47199 70.79810
#> [101] 163.86766 249.58251 116.47469 204.55687 65.54925
#> [106] 82.56397 79.27404 66.98101 73.04404 71.40696
#> [111] 154.46411 125.36555 66.49939 69.84952 69.16011
#> [116] 143.74835 244.83449 163.54088 83.88655 102.62380
#> [121] 75.58247 79.86012 205.00974 119.31826 93.75077
#> [126] 72.19327 74.73208 90.03096 118.46810 66.78613
#> [131] 96.99334 65.18144 64.02512 66.27139 149.59031
#> [136] 64.29314 106.75712 77.70031 128.64576 139.34196
#> [141] 110.68578 101.10946 85.94094 72.56803 75.84171
#> [146] 70.33664 78.30061 60.40272 83.54489 67.33000
#> [151] 284.25787 76.67058 61.52102 92.37880 127.17029
#> [156] 323.13617 471.09415 438.24117 72.70181 70.12060
#> [161] 240.03955 13775.37434 103.72909 70.71704 440.93570
#> [166] 154.44702 6727.26082 357.37286 62.82073 74.91784
#> [171] 62.73283 68.10007 138.94232 65.42964 317.84249
#> [176] 203.47917 234.57339 75.00726 79.89840 102.79602
#> [181] 82.31726 68.98146 144.42374 170.58799 161.14808
#> [186] 79.58501 91.15671 84.15348 234.67851 91.32904
#> [191] 93.26457 112.05874 73.41156 65.31186 66.86841
#> [196] 114.24426 155.17440 145.28373 79.94168 70.50849
#> [201] 75.93759 125.34963 165.79136 15598.96495 3868.02112
#> [206] 417.71831 69.06289 81.80921 72.37131 76.19250
#> [211] 69.57600 75.00859 79.91638 60.60269 64.21463
#> [216] 70.59369 104.75394 106.47340 152.73394 198.16350
#> [221] 80.38139 167.96884 72.94353 60.50025 82.00724
#> [226] 73.69266 94.16440 80.79327 95.35273 89.07167
#> [231] 105.41238 1240.29967 199.26018 107.87014 109.44294
#> [236] 100.39395 192.52276 387.97554 273.12256 275.01890
#> [241] 214.91937 156.40045 149.22934 66.35761 82.91028
#> [246] 140.88219 70.73816 95.70330 1751.53640 128.12089
#> [251] 69.22895 223.53861 64.52821 533.71265 103.53131
#> [256] 71.12731 120.94826 97.52018 79.31878 64.53827
#> [261] 76.94918 135.35107 193.76161 89.99029 69.62409
#> [266] 75.61608 76.40735 158.43810 72.66719 185.53340
#> [271] 92.14974 107.24008 118.64598 108.04960 80.96666
#> [276] 62.12373 88584.95226 86.84368 91.79838 6249.21515
#> [281] 80.71196 67.39811 72.50300 91.78939 75.89280
#> [286] 88.14210 69.16082 62.32276 194.00630 128.15627
#> [291] 191.65450 154.50435 84.15657 83.97467 78.07352
#> [296] 144.73808 92.87177 96.08915 166.68358 644.16958
#> [301] 367.73909 497.97570 95.93921 84.19389 136.93728
#> [306] 99.26569 92.31463 82.68316 117.99574 77.73734
#> [311] 82.49048 222.80178 113.93629 121.93720 90.31608
#> [316] 81.75422 63.80854 2367.37126 78.33737 75.72620
#> [321] 132.44242 99.96257 97.34953 398.22032 60.81579
#> [326] 78.59718 86.71068 68.29055 67.42987 536.73826
#> [331] 96.87668 182.75463 111.32310 152.77962 74.72671
#> [336] 74.89639 79.83327 72.75529 72.20853 320.82215
#> [341] 64.25194 77.25131 104.25966 74.20035 108.25829
#> [346] 215.58731 108.20478 82.50650 72.69823 116.13429
#> [351] 84.29504 131.62447 62.41751 91.66671 90.93031
#> [356] 71.62015 77.46459 67.72775 85.91167 92.47168
#> [361] 111.11341 63.46845 157.34329 95.18399 60.91203
#> [366] 83.92534 101.35070 124.33862 65.85369 122.23441
#> [371] 86.83674 346.53798 66.15415 73.70480 67.54324
#> [376] 99.08080 90.62302 74.07747 62.96842 67.53690
#> [381] 194.46923 114.59527 91.71066 94.44658 85.46631
#> [386] 89.80133 65.37979 123.74534 91.57081 434.61994
#> [391] 162.28619 68.11065 200.59337 70.23207 66.74012
#> [396] 86.54374 64.15592 82.83262 74.84663 73.75711
#> [401] 67.14343 69.60185 466.13459 58.46370 73.91705
#> [406] 69.61814 90.04694 91.01352 82.09338 66.78630
#> [411] 252.93952 124.08712 222.70784 198.13850 188.61651
#> [416] 66.58679 70.41863 89.15580 75.90160 97.09060
#> [421] 198.34953 108.58380 105.48124 77.90505 71.83692
#> [426] 76.31364 106.68863 1625.88742 117.60214 77.11319
#> [431] 76.83185 65.50323 86.62545 110.52615 78.57910
#> [436] 176.81094 119.54763 66.81996 72.89257 88.64310
#> [441] 91.78462 63.25501 132.00932 117.94765 91.11596
#> [446] 70.97073 70.28617 64.86115 63.98964 86.73913
#> [451] 360.39259 140.00034 233.85329 1217.57723 98.45169
#> [456] 108.74713 167.58496 9107.11466 371.46435 81.76429
#> [461] 78.07562 98.59984 77.88734 76.96672 133.72599
#> [466] 355.48327 74.10448 65.96442 82.71892 231.69964
#> [471] 68.30549 81.52806 62.64894 77.43660 91.21182
#> [476] 596.12166 110.16849 71.72536 69.34967 66.12392
#> [481] 71.38198 76.60737 75.82747 69.55884 84.68951
#> [486] 92.06705 65.95798 253.80384 173.86841 156.55990
#> [491] 79.00510 70.84516 71.29012 96.93634 95.17831
#> [496] 73.65644 68.86009 65.49457 74.23687 112.44038
#> [501] 122.21023 132.00549 104.70673 177.61968 70.83048
#> [506] 80.44383 244.01659 152.48394 62.08424 81.40204
#> [511] 208.93317 515.92108 232.86427 91.32508 103.40125
#> [516] 83.59339 68.78849 230.90752 130.94416 69.47811
#> [521] 192.43701 90.92145 69.58434 461.13548 252.77090
#> [526] 104.99608 76.52180 77.76665 68.35077 77.04989
#> [531] 87.67599 82.69424 73.04195 86.70908 86.60700
#> [536] 117.30118 108.61642 81.42649 101.29409 120.27039
#> [541] 267.48437 101.78996 87.26889 117.84765 77.29013
#> [546] 82.60150 87.57114 79.06559 86.83011 71.27230
#> [551] 81.91691 361.98537 80.22301 62.41715 121.39121
#> [556] 66.33161 92.95770 143.71161 90.26983 75.61842
#> [561] 65.24063 100.73408 62.45342 129.33578 108.37066
#> [566] 83.38896 74.77442 67.55494 67.99162 75.17366
#> [571] 68.24890 69.50017 3152.46193 1207.61385 107.41828
#> [576] 68.55592 81.00199 64.79889 68.78934 111.02918
#> [581] 98.68047 68.04996 94.97782 76.60970 244.51016
#> [586] 14615.87635 73.53016 83.29785 72.72749 70.02137
#> [591] 338.18768 6301.95932 86.84973 124.03593 445.64248
#> [596] 59.40858 68.84397 62.96731 71.69099 73.77970
#> [601] 171.92254 84.08096 5432.07065 69.52190 601.96063
#> [606] 92.09410 75.16988 98.73820 9024.21808 69.55952
#> [611] 139.69483 61.35433 238.50738 65.61194 78.62888
#> [616] 67.33968 237.54670 84.36551 73.69538 371.29678
#> [621] 100.72415 107.55487 72.36362 62.62495 611.60978
#> [626] 69.17156 216.17023 66.93627 62.73812 62.87232
#> [631] 107.92988 81.66834 77.00891 66.74674 67.06209
#> [636] 107.30977 97.70925 131.93074 268.51446 84.75081
#> [641] 143.57515 125.80992 245.45166 376.43703 70.66630
#> [646] 141.35213 92.30467 832.19820 563180.01029 58.73901
#> [651] 66.55761 63.40337 95.27069 75.48589 112.36557
#> [656] 63.31909 65.70439 90.94923 75.00967 68.74565
#> [661] 79.17157 72.72991 97.85504 90.85107 70.29910
#> [666] 111.75599 88.87272 113.56327 797.61152 161.27273
#> [671] 70.05213 231.89376 171.03793 70.28410 158.97345
#> [676] 66.98018 1110.66519 201.30390 94.70942 106.20767
#> [681] 112.81336 91.08096 132.45751 239.31257 1698.86969
#> [686] 852961.66282 2646.95550 560.74700 74.16879 101.38233
#> [691] 72.17387 70.62854 73.84507 79.26351 180.29942
#> [696] 71.65804 64.97531 601.79253 199.68507 223.59831
#> [701] 71.67589 66.05294 93.75988 86.21241 182.81139
#> [706] 152.50288 101.90468 102.72226 89.18931 2218.86459
#> [711] 1934.67767 116.61244 73.53288 80.84488 96.21513
#> [716] 138.33648 99.44121 105.32917 75.82753 159.97519
#> [721] 72.47744 101.80121 80.75119 71.76905 166.10974
#> [726] 157.48186 155.30966 103.79648 549.01849 99.40604
#> [731] 110.61267 112.21936 141.65508 116.30223 112.94305
#> [736] 66.69603 69.91262 123.73568 835.68133 90.76425
#> [741] 70.96036 73.74569 81.32729 94.52803 118.81464
#> [746] 75.68989 68.81629 96.12153 711.75015 78.56549
#> [751] 73.48569 108.56174 414.22851 212.01395 76.65573
#> [756] 64.90262 197.01841 100.30929 67.32056 97.23008
#> [761] 68.90504 137.68307 139.75181 131.38484 103.55383
#> [766] 71.86620 74.17898 75.85382 1377.01788 124.30337
#> [771] 81.67322 191.41182 807.17500 186.96784 135.51444
#> [776] 71.06977 594.66359 69.79100 84.43058 67.00293
#> [781] 85.17162 97.83210 78.21512 65.82439 399.16155
#> [786] 103.37959 162.10918 81.08682 110.35800 91.34171
#> [791] 91.01225 65.32992 68.38709 531.93185 73.49531
#> [796] 79.22412 135.85244 320.05027 72.32388 67.15913
#> [801] 66.47447 65.21084 76.47611 199.11517 79.57119
#> [806] 75.28144 83.78354 135.38128 70.85753 83.68132
#> [811] 84.22987 155.84661 64.88147 65.60911 209.04852
#> [816] 154.41001 63.51906 195.78889 84.65021 88.25581
#> [821] 117.75217 138.60166 115.65350 75.46866 73.34425
#> [826] 69.44493 8575.02664 238.84453 99.05833 65.70489
#> [831] 109.92259 63.75750 70.99184 106.44819 69.61061
#> [836] 92.69109 78.37580 71.96525 69.61399 71.27193
#> [841] 416.22601 77.79567 99.34051 212.22865 166.32286
#> [846] 78.79118 89.67337 73.27543 75.59139 72.72706
#> [851] 94.13064 674.80698 72.53654 96.97547 81.16880
#> [856] 76.91996 150.78357 66.16605 64.60852 74.64685
#> [861] 64.10390 123.59481 125.87424 327.16878 190.96970
#> [866] 79.38194 66.54357 241.24645 66.85302 1071.21444
#> [871] 726.98190 444.07377 155.40606 1802.33569 64.38569
#> [876] 96.76312 76.36165 89.78048 73.19837 70.65453
#> [881] 226.77419 309.39486 105.28018 85.55911 78.39444
#> [886] 72.57689 76.45782 84.24452 167.81950 101.75511
#> [891] 129.95400 67.06682 74.28873 130.49149 67.81211
#> [896] 935.55114 65.36861 66.48794 148.48389 110.78016
#> [901] 86.27052 77.75530 65.04924 77.89665 71.68465
#> [906] 72.30116 64.23332 65.53187 73.79525 140.64264
#> [911] 82.55027 6020.96566 67.02360 62.76597 121.87051
#> [916] 70.48768 84.75841 69.15312 75.10422 77.31581
#> [921] 121.62007 67.14714 70.87750 76.70958 81.24843
#> [926] 78.73668 158.97094 127.15290 78.17060 70.90795
#> [931] 129.91243 1517.42501 109.95234 72.82609 89.79352
#> [936] 2474.45607 63.72291 71.56531 134.05641 247.26110
#> [941] 117.00557 143.76257 78.42172 144.94041 78.26436
#> [946] 68.28765 136.58015 70.33832 137.61982 88.60269
#> [951] 267.31026 91.09819 235.05962 73.39794 73.38845
#> [956] 67.57935 69.29578 66.31610 80.02472 78.62896
#> [961] 69.74742 81.37886 123.94688 240.86518 138.41225
#> [966] 6472.55933 2470.78037 69.42932 93.00525 3930.34379
#> [971] 67.75068 84.71488 78.07691 81.88180 119.39790
#> [976] 65.80736 64.76640 67.03671 68.69964 708.23007
#> [981] 74.71848 173.19896 206.72401 61.38209 76.14799
#> [986] 69.40527 81.99074 100.99642 134.35750 64.57614
#> [991] 75.87246 91.72190 66.18222 175.22060 77.03751
#> [996] 88.47058 91.18605 78.84270 70.42452 107.90910
doseFunction(
model,
alpha2 = samples@data$alpha2,
beta = samples@data$beta,
grade = 2L
)(x = 0.25)
#> [1] 65.21197 843.06245 54.91180 59.46359 80.11738 61.44950
#> [7] 67.55348 67.60611 57.00299 60.71239 60.38822 764.96364
#> [13] 53.64358 146.96703 68.46140 59.71143 58.89278 73.48954
#> [19] 59.58844 69.54342 59.70604 78.67429 59.10322 61.53681
#> [25] 79.42656 69.21230 64.98366 66.14547 53.14927 81.93078
#> [31] 60.25901 63.46004 79.14862 96.92995 66.11873 210.63847
#> [37] 49.54640 64.10306 66.65168 57.26996 55.36207 58.95420
#> [43] 55.52689 60.47548 53.79073 57.74549 62.69988 59.17314
#> [49] 56.01726 60.09883 67.45800 71.15739 76.87765 60.87022
#> [55] 62.19660 61.50237 97.78747 102.88565 61.73756 59.04157
#> [61] 64.23260 64.47990 111.90525 59.07525 64.09058 90.98112
#> [67] 104.21112 66.46949 56.94749 51.12430 66.82370 61.68284
#> [73] 75.03068 57.47646 79.60737 83.65642 68.32628 92.64582
#> [79] 221.99517 78.27465 60.15457 143.02570 182.97028 51.94346
#> [85] 127.60184 52.74754 66.66731 60.68738 65.22932 56.15054
#> [91] 64.76270 67.64862 143.82188 58.09027 55.66261 56.43977
#> [97] 61.83280 110.75240 109.45938 58.09270 110.06483 84.26894
#> [103] 87.66806 119.64758 55.93406 60.29976 64.27560 58.72806
#> [109] 62.89868 64.29060 77.17199 66.38325 54.81986 61.28116
#> [115] 59.08940 96.61561 128.28416 94.69202 65.32885 74.32976
#> [121] 59.51003 60.74668 101.96579 69.35960 61.02878 57.29071
#> [127] 59.85037 73.44914 83.02359 56.96288 68.13469 56.73180
#> [133] 56.10470 53.33610 59.57936 55.39648 78.01127 65.37032
#> [139] 88.09829 81.17624 61.66144 58.93760 56.69110 61.88284
#> [145] 64.81576 61.44766 58.40505 53.79967 68.90412 54.85044
#> [151] 160.58538 61.76085 53.39155 69.37271 84.51807 107.30269
#> [157] 178.61356 155.16377 56.64936 60.06652 127.55844 1775.95889
#> [163] 76.36506 57.77378 142.45850 83.20542 899.42982 148.02508
#> [169] 55.23494 64.88196 58.45627 61.27424 88.81118 53.09075
#> [175] 142.41917 90.46095 104.42631 56.83312 59.15192 74.37874
#> [181] 61.95118 59.46677 84.52294 86.75787 80.32434 62.15489
#> [187] 73.48928 63.28020 127.56176 67.99565 62.64078 73.30005
#> [193] 53.75235 56.87802 58.90100 73.87741 52.31273 93.14477
#> [199] 56.17483 55.00129 62.90732 75.47145 111.65222 2812.26486
#> [205] 364.97613 142.88102 59.98503 66.10974 58.77216 60.33858
#> [211] 56.21221 63.91442 58.45518 55.90151 57.33601 61.26070
#> [217] 71.34865 77.28573 61.93334 86.10280 64.65163 73.71751
#> [223] 63.04605 53.17526 65.74944 59.93173 74.39018 58.98224
#> [229] 75.34164 53.69541 67.15014 178.71081 94.55627 68.57189
#> [235] 51.72776 59.12330 104.27587 104.06203 117.58918 98.63023
#> [241] 105.36884 88.37760 100.75190 53.23187 64.64926 81.34985
#> [247] 59.38970 73.66111 271.97709 91.76460 61.28015 64.96768
#> [253] 57.40782 215.60444 69.70189 52.80676 84.46464 66.98508
#> [259] 62.30272 54.32957 62.12966 74.37518 112.06749 54.32996
#> [265] 59.42185 62.14331 57.45612 92.75846 41.62810 108.34964
#> [271] 71.96436 79.70237 84.04205 77.76939 61.24577 50.98234
#> [277] 7590.30735 50.00575 64.92731 427.29014 53.66548 58.71724
#> [283] 61.01363 65.84034 57.14648 69.84903 60.73359 54.45547
#> [289] 114.51751 84.91381 76.95256 72.19959 62.17390 66.31209
#> [295] 58.66258 101.83763 73.74043 74.52933 82.73144 182.69644
#> [301] 145.17223 84.77223 55.51862 58.89734 92.30961 78.73495
#> [307] 66.96706 67.58776 85.53442 65.80546 55.61989 100.98559
#> [313] 59.93336 69.14006 53.28649 62.03232 52.46130 397.87973
#> [319] 63.43735 64.53308 99.93627 82.71069 70.62488 91.55689
#> [325] 54.80489 63.13770 53.54859 58.19141 59.00171 174.98876
#> [331] 71.35081 87.03285 82.56981 94.64784 64.59526 60.07734
#> [337] 65.65093 58.40425 58.19727 73.00517 56.59691 64.66039
#> [343] 66.86842 55.60167 72.73091 113.91986 57.56651 68.96473
#> [349] 57.89238 79.02053 66.86750 100.07728 55.96351 73.99656
#> [355] 75.16457 57.00135 65.90151 60.94348 67.59158 68.98690
#> [361] 68.83484 55.65565 97.61612 58.34286 56.63440 68.80611
#> [367] 67.26433 79.67196 58.75460 94.87581 74.02464 168.22321
#> [373] 50.37131 58.54335 58.36470 73.79301 62.73731 64.50380
#> [379] 56.14950 59.52272 103.10110 76.87414 66.10383 71.66419
#> [385] 56.33636 69.02998 51.75968 84.06037 69.43516 116.63566
#> [391] 67.85077 58.03704 88.60899 57.67914 60.41707 72.89041
#> [397] 59.26735 73.59616 67.87339 67.14696 62.58586 64.11889
#> [403] 177.79523 50.77625 65.66225 61.86008 59.23290 47.58923
#> [409] 68.29252 57.61699 86.34969 88.46765 119.98343 69.62955
#> [415] 83.79758 59.86655 59.51557 66.44406 61.85049 62.26498
#> [421] 92.21090 73.16013 73.58944 53.16291 60.76414 62.73328
#> [427] 76.07155 333.30200 78.54313 63.17834 62.46558 54.92410
#> [433] 65.07892 62.88554 56.41226 109.95063 83.11672 60.85754
#> [439] 66.12604 65.36085 64.15140 54.19582 75.42725 83.57813
#> [445] 63.57837 55.85433 60.21892 59.09302 59.67380 69.52452
#> [451] 68.70308 79.06803 134.41295 111.52645 70.51110 67.83844
#> [457] 87.09678 1320.20489 154.52619 63.42812 58.81843 53.66840
#> [463] 57.20669 57.95034 89.85542 136.90084 62.89812 58.62176
#> [469] 65.75220 106.02712 62.07520 58.08293 57.59503 66.14783
#> [475] 62.75108 199.47982 75.67509 63.38316 61.69266 57.35229
#> [481] 59.93274 65.33096 62.55858 61.17110 69.90776 67.48553
#> [487] 55.31754 135.55544 82.27697 51.84189 61.94416 63.14714
#> [493] 59.77391 70.56130 46.50670 47.73933 60.99701 59.83135
#> [499] 62.12676 82.60543 73.12232 90.51983 71.85622 90.70325
#> [505] 61.19057 67.70534 136.56198 95.97760 53.18237 54.77871
#> [511] 100.93948 152.69806 105.64590 69.15699 52.27316 60.94054
#> [517] 57.19451 99.46888 79.82100 55.70553 108.27733 67.35289
#> [523] 61.20142 135.82437 95.33139 63.70749 58.14168 50.30503
#> [529] 58.74928 61.34268 65.27114 66.30937 59.84011 71.23491
#> [535] 69.04662 88.16112 86.97031 68.68338 78.54133 79.92048
#> [541] 116.89160 78.19306 73.14952 85.94618 64.95335 66.25420
#> [547] 67.76494 65.84545 66.14404 60.44368 68.45839 124.85091
#> [553] 61.37500 54.90363 74.33480 54.46733 64.08249 85.17355
#> [559] 65.02625 65.22580 58.90600 77.08506 58.06179 79.85637
#> [565] 76.33290 67.92344 62.65556 57.57417 57.75693 66.15428
#> [571] 60.93034 62.48121 733.92204 188.77708 67.18191 60.68296
#> [577] 63.62103 56.12037 60.91334 67.42539 63.68381 53.11600
#> [583] 56.40561 51.41452 79.13094 781.85163 58.15924 67.49872
#> [589] 62.98968 58.44151 149.03131 655.33505 53.53497 83.56963
#> [595] 191.63303 51.23884 58.60586 55.54438 61.85952 62.88944
#> [601] 66.93743 64.06613 951.88578 51.80278 200.45370 61.83573
#> [607] 60.49242 74.89231 725.85749 58.22252 94.25718 55.61648
#> [613] 108.69834 54.14736 53.30707 56.17347 93.10175 65.91147
#> [619] 61.94425 172.24768 62.69456 74.35268 60.86794 56.21407
#> [625] 73.20216 53.64618 120.05730 60.08456 57.57647 58.30218
#> [631] 81.63680 70.00330 66.03300 57.85178 58.19203 73.66600
#> [637] 66.93374 78.07954 113.70061 53.92133 53.29888 70.62944
#> [643] 72.34250 141.63959 58.56046 89.31673 65.01436 284.56163
#> [649] 7837.11770 53.91643 58.78734 57.03624 60.20335 52.58766
#> [655] 74.70886 56.57534 58.58122 71.60373 62.08601 60.13751
#> [661] 65.79361 62.48046 65.47049 71.60587 57.28232 67.13266
#> [667] 66.91457 82.56213 274.65397 105.38585 59.44540 122.01513
#> [673] 103.19969 63.97987 111.61237 60.01777 231.80401 94.71459
#> [679] 71.42326 71.29456 53.46082 56.82121 69.60924 117.20354
#> [685] 100.18496 58429.11543 268.38373 182.17301 57.15872 77.29896
#> [691] 58.24209 58.36445 60.31145 61.55760 108.41674 62.79482
#> [697] 57.09223 198.48029 91.44184 104.05275 62.74892 60.03932
#> [703] 66.84301 60.30542 102.17956 88.48883 61.54564 71.44037
#> [709] 59.42040 152.29123 305.70770 76.10825 58.57978 61.86977
#> [715] 70.83403 72.44770 78.58217 71.60453 60.52860 95.05785
#> [721] 60.06495 67.00163 58.82318 59.50512 94.00709 73.99994
#> [727] 99.09406 77.13076 155.50922 72.34864 68.90702 65.32671
#> [733] 72.86630 60.91128 74.98038 57.96328 57.05537 65.82421
#> [739] 213.67436 56.41672 59.01678 62.99104 71.48186 79.12651
#> [745] 88.77440 56.10012 57.92688 66.53337 167.84663 56.42041
#> [751] 52.13931 60.93101 137.87307 96.34452 58.97467 54.78718
#> [757] 126.44154 72.85088 57.53379 60.47539 55.00877 75.08217
#> [763] 58.11285 89.61214 70.84047 58.49713 61.79245 58.75812
#> [769] 217.18843 67.72707 65.64922 102.30693 323.30375 120.74337
#> [775] 92.25223 57.12964 164.39237 56.57366 58.17921 56.91173
#> [781] 62.67683 65.68742 60.94167 50.49567 105.30745 57.38270
#> [787] 99.28755 61.68593 65.19038 67.42703 62.74027 56.82603
#> [793] 60.19518 201.27871 61.34597 56.06125 81.10834 124.58284
#> [799] 61.83366 60.11638 59.60228 59.68135 62.80681 88.68749
#> [805] 59.55147 62.46833 70.04278 89.68727 56.42841 50.07057
#> [811] 55.69456 111.49136 59.67867 59.81861 111.92893 91.69303
#> [817] 55.33716 117.80392 66.73429 70.25277 77.30216 79.08206
#> [823] 62.13756 58.83905 61.18873 54.59134 330.29489 72.25125
#> [829] 73.62941 55.38548 61.47054 56.15293 61.21913 76.30317
#> [835] 59.81242 67.43425 60.34672 61.42951 58.82948 60.77817
#> [841] 180.17494 63.21414 68.68328 110.69401 89.84237 59.34879
#> [847] 62.56905 57.20404 57.30951 54.96887 57.44278 215.59780
#> [853] 55.25236 61.54977 63.65688 60.99845 88.12142 58.92821
#> [859] 59.02695 62.93304 56.18859 76.88400 80.33608 153.67529
#> [865] 95.69489 60.29580 60.01635 130.96704 55.80461 320.13024
#> [871] 196.51621 104.20176 101.39174 268.79713 55.20373 66.84231
#> [877] 53.45878 62.32284 62.55501 59.07737 126.07555 170.57717
#> [883] 69.79372 66.68013 60.89525 48.68322 39.94208 72.27006
#> [889] 100.81271 70.11935 82.06684 60.17531 63.56123 77.31715
#> [895] 60.48821 217.74282 57.14600 57.72333 95.17836 84.56504
#> [901] 68.13669 66.02153 58.37415 65.11811 62.98225 64.04801
#> [907] 58.25733 58.52903 55.00054 84.32325 56.56071 199.61926
#> [913] 56.32508 53.87341 70.66922 57.98750 62.13644 60.78144
#> [919] 55.06956 53.99712 77.73262 55.42373 60.92302 63.29162
#> [925] 67.96424 65.76617 62.18289 64.83233 63.27259 57.22081
#> [931] 82.68846 411.17267 73.31382 60.86371 70.06223 575.14669
#> [937] 57.24891 63.02797 80.29135 127.98334 65.98440 77.38381
#> [943] 60.64306 78.42763 56.97633 59.28906 74.14120 49.57173
#> [949] 57.65968 65.28543 103.96871 55.97226 106.31148 59.84390
#> [955] 63.79613 57.13538 63.46750 59.89870 66.51259 63.21455
#> [961] 59.52999 62.44693 88.84838 101.13497 69.43033 700.29308
#> [967] 317.51687 59.99167 69.79348 311.56038 57.34655 69.11976
#> [973] 64.50135 66.48910 74.04854 59.26322 57.73889 59.89110
#> [979] 57.78428 193.93255 61.86580 100.78408 128.50793 56.46840
#> [985] 64.89316 61.94599 69.57016 64.36152 84.65711 58.07816
#> [991] 63.96216 71.79084 58.33739 97.90740 66.98015 74.64697
#> [997] 74.22033 65.06759 57.20611 73.67011