![[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):doseFunction(ModelPseudo):doseFunction(LogisticLogNormalOrdinal):
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] 8.265336e+01 4.373353e+02 7.446676e+01 1.193286e+04 6.530749e+02
#> [6] 1.843335e+02 2.338013e+02 9.168704e+01 7.383717e+01 1.125251e+02
#> [11] 6.910355e+01 6.334124e+01 1.211475e+02 8.471525e+01 6.275211e+01
#> [16] 7.900889e+01 8.117329e+02 9.909791e+01 2.318247e+02 6.256709e+02
#> [21] 2.058261e+02 2.054242e+03 1.044678e+02 7.321810e+01 1.230579e+02
#> [26] 1.011982e+02 2.238782e+02 6.608516e+01 6.721764e+01 1.043832e+02
#> [31] 2.346331e+02 2.184986e+02 1.487542e+02 6.393092e+01 8.425008e+01
#> [36] 2.104173e+02 1.055264e+02 3.317034e+02 8.005570e+01 7.752652e+01
#> [41] 1.837759e+02 7.561675e+01 6.907092e+01 2.911233e+02 7.455197e+01
#> [46] 8.032300e+01 8.046106e+01 9.998230e+01 6.428024e+01 6.667111e+01
#> [51] 7.525124e+01 7.319962e+01 1.231414e+02 6.284616e+01 1.819501e+02
#> [56] 7.712550e+01 6.583398e+01 1.382527e+02 1.930196e+02 6.780618e+01
#> [61] 7.282428e+01 6.853658e+01 1.887692e+02 8.239530e+01 1.874902e+02
#> [66] 6.684849e+01 1.339231e+02 7.851763e+01 7.304714e+01 7.421228e+01
#> [71] 7.214067e+01 8.122994e+01 1.066762e+02 1.115557e+02 2.339823e+02
#> [76] 9.387858e+01 6.576413e+01 6.328661e+01 7.638161e+01 6.820731e+01
#> [81] 7.026488e+01 8.825408e+01 8.642305e+01 3.113961e+02 6.553709e+01
#> [86] 7.946370e+01 7.006063e+01 1.929783e+02 8.758227e+01 7.893392e+01
#> [91] 1.799120e+02 1.278256e+02 8.258697e+01 7.471219e+01 7.370186e+01
#> [96] 7.641603e+01 6.414923e+01 7.459480e+01 8.639305e+01 2.797549e+02
#> [101] 8.081809e+01 2.289124e+02 1.182275e+02 1.686644e+02 1.134736e+02
#> [106] 1.072273e+02 7.514814e+01 1.180979e+02 1.162988e+02 1.557862e+02
#> [111] 1.449258e+02 6.177080e+01 5.986836e+01 7.102464e+01 7.471960e+01
#> [116] 7.107914e+01 1.158109e+02 6.864649e+01 9.814073e+01 2.652431e+02
#> [121] 6.491086e+02 8.590356e+01 3.833009e+02 6.822690e+01 8.592516e+01
#> [126] 6.612620e+01 1.963230e+02 7.274641e+01 1.671654e+02 9.036604e+01
#> [131] 3.820789e+02 2.084639e+02 6.667826e+01 8.421347e+01 1.003250e+02
#> [136] 1.013734e+02 1.028733e+02 7.990153e+02 8.950555e+01 7.117588e+01
#> [141] 7.771878e+01 7.861753e+01 2.932823e+02 6.028093e+02 7.780613e+01
#> [146] 1.259808e+02 3.684574e+03 8.362575e+01 6.703445e+01 7.339913e+01
#> [151] 1.443509e+02 6.633386e+01 1.002996e+02 7.340179e+01 7.533300e+01
#> [156] 2.164941e+04 2.775475e+04 1.798902e+02 8.857899e+01 6.379162e+01
#> [161] 6.545096e+01 6.986489e+01 6.922926e+01 1.641015e+02 6.382332e+01
#> [166] 6.664030e+01 2.396753e+02 8.254090e+01 8.671170e+01 8.742353e+01
#> [171] 8.128814e+01 6.445849e+01 7.884657e+01 1.039777e+02 6.355647e+01
#> [176] 1.840816e+03 8.544352e+01 1.628225e+02 6.816790e+01 8.479980e+01
#> [181] 8.117719e+01 6.422535e+01 7.405237e+01 1.296548e+02 7.701432e+01
#> [186] 6.537353e+01 6.829138e+01 2.864307e+04 7.265633e+01 9.894338e+01
#> [191] 7.016155e+01 7.231868e+01 1.428725e+02 8.524481e+01 8.457489e+01
#> [196] 1.046641e+02 1.241338e+02 6.293224e+01 8.659141e+01 6.677961e+01
#> [201] 8.600744e+03 3.960650e+02 1.033352e+02 6.154911e+01 6.845080e+01
#> [206] 6.519799e+01 8.066923e+01 6.669483e+01 6.272378e+01 6.306307e+01
#> [211] 6.927146e+01 7.359538e+01 8.618080e+01 1.344335e+02 1.945333e+02
#> [216] 1.119665e+02 9.246521e+01 7.066430e+01 7.925480e+01 9.175357e+01
#> [221] 8.843777e+01 8.127978e+01 7.078922e+01 1.799016e+02 9.826447e+01
#> [226] 7.085014e+01 1.017113e+02 8.567551e+01 6.716081e+01 8.397842e+01
#> [231] 6.566629e+02 1.410158e+02 6.777259e+01 7.712257e+01 7.855592e+01
#> [236] 7.226289e+01 6.592647e+01 7.148169e+01 6.455250e+01 7.081755e+01
#> [241] 7.632982e+01 8.634848e+01 2.942437e+02 2.650180e+02 1.880886e+02
#> [246] 5.833620e+02 7.970040e+01 1.047648e+02 1.523809e+02 8.121187e+01
#> [251] 9.759366e+01 7.062881e+01 7.612857e+01 1.000575e+02 9.354771e+01
#> [256] 1.192098e+02 3.618067e+02 6.817990e+01 8.991559e+01 8.174716e+01
#> [261] 7.034929e+01 1.780972e+02 1.170203e+02 1.434641e+02 8.237080e+01
#> [266] 1.176985e+02 9.699195e+01 7.984157e+01 9.545363e+01 8.219723e+01
#> [271] 2.619396e+02 8.828578e+01 8.108899e+01 7.257885e+01 4.537968e+02
#> [276] 6.777921e+02 8.024768e+01 8.863983e+01 7.013614e+01 6.420537e+01
#> [281] 3.391583e+04 1.167153e+02 7.467475e+01 7.368718e+01 8.927074e+01
#> [286] 8.648888e+01 1.048809e+02 1.273890e+02 1.009083e+02 1.160441e+02
#> [291] 1.886733e+02 1.404175e+02 9.459445e+01 6.564855e+01 8.142820e+01
#> [296] 1.102752e+02 6.489854e+01 1.909936e+02 8.742364e+01 1.479841e+02
#> [301] 1.532591e+02 1.123066e+02 2.176397e+02 1.458697e+02 8.597741e+01
#> [306] 1.029323e+02 6.560535e+01 7.898028e+01 7.262080e+01 7.458035e+01
#> [311] 8.054386e+01 1.126288e+02 5.060498e+02 1.116187e+02 1.108379e+02
#> [316] 8.187492e+01 6.376084e+01 1.037748e+02 6.579762e+01 7.176950e+01
#> [321] 1.134976e+02 1.270725e+02 1.987424e+02 1.398917e+02 9.198075e+01
#> [326] 6.782058e+02 1.122488e+04 7.010177e+01 1.332117e+02 1.072020e+03
#> [331] 9.782615e+01 1.213072e+02 5.764449e+01 6.226888e+01 1.002497e+02
#> [336] 7.626672e+01 8.402120e+01 1.983392e+02 2.104623e+02 6.574715e+01
#> [341] 9.117148e+01 6.826615e+01 9.306082e+01 6.563894e+01 5.441061e+02
#> [346] 8.895696e+01 6.590320e+01 7.441566e+01 1.180158e+02 7.314892e+01
#> [351] 8.605712e+01 6.672043e+01 6.459102e+01 5.310103e+02 5.136398e+04
#> [356] 1.242501e+02 1.390306e+02 9.446569e+01 9.607201e+01 7.246973e+01
#> [361] 7.350573e+01 8.478121e+01 1.066100e+02 3.353014e+02 1.068368e+02
#> [366] 6.511948e+01 9.637387e+01 7.775790e+01 7.936708e+01 1.667523e+02
#> [371] 6.397736e+01 6.804476e+01 1.068264e+02 8.827991e+01 9.046662e+01
#> [376] 9.972547e+01 6.322293e+01 1.285828e+02 1.390233e+02 9.425200e+01
#> [381] 7.110804e+01 3.233053e+02 8.241040e+01 9.076489e+01 7.408661e+01
#> [386] 6.000420e+01 1.306301e+04 6.890025e+03 6.624796e+01 9.682430e+01
#> [391] 2.526606e+05 6.035081e+01 6.997512e+01 9.702851e+01 7.168361e+01
#> [396] 6.279369e+01 6.310181e+01 6.678392e+01 6.643919e+01 7.434987e+01
#> [401] 6.796415e+01 6.143762e+01 6.957343e+01 9.999275e+01 6.447731e+01
#> [406] 1.436477e+02 1.129151e+02 1.062571e+02 6.655822e+01 8.706566e+01
#> [411] 7.335009e+01 2.786087e+02 1.077214e+02 1.194235e+02 7.167018e+01
#> [416] 1.047974e+02 8.492269e+01 8.036927e+01 6.979156e+01 6.365757e+01
#> [421] 1.437691e+05 2.053551e+05 8.681764e+01 6.871332e+01 7.398645e+01
#> [426] 1.277718e+02 1.248757e+02 1.355258e+02 7.113297e+01 7.697674e+01
#> [431] 6.638986e+01 9.033745e+01 7.121332e+01 8.504047e+01 6.781493e+01
#> [436] 7.215753e+01 2.845553e+02 7.436225e+01 8.503321e+01 3.959701e+04
#> [441] 2.000998e+02 8.932119e+01 7.258604e+01 8.906119e+01 8.456073e+01
#> [446] 9.183202e+01 6.172744e+01 7.749378e+01 7.745411e+01 8.111485e+01
#> [451] 1.193582e+02 6.693805e+01 1.171501e+02 1.117940e+02 1.123696e+02
#> [456] 7.635446e+01 1.111243e+02 9.573860e+01 1.085177e+02 3.868458e+02
#> [461] 5.994964e+02 1.067215e+02 7.247493e+01 6.110334e+01 2.094467e+02
#> [466] 6.488524e+01 9.454745e+01 7.820212e+01 1.136482e+02 1.054072e+02
#> [471] 9.486870e+01 1.357145e+02 1.266803e+02 8.718673e+01 1.773050e+02
#> [476] 1.008617e+02 6.648940e+01 7.104794e+01 8.783310e+01 3.504102e+02
#> [481] 1.525332e+03 5.543258e+02 3.207263e+02 1.648395e+03 6.554461e+01
#> [486] 7.071597e+01 1.085899e+02 1.449622e+02 8.958430e+01 7.749629e+01
#> [491] 8.707406e+01 8.927968e+01 8.149258e+01 1.122329e+02 8.470344e+01
#> [496] 6.682570e+01 3.703096e+02 3.989659e+02 9.676614e+01 1.636498e+02
#> [501] 8.588833e+01 3.532727e+02 8.648710e+01 9.762016e+01 2.916184e+02
#> [506] 6.909808e+01 1.115389e+02 1.957427e+02 2.537050e+02 1.750704e+02
#> [511] 8.044003e+01 8.606674e+01 2.038131e+02 7.100689e+01 7.995688e+01
#> [516] 7.964333e+01 7.893861e+01 1.213902e+02 8.184003e+01 6.634930e+01
#> [521] 1.815388e+02 2.074583e+02 7.112296e+01 9.543421e+01 2.706944e+02
#> [526] 6.336525e+01 2.023495e+02 6.659835e+01 8.620700e+01 1.198270e+02
#> [531] 1.481802e+02 7.490505e+01 7.771202e+01 3.880400e+02 6.720995e+01
#> [536] 6.909308e+01 7.958365e+01 7.415335e+01 8.200910e+01 7.611567e+01
#> [541] 7.660371e+01 6.226100e+01 6.489092e+01 7.374419e+01 8.836821e+01
#> [546] 6.940365e+01 8.394402e+01 7.049014e+06 2.488414e+02 6.955398e+01
#> [551] 7.287079e+01 6.392369e+01 6.773836e+01 7.178760e+01 9.103947e+01
#> [556] 4.149774e+03 8.861833e+01 6.788067e+01 7.010196e+01 7.927638e+01
#> [561] 9.396124e+01 9.479509e+01 1.121747e+02 8.535721e+01 7.265687e+01
#> [566] 6.543233e+01 6.370232e+01 1.482124e+02 6.832594e+01 9.257500e+01
#> [571] 6.557833e+01 1.051219e+03 2.041085e+04 2.449089e+02 9.114890e+01
#> [576] 7.654414e+01 1.348652e+02 7.864349e+01 1.213901e+02 7.304954e+01
#> [581] 7.011452e+01 8.723284e+01 7.717319e+01 6.481667e+01 6.573735e+01
#> [586] 9.768449e+01 9.409917e+01 1.159221e+02 1.369179e+02 7.150367e+01
#> [591] 8.761400e+01 6.923520e+01 7.727923e+01 7.815866e+01 6.887432e+01
#> [596] 1.042537e+02 6.938652e+01 6.914660e+01 7.009790e+01 7.452112e+01
#> [601] 1.140144e+02 7.472707e+01 7.219069e+01 1.619779e+02 1.110874e+02
#> [606] 6.820006e+01 9.591210e+01 1.349744e+02 8.356970e+01 1.133730e+02
#> [611] 4.333492e+02 6.791674e+01 7.299481e+01 1.365807e+02 6.066737e+01
#> [616] 7.075591e+01 1.040176e+02 9.290874e+01 7.621800e+01 7.144427e+01
#> [621] 9.446211e+01 7.327532e+01 1.120682e+02 1.169253e+02 9.325888e+01
#> [626] 6.542299e+01 2.488012e+02 9.482448e+01 1.283010e+02 8.455767e+01
#> [631] 7.649772e+01 7.047615e+01 7.799441e+01 7.643093e+01 8.655283e+02
#> [636] 1.430960e+02 1.510423e+02 1.057369e+03 1.932331e+03 1.136176e+02
#> [641] 6.714881e+01 7.485812e+01 1.510751e+02 1.023354e+03 1.491642e+02
#> [646] 3.350002e+02 8.275639e+01 1.334193e+02 6.365560e+14 2.175786e+02
#> [651] 2.176425e+02 8.583487e+01 6.570227e+01 1.059421e+02 1.267805e+02
#> [656] 1.099734e+02 6.648479e+01 7.325233e+01 1.136757e+02 6.116321e+02
#> [661] 7.917292e+01 3.705561e+02 5.870913e+01 9.620863e+01 8.530340e+01
#> [666] 2.564456e+02 6.776526e+01 7.131565e+01 1.375866e+02 9.333564e+01
#> [671] 1.525954e+02 1.881878e+02 7.375428e+01 1.413164e+02 5.541548e+02
#> [676] 6.511101e+01 9.267918e+01 7.084628e+01 1.032755e+02 7.383801e+01
#> [681] 1.008743e+02 1.133746e+02 6.813347e+01 1.434155e+02 6.982551e+01
#> [686] 5.140562e+02 9.676732e+01 3.829802e+02 7.252291e+01 4.379504e+02
#> [691] 5.964685e+01 6.722970e+02 1.875789e+02 2.469256e+03 1.123183e+02
#> [696] 3.367231e+02 1.717403e+02 6.987050e+01 8.929163e+01 5.909369e+02
#> [701] 1.095347e+02 7.446968e+01 1.793708e+02 2.049428e+02 8.047581e+01
#> [706] 7.298158e+01 6.836821e+01 7.378225e+01 7.542144e+01 6.584642e+02
#> [711] 7.488896e+01 7.409175e+01 6.687601e+01 6.526492e+01 6.999885e+01
#> [716] 9.668825e+01 1.924011e+02 6.631592e+01 1.927178e+02 8.854585e+01
#> [721] 1.036564e+02 7.599173e+01 2.273153e+02 1.151291e+02 2.043605e+02
#> [726] 7.616122e+01 1.403342e+02 1.073838e+02 8.232758e+01 8.950544e+01
#> [731] 6.650241e+01 1.176155e+02 1.393223e+02 1.997919e+02 2.197927e+03
#> [736] 8.185790e+01 8.898520e+01 8.172400e+01 7.326531e+01 1.396024e+02
#> [741] 7.855106e+01 1.050153e+02 1.405160e+02 7.822856e+01 2.589979e+02
#> [746] 5.739583e+03 8.131769e+01 1.075021e+02 3.187400e+02 8.185986e+01
#> [751] 6.499444e+01 1.219116e+02 7.639500e+01 4.911975e+02 6.469723e+01
#> [756] 3.966820e+02 7.745113e+02 7.716770e+01 2.034380e+02 5.597968e+02
#> [761] 2.372317e+03 1.692539e+02 1.790528e+02 2.781232e+02 9.098867e+01
#> [766] 7.312907e+01 6.648442e+01 9.699935e+01 1.638098e+02 9.826368e+01
#> [771] 3.085388e+02 8.179669e+01 8.390368e+01 6.391822e+01 2.719887e+03
#> [776] 7.167778e+01 7.201794e+01 8.362747e+01 8.194474e+01 8.313007e+01
#> [781] 7.264018e+01 1.049917e+02 1.294563e+03 1.160591e+02 9.053761e+01
#> [786] 7.934206e+01 8.125155e+01 6.331455e+01 5.461751e+02 2.180947e+04
#> [791] 2.390873e+05 3.526797e+02 9.958475e+01 1.185059e+02 6.179048e+01
#> [796] 9.454701e+01 6.346567e+01 2.369417e+02 7.687314e+01 8.218615e+02
#> [801] 1.292955e+02 8.108004e+01 9.051093e+01 1.988272e+04 6.210152e+01
#> [806] 7.915829e+01 7.944382e+01 8.646008e+01 1.066326e+02 7.622727e+01
#> [811] 7.085050e+01 8.677253e+01 6.756836e+01 7.864817e+01 9.401419e+01
#> [816] 6.818264e+01 5.449587e+02 9.692765e+01 1.971865e+02 1.067313e+02
#> [821] 8.136841e+01 7.560886e+01 7.204232e+01 1.739020e+02 1.891364e+02
#> [826] 9.017260e+01 6.149439e+01 2.290925e+02 1.640381e+02 6.552404e+01
#> [831] 1.270828e+02 6.388283e+01 6.689087e+01 6.114690e+01 6.453406e+01
#> [836] 1.350660e+02 8.186972e+01 1.342676e+03 1.884633e+02 1.769228e+02
#> [841] 6.712611e+01 9.489666e+01 9.338620e+01 7.322667e+01 1.126727e+02
#> [846] 6.898274e+01 6.581552e+01 7.931251e+01 7.916628e+01 9.194439e+01
#> [851] 6.708373e+02 7.602036e+01 1.831606e+02 7.163704e+01 1.614809e+02
#> [856] 1.401308e+02 2.143748e+02 6.194132e+01 1.081715e+02 3.065419e+02
#> [861] 7.972903e+01 3.610787e+04 7.587504e+01 7.079447e+01 1.026889e+02
#> [866] 7.383973e+01 7.779190e+01 7.727374e+01 6.561282e+01 1.207004e+02
#> [871] 7.679826e+01 8.277711e+01 1.109706e+02 1.496165e+02 7.455632e+01
#> [876] 1.444653e+02 1.458612e+05 2.479686e+02 6.977060e+01 7.190748e+01
#> [881] 7.942743e+01 8.256476e+01 7.243279e+01 8.827995e+01 1.104836e+02
#> [886] 7.540153e+01 7.577309e+01 1.038691e+02 8.441591e+01 1.242996e+02
#> [891] 2.348594e+02 7.206122e+01 7.062523e+01 7.154905e+01 8.640523e+01
#> [896] 3.691567e+02 1.849517e+02 8.838988e+01 6.843694e+01 9.172286e+01
#> [901] 9.499622e+01 7.090845e+01 8.313795e+01 8.654483e+01 6.874792e+01
#> [906] 8.256345e+01 3.123763e+02 7.037769e+01 9.251913e+01 9.063126e+01
#> [911] 6.529718e+01 7.763725e+01 7.192378e+01 6.595504e+01 7.265991e+01
#> [916] 8.261354e+01 7.606897e+01 7.579784e+01 2.464078e+02 6.866342e+01
#> [921] 7.127633e+01 7.153705e+01 6.848714e+01 7.939975e+01 4.898835e+03
#> [926] 7.600691e+01 6.976933e+01 7.825347e+01 6.274450e+01 8.059074e+01
#> [931] 9.394203e+01 6.816749e+01 1.055213e+02 9.294234e+01 1.039821e+02
#> [936] 1.325784e+02 1.150805e+02 1.404908e+02 6.877556e+01 6.416925e+01
#> [941] 8.452664e+01 8.575512e+01 1.745468e+02 7.835886e+01 6.798030e+01
#> [946] 6.712832e+01 8.105488e+02 1.660619e+02 9.792193e+01 7.668657e+01
#> [951] 7.271785e+01 8.471606e+01 7.130070e+01 7.264483e+01 2.876759e+02
#> [956] 7.571665e+01 7.315582e+01 1.200260e+02 6.835561e+01 6.989760e+01
#> [961] 7.955451e+01 1.009174e+02 8.453029e+01 1.345520e+02 8.171360e+01
#> [966] 1.091152e+02 6.091758e+01 6.476661e+01 7.744898e+01 6.702291e+01
#> [971] 2.488895e+02 1.656385e+02 7.722116e+01 1.168539e+02 6.619603e+01
#> [976] 1.703761e+03 1.724146e+03 5.600782e+03 1.614745e+03 7.626680e+01
#> [981] 2.097435e+02 1.464701e+02 1.128155e+02 8.908379e+01 6.098406e+01
#> [986] 6.345085e+01 7.306908e+03 7.742116e+01 5.557791e+02 4.523514e+02
#> [991] 7.535520e+01 3.435549e+02 8.460989e+01 8.978909e+01 1.909683e+02
#> [996] 6.143049e+01 8.211561e+01 9.463174e+01 1.152613e+02 1.232190e+02
doseFunction(model, alpha2 = samples@data$alpha2, beta = samples@data$beta, grade = 2L)(x = 0.25)
#> [1] 6.960195e+01 9.932385e+01 6.075968e+01 1.363302e+03 2.063804e+02
#> [6] 1.023067e+02 9.699768e+01 6.554945e+01 5.741499e+01 8.072601e+01
#> [11] 5.768438e+01 5.716629e+01 9.032426e+01 6.833073e+01 5.749425e+01
#> [16] 6.576595e+01 2.200031e+02 6.286792e+01 9.834050e+01 7.855329e+01
#> [21] 1.108311e+02 2.068332e+02 6.965460e+01 5.607804e+01 7.720218e+01
#> [26] 7.081063e+01 1.276347e+02 5.493108e+01 6.000848e+01 7.844908e+01
#> [31] 9.848487e+01 8.751915e+01 7.765789e+01 5.383185e+01 5.988683e+01
#> [36] 1.061693e+02 7.226730e+01 1.282435e+02 6.092016e+01 6.413435e+01
#> [41] 9.374873e+01 6.265631e+01 5.805641e+01 1.571965e+02 5.071083e+01
#> [46] 5.193307e+01 6.244632e+01 7.789922e+01 5.927111e+01 5.964109e+01
#> [51] 6.204930e+01 6.249577e+01 8.192148e+01 5.396933e+01 1.005981e+02
#> [56] 5.964264e+01 5.772927e+01 9.057003e+01 8.004580e+01 5.833457e+01
#> [61] 5.952143e+01 6.132379e+01 1.147784e+02 6.384963e+01 1.245934e+02
#> [66] 5.863608e+01 8.472976e+01 6.819545e+01 6.064205e+01 6.327739e+01
#> [71] 5.543249e+01 7.025545e+01 7.769011e+01 7.818204e+01 1.021211e+02
#> [76] 6.177004e+01 5.835631e+01 5.513306e+01 6.673265e+01 5.915305e+01
#> [81] 6.165698e+01 6.844404e+01 6.485812e+01 1.289597e+02 5.963141e+01
#> [86] 6.507195e+01 5.845564e+01 9.510875e+01 6.593403e+01 6.061642e+01
#> [91] 8.805082e+01 7.549977e+01 6.291239e+01 6.395587e+01 6.177158e+01
#> [96] 6.203529e+01 5.748841e+01 6.579903e+01 7.139826e+01 1.650784e+02
#> [101] 6.736191e+01 1.334447e+02 8.606890e+01 1.068758e+02 6.635673e+01
#> [106] 7.656527e+01 5.732241e+01 8.524603e+01 8.064744e+01 1.109960e+02
#> [111] 7.595839e+01 5.640794e+01 5.610494e+01 6.308539e+01 6.302486e+01
#> [116] 5.969409e+01 7.712926e+01 5.397264e+01 6.854337e+01 8.993986e+01
#> [121] 1.283920e+02 6.796375e+01 1.178116e+02 5.866937e+01 6.629484e+01
#> [126] 5.282531e+01 1.070056e+02 5.420776e+01 6.551338e+01 6.915750e+01
#> [131] 1.212641e+02 1.295151e+02 4.988383e+01 5.612010e+01 6.983271e+01
#> [136] 6.797538e+01 7.407597e+01 3.291898e+02 5.828703e+01 6.339855e+01
#> [141] 5.916244e+01 6.151595e+01 1.427384e+02 2.133686e+02 6.105873e+01
#> [146] 7.897619e+01 6.527299e+02 5.934124e+01 5.725389e+01 6.008308e+01
#> [151] 8.714625e+01 5.629807e+01 7.566406e+01 5.796902e+01 6.007832e+01
#> [156] 5.731570e+02 2.043245e+03 8.254409e+01 6.992410e+01 5.549081e+01
#> [161] 5.772428e+01 5.776133e+01 5.933899e+01 7.945155e+01 5.862385e+01
#> [166] 5.715344e+01 1.265028e+02 6.830013e+01 6.931220e+01 6.576832e+01
#> [171] 6.323384e+01 5.706080e+01 6.706930e+01 7.328447e+01 5.443117e+01
#> [176] 2.152416e+02 5.472260e+01 6.333573e+01 5.767293e+01 6.659539e+01
#> [181] 6.429912e+01 5.661194e+01 5.856266e+01 7.693605e+01 6.386797e+01
#> [186] 5.696864e+01 5.707866e+01 9.524720e+02 5.989766e+01 6.197641e+01
#> [191] 5.780280e+01 5.932716e+01 8.165878e+01 6.646840e+01 6.198138e+01
#> [196] 6.897582e+01 6.850689e+01 5.440974e+01 6.939654e+01 5.884111e+01
#> [201] 1.258957e+03 8.308524e+01 6.501927e+01 5.653631e+01 6.199118e+01
#> [206] 5.936433e+01 6.474505e+01 5.834980e+01 5.870520e+01 5.818982e+01
#> [211] 6.232707e+01 6.364410e+01 6.271267e+01 7.091652e+01 1.135640e+02
#> [216] 6.529727e+01 6.744471e+01 6.148750e+01 6.624010e+01 7.025595e+01
#> [221] 7.371731e+01 6.824030e+01 5.743670e+01 8.708271e+01 6.725602e+01
#> [226] 5.672758e+01 6.553699e+01 6.777234e+01 5.315265e+01 6.560300e+01
#> [231] 1.412315e+02 1.016018e+02 5.581750e+01 6.712182e+01 6.746360e+01
#> [236] 6.345075e+01 5.993298e+01 6.206299e+01 5.737524e+01 5.899358e+01
#> [241] 6.155752e+01 6.154450e+01 1.034261e+02 1.205828e+02 1.003084e+02
#> [246] 2.259147e+02 5.934394e+01 7.659352e+01 7.125726e+01 5.824998e+01
#> [251] 6.533745e+01 6.128644e+01 6.794293e+01 7.033265e+01 5.772283e+01
#> [256] 7.817471e+01 1.280331e+02 5.855960e+01 6.343813e+01 6.334080e+01
#> [261] 5.975811e+01 8.190326e+01 9.061032e+01 9.186628e+01 5.731573e+01
#> [266] 7.927543e+01 6.087985e+01 6.257770e+01 6.826414e+01 6.906888e+01
#> [271] 1.502689e+02 6.651416e+01 5.798676e+01 6.207753e+01 1.425137e+02
#> [276] 1.788004e+02 6.394900e+01 6.979504e+01 5.719004e+01 5.626461e+01
#> [281] 1.702757e+03 7.554127e+01 6.508850e+01 6.462624e+01 6.904797e+01
#> [286] 6.784560e+01 7.569678e+01 8.885075e+01 7.668318e+01 5.994655e+01
#> [291] 7.140014e+01 7.951840e+01 6.525351e+01 5.651139e+01 5.996219e+01
#> [296] 7.582668e+01 5.680298e+01 9.418163e+01 4.963553e+01 8.264457e+01
#> [301] 8.922955e+01 7.268923e+01 5.780241e+01 8.368527e+01 6.813375e+01
#> [306] 6.804083e+01 5.667870e+01 6.884729e+01 6.141163e+01 6.156960e+01
#> [311] 5.927357e+01 7.535642e+01 1.293643e+02 6.003166e+01 7.863511e+01
#> [316] 5.565045e+01 5.508700e+01 7.498939e+01 5.445589e+01 6.236884e+01
#> [321] 7.482360e+01 7.849467e+01 1.160519e+02 8.019347e+01 6.124657e+01
#> [326] 2.127863e+02 5.289913e+02 4.763392e+01 6.772906e+01 2.102068e+02
#> [331] 7.173230e+01 6.371736e+01 5.163378e+01 5.647185e+01 8.413663e+01
#> [336] 6.700088e+01 5.168663e+01 9.420989e+01 1.340394e+02 5.911840e+01
#> [341] 6.693669e+01 5.085109e+01 6.282526e+01 5.544802e+01 2.086792e+02
#> [346] 6.340083e+01 5.682721e+01 6.069698e+01 7.341382e+01 5.661642e+01
#> [351] 6.179999e+01 6.024250e+01 5.926935e+01 2.428068e+02 7.247798e+02
#> [356] 6.695149e+01 8.704521e+01 7.250913e+01 7.517056e+01 6.517020e+01
#> [361] 6.127646e+01 6.376320e+01 8.192787e+01 8.442175e+01 6.543205e+01
#> [366] 5.684427e+01 6.995219e+01 6.397352e+01 6.657363e+01 8.489316e+01
#> [371] 5.799904e+01 6.176093e+01 7.064761e+01 7.150976e+01 6.355477e+01
#> [376] 6.829331e+01 5.515827e+01 8.495511e+01 6.745913e+01 6.957235e+01
#> [381] 5.569203e+01 9.927501e+01 5.852686e+01 7.402613e+01 6.460539e+01
#> [386] 5.180900e+01 2.327759e+03 1.598169e+03 6.093477e+01 6.502555e+01
#> [391] 7.916924e+02 5.225829e+01 5.911721e+01 6.594385e+01 6.345540e+01
#> [396] 5.839866e+01 5.849700e+01 5.918335e+01 5.824829e+01 6.319772e+01
#> [401] 5.818904e+01 5.452933e+01 6.078474e+01 7.425918e+01 5.715557e+01
#> [406] 8.320384e+01 7.704780e+01 8.186927e+01 6.050178e+01 6.988066e+01
#> [411] 6.206780e+01 1.331518e+02 7.459926e+01 8.389361e+01 6.128061e+01
#> [416] 7.922044e+01 6.824796e+01 6.773516e+01 6.064324e+01 5.791951e+01
#> [421] 4.943047e+03 3.495144e+03 6.869076e+01 5.989655e+01 6.078939e+01
#> [426] 5.528327e+01 6.890768e+01 8.662905e+01 5.964057e+01 6.550480e+01
#> [431] 5.769352e+01 7.443531e+01 5.934308e+01 6.292008e+01 5.627202e+01
#> [436] 5.971470e+01 1.381341e+02 4.942395e+01 5.327343e+01 1.765066e+03
#> [441] 1.093029e+02 7.236595e+01 6.186529e+01 7.180886e+01 7.008605e+01
#> [446] 7.241036e+01 5.475227e+01 6.639019e+01 6.535266e+01 5.795366e+01
#> [451] 7.803551e+01 5.994782e+01 8.468753e+01 7.755753e+01 6.371265e+01
#> [456] 5.660303e+01 6.742967e+01 7.543115e+01 5.844042e+01 1.638376e+02
#> [461] 1.573334e+02 8.268885e+01 5.644422e+01 5.606153e+01 1.355746e+02
#> [466] 5.485758e+01 7.686486e+01 5.914653e+01 8.470989e+01 7.431073e+01
#> [471] 6.685676e+01 9.404111e+01 7.018958e+01 6.568333e+01 1.105693e+02
#> [476] 7.873233e+01 5.716953e+01 5.918103e+01 6.465580e+01 9.308026e+01
#> [481] 1.747617e+02 1.010833e+02 1.627431e+02 3.141875e+02 5.625843e+01
#> [486] 6.233737e+01 6.109326e+01 7.493235e+01 6.846114e+01 6.174029e+01
#> [491] 6.330541e+01 6.746868e+01 6.158146e+01 8.642071e+01 7.058997e+01
#> [496] 5.917211e+01 1.845172e+02 7.557073e+01 5.955060e+01 6.582712e+01
#> [501] 5.985098e+01 1.644563e+02 6.071549e+01 7.024388e+01 1.289366e+02
#> [506] 5.831346e+01 8.072043e+01 1.100883e+02 1.103537e+02 8.547176e+01
#> [511] 6.536189e+01 6.929370e+01 1.292805e+02 6.281541e+01 6.662506e+01
#> [516] 6.584338e+01 5.651960e+01 4.904607e+01 6.221919e+01 5.825422e+01
#> [521] 1.226587e+02 1.328959e+02 5.931688e+01 6.779382e+01 6.255154e+01
#> [526] 5.641305e+01 9.165461e+01 5.803362e+01 6.318814e+01 8.057568e+01
#> [531] 6.464973e+01 5.264134e+01 6.056002e+01 1.517277e+02 6.044070e+01
#> [536] 6.183602e+01 5.962658e+01 5.717515e+01 6.360683e+01 6.189486e+01
#> [541] 6.128479e+01 5.509931e+01 5.931083e+01 6.482588e+01 7.407780e+01
#> [546] 6.080548e+01 6.432944e+01 1.052290e+05 1.424745e+02 5.865731e+01
#> [551] 6.306736e+01 5.667254e+01 6.014464e+01 6.232251e+01 7.503559e+01
#> [556] 8.659758e+02 4.917726e+01 5.375492e+01 6.160192e+01 6.407032e+01
#> [561] 6.934775e+01 7.155979e+01 7.377067e+01 7.017011e+01 6.164293e+01
#> [566] 5.474194e+01 5.612159e+01 8.328620e+01 5.979145e+01 6.723569e+01
#> [571] 5.646186e+01 1.442159e+02 1.568084e+02 9.896066e+01 6.361382e+01
#> [576] 6.101379e+01 8.954685e+01 5.465731e+01 6.870365e+01 5.494613e+01
#> [581] 5.828471e+01 6.020131e+01 6.508278e+01 5.600057e+01 5.890790e+01
#> [586] 6.659659e+01 6.033324e+01 7.054418e+01 9.236940e+01 5.896502e+01
#> [591] 6.988714e+01 5.938336e+01 6.246989e+01 6.285711e+01 5.767335e+01
#> [596] 8.071460e+01 5.949345e+01 6.172668e+01 5.929593e+01 5.344095e+01
#> [601] 8.179520e+01 6.585106e+01 6.298873e+01 8.366478e+01 7.031894e+01
#> [606] 5.962062e+01 7.799485e+01 8.399711e+01 6.471320e+01 4.886406e+01
#> [611] 9.659914e+01 5.643272e+01 6.053960e+01 8.753005e+01 5.209168e+01
#> [616] 6.310391e+01 7.928739e+01 7.021912e+01 5.795857e+01 6.113754e+01
#> [621] 7.345051e+01 6.174696e+01 8.504628e+01 8.393045e+01 7.094906e+01
#> [626] 5.725182e+01 1.126373e+02 5.961476e+01 6.357958e+01 5.664090e+01
#> [631] 6.391992e+01 5.873078e+01 6.181751e+01 5.751120e+01 2.664483e+02
#> [636] 6.911929e+01 7.163757e+01 2.962648e+02 3.535784e+02 7.215096e+01
#> [641] 6.020425e+01 6.342392e+01 9.200636e+01 3.044217e+02 8.543952e+01
#> [646] 1.053713e+02 6.321070e+01 8.386563e+01 5.644736e+08 8.121385e+01
#> [651] 8.380624e+01 5.777652e+01 5.564318e+01 7.080861e+01 8.252471e+01
#> [656] 6.846808e+01 5.381853e+01 6.160022e+01 6.852592e+01 1.826578e+02
#> [661] 5.503938e+01 1.284865e+02 5.120304e+01 5.118384e+01 6.486578e+01
#> [666] 1.096368e+02 5.553580e+01 6.063481e+01 9.796559e+01 5.537312e+01
#> [671] 8.866651e+01 1.046394e+02 6.325244e+01 8.457621e+01 1.667864e+02
#> [676] 5.341799e+01 6.688403e+01 6.136456e+01 7.481712e+01 6.191601e+01
#> [681] 7.788408e+01 6.594682e+01 6.091422e+01 9.320971e+01 5.938894e+01
#> [686] 2.207901e+02 6.182001e+01 1.225232e+02 5.823966e+01 1.640786e+02
#> [691] 5.422788e+01 2.485396e+02 5.426396e+01 4.501936e+02 5.280441e+01
#> [696] 1.493137e+02 1.055145e+02 6.255234e+01 7.604086e+01 2.262499e+02
#> [701] 7.928527e+01 6.236226e+01 1.045538e+02 7.393585e+01 5.581888e+01
#> [706] 6.353656e+01 5.322931e+01 5.819228e+01 6.376747e+01 1.764874e+02
#> [711] 6.253853e+01 6.319984e+01 5.726711e+01 5.798973e+01 5.748636e+01
#> [716] 7.028530e+01 1.040584e+02 5.794596e+01 1.308496e+02 7.356920e+01
#> [721] 6.971341e+01 6.302002e+01 1.216973e+02 7.856116e+01 7.259432e+01
#> [726] 6.341878e+01 7.205269e+01 5.890445e+01 5.212117e+01 6.568908e+01
#> [731] 5.092692e+01 6.493449e+01 7.408539e+01 1.067832e+02 5.447812e+02
#> [736] 6.676704e+01 6.512341e+01 6.166439e+01 6.136718e+01 6.073720e+01
#> [741] 6.260402e+01 7.952663e+01 8.924433e+01 5.444114e+01 1.195194e+02
#> [746] 8.453371e+02 6.072125e+01 7.247066e+01 1.610447e+02 6.785159e+01
#> [751] 5.524703e+01 8.863358e+01 6.569786e+01 2.095557e+02 5.665025e+01
#> [756] 1.720252e+02 9.695972e+01 6.377035e+01 1.255840e+02 1.303955e+02
#> [761] 3.159939e+02 8.432256e+01 1.055164e+02 1.239763e+02 7.110999e+01
#> [766] 6.087910e+01 5.857022e+01 7.268812e+01 8.392330e+01 7.626777e+01
#> [771] 1.002238e+02 6.111042e+01 5.585214e+01 5.588747e+01 7.010046e+02
#> [776] 4.514325e+01 4.440547e+01 6.227101e+01 6.158768e+01 6.770211e+01
#> [781] 5.913647e+01 7.428752e+01 2.410648e+02 6.551646e+01 6.619601e+01
#> [786] 6.746822e+01 6.864131e+01 5.706951e+01 1.466633e+02 8.931173e+02
#> [791] 2.700004e+03 1.219173e+02 7.426541e+01 7.828114e+01 5.433387e+01
#> [796] 7.299432e+01 5.689797e+01 1.266491e+02 6.114732e+01 2.255390e+02
#> [801] 7.870103e+01 5.502934e+01 6.654846e+01 1.243999e+03 5.763456e+01
#> [806] 6.720542e+01 6.108122e+01 6.371183e+01 7.387542e+01 6.082468e+01
#> [811] 5.799490e+01 6.788905e+01 5.762774e+01 6.124572e+01 6.634039e+01
#> [816] 5.914793e+01 1.385683e+02 6.561939e+01 1.083811e+02 7.689474e+01
#> [821] 6.520329e+01 5.994342e+01 6.176473e+01 8.807655e+01 8.844213e+01
#> [826] 5.919832e+01 4.956856e+01 1.163238e+02 5.875551e+01 5.683728e+01
#> [831] 8.399039e+01 5.597233e+01 6.211498e+01 5.805473e+01 5.964011e+01
#> [836] 7.657906e+01 5.503626e+01 3.212651e+02 1.198160e+02 1.171330e+02
#> [841] 6.047181e+01 6.573531e+01 6.299634e+01 6.121750e+01 7.922590e+01
#> [846] 6.169031e+01 6.000873e+01 6.581746e+01 5.734013e+01 6.038889e+01
#> [851] 1.760469e+02 5.092729e+01 6.803270e+01 5.212302e+01 1.023621e+02
#> [856] 8.576400e+01 1.172089e+02 5.622205e+01 7.435950e+01 9.275600e+01
#> [861] 6.452296e+01 1.838273e+03 6.650698e+01 6.242742e+01 8.182379e+01
#> [866] 5.596237e+01 5.943589e+01 6.418666e+01 5.859839e+01 6.336085e+01
#> [871] 5.908199e+01 6.898264e+01 8.223411e+01 9.933074e+01 6.256079e+01
#> [876] 6.752438e+01 7.114333e+02 1.137549e+02 5.552834e+01 5.722097e+01
#> [881] 5.699077e+01 5.982850e+01 6.090215e+01 6.849982e+01 6.361671e+01
#> [886] 5.991215e+01 6.235310e+01 5.855430e+01 6.767235e+01 7.324243e+01
#> [891] 8.905046e+01 6.122469e+01 5.956901e+01 6.068218e+01 6.567090e+01
#> [896] 1.197762e+02 8.205905e+01 6.714415e+01 5.683610e+01 4.296761e+01
#> [901] 6.781343e+01 6.019734e+01 6.439158e+01 6.558382e+01 5.719418e+01
#> [906] 6.239933e+01 1.276342e+02 5.887776e+01 7.250553e+01 6.785795e+01
#> [911] 5.734965e+01 6.619378e+01 5.817362e+01 5.920200e+01 6.002201e+01
#> [916] 6.671975e+01 6.111754e+01 6.394558e+01 7.902163e+01 6.098560e+01
#> [921] 6.324736e+01 6.213111e+01 6.190219e+01 6.676204e+01 9.556848e+02
#> [926] 6.679178e+01 5.682946e+01 6.319739e+01 5.601855e+01 6.476892e+01
#> [931] 6.625479e+01 5.607709e+01 7.867604e+01 7.391186e+01 8.355441e+01
#> [936] 7.935551e+01 6.714838e+01 7.018258e+01 6.015604e+01 5.712610e+01
#> [941] 6.247613e+01 6.435590e+01 1.040989e+02 5.842402e+01 6.110708e+01
#> [946] 5.437714e+01 2.129749e+02 9.864414e+01 7.821238e+01 6.802751e+01
#> [951] 6.374840e+01 7.115129e+01 6.162040e+01 6.148024e+01 1.105295e+02
#> [956] 5.994336e+01 5.846067e+01 8.161215e+01 6.111759e+01 5.950099e+01
#> [961] 6.284138e+01 7.375617e+01 6.347391e+01 5.979000e+01 5.053043e+01
#> [966] 7.988565e+01 4.973386e+01 5.179315e+01 6.111058e+01 5.734548e+01
#> [971] 1.406432e+02 1.082246e+02 6.380108e+01 7.834880e+01 6.053050e+01
#> [976] 3.581192e+02 2.420185e+02 8.342451e+02 3.540699e+02 5.474027e+01
#> [981] 1.163364e+02 7.666808e+01 7.042008e+01 6.767070e+01 5.504665e+01
#> [986] 5.644779e+01 9.448273e+02 5.979995e+01 1.407716e+02 1.806418e+02
#> [991] 6.325622e+01 1.054409e+02 6.368250e+01 6.687321e+01 9.545095e+01
#> [996] 5.492276e+01 6.177970e+01 6.178334e+01 5.727741e+01 7.621633e+01