![[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] 2.917765e+02 9.444544e+03 8.671283e+01 7.108075e+01 1.041848e+02
#> [6] 8.271564e+01 1.349324e+02 7.318094e+01 7.829202e+01 1.001352e+02
#> [11] 6.641409e+01 6.897200e+01 7.039531e+01 1.041486e+02 2.481004e+02
#> [16] 2.209378e+02 7.590687e+01 6.703321e+01 2.056634e+03 1.590168e+02
#> [21] 6.379258e+01 9.957523e+01 2.508425e+02 7.846792e+01 6.019107e+02
#> [26] 1.551138e+02 1.968434e+02 2.432794e+02 9.811087e+04 7.736859e+01
#> [31] 2.626222e+02 1.154657e+02 9.819949e+01 6.249913e+01 6.585674e+01
#> [36] 6.492914e+01 1.970404e+03 1.348859e+02 3.496461e+02 8.987898e+01
#> [41] 6.781684e+01 6.657597e+01 1.153914e+02 7.637692e+01 8.885663e+01
#> [46] 7.269483e+01 6.651691e+01 8.508750e+01 9.764885e+01 9.862063e+01
#> [51] 1.869572e+02 7.814174e+01 7.954565e+01 7.066664e+01 8.028842e+01
#> [56] 7.497925e+01 1.116378e+02 7.286222e+01 6.954606e+01 7.406940e+01
#> [61] 1.686935e+02 1.791445e+02 7.773008e+01 7.471311e+01 7.103840e+01
#> [66] 3.074887e+02 7.183569e+01 7.575691e+01 1.008079e+02 2.558874e+04
#> [71] 7.498364e+01 7.047163e+01 9.216627e+01 8.451054e+01 2.498029e+02
#> [76] 7.448082e+01 7.028867e+02 1.853641e+02 9.029790e+01 6.804519e+01
#> [81] 7.225798e+01 6.546167e+01 6.887695e+01 6.667628e+01 6.817120e+01
#> [86] 1.103333e+02 8.172278e+01 6.397941e+01 1.066293e+02 7.961405e+01
#> [91] 1.190516e+02 1.331400e+02 1.136592e+05 1.898229e+02 9.169259e+01
#> [96] 8.457352e+01 6.049968e+01 4.312231e+02 4.088909e+04 1.933160e+02
#> [101] 1.638948e+02 7.106660e+01 1.756389e+02 6.764786e+01 8.306161e+01
#> [106] 1.161360e+02 1.885227e+02 6.757121e+01 2.348614e+02 1.166774e+02
#> [111] 1.190772e+02 9.295884e+01 5.992691e+01 8.349067e+02 9.533393e+01
#> [116] 7.455932e+01 6.912797e+01 9.333011e+01 8.003669e+01 6.376328e+01
#> [121] 4.389687e+02 1.876844e+03 1.042915e+02 8.408941e+02 2.672965e+02
#> [126] 6.280620e+01 9.499486e+01 6.282137e+01 7.098744e+01 8.735660e+01
#> [131] 8.179231e+01 6.976141e+01 1.402728e+10 3.575148e+04 6.416431e+01
#> [136] 8.118309e+01 7.054548e+01 8.348429e+01 7.655030e+01 6.917874e+01
#> [141] 8.509334e+01 8.114612e+01 7.106371e+01 8.869393e+01 6.985952e+01
#> [146] 1.030456e+02 1.100270e+02 9.574633e+01 2.106184e+02 6.375682e+01
#> [151] 1.107017e+02 8.050368e+01 1.914540e+02 9.008195e+01 2.867748e+02
#> [156] 7.956545e+01 7.384807e+01 8.598250e+01 7.561312e+01 6.914415e+01
#> [161] 8.971416e+01 6.486546e+01 8.830819e+01 7.566134e+01 8.728978e+01
#> [166] 1.476104e+02 1.830576e+02 5.809691e+02 1.612673e+02 2.968640e+03
#> [171] 9.423382e+02 2.105838e+04 2.953045e+03 7.432047e+03 9.523514e+03
#> [176] 4.248740e+02 6.206186e+01 6.994511e+01 6.438374e+01 7.920362e+01
#> [181] 7.676306e+01 1.731206e+02 1.782831e+02 2.649493e+02 6.918462e+01
#> [186] 8.299826e+01 6.325208e+01 1.359608e+02 1.194991e+03 9.147208e+01
#> [191] 2.043746e+02 8.711590e+01 1.488842e+03 1.069141e+02 8.935333e+01
#> [196] 1.404197e+02 2.156856e+02 7.808674e+01 1.415592e+02 9.173631e+01
#> [201] 1.452889e+02 1.419822e+02 1.252187e+02 1.884930e+03 1.532146e+02
#> [206] 6.588292e+01 8.130858e+01 8.390487e+01 1.040367e+02 1.749558e+02
#> [211] 2.600387e+02 6.794557e+01 8.994303e+01 2.132900e+02 1.319046e+02
#> [216] 6.453464e+01 7.186046e+01 6.514321e+01 7.242822e+01 6.459934e+01
#> [221] 1.456020e+02 1.897857e+02 4.083143e+02 2.178303e+02 1.728617e+02
#> [226] 1.533945e+02 1.248212e+02 6.861162e+01 7.456868e+01 1.830700e+03
#> [231] 6.120884e+01 7.391286e+01 2.227624e+02 3.692719e+02 8.333271e+01
#> [236] 7.663060e+01 7.310681e+01 8.195156e+01 7.197603e+01 2.499030e+02
#> [241] 8.648790e+02 8.006460e+01 6.757044e+01 7.049862e+01 1.340928e+02
#> [246] 7.390694e+01 2.366405e+02 3.320562e+02 1.217843e+02 1.063365e+02
#> [251] 7.110696e+01 1.176625e+02 1.200058e+02 8.191916e+01 2.110143e+02
#> [256] 1.118096e+02 7.907651e+01 6.445145e+01 7.091793e+01 6.614041e+01
#> [261] 7.159704e+01 6.718227e+01 1.139307e+02 6.336351e+01 8.056613e+01
#> [266] 1.346357e+02 3.833028e+02 9.152557e+01 1.620516e+02 1.132221e+02
#> [271] 7.540468e+02 8.095662e+01 6.762017e+01 1.885992e+02 3.744328e+02
#> [276] 7.654241e+01 1.116488e+02 7.559896e+01 3.563872e+02 6.424921e+01
#> [281] 3.544116e+02 6.212624e+01 1.123542e+02 1.061958e+02 6.663443e+01
#> [286] 8.139648e+01 7.314091e+01 8.037548e+01 1.005209e+02 9.696620e+01
#> [291] 1.125700e+02 9.094851e+01 9.263082e+01 2.053450e+02 6.503179e+01
#> [296] 7.486540e+01 9.712723e+01 3.380448e+02 1.334911e+02 9.928067e+01
#> [301] 8.157810e+01 7.752559e+01 7.858096e+01 7.127028e+01 1.506408e+02
#> [306] 8.059039e+01 1.522490e+02 1.097786e+02 7.158250e+01 8.892759e+01
#> [311] 9.484163e+01 8.522883e+01 7.713023e+01 1.167391e+02 6.967892e+01
#> [316] 8.606208e+01 8.603614e+01 1.426118e+02 7.333442e+01 9.884995e+01
#> [321] 1.225898e+02 8.536877e+01 7.740009e+01 7.486546e+01 1.352430e+02
#> [326] 7.902331e+01 9.347030e+01 7.598345e+01 3.714488e+02 8.632816e+01
#> [331] 8.043497e+01 1.378413e+02 1.070794e+02 1.515351e+02 9.285804e+01
#> [336] 9.306132e+01 6.981776e+01 9.355621e+02 7.848942e+01 7.389167e+01
#> [341] 6.630251e+01 9.411940e+01 6.543051e+01 6.517880e+01 6.811693e+01
#> [346] 7.887649e+01 6.842298e+01 7.140054e+01 1.583679e+02 8.152427e+01
#> [351] 1.322117e+02 1.501026e+02 1.876673e+02 1.424800e+02 1.145219e+02
#> [356] 8.835394e+01 3.312562e+02 7.802423e+03 6.245059e+01 1.126478e+02
#> [361] 6.752824e+01 6.811586e+01 9.308554e+01 8.125876e+01 1.181240e+02
#> [366] 8.876836e+01 1.209516e+02 9.393210e+01 9.025305e+02 8.182790e+01
#> [371] 1.223011e+02 6.546359e+01 7.328601e+01 6.671257e+01 7.267714e+01
#> [376] 7.276781e+01 7.940043e+01 7.544827e+01 3.553762e+02 7.821231e+01
#> [381] 2.838671e+02 7.000126e+01 5.974208e+01 1.760748e+02 7.266847e+01
#> [386] 4.536822e+02 9.586913e+01 6.275525e+01 6.586689e+01 1.398454e+02
#> [391] 2.501324e+02 1.843635e+02 9.460668e+01 6.570171e+01 7.474748e+01
#> [396] 8.255032e+01 2.230026e+02 6.442289e+01 6.654240e+01 6.730087e+01
#> [401] 1.526962e+02 6.770618e+01 1.870169e+02 2.316913e+02 6.818002e+01
#> [406] 6.684078e+01 6.353626e+01 8.243509e+01 7.450844e+01 1.196165e+02
#> [411] 8.791380e+01 6.270758e+01 6.938203e+01 2.715751e+02 8.982832e+01
#> [416] 1.460150e+02 2.445167e+02 7.096853e+01 2.960844e+02 7.796312e+01
#> [421] 8.325004e+01 5.489530e+02 8.012380e+01 7.492089e+01 6.389115e+01
#> [426] 6.870360e+01 5.836590e+02 7.967514e+01 6.910953e+01 2.381639e+02
#> [431] 8.867740e+01 1.057624e+02 7.019969e+01 7.915736e+01 1.345648e+02
#> [436] 8.988211e+01 7.013619e+01 1.114530e+02 1.370365e+02 6.693027e+01
#> [441] 7.501970e+01 7.255186e+01 1.272761e+02 7.929117e+01 8.873235e+02
#> [446] 7.428527e+01 1.371382e+02 8.082856e+01 1.925943e+03 7.525251e+01
#> [451] 9.878602e+01 1.297858e+02 7.679457e+01 6.502105e+01 6.268173e+01
#> [456] 7.387726e+01 6.861112e+01 7.096317e+01 1.137696e+02 8.546219e+01
#> [461] 7.926590e+01 2.003799e+02 7.576784e+01 6.797412e+01 1.448656e+02
#> [466] 1.598236e+02 7.457949e+01 1.329072e+02 1.031416e+02 3.737908e+02
#> [471] 7.237051e+01 7.989141e+01 8.858445e+01 1.105851e+02 1.562498e+02
#> [476] 8.702432e+02 4.957281e+02 2.271460e+02 1.055023e+02 1.070838e+02
#> [481] 6.337690e+01 7.088720e+01 1.432447e+02 1.501400e+02 7.211511e+01
#> [486] 7.515486e+01 1.921069e+02 9.458458e+01 2.208853e+02 8.444424e+01
#> [491] 7.262623e+01 1.056253e+02 1.145915e+02 6.630562e+01 1.162634e+02
#> [496] 4.066801e+06 1.641563e+05 8.241364e+05 6.266757e+01 6.868315e+01
#> [501] 1.502675e+02 1.038174e+02 4.382257e+02 7.339561e+02 1.313365e+02
#> [506] 6.748040e+01 1.729126e+02 2.864253e+02 5.984949e+01 1.150359e+02
#> [511] 6.405288e+01 6.461201e+01 7.018128e+01 6.797978e+01 2.064748e+02
#> [516] 9.034519e+01 7.158077e+01 1.655340e+04 8.627539e+01 7.740316e+02
#> [521] 8.810418e+02 1.067639e+02 1.085719e+02 6.455592e+01 7.362882e+01
#> [526] 9.036313e+01 1.307377e+02 1.702592e+02 7.553378e+01 6.432538e+01
#> [531] 7.364606e+01 1.424840e+02 4.793796e+02 8.090648e+01 1.540343e+02
#> [536] 7.796958e+01 4.534228e+02 7.765206e+01 1.272759e+02 9.834886e+01
#> [541] 6.584193e+01 8.790056e+01 1.075545e+02 7.221992e+01 7.311144e+01
#> [546] 5.510064e+02 6.924542e+01 1.289516e+02 1.779277e+02 1.363609e+02
#> [551] 8.312681e+01 7.624027e+01 7.167795e+01 6.667586e+01 6.600591e+01
#> [556] 7.099056e+01 6.528484e+01 6.413882e+01 9.381218e+01 6.836367e+01
#> [561] 7.180353e+01 1.719750e+02 6.243526e+01 8.659269e+01 8.727129e+01
#> [566] 6.681802e+01 7.408438e+01 6.126557e+01 6.789774e+01 9.953367e+01
#> [571] 9.308643e+01 8.792228e+01 9.780943e+01 8.061300e+01 8.924827e+01
#> [576] 6.358572e+01 8.645715e+01 6.653804e+01 7.136457e+01 1.348064e+02
#> [581] 1.197095e+02 1.625505e+03 1.688931e+02 8.863496e+01 7.377157e+01
#> [586] 8.411570e+01 1.358861e+02 1.353401e+02 1.061036e+02 2.451105e+02
#> [591] 4.924406e+02 8.584085e+01 6.757582e+01 7.863650e+01 9.034047e+02
#> [596] 1.025415e+02 8.786238e+01 9.424815e+01 9.163022e+01 9.509363e+01
#> [601] 9.008646e+01 6.564744e+01 7.057663e+01 7.386059e+01 1.072360e+02
#> [606] 7.422720e+03 5.937798e+03 8.987442e+01 8.379614e+01 1.226324e+02
#> [611] 7.118540e+01 1.110011e+02 1.069084e+02 7.548745e+01 8.578746e+01
#> [616] 7.136020e+01 7.287715e+01 1.156087e+02 2.791167e+02 2.753434e+02
#> [621] 1.830155e+02 9.090926e+01 1.646330e+02 7.191617e+01 7.872867e+01
#> [626] 3.199232e+02 7.197696e+01 1.939481e+02 2.711542e+02 6.629130e+01
#> [631] 1.527603e+02 5.867039e+02 9.876551e+01 7.258430e+01 1.435969e+02
#> [636] 4.660390e+02 2.821174e+02 9.568665e+01 1.802324e+02 8.222141e+01
#> [641] 6.753504e+01 6.696080e+01 1.149174e+02 7.755839e+01 8.588447e+01
#> [646] 1.059365e+02 1.208195e+02 7.102296e+01 6.834881e+01 6.474900e+01
#> [651] 1.866170e+02 1.296419e+02 6.301736e+01 8.780450e+01 1.204047e+02
#> [656] 9.305011e+01 1.695733e+02 9.506382e+01 2.005105e+02 7.118228e+01
#> [661] 8.728821e+01 1.458919e+02 1.106879e+02 1.274970e+02 2.330951e+02
#> [666] 9.849359e+01 9.295738e+02 2.090848e+02 1.855869e+02 9.383628e+01
#> [671] 1.188556e+02 1.028912e+02 1.449265e+02 7.529179e+02 7.412151e+01
#> [676] 6.811995e+01 8.818433e+01 6.562436e+01 6.510201e+01 1.273700e+02
#> [681] 9.143620e+01 6.331965e+01 1.816817e+03 9.174120e+01 7.825270e+01
#> [686] 8.957301e+01 8.126039e+01 7.270474e+01 9.706197e+01 9.089610e+01
#> [691] 8.746595e+01 1.104934e+02 9.329116e+01 2.309302e+02 8.414951e+02
#> [696] 4.125310e+02 7.000226e+03 2.112786e+02 6.345856e+01 6.777989e+01
#> [701] 7.340756e+01 1.612722e+02 1.424427e+02 8.142002e+01 2.281252e+02
#> [706] 1.449485e+02 4.376016e+02 4.383704e+03 1.583275e+03 1.306833e+02
#> [711] 1.074497e+02 1.064863e+02 7.180542e+01 1.874879e+02 4.531136e+02
#> [716] 1.052223e+03 1.114223e+02 1.319944e+02 1.906361e+02 7.540315e+01
#> [721] 9.745343e+01 3.217469e+05 1.147019e+03 3.261551e+03 6.757699e+01
#> [726] 9.451107e+01 1.164987e+02 8.311499e+01 2.034025e+02 1.503138e+02
#> [731] 8.932448e+01 1.545159e+02 6.132811e+01 6.784187e+01 1.226919e+02
#> [736] 1.081765e+02 8.590549e+01 8.581258e+01 1.668217e+02 7.363651e+01
#> [741] 7.649707e+01 2.702631e+02 1.506076e+02 6.807260e+01 6.027522e+01
#> [746] 7.718155e+01 7.111250e+01 7.225548e+01 7.839939e+01 6.523650e+01
#> [751] 1.906288e+02 1.325515e+02 2.094040e+02 8.823884e+01 4.135366e+02
#> [756] 1.094600e+02 2.396339e+02 6.590650e+01 7.505847e+01 4.019928e+02
#> [761] 1.165850e+02 1.277406e+02 9.420180e+01 7.042815e+01 6.874286e+01
#> [766] 6.178106e+01 1.353343e+02 4.875443e+02 9.662715e+01 8.421242e+01
#> [771] 1.549465e+02 1.993981e+03 1.699041e+02 1.369244e+02 6.613800e+01
#> [776] 9.639300e+01 9.675607e+01 7.191820e+01 9.649458e+01 7.211039e+01
#> [781] 7.155064e+01 6.876304e+01 3.638491e+02 7.500368e+01 1.039946e+02
#> [786] 6.840760e+01 8.352964e+01 2.597260e+02 7.678509e+01 7.332868e+01
#> [791] 1.482852e+02 1.673036e+02 8.120527e+02 2.227283e+02 1.129734e+03
#> [796] 8.110639e+01 8.937201e+01 1.465955e+02 8.864453e+01 1.232436e+02
#> [801] 1.301129e+08 7.760387e+01 8.001705e+01 5.808595e+02 9.423823e+01
#> [806] 7.807330e+01 7.936730e+01 8.474904e+01 1.009171e+02 8.216811e+03
#> [811] 3.824493e+03 1.652213e+02 1.205115e+02 3.723470e+03 7.436007e+01
#> [816] 1.100290e+02 8.679681e+01 3.171516e+02 6.831270e+01 1.161341e+02
#> [821] 7.606444e+03 6.577335e+01 8.297289e+01 1.942084e+03 9.352682e+01
#> [826] 7.976409e+01 6.634164e+01 2.504993e+02 2.088523e+02 1.000460e+02
#> [831] 9.263330e+01 1.574335e+02 1.038310e+02 1.053930e+02 1.232447e+02
#> [836] 7.081470e+01 1.706456e+02 6.868263e+01 8.709119e+01 6.912311e+01
#> [841] 8.080471e+01 8.411188e+01 2.103752e+02 7.917405e+01 1.593078e+02
#> [846] 9.964892e+01 6.426415e+01 7.079388e+01 2.704320e+02 1.278904e+02
#> [851] 6.110165e+01 6.494979e+01 7.243088e+01 9.627211e+01 7.327329e+01
#> [856] 6.382925e+01 7.065451e+01 9.343972e+01 6.254460e+01 6.196599e+01
#> [861] 7.576029e+01 8.094233e+01 7.770856e+01 9.209015e+01 7.767701e+01
#> [866] 1.105413e+02 7.098648e+01 9.277607e+02 8.993819e+01 1.544419e+03
#> [871] 8.224256e+01 7.096748e+01 3.146626e+02 8.604951e+01 1.270485e+02
#> [876] 6.352877e+02 1.287671e+03 8.744555e+01 1.658865e+02 3.716796e+03
#> [881] 1.296195e+02 1.175251e+02 1.063636e+02 1.234904e+02 8.530293e+01
#> [886] 9.332946e+01 1.208231e+02 6.941333e+01 6.956836e+01 7.947530e+01
#> [891] 7.558084e+01 7.044039e+01 9.128352e+01 8.921993e+01 1.790258e+02
#> [896] 1.033982e+02 8.486570e+01 7.334896e+01 1.035236e+02 1.148908e+02
#> [901] 9.035731e+01 7.409019e+01 7.695359e+01 3.526243e+02 1.006031e+04
#> [906] 6.884600e+02 9.415892e+01 6.780317e+01 1.795530e+02 2.074905e+02
#> [911] 7.586867e+01 1.248120e+02 9.558834e+01 2.557169e+02 1.167042e+02
#> [916] 1.064415e+02 6.764447e+01 1.126247e+02 9.072593e+01 7.241234e+01
#> [921] 1.632971e+02 9.690673e+01 1.043952e+02 9.641869e+01 1.263069e+02
#> [926] 8.790795e+01 1.142030e+02 8.785220e+01 7.362105e+02 1.386019e+02
#> [931] 8.040558e+01 2.461278e+02 1.819853e+02 6.665334e+01 8.321186e+01
#> [936] 1.929866e+02 9.000374e+01 2.816988e+03 3.130925e+02 7.226535e+01
#> [941] 1.077045e+02 8.795521e+01 1.197024e+02 1.088096e+02 7.512077e+01
#> [946] 7.952182e+01 6.991666e+01 6.694740e+01 1.561643e+02 1.782365e+02
#> [951] 9.570582e+01 6.875543e+01 7.126490e+01 8.172743e+01 9.914287e+01
#> [956] 6.849486e+01 9.333854e+01 7.236769e+01 6.259164e+01 1.594144e+02
#> [961] 9.008780e+01 1.126477e+02 8.239109e+01 7.031025e+02 7.795893e+01
#> [966] 8.469540e+01 1.519972e+02 6.732173e+01 7.805280e+01 7.396677e+01
#> [971] 2.693857e+02 1.735810e+02 1.225497e+02 6.703287e+01 6.526968e+01
#> [976] 4.339262e+02 8.166890e+01 9.781861e+01 1.465677e+02 1.487856e+02
#> [981] 6.656202e+01 7.368101e+01 1.077856e+02 6.656175e+01 1.019136e+02
#> [986] 7.511591e+01 8.511799e+01 1.461366e+02 1.376424e+02 9.064070e+01
#> [991] 1.169522e+02 8.319187e+01 7.974323e+01 1.811345e+02 1.482072e+02
#> [996] 1.088542e+02 1.694896e+02 9.621954e+01 6.905992e+01 8.395692e+01
doseFunction(
model,
alpha2 = samples@data$alpha2,
beta = samples@data$beta,
grade = 2L
)(x = 0.25)
#> [1] 8.643145e+01 1.050257e+03 6.303225e+01 5.624191e+01 8.140205e+01
#> [6] 6.537750e+01 7.112732e+01 6.443470e+01 6.854758e+01 7.815736e+01
#> [11] 5.908495e+01 6.228435e+01 6.334586e+01 7.567062e+01 1.351571e+02
#> [16] 6.355463e+01 6.498640e+01 5.951363e+01 4.382489e+02 8.527291e+01
#> [21] 5.656611e+01 7.833092e+01 1.091046e+02 6.450434e+01 1.464226e+02
#> [26] 9.445557e+01 1.083286e+02 6.936623e+01 1.062915e+03 6.912279e+01
#> [31] 1.575920e+02 6.482507e+01 6.031971e+01 5.607948e+01 6.006636e+01
#> [36] 5.867809e+01 4.110822e+02 6.433694e+01 1.163021e+02 4.710156e+01
#> [41] 4.830371e+01 5.853151e+01 8.495050e+01 5.574121e+01 6.815409e+01
#> [46] 6.200064e+01 5.937173e+01 6.467409e+01 6.599140e+01 7.063264e+01
#> [51] 7.460381e+01 5.980058e+01 6.430431e+01 5.631401e+01 6.654091e+01
#> [56] 6.344564e+01 8.223774e+01 5.765310e+01 6.282823e+01 6.487811e+01
#> [61] 1.196953e+02 1.220861e+02 6.175770e+01 6.103434e+01 5.231894e+01
#> [66] 1.133920e+02 5.757493e+01 6.162047e+01 6.504366e+01 9.149936e+02
#> [71] 6.379844e+01 6.086457e+01 6.859434e+01 5.152503e+01 5.723472e+01
#> [76] 5.891301e+01 2.680023e+02 1.080044e+02 6.963100e+01 5.869711e+01
#> [81] 6.115403e+01 6.013216e+01 6.288730e+01 6.038912e+01 6.117192e+01
#> [86] 8.135992e+01 6.737701e+01 5.449231e+01 7.207697e+01 6.228264e+01
#> [91] 6.765892e+01 7.241901e+01 8.708449e+02 8.165361e+01 6.784613e+01
#> [96] 6.432842e+01 5.561774e+01 1.879023e+02 2.507205e+03 7.608276e+01
#> [101] 8.839548e+01 5.806639e+01 9.440498e+01 5.795322e+01 7.004481e+01
#> [106] 7.893069e+01 6.543684e+01 5.785710e+01 1.207256e+02 6.130989e+01
#> [111] 6.745086e+01 6.206607e+01 5.094181e+01 2.110074e+02 5.969576e+01
#> [116] 5.496442e+01 5.220612e+01 7.264203e+01 6.485860e+01 5.538515e+01
#> [121] 1.266693e+02 5.725618e+02 7.313483e+01 9.980494e+01 8.186647e+01
#> [126] 5.400315e+01 5.674237e+01 5.573570e+01 6.423821e+01 7.302765e+01
#> [131] 6.909091e+01 6.181998e+01 1.690459e+07 2.719287e+03 5.409909e+01
#> [136] 6.105749e+01 5.895041e+01 6.499455e+01 6.185551e+01 5.532888e+01
#> [141] 5.692777e+01 6.409822e+01 6.168407e+01 6.871509e+01 6.094448e+01
#> [146] 7.090475e+01 6.829679e+01 7.007191e+01 1.254627e+02 5.817534e+01
#> [151] 8.476843e+01 6.469773e+01 9.544456e+01 6.235813e+01 1.366631e+02
#> [156] 6.233295e+01 5.965968e+01 6.789919e+01 6.604362e+01 6.069069e+01
#> [161] 7.015220e+01 5.761206e+01 7.600409e+01 6.483052e+01 6.333917e+01
#> [166] 7.471304e+01 9.087139e+01 1.742840e+02 6.943509e+01 3.902207e+02
#> [171] 2.460503e+02 1.596266e+03 3.159401e+02 5.175190e+02 7.028182e+02
#> [176] 1.482160e+02 5.446745e+01 6.066958e+01 6.026090e+01 7.039281e+01
#> [181] 6.455715e+01 1.007985e+02 8.473676e+01 6.992018e+01 5.678444e+01
#> [186] 5.789738e+01 5.315579e+01 9.292654e+01 3.336090e+02 7.135495e+01
#> [191] 8.577073e+01 6.576171e+01 1.993496e+02 5.974241e+01 5.868281e+01
#> [196] 8.421734e+01 8.323527e+01 6.325579e+01 3.323777e+01 4.965108e+01
#> [201] 8.108700e+01 8.147556e+01 8.887056e+01 5.113399e+02 9.084368e+01
#> [206] 5.766517e+01 6.570429e+01 6.385642e+01 5.175581e+01 7.780446e+01
#> [211] 1.257290e+02 6.082124e+01 7.317280e+01 1.007726e+02 7.314586e+01
#> [216] 5.610556e+01 6.249813e+01 5.840815e+01 6.541904e+01 5.846366e+01
#> [221] 1.034966e+02 1.266246e+02 1.910668e+02 1.071480e+02 7.645685e+01
#> [226] 8.771320e+01 6.142249e+01 5.786243e+01 6.129631e+01 1.727978e+02
#> [231] 5.401629e+01 6.322436e+01 1.243420e+02 1.276810e+02 6.566063e+01
#> [236] 6.247903e+01 6.216937e+01 6.619571e+01 6.192688e+01 7.608988e+01
#> [241] 1.259392e+02 6.449159e+01 5.569393e+01 6.125509e+01 7.350661e+01
#> [246] 6.374378e+01 1.225426e+02 1.403140e+02 7.293975e+01 7.600431e+01
#> [251] 5.856973e+01 6.916219e+01 7.646095e+01 5.985392e+01 5.720755e+01
#> [256] 6.079628e+01 5.107224e+01 5.765383e+01 6.310343e+01 6.063038e+01
#> [261] 6.201905e+01 5.241432e+01 8.079162e+01 5.889922e+01 5.981696e+01
#> [266] 7.698694e+01 9.903450e+01 6.048926e+01 8.247634e+01 7.346194e+01
#> [271] 1.375843e+02 6.105540e+01 5.990692e+01 7.238090e+01 9.282521e+01
#> [276] 6.301836e+01 8.474273e+01 5.229423e+01 1.419188e+02 4.986032e+01
#> [281] 7.771681e+01 5.555321e+01 7.205413e+01 7.075083e+01 5.878145e+01
#> [286] 6.402202e+01 5.986173e+01 6.074555e+01 7.183546e+01 7.021898e+01
#> [291] 6.207372e+01 6.367366e+01 6.880987e+01 8.150012e+01 5.774120e+01
#> [296] 6.251094e+01 6.969798e+01 1.680847e+02 8.951116e+01 6.574690e+01
#> [301] 6.607412e+01 5.867112e+01 5.495263e+01 6.178309e+01 8.454485e+01
#> [306] 6.195276e+01 1.119881e+02 7.839606e+01 4.980299e+01 6.578532e+01
#> [311] 6.984929e+01 6.796656e+01 6.560771e+01 8.869859e+01 6.098200e+01
#> [316] 6.422098e+01 6.035687e+01 7.279030e+01 5.432531e+01 7.725374e+01
#> [321] 7.101646e+01 6.639673e+01 6.373736e+01 6.324182e+01 7.341284e+01
#> [326] 5.956313e+01 7.365921e+01 5.654020e+01 1.497836e+02 7.442541e+01
#> [331] 7.204825e+01 9.505651e+01 7.999479e+01 5.718501e+01 5.217655e+01
#> [336] 6.675252e+01 5.271188e+01 1.407300e+02 6.443851e+01 6.262848e+01
#> [341] 5.817188e+01 7.596355e+01 5.888804e+01 6.092122e+01 6.158807e+01
#> [346] 6.512982e+01 5.817666e+01 6.247841e+01 9.789317e+01 4.904510e+01
#> [351] 9.016587e+01 5.664060e+01 7.913952e+01 7.668773e+01 6.281610e+01
#> [356] 7.084905e+01 1.472743e+02 7.283288e+02 5.534244e+01 8.274527e+01
#> [361] 5.879462e+01 6.053837e+01 7.635375e+01 6.932647e+01 7.870580e+01
#> [366] 6.349909e+01 8.439730e+01 7.311118e+01 3.325286e+02 6.259565e+01
#> [371] 8.187459e+01 5.554363e+01 5.786031e+01 5.593397e+01 6.356140e+01
#> [376] 6.085741e+01 6.509486e+01 6.535535e+01 8.517671e+01 5.120146e+01
#> [381] 1.597538e+02 5.284982e+01 5.234383e+01 9.838282e+01 5.768654e+01
#> [386] 1.754321e+02 7.287958e+01 5.669669e+01 5.840460e+01 7.697617e+01
#> [391] 1.145149e+02 9.642871e+01 6.781788e+01 5.887738e+01 6.106598e+01
#> [396] 6.298117e+01 8.528520e+01 5.178512e+01 5.609558e+01 5.854454e+01
#> [401] 9.621878e+01 5.564778e+01 9.473364e+01 1.064802e+02 6.063804e+01
#> [406] 5.928106e+01 5.608303e+01 6.952989e+01 6.339690e+01 8.307489e+01
#> [411] 7.126211e+01 5.508984e+01 6.175804e+01 1.167253e+02 6.867671e+01
#> [416] 8.313295e+01 1.311543e+02 5.413945e+01 1.314089e+02 6.124112e+01
#> [421] 6.289803e+01 7.691037e+01 6.440882e+01 6.165290e+01 5.659148e+01
#> [426] 6.232223e+01 1.259930e+02 6.360783e+01 6.093489e+01 9.733082e+01
#> [431] 5.894980e+01 6.119302e+01 5.713560e+01 5.925735e+01 8.154352e+01
#> [436] 7.102220e+01 5.703555e+01 7.270208e+01 7.266728e+01 5.665386e+01
#> [441] 6.140032e+01 6.323074e+01 8.452100e+01 6.261679e+01 1.099444e+02
#> [446] 6.275535e+01 7.629845e+01 5.699669e+01 3.371531e+02 6.115601e+01
#> [451] 7.716194e+01 9.099260e+01 6.372531e+01 5.952707e+01 5.773943e+01
#> [456] 6.203466e+01 5.685612e+01 6.048681e+01 8.002995e+01 6.653803e+01
#> [461] 6.314330e+01 1.206970e+02 6.314344e+01 5.861474e+01 8.884404e+01
#> [466] 9.810329e+01 6.245044e+01 6.564406e+01 6.159329e+01 1.364509e+02
#> [471] 5.902018e+01 6.484866e+01 7.388891e+01 8.269524e+01 8.580998e+01
#> [476] 1.772245e+02 2.045610e+02 1.253237e+02 7.383830e+01 6.267233e+01
#> [481] 5.804465e+01 6.231395e+01 7.865590e+01 8.348111e+01 5.937593e+01
#> [486] 5.976588e+01 9.295868e+01 5.261376e+01 1.089749e+02 6.164427e+01
#> [491] 6.033073e+01 6.176533e+01 7.351948e+01 5.817809e+01 7.062114e+01
#> [496] 4.708920e+03 2.435656e+03 4.338470e+03 5.463819e+01 5.902527e+01
#> [501] 9.683891e+01 6.602539e+01 8.059699e+01 2.975422e+02 8.095237e+01
#> [506] 5.964207e+01 1.005447e+02 1.118896e+02 5.399444e+01 7.713120e+01
#> [511] 5.550477e+01 5.827513e+01 5.955294e+01 6.052046e+01 1.299198e+02
#> [516] 6.903320e+01 5.799012e+01 1.269047e+03 5.923950e+01 1.575404e+02
#> [521] 1.619978e+02 6.635012e+01 7.296511e+01 5.403407e+01 5.375995e+01
#> [526] 7.243881e+01 7.547860e+01 8.418113e+01 6.184411e+01 5.462739e+01
#> [531] 6.231514e+01 8.192621e+01 1.868041e+02 6.129422e+01 1.016114e+02
#> [536] 5.935900e+01 1.589520e+02 5.780053e+01 8.387367e+01 6.960870e+01
#> [541] 5.618643e+01 6.417307e+01 7.488201e+01 6.175044e+01 6.011079e+01
#> [546] 1.692604e+02 5.445712e+01 8.781053e+01 5.768356e+01 8.824468e+01
#> [551] 6.703498e+01 6.421478e+01 6.130341e+01 5.836267e+01 5.977030e+01
#> [556] 6.229896e+01 5.644176e+01 5.654110e+01 7.524576e+01 6.255631e+01
#> [561] 6.333558e+01 8.911167e+01 5.595241e+01 6.787631e+01 5.941252e+01
#> [566] 5.751644e+01 6.239238e+01 5.594640e+01 5.039725e+01 5.124247e+01
#> [571] 6.450678e+01 6.399089e+01 5.816856e+01 6.373943e+01 6.159892e+01
#> [576] 5.721130e+01 7.301232e+01 5.993514e+01 5.532422e+01 7.744916e+01
#> [581] 5.838446e+01 3.651448e+02 9.260900e+01 6.421003e+01 6.148094e+01
#> [586] 6.016757e+01 7.741676e+01 8.661519e+01 8.377861e+01 1.443456e+02
#> [591] 2.122531e+02 6.570647e+01 5.996281e+01 6.606533e+01 9.649342e+01
#> [596] 7.701022e+01 6.885436e+01 6.716069e+01 6.655622e+01 6.320634e+01
#> [601] 6.488043e+01 5.529271e+01 6.213437e+01 5.511878e+01 6.744328e+01
#> [606] 1.029429e+03 9.005444e+02 7.153745e+01 5.942001e+01 8.575059e+01
#> [611] 5.993787e+01 7.571826e+01 7.938162e+01 6.389608e+01 7.283192e+01
#> [616] 6.305679e+01 6.194509e+01 8.654514e+01 1.414472e+02 1.348827e+02
#> [621] 5.546417e+01 6.405844e+01 1.007040e+02 5.414035e+01 6.943993e+01
#> [626] 1.518741e+02 5.870408e+01 1.118281e+02 1.411662e+02 5.890756e+01
#> [631] 7.183311e+01 1.083497e+02 7.904390e+01 6.135679e+01 7.987979e+01
#> [636] 1.475292e+02 1.698476e+02 5.723659e+01 7.498971e+01 6.771326e+01
#> [641] 5.720116e+01 5.240728e+01 7.779629e+01 6.419084e+01 6.824676e+01
#> [646] 6.685341e+01 6.205942e+01 5.887955e+01 5.883422e+01 5.509636e+01
#> [651] 1.149412e+02 6.514149e+01 5.570501e+01 6.734630e+01 7.593593e+01
#> [656] 6.572335e+01 7.419754e+01 7.086590e+01 7.881827e+01 5.933371e+01
#> [661] 6.458974e+01 9.124058e+01 7.303763e+01 7.600102e+01 9.963326e+01
#> [666] 5.213335e+01 2.093172e+02 9.823743e+01 6.343762e+01 6.188745e+01
#> [671] 7.386001e+01 7.926089e+01 8.978195e+01 2.519913e+02 6.045429e+01
#> [676] 5.241129e+01 7.074070e+01 5.874297e+01 5.715657e+01 6.673358e+01
#> [681] 5.988184e+01 5.673779e+01 3.516517e+02 7.098027e+01 6.486500e+01
#> [686] 7.010416e+01 6.696687e+01 6.058101e+01 7.082456e+01 5.154058e+01
#> [691] 7.137317e+01 6.747160e+01 7.163157e+01 1.114741e+02 2.491100e+02
#> [696] 1.871033e+02 8.533834e+02 8.483365e+01 5.508612e+01 5.998753e+01
#> [701] 6.171468e+01 8.084566e+01 6.057047e+01 6.050623e+01 6.287623e+01
#> [706] 7.318504e+01 1.128533e+02 3.988855e+02 4.076955e+02 8.698514e+01
#> [711] 5.675808e+01 7.734486e+01 6.163047e+01 8.032956e+01 2.353836e+02
#> [716] 2.844122e+02 7.117655e+01 7.290079e+01 9.810545e+01 6.146801e+01
#> [721] 6.911532e+01 1.351722e+04 3.197020e+01 5.365740e+02 5.912495e+01
#> [726] 7.456201e+01 7.189409e+01 5.948094e+01 1.066114e+02 8.974717e+01
#> [731] 5.736076e+01 9.528595e+01 5.347142e+01 6.122159e+01 8.253157e+01
#> [736] 7.287448e+01 6.480032e+01 7.273006e+01 1.105214e+02 6.287921e+01
#> [741] 6.620174e+01 1.618760e+02 9.434369e+01 5.480508e+01 5.585835e+01
#> [746] 6.647435e+01 5.959117e+01 6.136419e+01 6.574458e+01 5.652303e+01
#> [751] 1.101740e+02 8.706966e+01 1.341319e+02 6.960212e+01 1.747289e+02
#> [756] 5.122092e+01 9.567962e+01 5.929205e+01 6.163324e+01 1.444231e+02
#> [761] 8.002790e+01 8.798377e+01 7.069637e+01 6.070850e+01 6.057419e+01
#> [766] 5.573056e+01 7.912962e+01 1.345687e+02 6.594247e+01 6.123875e+01
#> [771] 1.064676e+02 2.277883e+02 9.059058e+01 8.387108e+01 5.758521e+01
#> [776] 7.305131e+01 5.931595e+01 6.136206e+01 6.736029e+01 6.286977e+01
#> [781] 6.214898e+01 5.849415e+01 1.425659e+02 5.901372e+01 7.881480e+01
#> [786] 4.760285e+01 5.979075e+01 4.962310e+01 5.277465e+01 5.792336e+01
#> [791] 7.119383e+01 4.596367e+01 1.814386e+02 9.338631e+01 3.168678e+02
#> [796] 6.164715e+01 6.496211e+01 9.480358e+01 6.037548e+01 7.631149e+01
#> [801] 4.901900e+04 6.160066e+01 6.032471e+01 2.099653e+02 6.822679e+01
#> [806] 6.637006e+01 5.795008e+01 6.204520e+01 6.283589e+01 1.450195e+03
#> [811] 8.315238e+02 9.833367e+01 6.111881e+01 2.558529e+02 5.961090e+01
#> [816] 6.862099e+01 7.062613e+01 1.414927e+02 5.768371e+01 8.200038e+01
#> [821] 1.048464e+03 5.857459e+01 6.441491e+01 3.495532e+02 6.413400e+01
#> [826] 6.791335e+01 5.761136e+01 1.150261e+02 9.783103e+01 6.785414e+01
#> [831] 6.587289e+01 7.976776e+01 7.780424e+01 8.020049e+01 8.420820e+01
#> [836] 5.891111e+01 9.256684e+01 5.970597e+01 7.155329e+01 4.748319e+01
#> [841] 6.643247e+01 6.899154e+01 1.050365e+02 6.400591e+01 1.030288e+02
#> [846] 7.576328e+01 5.709638e+01 6.105953e+01 9.063576e+01 7.459514e+01
#> [851] 5.554582e+01 6.015888e+01 6.486503e+01 7.813653e+01 6.509394e+01
#> [856] 5.688628e+01 5.862225e+01 6.509540e+01 5.621669e+01 5.601035e+01
#> [861] 6.842063e+01 6.702753e+01 6.509193e+01 5.794017e+01 5.775849e+01
#> [866] 7.652327e+01 5.823645e+01 1.228037e+02 4.593174e+01 2.205000e+02
#> [871] 5.648245e+01 6.243569e+01 1.971238e+02 6.332403e+01 8.451560e+01
#> [876] 1.313231e+02 3.347843e+02 6.230567e+01 1.015223e+02 1.628365e+02
#> [881] 5.413262e+01 7.661336e+01 7.847360e+01 5.966578e+01 6.057354e+01
#> [886] 6.134273e+01 8.165148e+01 5.616842e+01 5.897077e+01 5.790702e+01
#> [891] 5.843520e+01 5.745441e+01 6.585141e+01 6.379222e+01 8.553052e+01
#> [896] 6.977077e+01 6.951151e+01 6.078822e+01 6.821102e+01 9.052645e+01
#> [901] 7.616603e+01 6.516727e+01 6.231692e+01 1.174469e+02 8.751639e+02
#> [906] 1.890452e+02 6.853905e+01 5.307525e+01 8.359832e+01 1.087927e+02
#> [911] 5.733600e+01 7.413282e+01 6.812047e+01 1.521075e+02 8.832374e+01
#> [916] 7.559007e+01 5.641374e+01 8.670714e+01 7.716254e+01 6.507049e+01
#> [921] 6.093538e+01 6.766269e+01 7.349775e+01 7.056003e+01 7.372770e+01
#> [926] 5.693976e+01 7.593417e+01 6.414667e+01 1.880627e+02 9.006995e+01
#> [931] 6.251149e+01 1.336812e+02 6.205972e+01 5.749610e+01 5.349768e+01
#> [936] 9.173024e+01 6.478905e+01 7.668363e+02 7.984305e+01 5.813707e+01
#> [941] 7.483788e+01 6.880984e+01 8.125175e+01 7.496507e+01 6.011738e+01
#> [946] 5.758444e+01 5.459512e+01 5.452985e+01 5.508779e+01 1.141991e+02
#> [951] 7.711544e+01 5.845538e+01 6.181069e+01 6.565423e+01 6.461329e+01
#> [956] 5.482467e+01 7.247379e+01 6.429342e+01 5.480560e+01 8.801688e+01
#> [961] 5.572390e+01 6.169601e+01 5.577257e+01 1.194846e+02 5.437160e+01
#> [966] 6.660950e+01 1.056445e+02 6.055921e+01 6.261109e+01 6.281931e+01
#> [971] 1.222260e+02 9.767586e+01 5.221333e+01 5.692838e+01 5.789856e+01
#> [976] 1.621645e+02 6.141506e+01 6.888366e+01 6.173896e+01 1.062768e+02
#> [981] 5.937631e+01 6.134590e+01 8.233342e+01 5.885546e+01 7.810274e+01
#> [986] 6.062978e+01 6.180225e+01 9.938376e+01 7.962578e+01 5.772374e+01
#> [991] 8.403815e+01 7.223131e+01 5.767946e+01 1.028213e+02 7.940749e+01
#> [996] 7.663437e+01 1.134844e+02 7.591718e+01 5.914011e+01 6.530218e+01