![[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] 1.899253e+02 7.210316e+01 3.535938e+02 1.562597e+02 8.536496e+01
#> [6] 1.093289e+02 8.503319e+01 7.642143e+01 1.705831e+02 7.075211e+01
#> [11] 1.132755e+02 6.347500e+01 8.512269e+01 8.648628e+01 2.395530e+02
#> [16] 7.009757e+01 7.566844e+01 1.822715e+02 7.328008e+01 1.036313e+02
#> [21] 9.326659e+01 6.706644e+01 1.038256e+02 5.766053e+02 8.288872e+01
#> [26] 7.874573e+01 6.382801e+01 8.608218e+01 8.545663e+01 1.818461e+02
#> [31] 7.409687e+01 1.946884e+02 8.705562e+01 1.260731e+02 8.094242e+01
#> [36] 1.431826e+02 6.680950e+01 6.270038e+02 1.012503e+03 8.308751e+01
#> [41] 6.653308e+01 4.317225e+02 1.015065e+02 7.701153e+01 1.322273e+02
#> [46] 7.465916e+01 2.008151e+02 8.408163e+01 1.516111e+02 1.081894e+02
#> [51] 7.850842e+01 6.263104e+01 7.322825e+01 3.248012e+02 1.153021e+02
#> [56] 3.441049e+02 7.441062e+01 1.038674e+02 6.857497e+01 6.450786e+01
#> [61] 6.766659e+01 8.866065e+01 9.949463e+02 1.054098e+02 6.939841e+01
#> [66] 1.027427e+02 9.797306e+01 2.446267e+02 9.592142e+01 8.378427e+01
#> [71] 6.985571e+01 9.050053e+01 9.309773e+01 7.145811e+01 6.264178e+01
#> [76] 1.714859e+02 6.150944e+01 6.284055e+02 6.855418e+01 7.036518e+01
#> [81] 8.018313e+01 9.555173e+01 7.028509e+01 7.522621e+01 7.046674e+01
#> [86] 6.971762e+01 1.358024e+02 7.048523e+01 1.522838e+02 1.028656e+02
#> [91] 6.805747e+01 1.957116e+02 7.473484e+01 3.659161e+02 7.684033e+01
#> [96] 6.807183e+01 7.741919e+01 8.324544e+01 6.677841e+01 7.734601e+01
#> [101] 8.915841e+01 1.027114e+02 1.145277e+02 7.384320e+01 7.993452e+01
#> [106] 2.009096e+02 9.619184e+01 6.392751e+01 6.926916e+01 5.803814e+01
#> [111] 5.485875e+02 8.263789e+01 8.270602e+01 8.062282e+01 7.413633e+01
#> [116] 7.412647e+01 8.250964e+01 2.344579e+02 1.760883e+02 7.981896e+01
#> [121] 8.863541e+02 1.266018e+02 1.070238e+02 6.595397e+01 9.975668e+01
#> [126] 6.619406e+01 1.320744e+02 1.790120e+02 7.485876e+01 1.605664e+02
#> [131] 1.044409e+02 1.206140e+02 1.064786e+02 2.322767e+02 9.171316e+01
#> [136] 9.808180e+01 6.614407e+01 1.329810e+02 3.107622e+02 1.012657e+02
#> [141] 2.463954e+02 8.296627e+02 8.137269e+01 1.043190e+04 6.248081e+01
#> [146] 2.059841e+02 9.295631e+01 2.603766e+04 3.207198e+02 6.520943e+01
#> [151] 6.310889e+01 6.589242e+01 7.479556e+01 1.041094e+02 1.094253e+02
#> [156] 6.998186e+01 1.121679e+02 9.781946e+01 7.204423e+01 1.067015e+02
#> [161] 1.304876e+02 7.031037e+01 8.665130e+01 1.171584e+02 6.817059e+01
#> [166] 7.780835e+01 6.202361e+01 8.542107e+02 1.998956e+02 1.131121e+02
#> [171] 8.712225e+01 1.088735e+02 8.050726e+01 1.395706e+02 9.209887e+01
#> [176] 1.013009e+02 1.135640e+02 1.523077e+02 1.020677e+02 6.693271e+01
#> [181] 8.124513e+01 6.735178e+01 7.074675e+01 7.844078e+01 7.198800e+01
#> [186] 2.263116e+02 4.726018e+02 6.669855e+01 1.258192e+02 6.042581e+01
#> [191] 1.260307e+02 7.117489e+01 7.490956e+01 7.037894e+01 6.613427e+01
#> [196] 7.651111e+01 8.268786e+01 2.594355e+02 1.371482e+02 7.601833e+01
#> [201] 6.836281e+02 2.908699e+04 1.100684e+02 1.311768e+02 7.203276e+01
#> [206] 1.388727e+02 1.850207e+02 1.386313e+02 2.264633e+02 4.913417e+02
#> [211] 6.980657e+02 1.034512e+02 9.843338e+01 9.104433e+01 1.312307e+02
#> [216] 8.081328e+01 9.476826e+01 1.083663e+02 1.405107e+02 6.570748e+02
#> [221] 8.376674e+01 6.616114e+01 6.960264e+01 1.615025e+02 6.591238e+01
#> [226] 8.263396e+02 6.540856e+01 1.114584e+02 3.276422e+02 1.361879e+02
#> [231] 1.870375e+02 1.624742e+02 1.035974e+02 8.827328e+01 2.528129e+02
#> [236] 1.402930e+02 2.190356e+02 6.037922e+01 8.997943e+01 6.463610e+02
#> [241] 7.677983e+01 8.290103e+01 7.671520e+01 1.352775e+02 7.466288e+01
#> [246] 7.952236e+01 9.783116e+01 8.505560e+01 9.899747e+01 7.971645e+01
#> [251] 1.033824e+02 6.470677e+01 6.492563e+01 1.011926e+02 8.194359e+01
#> [256] 9.620633e+01 6.699771e+01 1.234555e+02 9.783992e+01 7.549969e+01
#> [261] 1.189548e+02 5.484790e+02 7.629050e+01 8.340208e+01 6.678892e+01
#> [266] 7.890130e+01 8.611209e+01 8.859279e+01 8.696694e+01 7.580667e+01
#> [271] 3.574069e+02 8.522384e+01 9.685377e+01 1.221457e+02 6.666616e+01
#> [276] 3.933437e+03 1.360886e+02 6.502051e+01 1.471830e+02 6.864889e+01
#> [281] 7.986447e+01 8.460277e+01 9.565554e+01 8.918841e+01 8.611870e+01
#> [286] 1.685151e+02 7.215940e+01 7.939348e+01 9.486758e+01 1.248336e+02
#> [291] 7.836792e+01 7.242830e+01 5.093625e+02 1.586887e+02 9.619293e+01
#> [296] 1.270757e+02 1.082284e+02 1.664281e+02 1.086518e+02 6.572478e+01
#> [301] 6.262203e+01 7.239095e+01 6.324825e+01 7.220589e+01 6.715313e+01
#> [306] 3.704896e+02 7.688321e+01 8.834417e+01 6.364845e+01 7.213206e+01
#> [311] 9.548360e+01 4.160982e+02 1.059225e+02 8.049105e+01 1.004211e+02
#> [316] 1.052279e+02 2.151915e+02 6.795727e+01 9.652330e+01 9.446693e+01
#> [321] 6.697544e+01 8.366874e+01 1.380294e+02 6.413371e+01 2.909363e+02
#> [326] 1.262064e+02 8.107558e+01 4.491564e+02 2.458672e+02 7.112515e+01
#> [331] 9.668809e+01 9.522261e+01 9.159064e+01 7.496928e+01 7.695279e+01
#> [336] 1.151138e+02 1.666028e+02 8.036996e+01 1.503690e+02 1.322241e+02
#> [341] 7.737383e+01 6.473642e+01 9.227476e+01 8.930138e+01 1.553775e+02
#> [346] 8.290936e+01 1.533572e+03 7.033723e+01 7.566502e+01 2.527476e+02
#> [351] 6.850629e+01 7.658924e+01 1.288616e+02 1.337214e+03 9.570968e+01
#> [356] 8.662736e+01 7.200773e+01 1.301948e+02 6.899182e+01 9.012244e+01
#> [361] 8.005913e+01 1.033318e+02 1.711618e+02 1.246647e+02 2.007221e+02
#> [366] 7.094974e+02 1.103332e+03 6.747210e+01 7.155033e+01 6.054751e+01
#> [371] 1.070992e+02 7.339593e+01 8.502121e+01 1.703693e+02 2.826856e+04
#> [376] 8.685945e+01 6.141937e+01 7.237666e+01 7.211551e+01 8.083836e+01
#> [381] 7.772180e+01 1.276339e+02 1.271241e+02 8.825897e+01 8.100721e+01
#> [386] 1.754117e+02 7.774547e+01 2.719784e+02 6.672995e+01 6.784379e+01
#> [391] 1.167656e+02 1.104176e+02 1.509272e+02 9.563480e+01 1.053499e+02
#> [396] 1.375935e+02 1.257711e+03 8.820953e+01 7.056613e+01 1.536017e+02
#> [401] 9.094013e+01 8.503292e+01 7.196735e+01 1.121101e+02 1.576330e+02
#> [406] 8.338428e+01 1.023494e+02 9.217789e+01 6.467819e+01 2.436713e+02
#> [411] 7.261710e+01 1.384152e+02 1.809677e+06 1.075622e+02 1.284460e+02
#> [416] 5.512033e+02 8.246361e+02 2.810532e+09 1.844409e+02 8.164320e+01
#> [421] 2.183826e+02 7.618722e+01 8.484895e+01 8.775193e+01 1.552692e+02
#> [426] 8.326317e+01 9.574467e+01 4.052963e+02 2.095775e+02 7.452703e+01
#> [431] 1.493625e+02 1.182754e+02 1.643854e+02 1.238164e+02 6.498251e+01
#> [436] 2.505754e+02 7.177916e+01 7.661930e+01 7.373007e+01 8.162911e+01
#> [441] 7.245966e+01 8.091598e+01 9.315712e+01 8.514951e+01 1.050173e+02
#> [446] 8.202230e+01 6.393042e+01 8.456569e+01 6.482849e+01 8.038929e+01
#> [451] 2.369267e+02 1.123485e+03 1.383257e+02 7.942424e+01 6.437854e+03
#> [456] 7.317149e+01 7.453411e+01 8.530818e+01 9.904689e+01 1.101535e+02
#> [461] 5.529052e+02 1.006942e+02 8.943000e+01 8.964127e+01 7.531381e+01
#> [466] 7.047637e+01 4.039286e+02 9.195183e+01 8.505524e+01 8.692735e+01
#> [471] 7.576442e+01 8.786990e+01 2.675560e+02 6.852511e+01 3.537584e+03
#> [476] 4.171369e+02 9.066856e+01 1.586409e+02 5.061987e+02 1.016028e+03
#> [481] 1.721506e+02 7.795029e+01 9.453472e+01 1.802546e+02 8.604497e+01
#> [486] 6.940092e+01 6.552385e+01 8.999454e+01 1.634949e+02 7.877126e+01
#> [491] 9.717815e+01 8.863472e+01 1.468885e+02 6.761337e+01 9.604090e+01
#> [496] 9.040411e+01 2.835185e+02 3.551161e+02 1.889999e+03 2.029883e+04
#> [501] 7.125440e+01 1.879722e+02 1.081571e+02 5.258949e+02 6.488295e+01
#> [506] 4.363777e+04 7.735405e+01 1.634552e+02 8.512614e+01 7.024665e+01
#> [511] 9.351012e+01 1.366337e+02 8.481524e+01 1.855466e+02 2.539782e+02
#> [516] 8.467963e+01 2.388859e+02 6.995088e+01 9.665669e+01 1.777020e+02
#> [521] 1.148712e+02 8.614416e+01 1.548816e+02 8.590737e+01 8.845907e+01
#> [526] 1.857610e+03 1.312386e+02 6.140989e+01 6.833094e+01 6.632971e+01
#> [531] 8.532495e+01 1.030074e+02 6.551596e+01 6.782157e+01 7.726828e+01
#> [536] 7.042751e+01 9.119885e+01 9.802024e+01 1.178666e+02 7.121327e+01
#> [541] 1.210249e+02 3.038331e+02 3.183463e+03 1.087681e+02 8.914459e+01
#> [546] 1.818551e+02 2.672663e+02 6.375644e+01 6.846447e+01 6.759599e+01
#> [551] 9.625713e+01 1.632439e+02 8.474054e+01 1.435261e+03 8.613178e+01
#> [556] 2.820842e+02 8.816483e+01 1.676126e+02 9.202197e+01 4.296671e+03
#> [561] 8.043790e+01 2.095634e+02 1.999632e+02 7.334693e+01 6.152341e+01
#> [566] 9.049155e+01 7.310625e+01 1.328012e+02 2.184808e+02 9.242487e+01
#> [571] 1.185482e+02 2.691099e+02 1.461094e+02 1.156319e+02 1.408914e+02
#> [576] 1.549403e+02 6.391534e+01 6.483654e+01 7.949958e+01 8.015917e+01
#> [581] 8.538828e+01 7.879883e+01 7.483421e+01 6.649867e+01 7.925305e+01
#> [586] 6.899562e+01 8.810203e+01 7.428122e+01 3.682111e+03 8.228187e+01
#> [591] 1.910382e+02 1.235673e+02 9.778071e+01 7.664352e+01 1.448531e+02
#> [596] 6.706621e+01 8.520936e+01 1.252789e+02 5.978190e+01 8.612887e+01
#> [601] 6.880003e+01 1.003163e+02 9.689072e+01 6.115222e+02 7.923943e+01
#> [606] 6.272863e+01 7.626943e+01 1.305677e+02 1.382240e+03 1.440268e+02
#> [611] 8.836703e+01 8.204740e+01 7.129581e+01 1.639076e+02 4.179735e+02
#> [616] 8.121882e+01 8.032346e+01 1.047780e+02 6.808084e+01 1.016195e+02
#> [621] 8.269291e+01 1.361976e+02 8.304840e+01 6.758053e+01 1.106337e+02
#> [626] 8.736018e+01 8.698653e+01 9.091273e+01 7.502767e+01 8.030143e+01
#> [631] 1.564888e+02 7.522294e+01 9.147997e+01 1.105888e+02 1.038078e+02
#> [636] 2.142664e+02 8.434765e+01 8.129562e+01 1.994080e+02 2.228047e+02
#> [641] 9.650814e+01 7.646683e+01 8.165241e+01 1.314234e+02 2.132922e+02
#> [646] 7.715325e+01 7.446385e+01 6.719170e+01 9.215095e+01 7.391271e+01
#> [651] 8.872259e+01 1.011552e+02 7.460419e+01 1.359275e+03 9.952033e+01
#> [656] 1.035854e+02 8.287560e+01 7.745809e+01 8.496017e+01 1.225570e+02
#> [661] 6.201457e+01 1.276392e+02 1.675472e+02 1.107240e+02 1.298349e+02
#> [666] 2.731140e+02 1.162032e+02 1.148022e+02 7.613719e+01 6.294735e+01
#> [671] 8.863275e+01 1.054856e+02 7.963929e+01 6.863776e+01 1.397305e+02
#> [676] 7.464834e+01 7.336716e+01 7.237895e+01 6.943968e+01 1.681960e+02
#> [681] 9.437560e+01 1.516479e+02 2.239832e+02 6.114958e+01 7.727890e+01
#> [686] 7.736109e+01 1.395083e+02 3.518145e+02 1.596659e+02 7.500372e+01
#> [691] 1.481925e+02 1.095654e+02 6.157744e+01 9.615174e+01 9.309486e+01
#> [696] 8.358438e+01 8.287619e+01 9.297559e+01 1.822796e+02 1.000575e+02
#> [701] 1.436024e+02 6.681088e+01 6.289630e+01 6.311801e+01 6.211738e+01
#> [706] 6.566158e+01 7.628378e+01 5.501438e+02 1.492863e+02 1.614702e+02
#> [711] 1.220170e+02 6.771976e+01 7.277755e+01 8.082156e+01 1.415625e+02
#> [716] 7.024896e+01 7.805185e+01 7.560905e+01 1.519449e+02 6.805282e+01
#> [721] 1.227221e+02 5.684177e+02 9.386365e+01 7.210791e+01 1.257768e+02
#> [726] 7.513157e+01 7.502927e+01 6.394285e+01 8.204694e+01 9.365892e+01
#> [731] 9.458210e+01 1.504103e+03 9.974000e+01 3.054155e+02 7.551109e+01
#> [736] 1.763059e+02 8.390206e+01 8.370082e+01 6.538839e+01 8.168563e+01
#> [741] 2.365048e+02 9.019243e+01 6.087565e+01 2.813296e+02 1.015841e+02
#> [746] 9.608011e+01 6.944062e+01 9.238129e+01 7.509809e+01 6.250212e+01
#> [751] 6.440981e+01 6.277476e+01 9.566519e+01 4.667032e+02 8.865700e+01
#> [756] 8.407849e+01 1.608751e+02 8.267786e+01 6.678206e+01 8.503097e+01
#> [761] 7.301790e+01 6.597388e+01 6.988089e+01 7.409505e+01 6.402489e+01
#> [766] 2.650187e+02 9.317179e+01 2.128907e+02 1.306620e+02 8.090369e+01
#> [771] 1.004584e+02 7.370794e+01 1.149097e+02 7.869787e+01 1.157949e+02
#> [776] 6.646760e+01 7.768482e+01 1.124262e+02 7.440082e+01 8.251954e+01
#> [781] 4.232117e+02 1.013318e+02 7.529818e+01 7.102529e+01 9.636006e+01
#> [786] 6.211210e+01 1.392214e+02 1.133618e+02 6.797096e+01 1.173762e+02
#> [791] 6.378547e+01 7.356101e+01 7.922991e+01 8.405587e+01 6.648895e+02
#> [796] 5.162163e+02 7.716743e+01 6.682437e+01 9.580892e+01 6.481411e+01
#> [801] 7.546374e+01 6.744530e+01 1.589866e+02 7.316685e+01 8.012531e+01
#> [806] 1.298808e+03 2.879390e+02 2.619286e+02 1.706151e+02 1.030047e+02
#> [811] 9.380941e+01 7.795707e+01 7.813825e+01 7.953513e+01 1.316127e+02
#> [816] 1.807081e+02 1.091858e+02 6.672561e+01 6.572650e+01 6.580774e+01
#> [821] 8.517721e+01 1.033448e+02 8.355914e+01 7.867995e+01 8.538970e+01
#> [826] 6.175259e+01 1.132430e+02 1.222808e+02 7.290201e+01 1.339516e+02
#> [831] 1.092147e+02 7.733599e+02 8.087841e+01 6.818079e+01 1.576469e+02
#> [836] 1.199823e+02 7.035733e+01 7.331103e+01 8.373456e+01 2.099336e+02
#> [841] 1.198203e+02 1.341559e+03 1.082638e+02 8.024230e+01 6.831707e+01
#> [846] 1.187498e+02 7.227425e+01 1.172565e+02 9.390128e+01 1.121229e+02
#> [851] 1.743838e+02 9.271570e+01 6.197891e+01 8.881067e+01 9.971996e+01
#> [856] 6.671026e+01 1.144589e+02 8.236123e+01 1.405471e+02 2.119886e+02
#> [861] 1.106702e+03 6.858562e+01 7.218689e+01 1.115030e+02 2.280712e+02
#> [866] 8.330407e+01 1.517192e+02 8.875730e+01 1.565014e+02 7.390734e+01
#> [871] 9.697715e+01 3.207983e+02 6.891264e+01 6.907791e+01 6.779972e+01
#> [876] 2.247422e+02 9.173925e+01 7.744109e+01 6.473107e+01 7.280026e+01
#> [881] 7.594861e+01 9.155790e+01 3.323651e+02 1.617464e+04 7.362698e+05
#> [886] 6.803863e+20 3.356166e+02 1.178247e+02 7.825418e+01 2.077981e+02
#> [891] 8.357763e+01 1.759966e+02 1.186700e+02 7.251192e+01 7.885181e+01
#> [896] 7.168923e+01 6.283884e+01 6.671753e+01 6.968260e+01 8.570932e+01
#> [901] 5.004657e+02 7.703419e+01 1.084568e+02 2.556262e+02 9.105021e+01
#> [906] 1.072442e+02 9.090012e+01 8.583380e+01 6.604920e+01 6.389859e+01
#> [911] 7.819470e+01 1.492089e+02 7.032579e+01 1.426525e+02 1.137306e+02
#> [916] 8.726924e+01 6.203327e+01 1.031301e+03 6.577740e+01 1.695286e+02
#> [921] 1.314985e+02 1.114262e+02 5.034173e+02 2.213630e+02 1.332245e+02
#> [926] 1.013463e+02 6.740751e+01 7.917646e+01 7.233160e+01 2.029396e+02
#> [931] 1.111420e+02 7.838018e+01 1.216629e+02 1.555135e+02 8.810530e+01
#> [936] 7.274102e+01 6.580811e+02 1.020580e+02 7.433056e+01 1.134956e+02
#> [941] 1.159848e+02 1.434319e+02 1.168787e+02 2.324638e+02 9.132884e+01
#> [946] 5.928998e+01 1.168093e+02 1.122004e+02 6.916570e+01 6.321845e+01
#> [951] 6.636059e+01 7.978052e+01 8.500105e+01 8.390123e+01 1.714716e+02
#> [956] 4.762262e+02 7.700163e+01 8.363786e+01 6.319106e+01 2.165281e+02
#> [961] 8.917237e+01 1.101734e+02 6.222937e+01 7.844804e+01 8.183313e+01
#> [966] 1.027565e+02 9.712427e+01 8.786106e+01 8.069461e+01 9.846549e+01
#> [971] 2.954797e+04 1.134608e+02 1.107632e+02 6.424490e+01 8.517483e+01
#> [976] 1.885548e+02 7.970198e+01 6.770662e+01 9.412767e+01 6.474098e+01
#> [981] 8.667842e+01 2.620892e+02 6.958593e+01 7.570569e+01 1.035408e+03
#> [986] 9.712314e+01 5.710322e+02 4.747393e+02 1.201580e+03 1.427649e+02
#> [991] 1.182414e+03 5.078037e+02 1.331779e+02 7.921950e+01 1.276403e+02
#> [996] 6.525266e+01 1.250348e+02 6.914919e+01 8.125314e+01 7.612723e+01
doseFunction(
model,
alpha2 = samples@data$alpha2,
beta = samples@data$beta,
grade = 2L
)(x = 0.25)
#> [1] 7.831900e+01 5.865592e+01 1.552751e+02 7.157685e+01 6.345620e+01
#> [6] 5.937880e+01 6.713107e+01 6.134320e+01 9.871627e+01 5.176306e+01
#> [11] 6.927787e+01 5.509622e+01 6.638689e+01 6.222174e+01 1.280995e+02
#> [16] 5.755882e+01 5.634983e+01 7.662849e+01 6.037596e+01 7.569717e+01
#> [21] 6.409864e+01 5.514268e+01 7.252506e+01 2.202942e+02 5.646370e+01
#> [26] 6.336712e+01 5.519964e+01 6.859134e+01 5.274279e+01 1.083874e+02
#> [31] 5.506974e+01 1.041202e+02 5.916204e+01 6.573278e+01 5.619506e+01
#> [36] 7.341366e+01 5.859499e+01 1.255572e+02 1.419050e+02 6.739821e+01
#> [41] 5.880897e+01 1.556186e+02 7.144384e+01 6.499937e+01 9.349573e+01
#> [46] 6.227815e+01 1.138243e+02 5.322616e+01 9.359818e+01 7.927023e+01
#> [51] 6.734285e+01 5.813800e+01 6.119568e+01 1.227949e+02 7.138708e+01
#> [56] 1.097332e+02 6.177177e+01 6.473160e+01 6.218165e+01 5.645824e+01
#> [61] 5.367612e+01 6.687581e+01 1.652898e+02 7.769052e+01 5.618346e+01
#> [66] 8.062847e+01 7.643418e+01 1.159730e+02 7.191594e+01 6.659562e+01
#> [71] 5.841046e+01 6.751283e+01 6.677951e+01 5.875742e+01 5.229812e+01
#> [76] 1.005197e+02 5.596770e+01 3.031520e+02 6.064718e+01 6.065095e+01
#> [81] 5.984701e+01 7.048821e+01 6.196008e+01 6.072411e+01 5.847427e+01
#> [86] 5.612132e+01 7.648043e+01 5.968563e+01 9.036562e+01 7.624348e+01
#> [91] 6.038949e+01 1.187035e+02 6.286981e+01 1.783468e+02 6.568036e+01
#> [96] 5.807481e+01 6.283103e+01 6.683526e+01 5.968011e+01 6.445841e+01
#> [101] 6.504334e+01 6.838532e+01 7.388710e+01 6.153804e+01 5.605293e+01
#> [106] 9.873789e+01 5.845326e+01 5.291357e+01 5.954185e+01 5.280498e+01
#> [111] 2.868489e+02 5.749630e+01 4.852317e+01 6.135996e+01 6.008311e+01
#> [116] 5.845172e+01 6.290431e+01 9.430754e+01 9.764581e+01 6.495282e+01
#> [121] 2.208024e+02 7.166937e+01 8.035768e+01 6.000057e+01 7.902105e+01
#> [126] 6.054610e+01 9.243016e+01 7.809759e+01 6.133544e+01 9.050694e+01
#> [131] 4.707903e+01 7.098870e+01 8.149207e+01 1.404108e+02 6.612515e+01
#> [136] 7.357817e+01 6.001122e+01 9.057768e+01 1.390903e+02 7.652879e+01
#> [141] 1.104515e+02 1.937285e+02 6.644237e+01 1.049204e+03 5.550148e+01
#> [146] 1.164507e+02 6.872521e+01 9.842470e+02 1.456423e+02 5.274095e+01
#> [151] 5.718089e+01 5.938776e+01 6.456672e+01 7.877986e+01 7.974095e+01
#> [156] 5.915147e+01 8.169556e+01 7.279351e+01 6.185057e+01 7.532151e+01
#> [161] 4.580906e+01 5.787707e+01 6.318705e+01 7.084041e+01 5.558806e+01
#> [166] 5.770383e+01 5.792564e+01 3.941551e+02 1.368142e+02 8.475186e+01
#> [171] 6.762780e+01 8.051755e+01 6.236160e+01 9.279922e+01 7.357498e+01
#> [176] 7.078601e+01 7.095933e+01 9.451047e+01 6.452765e+01 5.384564e+01
#> [181] 6.853416e+01 6.140317e+01 6.218990e+01 6.859669e+01 5.751037e+01
#> [186] 6.366841e+01 1.786750e+02 5.199067e+01 7.751395e+01 5.553814e+01
#> [191] 9.173939e+01 6.201570e+01 5.963737e+01 6.000083e+01 6.013302e+01
#> [196] 6.615598e+01 7.233644e+01 1.408574e+02 8.744859e+01 6.045355e+01
#> [201] 1.068883e+02 2.262682e+03 7.563365e+01 8.504752e+01 6.226900e+01
#> [206] 7.935777e+01 1.085136e+02 7.432316e+01 1.043466e+02 1.832331e+02
#> [211] 1.230929e+02 6.302015e+01 7.057012e+01 6.178420e+01 7.000323e+01
#> [216] 5.944700e+01 6.764894e+01 6.094209e+01 9.056423e+01 1.839477e+02
#> [221] 6.686862e+01 5.986315e+01 6.182263e+01 1.095878e+02 5.907035e+01
#> [226] 3.473549e+02 5.669591e+01 8.801496e+01 2.057700e+02 8.747842e+01
#> [231] 1.045696e+02 6.394844e+01 7.419870e+01 6.579212e+01 9.264394e+01
#> [236] 9.237260e+01 8.260256e+01 5.493386e+01 7.131508e+01 2.020859e+02
#> [241] 6.261180e+01 6.815701e+01 6.503340e+01 8.848247e+01 6.359126e+01
#> [246] 6.571871e+01 7.951470e+01 6.707317e+01 6.393432e+01 5.990971e+01
#> [251] 8.008289e+01 5.836176e+01 5.926505e+01 7.556074e+01 6.422018e+01
#> [256] 7.382475e+01 5.942739e+01 8.766786e+01 7.677792e+01 5.894648e+01
#> [261] 5.812103e+01 1.806221e+02 6.044216e+01 6.766564e+01 5.514500e+01
#> [266] 6.225204e+01 7.482054e+01 7.692739e+01 7.451071e+01 6.173730e+01
#> [271] 1.886032e+02 5.710390e+01 7.639401e+01 7.727530e+01 5.778703e+01
#> [276] 3.347352e+02 8.312517e+01 5.074670e+01 8.296124e+01 6.063093e+01
#> [281] 6.446617e+01 6.985799e+01 6.509895e+01 6.328284e+01 5.321257e+01
#> [286] 1.106697e+02 5.812265e+01 6.113895e+01 6.881804e+01 6.312707e+01
#> [291] 6.021429e+01 5.937683e+01 1.341131e+02 9.877698e+01 6.755528e+01
#> [296] 7.724438e+01 6.483238e+01 9.297859e+01 6.367708e+01 5.451385e+01
#> [301] 5.631891e+01 6.051496e+01 5.597141e+01 6.192406e+01 5.941807e+01
#> [306] 1.382331e+02 6.430939e+01 6.694282e+01 5.410960e+01 5.943264e+01
#> [311] 7.207413e+01 1.447711e+02 6.369373e+01 6.033753e+01 6.737709e+01
#> [316] 7.392013e+01 9.904185e+01 5.763244e+01 8.123474e+01 6.988819e+01
#> [321] 5.809506e+01 6.646797e+01 4.895904e+01 5.533329e+01 3.236373e+01
#> [326] 6.361145e+01 6.708160e+01 1.951742e+02 1.616702e+02 6.250749e+01
#> [331] 7.021002e+01 6.660266e+01 7.385041e+01 6.434043e+01 5.934060e+01
#> [336] 6.935954e+01 8.598840e+01 6.262878e+01 8.652620e+01 8.191427e+01
#> [341] 6.739715e+01 5.737689e+01 7.490065e+01 6.689668e+01 1.020497e+02
#> [346] 6.746071e+01 2.197235e+02 5.664061e+01 6.257654e+01 1.158979e+02
#> [351] 5.700378e+01 5.800477e+01 7.096279e+01 1.647754e+02 7.232976e+01
#> [356] 5.940382e+01 5.726033e+01 7.712021e+01 5.920841e+01 6.562412e+01
#> [361] 6.560992e+01 8.067600e+01 1.085686e+02 6.241488e+01 7.538343e+01
#> [366] 3.182856e+02 4.926461e+02 6.168316e+01 6.037592e+01 5.528308e+01
#> [371] 7.905879e+01 6.147129e+01 6.814184e+01 1.027210e+02 2.439242e+03
#> [376] 6.891589e+01 5.343596e+01 5.992907e+01 6.065784e+01 6.392954e+01
#> [381] 6.084716e+01 7.539522e+01 7.215416e+01 6.775256e+01 6.980700e+01
#> [386] 1.214011e+02 6.839589e+01 1.574745e+02 5.320333e+01 5.670127e+01
#> [391] 8.267948e+01 6.888810e+01 8.305421e+01 6.852832e+01 6.754652e+01
#> [396] 8.184239e+01 3.333511e+02 6.665175e+01 5.998050e+01 8.034173e+01
#> [401] 6.086035e+01 6.625516e+01 6.119099e+01 7.515097e+01 9.727839e+01
#> [406] 5.951409e+01 7.648947e+01 6.515256e+01 5.675247e+01 1.213940e+02
#> [411] 5.458961e+01 7.294458e+01 1.221346e+04 6.704032e+01 6.219967e+01
#> [416] 2.429549e+02 1.794710e+02 8.466219e+05 9.480459e+01 6.054241e+01
#> [421] 1.248760e+02 5.477303e+01 6.643859e+01 6.378219e+01 9.503168e+01
#> [426] 6.007755e+01 7.368771e+01 8.714737e+01 7.813565e+01 6.098532e+01
#> [431] 6.692645e+01 6.033366e+01 5.291238e+01 9.113045e+01 5.848525e+01
#> [436] 1.145793e+02 5.136807e+01 6.550814e+01 6.002486e+01 6.665901e+01
#> [441] 6.014810e+01 6.211825e+01 6.413395e+01 6.770071e+01 6.441616e+01
#> [446] 7.175465e+01 5.940537e+01 7.089459e+01 5.596049e+01 6.705316e+01
#> [451] 1.231534e+02 3.275035e+02 6.664500e+01 5.889393e+01 6.915429e+02
#> [456] 6.221309e+01 6.128434e+01 6.701625e+01 7.135626e+01 7.085812e+01
#> [461] 1.931720e+02 7.001929e+01 5.621316e+01 6.956911e+01 6.420736e+01
#> [466] 6.040287e+01 1.360664e+02 6.261396e+01 5.885647e+01 6.576325e+01
#> [471] 6.402559e+01 6.345243e+01 1.118763e+02 5.715054e+01 8.553637e+02
#> [476] 1.597710e+02 5.478078e+01 8.769782e+01 1.417776e+02 1.860130e+02
#> [481] 9.540714e+01 6.220280e+01 5.512322e+01 9.102426e+01 6.856559e+01
#> [486] 6.114127e+01 5.857681e+01 7.121355e+01 9.333087e+01 6.919171e+01
#> [491] 7.109128e+01 7.243050e+01 7.111888e+01 5.857813e+01 6.927813e+01
#> [496] 6.195116e+01 7.205176e+01 1.099393e+02 1.847200e+02 5.543074e+02
#> [501] 5.211829e+01 8.551336e+01 6.638786e+01 1.584026e+02 5.009680e+01
#> [506] 8.178561e+02 5.837881e+01 9.768655e+01 6.876763e+01 5.773307e+01
#> [511] 6.783802e+01 8.248936e+01 6.737223e+01 1.085587e+02 8.873312e+01
#> [516] 6.655308e+01 1.008224e+02 5.680513e+01 7.363584e+01 1.193294e+02
#> [521] 8.207858e+01 6.525341e+01 9.319225e+01 6.321181e+01 6.415573e+01
#> [526] 5.876341e+02 7.097297e+01 5.545617e+01 5.965793e+01 5.833119e+01
#> [531] 7.119132e+01 6.707229e+01 5.773262e+01 6.104723e+01 6.437631e+01
#> [536] 5.948380e+01 6.849521e+01 7.500769e+01 8.661335e+01 5.740213e+01
#> [541] 8.177520e+01 1.464200e+02 3.830990e+02 8.177737e+01 5.625063e+01
#> [546] 1.023700e+02 1.354817e+02 5.841632e+01 6.150280e+01 6.104123e+01
#> [551] 7.935209e+01 1.045184e+02 6.378202e+01 3.267711e+02 6.867325e+01
#> [556] 1.359396e+02 5.734437e+01 1.013069e+02 4.491919e+01 4.434019e+02
#> [561] 6.010798e+01 1.134621e+02 8.711944e+01 6.298908e+01 5.255958e+01
#> [566] 6.781789e+01 5.354299e+01 6.782958e+01 5.672046e+01 5.789563e+01
#> [571] 8.199495e+01 1.445012e+02 1.015804e+02 7.835286e+01 8.809927e+01
#> [576] 9.245381e+01 5.830264e+01 5.932615e+01 6.620172e+01 6.108679e+01
#> [581] 6.629079e+01 5.939870e+01 5.932370e+01 5.785520e+01 5.867163e+01
#> [586] 5.769089e+01 6.788043e+01 5.959771e+01 5.411625e+02 5.222999e+01
#> [591] 9.537235e+01 7.659231e+01 6.579708e+01 5.922673e+01 7.141751e+01
#> [596] 5.645547e+01 6.390354e+01 6.197488e+01 5.269519e+01 7.171784e+01
#> [601] 5.747040e+01 7.593154e+01 6.099183e+01 2.290362e+02 6.646706e+01
#> [606] 5.719446e+01 5.625118e+01 7.254263e+01 2.995996e+02 8.664610e+01
#> [611] 6.895982e+01 6.460317e+01 5.591399e+01 9.227407e+01 1.306315e+02
#> [616] 6.287472e+01 6.493652e+01 7.104580e+01 6.035777e+01 8.014787e+01
#> [621] 5.772461e+01 7.717844e+01 6.638119e+01 5.830381e+01 8.358719e+01
#> [626] 6.379743e+01 7.045518e+01 7.343681e+01 5.933890e+01 6.693282e+01
#> [631] 4.895603e+01 6.169934e+01 6.439642e+01 7.331299e+01 7.577208e+01
#> [636] 8.899347e+01 5.989040e+01 6.548539e+01 1.251189e+02 9.133802e+01
#> [641] 6.933985e+01 6.459164e+01 6.103635e+01 9.465176e+01 1.315487e+02
#> [646] 5.884353e+01 6.185216e+01 5.655666e+01 7.402008e+01 6.386817e+01
#> [651] 7.261688e+01 5.625880e+01 6.237833e+01 5.388798e+02 4.840283e+01
#> [656] 7.003132e+01 6.015733e+01 6.463704e+01 6.883643e+01 8.371958e+01
#> [661] 5.617539e+01 8.654151e+01 7.905904e+01 7.445693e+01 6.678785e+01
#> [666] 1.348082e+02 7.391320e+01 8.382180e+01 6.085871e+01 4.663464e+01
#> [671] 6.267714e+01 6.834241e+01 6.263985e+01 5.827049e+01 7.864597e+01
#> [676] 6.109338e+01 5.615023e+01 5.863906e+01 5.718293e+01 6.591104e+01
#> [681] 6.021107e+01 8.160730e+01 1.197012e+02 5.409526e+01 6.367112e+01
#> [686] 6.588870e+01 8.493013e+01 1.024298e+02 8.992103e+01 6.201064e+01
#> [691] 9.187072e+01 7.819965e+01 5.550306e+01 7.183707e+01 6.846780e+01
#> [696] 5.913032e+01 6.675439e+01 7.402648e+01 1.007253e+02 7.058076e+01
#> [701] 9.259576e+01 5.781646e+01 5.701679e+01 5.803001e+01 5.721186e+01
#> [706] 5.851958e+01 5.570680e+01 1.177287e+02 8.927496e+01 8.489905e+01
#> [711] 7.284837e+01 6.056061e+01 5.743546e+01 5.803160e+01 8.828458e+01
#> [716] 6.016949e+01 5.941858e+01 6.672438e+01 9.028332e+01 5.050338e+01
#> [721] 7.713406e+01 2.043575e+02 7.628641e+01 6.374450e+01 9.043596e+01
#> [726] 6.453725e+01 6.471715e+01 5.878327e+01 6.945584e+01 7.284012e+01
#> [731] 6.959255e+01 3.265521e+02 7.385747e+01 1.523070e+02 6.448290e+01
#> [736] 8.938205e+01 6.181426e+01 6.701001e+01 5.656662e+01 6.125332e+01
#> [741] 9.686853e+01 6.303020e+01 5.201478e+01 7.486419e+01 6.804314e+01
#> [746] 6.718856e+01 6.195549e+01 7.643678e+01 6.368560e+01 5.593054e+01
#> [751] 5.894656e+01 5.743796e+01 7.699636e+01 1.426979e+02 6.659547e+01
#> [756] 6.278549e+01 9.952636e+01 6.350466e+01 5.731001e+01 6.132388e+01
#> [761] 5.622464e+01 5.625356e+01 6.026998e+01 5.855305e+01 5.532580e+01
#> [766] 1.209140e+02 7.536759e+01 1.091051e+02 8.120488e+01 6.101366e+01
#> [771] 7.076127e+01 5.768227e+01 7.729665e+01 5.756105e+01 7.743922e+01
#> [776] 5.827035e+01 6.532879e+01 8.510291e+01 5.973309e+01 5.835702e+01
#> [781] 1.374171e+02 6.559635e+01 5.397599e+01 5.510372e+01 7.498054e+01
#> [786] 5.484133e+01 9.044920e+01 8.056163e+01 5.294652e+01 8.026208e+01
#> [791] 5.489261e+01 6.105271e+01 6.471414e+01 5.960390e+01 1.681525e+02
#> [796] 2.085491e+02 6.452310e+01 5.683128e+01 7.415947e+01 5.592335e+01
#> [801] 6.283755e+01 5.607663e+01 6.576639e+01 5.659347e+01 6.522086e+01
#> [806] 2.915691e+02 1.076861e+02 1.395891e+02 9.312506e+01 7.056201e+01
#> [811] 5.977253e+01 5.959229e+01 5.773889e+01 6.475067e+01 8.321562e+01
#> [816] 6.492458e+01 8.276992e+01 6.160046e+01 5.983180e+01 5.629270e+01
#> [821] 6.227582e+01 6.950073e+01 5.389942e+01 6.049533e+01 6.722479e+01
#> [826] 5.576221e+01 6.673735e+01 7.283943e+01 5.760647e+01 7.612133e+01
#> [831] 5.795867e+01 1.741873e+02 5.546117e+01 5.172248e+01 8.034342e+01
#> [836] 7.296674e+01 5.606640e+01 6.105437e+01 6.423052e+01 9.970390e+01
#> [841] 7.910719e+01 2.460574e+02 7.732616e+01 6.521858e+01 5.502611e+01
#> [846] 6.412102e+01 6.197149e+01 7.233473e+01 7.389763e+01 7.311360e+01
#> [851] 8.991205e+01 5.957593e+01 5.345191e+01 6.531910e+01 6.543903e+01
#> [856] 4.982564e+01 7.075114e+01 6.354013e+01 8.569161e+01 1.197456e+02
#> [861] 3.143219e+02 6.060091e+01 6.396064e+01 8.066181e+01 8.513369e+01
#> [866] 6.034873e+01 8.452935e+01 6.625616e+01 9.464697e+01 6.429347e+01
#> [871] 6.620508e+01 1.465878e+02 5.736049e+01 5.990605e+01 5.823056e+01
#> [876] 8.606757e+01 7.161657e+01 6.223888e+01 5.580827e+01 6.263051e+01
#> [881] 6.504261e+01 7.193015e+01 1.579513e+02 1.013175e+03 1.506206e+02
#> [886] 6.063119e+14 1.822830e+02 9.083462e+01 6.054005e+01 1.006943e+02
#> [891] 6.347413e+01 4.880804e+01 7.242059e+01 5.823623e+01 6.189242e+01
#> [896] 6.103087e+01 5.623617e+01 6.002697e+01 6.144719e+01 7.121686e+01
#> [901] 1.656432e+02 5.558824e+01 7.542789e+01 1.018858e+02 6.825827e+01
#> [906] 7.269481e+01 6.179435e+01 6.844194e+01 5.896422e+01 5.598541e+01
#> [911] 6.603896e+01 1.050593e+02 6.282059e+01 8.511672e+01 8.453201e+01
#> [916] 6.804382e+01 5.450759e+01 3.215600e+02 5.710444e+01 1.012942e+02
#> [921] 8.626747e+01 7.830776e+01 2.245904e+02 1.282532e+02 6.386259e+01
#> [926] 6.151580e+01 6.093575e+01 6.638916e+01 6.110509e+01 1.000326e+02
#> [931] 7.041866e+01 6.209841e+01 7.288467e+01 7.320309e+01 6.603051e+01
#> [936] 6.311715e+01 1.467931e+02 7.173646e+01 6.037093e+01 7.803892e+01
#> [941] 8.035903e+01 8.365866e+01 4.856339e+01 1.027338e+02 6.039040e+01
#> [946] 5.262606e+01 6.828729e+01 8.216888e+01 6.070784e+01 5.291920e+01
#> [951] 5.662717e+01 6.046026e+01 6.249911e+01 5.589841e+01 5.310470e+01
#> [956] 1.967428e+02 6.145341e+01 6.616336e+01 5.728231e+01 1.145579e+02
#> [961] 6.217823e+01 7.398618e+01 5.682351e+01 6.756629e+01 6.530105e+01
#> [966] 7.305251e+01 7.860132e+01 7.036778e+01 6.594132e+01 7.522406e+01
#> [971] 1.644040e+03 8.465161e+01 5.830207e+01 5.730967e+01 7.019841e+01
#> [976] 8.643980e+01 6.238895e+01 5.538896e+01 6.373050e+01 5.423406e+01
#> [981] 7.126871e+01 1.483264e+02 6.134345e+01 6.431989e+01 2.677469e+02
#> [986] 6.229860e+01 2.468789e+02 1.485933e+02 1.755721e+02 6.721798e+01
#> [991] 1.928857e+02 1.902343e+02 9.155413e+01 6.466272e+01 8.282751e+01
#> [996] 5.340433e+01 5.977147e+01 5.636390e+01 6.639420e+01 6.015183e+01