![[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] 2.436625e+09 2.655774e+02 2.886738e+03 2.905966e+03 6.825398e+01
#> [6] 6.801262e+01 8.145520e+01 1.457973e+02 1.189426e+02 1.958782e+02
#> [11] 1.723085e+02 6.714849e+01 9.941373e+01 7.927979e+01 1.061992e+03
#> [16] 7.795374e+01 7.039830e+01 7.292439e+01 3.376365e+02 4.261643e+02
#> [21] 7.017428e+01 1.308992e+02 8.661852e+01 6.649135e+01 1.223478e+03
#> [26] 1.856582e+03 1.761468e+02 7.363789e+01 8.385442e+01 8.343852e+01
#> [31] 3.395170e+02 1.078982e+02 2.209468e+02 1.472909e+02 3.170984e+02
#> [36] 1.537143e+02 5.919506e+02 2.120463e+02 9.990339e+01 1.708131e+02
#> [41] 7.361458e+01 6.656649e+01 7.626503e+01 6.410798e+01 7.297785e+01
#> [46] 6.963516e+01 6.793527e+01 2.456818e+02 6.672437e+01 8.617022e+01
#> [51] 8.230806e+01 7.049840e+01 6.312532e+01 6.439511e+01 1.110934e+02
#> [56] 1.254231e+02 2.686480e+02 7.322721e+01 7.235622e+01 8.055525e+01
#> [61] 7.561807e+01 2.490845e+02 1.016476e+02 1.343733e+02 1.001846e+02
#> [66] 4.842361e+02 2.899506e+02 7.822654e+01 1.023270e+02 7.653105e+01
#> [71] 5.920684e+02 1.325780e+02 7.222303e+01 6.547300e+01 1.087628e+03
#> [76] 6.716853e+01 1.305125e+02 3.613441e+02 7.034288e+01 7.362674e+01
#> [81] 9.288662e+01 9.122834e+01 7.817185e+01 7.112678e+01 8.784054e+01
#> [86] 7.042994e+01 8.624718e+01 9.827487e+01 1.486809e+03 1.131031e+02
#> [91] 7.332143e+01 2.477410e+02 9.457810e+01 1.184601e+02 1.184203e+02
#> [96] 1.008081e+02 8.353546e+01 1.132515e+02 1.103251e+02 6.453601e+01
#> [101] 7.828479e+01 1.375136e+02 8.044610e+01 6.806097e+01 9.452016e+01
#> [106] 1.253437e+02 6.765072e+01 7.594482e+01 3.347376e+02 7.131892e+03
#> [111] 1.850873e+02 9.956782e+01 7.916968e+01 7.306294e+02 2.316822e+02
#> [116] 7.420387e+01 6.402318e+01 5.788517e+02 3.806112e+02 1.957253e+02
#> [121] 7.206999e+01 5.439547e+02 2.144051e+02 8.500077e+01 1.496048e+02
#> [126] 1.644702e+02 7.067340e+01 1.012362e+02 9.368702e+01 7.804958e+01
#> [131] 8.996067e+01 7.938784e+01 8.078724e+01 8.037072e+01 7.414912e+01
#> [136] 9.042887e+01 1.276838e+02 9.338598e+01 8.962577e+01 1.219781e+02
#> [141] 6.729237e+01 7.379449e+01 9.387151e+02 6.632880e+01 1.737782e+02
#> [146] 6.834001e+01 1.981230e+02 1.090020e+03 7.661081e+01 8.179813e+01
#> [151] 1.214182e+02 1.293993e+02 1.520376e+02 7.319436e+01 6.411137e+01
#> [156] 7.196700e+01 1.911434e+02 1.092214e+02 1.842336e+02 3.128224e+02
#> [161] 1.581724e+03 6.382194e+04 1.352873e+02 7.625331e+01 3.762091e+09
#> [166] 1.155359e+08 3.888512e+07 1.368418e+02 6.345694e+01 9.219965e+01
#> [171] 8.995384e+01 6.897489e+01 8.046019e+01 8.848143e+01 2.180176e+02
#> [176] 6.756247e+01 7.404899e+01 1.009951e+02 6.842552e+01 7.065276e+01
#> [181] 8.507258e+01 9.272932e+01 7.255982e+01 7.489015e+01 1.883686e+02
#> [186] 2.162138e+02 2.870983e+02 1.113344e+02 3.624936e+02 6.936948e+01
#> [191] 6.512749e+01 6.188930e+01 7.656304e+01 7.104299e+01 6.628830e+01
#> [196] 8.423400e+01 6.163560e+01 8.951208e+01 9.664821e+02 1.072642e+02
#> [201] 5.915759e+02 1.149523e+02 8.684897e+01 3.135217e+02 1.012716e+02
#> [206] 7.548391e+01 1.208608e+02 8.842235e+01 9.630986e+01 1.107065e+02
#> [211] 1.219784e+03 1.166423e+02 3.815984e+03 2.165581e+03 7.062491e+01
#> [216] 1.035822e+02 9.336669e+02 1.151655e+02 7.496991e+02 7.256868e+01
#> [221] 2.053575e+02 9.225114e+01 4.248384e+03 1.273177e+02 2.031016e+02
#> [226] 7.130838e+01 1.685134e+02 7.431073e+01 3.716586e+02 1.255587e+03
#> [231] 7.603700e+01 7.178249e+01 7.967191e+01 7.776528e+01 1.066040e+02
#> [236] 8.843416e+01 6.584980e+01 8.798805e+01 6.941887e+01 8.277197e+01
#> [241] 6.992763e+01 8.128798e+01 7.912631e+01 6.367905e+01 6.702473e+01
#> [246] 1.524817e+02 8.982634e+01 1.360349e+02 1.277842e+02 2.031829e+02
#> [251] 1.141387e+02 9.282635e+01 9.489663e+01 8.212520e+01 7.608393e+01
#> [256] 8.575635e+01 9.165220e+01 1.346915e+02 6.804642e+01 5.912261e+03
#> [261] 1.028236e+02 3.651355e+03 6.667728e+01 9.952553e+02 8.053733e+01
#> [266] 1.000701e+02 6.973476e+01 7.409182e+01 6.691420e+01 1.322091e+02
#> [271] 7.070597e+01 8.063491e+01 9.371013e+01 6.828160e+01 6.716026e+01
#> [276] 7.836687e+01 7.219857e+01 7.048510e+01 1.475003e+02 7.391300e+01
#> [281] 7.383746e+01 7.691400e+01 8.768865e+01 8.284233e+01 5.026145e+02
#> [286] 7.288076e+01 9.909354e+01 7.852620e+01 4.391363e+02 1.118761e+03
#> [291] 9.278899e+01 7.126438e+01 7.429976e+01 6.715156e+01 8.573567e+01
#> [296] 7.452954e+01 5.662867e+02 1.092522e+02 9.235830e+01 6.887346e+01
#> [301] 2.373101e+02 1.658464e+02 1.235264e+02 1.854992e+02 2.068501e+04
#> [306] 1.012790e+02 1.502037e+02 1.219817e+02 7.421780e+01 8.131789e+01
#> [311] 6.356912e+01 6.270186e+01 6.713160e+01 7.655393e+01 1.072738e+02
#> [316] 2.211866e+03 1.157787e+02 8.922319e+01 6.722170e+01 5.680452e+02
#> [321] 8.802208e+01 7.378813e+01 6.878515e+01 6.328045e+01 8.135930e+01
#> [326] 9.858704e+01 5.790444e+02 7.924544e+01 6.801713e+01 9.775592e+01
#> [331] 7.717405e+01 6.750956e+01 3.652261e+02 1.472407e+08 7.555325e+01
#> [336] 9.377186e+01 1.309697e+02 6.695609e+01 7.546372e+01 9.480322e+01
#> [341] 1.107273e+05 1.098577e+02 1.071846e+02 1.029588e+02 8.297335e+01
#> [346] 1.904813e+02 3.382134e+02 8.967963e+01 1.733887e+03 6.719109e+01
#> [351] 7.288284e+01 1.708132e+02 1.345401e+02 6.547541e+01 1.125212e+02
#> [356] 7.688190e+01 9.587326e+01 1.235963e+02 1.544822e+02 2.941585e+02
#> [361] 1.121442e+02 7.508090e+01 7.508830e+01 6.961992e+01 2.042136e+02
#> [366] 7.744555e+11 7.695324e+01 9.707855e+01 7.564563e+01 1.037368e+02
#> [371] 7.029882e+01 1.257751e+02 7.797427e+01 8.539904e+01 2.325337e+02
#> [376] 8.867011e+01 7.475310e+01 9.621077e+01 8.480582e+01 8.966706e+01
#> [381] 6.990047e+01 6.668920e+01 6.575244e+01 2.817269e+02 7.055155e+01
#> [386] 3.476734e+02 7.547409e+01 7.699073e+01 6.610310e+01 1.395592e+02
#> [391] 1.064587e+02 5.863829e+01 1.209556e+02 6.698435e+01 7.450012e+01
#> [396] 1.250677e+02 7.389011e+01 1.083517e+02 6.434741e+01 1.328275e+02
#> [401] 3.130159e+02 9.517170e+01 7.131256e+01 1.095194e+02 7.223596e+01
#> [406] 3.006148e+02 1.098200e+02 7.824879e+01 1.092923e+02 1.198308e+02
#> [411] 1.808473e+02 2.569187e+02 8.378728e+01 1.169628e+02 7.745265e+01
#> [416] 7.397594e+01 9.211668e+01 8.353100e+01 7.985177e+01 6.544807e+01
#> [421] 6.427155e+01 9.127519e+01 6.248548e+01 6.368214e+01 7.344693e+01
#> [426] 8.919562e+01 3.292218e+02 9.267727e+01 6.808224e+01 6.840983e+01
#> [431] 8.445690e+01 6.561324e+01 6.665488e+01 7.116289e+01 2.053338e+02
#> [436] 1.573236e+02 9.113672e+01 6.127232e+01 6.583408e+01 3.484179e+02
#> [441] 6.324647e+01 6.221828e+01 1.104400e+03 7.048384e+01 1.536919e+02
#> [446] 1.115033e+02 7.091011e+01 1.191527e+02 7.244071e+01 5.000996e+02
#> [451] 6.719627e+01 6.488637e+01 1.173235e+02 1.197312e+02 7.505804e+02
#> [456] 6.670849e+01 7.660821e+01 6.688776e+01 7.335631e+01 1.053705e+02
#> [461] 1.072558e+02 2.207181e+04 1.166120e+02 7.989662e+02 7.293956e+01
#> [466] 4.127250e+02 9.787324e+01 1.831142e+02 5.999285e+01 6.898228e+01
#> [471] 9.864609e+01 8.399696e+01 9.685706e+01 9.162436e+01 1.606082e+02
#> [476] 8.354563e+01 1.308839e+02 6.461284e+01 6.583651e+01 6.839396e+01
#> [481] 7.120089e+01 8.187449e+01 8.308581e+01 1.708770e+02 6.756756e+01
#> [486] 8.003531e+01 1.369379e+02 8.298008e+01 6.928393e+01 9.781515e+02
#> [491] 2.465251e+04 4.137358e+02 1.515547e+02 4.210566e+02 7.085867e+01
#> [496] 9.670456e+01 1.875502e+02 1.249976e+02 6.888786e+01 6.391086e+01
#> [501] 1.462471e+02 1.580528e+02 1.124158e+02 1.023412e+02 8.237583e+01
#> [506] 2.609048e+02 6.979054e+01 6.744845e+01 9.175837e+01 7.550944e+01
#> [511] 7.875161e+01 7.855625e+01 8.751988e+01 8.535048e+01 7.238733e+01
#> [516] 7.409177e+01 6.663427e+01 7.468339e+01 7.053577e+01 9.108564e+01
#> [521] 1.205688e+02 2.021260e+02 4.035846e+03 1.090811e+02 1.741501e+02
#> [526] 1.758900e+02 8.020596e+01 1.056278e+02 8.633897e+01 1.839306e+02
#> [531] 5.158104e+02 1.515767e+02 6.730060e+01 6.941199e+01 1.056809e+02
#> [536] 1.303840e+02 7.989223e+01 7.466255e+01 7.897493e+01 9.100523e+01
#> [541] 6.706002e+01 2.192044e+02 7.251519e+01 6.720545e+01 8.989851e+01
#> [546] 7.459082e+01 9.125794e+01 2.324718e+02 1.201616e+02 2.031199e+02
#> [551] 1.632092e+02 7.564753e+01 8.193298e+01 1.616048e+02 7.523349e+01
#> [556] 7.195446e+01 7.237358e+01 7.375111e+01 1.129348e+02 6.689696e+01
#> [561] 6.916959e+01 5.234616e+02 2.465197e+02 7.948167e+01 7.595758e+01
#> [566] 9.744982e+01 1.132474e+02 8.596771e+01 9.196322e+01 6.744101e+01
#> [571] 7.451700e+01 6.438285e+01 7.212286e+01 1.856109e+02 2.019440e+02
#> [576] 7.683287e+01 1.209249e+02 6.547437e+01 8.556369e+01 1.809764e+03
#> [581] 1.153907e+02 7.735709e+01 8.000823e+01 1.156724e+02 1.144153e+02
#> [586] 8.772074e+01 2.170836e+02 7.084055e+01 6.527414e+01 1.261859e+02
#> [591] 8.886126e+01 7.664432e+01 1.143455e+02 7.570902e+01 1.442578e+02
#> [596] 9.660851e+01 6.391528e+01 3.525607e+02 2.910592e+02 7.235956e+01
#> [601] 6.806685e+01 6.713665e+01 8.551483e+01 6.981182e+01 9.163834e+01
#> [606] 7.997890e+01 8.591471e+01 1.002289e+02 2.672882e+02 1.227084e+03
#> [611] 1.115409e+02 9.218380e+01 1.437589e+02 6.513548e+01 7.508505e+01
#> [616] 1.323437e+03 1.020849e+02 8.926576e+01 9.796409e+01 7.126015e+01
#> [621] 1.607451e+02 1.252074e+02 9.314665e+01 1.582276e+02 6.464432e+01
#> [626] 7.739549e+01 7.877237e+01 7.507219e+01 7.447831e+01 8.536100e+01
#> [631] 8.921642e+01 7.943839e+01 6.237600e+01 6.972810e+01 7.524252e+01
#> [636] 8.873789e+01 4.855908e+05 3.991597e+03 3.213481e+02 3.612262e+02
#> [641] 1.788215e+03 5.990090e+01 1.139348e+02 2.380860e+02 9.803459e+01
#> [646] 8.514470e+01 8.394247e+01 7.166238e+01 7.623841e+01 6.528485e+01
#> [651] 8.321496e+01 8.231818e+01 1.364402e+02 1.127559e+02 2.098806e+02
#> [656] 1.166805e+03 8.919762e+01 8.751316e+01 1.810110e+02 8.600099e+01
#> [661] 9.134030e+01 1.831061e+02 7.292736e+02 8.898127e+01 7.394713e+01
#> [666] 1.606636e+02 9.729386e+01 1.422240e+02 6.866213e+01 7.820121e+02
#> [671] 1.617610e+02 6.634183e+01 7.159261e+01 1.300340e+02 7.306742e+01
#> [676] 1.044114e+04 1.177919e+02 8.272987e+01 7.825166e+01 7.151078e+01
#> [681] 6.882388e+01 7.981432e+01 7.964302e+01 1.026630e+02 9.187948e+01
#> [686] 8.279169e+01 1.407909e+03 6.879096e+01 9.036876e+01 7.200165e+01
#> [691] 7.633145e+01 3.870355e+05 1.269991e+02 8.703728e+01 6.770692e+01
#> [696] 1.457109e+02 6.878852e+01 6.394627e+01 7.060079e+01 7.028284e+01
#> [701] 8.632508e+01 9.481230e+01 3.282710e+02 1.176502e+02 1.905265e+02
#> [706] 2.102244e+02 7.186022e+01 7.363411e+01 1.146841e+02 6.704666e+01
#> [711] 1.461622e+02 1.140327e+02 1.085387e+02 1.238029e+02 8.211497e+01
#> [716] 1.478044e+02 1.105876e+02 1.357321e+02 6.553881e+01 6.870070e+01
#> [721] 9.959880e+01 8.890154e+01 1.623911e+02 6.542121e+01 7.542184e+01
#> [726] 1.314784e+02 1.134037e+02 7.586932e+01 6.373297e+01 6.882074e+01
#> [731] 9.377482e+01 6.903944e+01 1.138334e+02 6.835242e+01 8.006814e+01
#> [736] 7.593436e+01 7.537611e+01 8.780127e+03 1.861100e+02 1.070466e+02
#> [741] 6.961570e+01 7.068021e+01 9.760843e+01 1.326639e+02 2.705025e+02
#> [746] 1.026049e+02 6.919780e+01 6.568787e+01 6.292747e+01 3.488074e+02
#> [751] 1.173503e+02 1.163057e+02 1.015651e+04 1.921131e+02 1.105412e+02
#> [756] 1.039982e+02 2.834952e+02 6.481855e+01 7.418723e+01 9.926671e+01
#> [761] 7.555814e+01 7.612084e+01 9.371679e+01 4.932100e+03 6.952848e+01
#> [766] 1.812402e+02 8.794948e+01 2.333974e+03 9.750448e+01 7.114455e+01
#> [771] 1.504705e+02 7.709267e+01 8.059194e+01 1.865757e+02 6.869445e+01
#> [776] 9.402495e+01 6.997929e+01 7.815152e+01 9.070200e+01 8.832504e+01
#> [781] 9.191409e+01 6.252285e+01 1.241353e+02 8.848285e+02 3.507674e+04
#> [786] 7.898254e+01 6.717571e+01 1.670388e+02 6.755418e+01 6.241350e+02
#> [791] 1.491287e+02 6.580884e+01 6.732159e+01 6.711205e+01 8.217686e+01
#> [796] 6.448309e+01 7.307365e+01 2.232591e+02 1.745499e+02 6.615939e+02
#> [801] 1.251475e+04 1.642321e+02 6.386348e+01 7.227089e+01 7.001619e+01
#> [806] 8.899822e+01 8.899483e+01 1.250895e+02 1.254232e+02 4.275212e+02
#> [811] 7.999675e+01 6.303217e+01 7.702600e+01 7.233577e+01 8.547734e+01
#> [816] 7.799041e+01 6.878422e+01 7.627892e+01 4.391326e+04 1.679487e+02
#> [821] 1.679682e+02 6.952248e+01 2.100989e+02 1.505052e+02 3.922321e+02
#> [826] 1.395678e+02 6.879577e+01 1.025935e+02 7.780419e+01 2.996182e+02
#> [831] 7.515447e+01 7.321643e+01 7.137588e+01 6.337892e+01 2.252501e+02
#> [836] 2.319990e+03 2.617805e+02 4.133164e+04 9.372788e+01 7.726869e+01
#> [841] 7.473600e+01 1.811667e+02 7.361870e+01 1.369371e+02 9.868634e+01
#> [846] 1.014356e+02 6.147858e+02 2.418860e+02 6.470477e+01 1.197714e+02
#> [851] 1.506307e+02 1.756653e+02 8.181498e+01 8.800544e+01 7.459428e+01
#> [856] 7.567330e+01 1.263014e+02 2.451597e+04 4.866865e+02 6.485495e+01
#> [861] 8.838477e+01 7.026039e+01 9.062273e+01 6.556088e+01 6.403473e+01
#> [866] 2.700998e+02 7.911891e+01 9.665168e+01 4.138182e+06 1.119570e+05
#> [871] 9.941660e+01 1.062703e+02 7.643860e+01 9.336610e+01 1.066506e+02
#> [876] 1.059983e+02 7.626239e+01 7.910415e+02 1.177132e+02 1.635615e+02
#> [881] 8.011304e+01 8.598674e+01 6.467025e+01 7.745594e+01 7.039532e+01
#> [886] 7.545259e+01 8.250195e+01 8.400794e+01 3.390841e+02 1.551088e+02
#> [891] 6.546732e+01 1.645600e+02 7.211271e+01 1.233400e+02 5.777572e+02
#> [896] 7.266401e+01 7.357305e+01 1.040546e+02 1.352558e+03 1.010677e+02
#> [901] 8.742532e+01 4.643622e+02 2.381967e+02 1.429726e+02 1.089301e+02
#> [906] 1.139579e+02 8.358617e+01 7.037728e+01 1.220570e+02 7.780766e+01
#> [911] 1.291937e+02 1.161754e+02 7.276940e+01 2.078181e+02 6.755359e+01
#> [916] 1.123922e+02 6.050147e+01 6.068661e+01 6.181063e+01 3.272063e+02
#> [921] 6.381742e+01 9.525286e+01 7.756382e+01 1.286529e+02 6.230034e+01
#> [926] 6.702722e+01 8.684585e+01 7.009501e+01 8.189363e+01 7.333696e+01
#> [931] 7.446244e+01 8.741742e+01 8.728662e+01 1.267753e+02 1.142578e+02
#> [936] 7.128662e+01 7.495903e+01 7.231153e+01 6.956964e+01 9.780326e+02
#> [941] 9.602808e+01 8.128579e+01 2.949292e+02 5.620475e+02 8.228989e+01
#> [946] 9.425967e+01 1.020636e+02 8.652656e+01 1.835339e+02 7.532456e+01
#> [951] 8.019890e+01 8.858181e+01 1.303813e+02 7.254478e+01 6.647431e+01
#> [956] 3.138752e+02 8.979351e+01 7.583383e+01 7.293587e+01 8.354520e+01
#> [961] 7.511791e+01 8.108720e+01 7.427921e+01 7.012914e+01 7.259914e+01
#> [966] 6.885846e+01 7.829989e+01 6.873052e+01 7.511149e+01 7.843561e+01
#> [971] 9.056782e+01 1.199712e+02 5.905213e+01 6.245395e+01 6.496077e+01
#> [976] 6.941488e+01 7.259716e+01 7.122996e+01 7.459109e+01 7.915826e+01
#> [981] 2.180100e+02 7.777572e+01 6.845907e+01 8.867043e+01 1.286216e+02
#> [986] 7.798795e+01 6.364073e+01 9.680225e+01 1.088981e+02 9.917795e+01
#> [991] 4.140843e+02 9.969872e+01 1.904908e+02 6.761678e+01 9.005532e+01
#> [996] 8.318257e+01 7.284497e+01 1.022245e+02 1.093596e+02 1.311799e+02
doseFunction(
model,
alpha2 = samples@data$alpha2,
beta = samples@data$beta,
grade = 2L
)(x = 0.25)
#> [1] 1.750826e+06 1.140200e+02 2.830166e+02 4.798498e+02 5.915681e+01
#> [6] 5.859984e+01 6.433581e+01 8.326170e+01 7.675165e+01 1.139845e+02
#> [11] 1.005380e+02 6.049247e+01 7.934413e+01 6.722656e+01 3.697254e+02
#> [16] 6.496790e+01 6.306447e+01 6.217163e+01 1.535798e+02 1.219032e+02
#> [21] 6.148686e+01 7.842686e+01 6.629929e+01 5.794765e+01 2.140532e+02
#> [26] 1.713573e+02 1.039753e+02 6.230172e+01 6.846018e+01 5.565871e+01
#> [31] 1.124388e+02 8.039983e+01 1.089645e+02 9.208207e+01 1.399082e+02
#> [36] 6.079775e+01 1.842170e+02 9.179894e+01 7.420715e+01 1.049973e+02
#> [41] 6.379846e+01 6.013285e+01 6.724858e+01 5.804097e+01 6.414232e+01
#> [46] 5.987935e+01 5.911685e+01 1.285497e+02 5.952522e+01 5.874340e+01
#> [51] 6.328526e+01 5.891702e+01 5.582118e+01 5.725328e+01 7.530949e+01
#> [56] 7.105804e+01 1.309635e+02 6.227162e+01 5.995251e+01 6.599650e+01
#> [61] 6.451974e+01 9.531734e+01 6.916644e+01 7.424861e+01 5.237617e+01
#> [66] 1.167421e+02 1.343474e+02 5.844791e+01 6.929140e+01 6.126776e+01
#> [71] 1.474542e+02 7.257494e+01 5.930590e+01 5.566788e+01 2.827663e+02
#> [76] 5.663923e+01 8.664225e+01 9.735533e+01 5.706631e+01 5.601888e+01
#> [81] 7.587150e+01 5.791456e+01 6.639175e+01 6.313164e+01 6.962265e+01
#> [86] 6.069531e+01 6.600654e+01 7.566892e+01 8.082916e+01 8.379883e+01
#> [91] 5.994472e+01 1.075564e+02 6.215633e+01 7.735786e+01 7.359631e+01
#> [96] 6.317160e+01 5.399882e+01 7.656950e+01 7.399559e+01 5.688218e+01
#> [101] 6.393468e+01 9.517676e+01 6.387692e+01 5.704150e+01 6.925297e+01
#> [106] 8.960211e+01 5.770931e+01 6.353387e+01 7.059056e+01 3.948493e+02
#> [111] 9.214993e+01 6.480624e+01 5.606328e+01 3.075309e+02 1.269730e+02
#> [116] 5.982988e+01 5.535657e+01 1.484622e+02 1.154905e+02 1.043150e+02
#> [121] 6.004096e+01 1.760841e+02 8.955064e+01 5.629088e+01 7.666434e+01
#> [126] 8.640774e+01 5.451516e+01 6.603862e+01 6.086628e+01 6.525163e+01
#> [131] 6.839621e+01 6.561838e+01 6.689175e+01 6.370294e+01 6.180127e+01
#> [136] 6.156568e+01 7.753166e+01 6.172390e+01 6.845108e+01 8.481162e+01
#> [141] 6.113282e+01 6.432827e+01 2.096594e+02 5.490756e+01 8.373387e+01
#> [146] 5.353480e+01 1.107688e+02 2.938961e+02 5.730219e+01 6.553728e+01
#> [151] 7.266338e+01 7.000646e+01 7.542954e+01 5.687630e+01 5.699750e+01
#> [156] 5.797441e+01 9.833288e+01 7.006402e+01 9.444568e+01 1.815060e+02
#> [161] 3.906329e+02 3.153704e+03 7.969473e+01 5.645038e+01 1.803735e+06
#> [166] 4.561744e+04 3.335402e+04 8.185886e+01 5.359030e+01 6.688755e+01
#> [171] 6.537643e+01 5.857924e+01 6.223674e+01 6.379973e+01 1.023725e+02
#> [176] 6.012239e+01 6.060752e+01 7.354838e+01 5.881894e+01 6.042237e+01
#> [181] 6.660705e+01 6.785294e+01 6.217616e+01 5.832552e+01 6.954861e+01
#> [186] 9.274524e+01 1.343467e+02 7.459304e+01 1.525765e+02 6.031260e+01
#> [191] 5.359982e+01 5.566812e+01 6.438167e+01 6.180769e+01 5.743567e+01
#> [196] 5.647479e+01 5.451452e+01 6.902695e+01 2.320344e+02 7.420596e+01
#> [201] 9.018363e+01 7.418140e+01 6.417761e+01 1.358599e+02 7.395496e+01
#> [206] 5.613648e+01 8.526077e+01 7.064580e+01 7.018543e+01 8.399683e+01
#> [211] 4.050640e+02 6.447088e+01 2.000751e+02 3.569899e+02 5.744346e+01
#> [216] 7.208317e+01 2.233840e+02 7.296114e+01 1.233722e+02 6.133397e+01
#> [221] 1.244771e+02 6.511037e+01 4.246327e+02 8.119920e+01 1.113732e+02
#> [226] 5.942880e+01 1.013040e+02 5.898568e+01 1.477294e+02 3.478104e+02
#> [231] 6.528444e+01 6.191613e+01 6.529964e+01 6.308337e+01 7.380705e+01
#> [236] 5.905292e+01 5.479860e+01 6.719492e+01 5.637587e+01 6.856792e+01
#> [241] 5.967045e+01 7.091572e+01 6.699800e+01 5.905246e+01 6.117763e+01
#> [246] 1.044908e+02 6.120686e+01 7.435984e+01 7.668172e+01 1.105078e+02
#> [251] 8.710925e+01 6.980088e+01 7.601085e+01 6.541760e+01 6.327341e+01
#> [256] 5.582614e+01 7.374622e+01 9.678783e+01 6.176284e+01 6.070612e+02
#> [261] 5.999817e+01 2.660995e+02 5.734160e+01 1.707223e+02 5.797089e+01
#> [266] 6.378030e+01 5.741539e+01 6.512077e+01 5.748972e+01 9.211980e+01
#> [271] 5.475681e+01 5.907030e+01 6.924429e+01 5.867681e+01 5.513871e+01
#> [276] 6.111302e+01 6.061865e+01 6.164279e+01 8.373805e+01 6.111788e+01
#> [281] 6.115899e+01 6.637370e+01 6.801793e+01 6.564351e+01 1.377519e+02
#> [286] 6.002996e+01 6.776829e+01 6.176490e+01 1.722522e+02 3.972758e+02
#> [291] 7.015755e+01 5.565307e+01 6.410213e+01 5.769941e+01 6.411773e+01
#> [296] 5.876142e+01 1.355241e+02 7.305511e+01 6.796656e+01 5.872118e+01
#> [301] 1.363804e+02 6.126360e+01 6.129857e+01 9.755653e+01 1.422109e+03
#> [306] 7.552305e+01 8.139174e+01 8.157292e+01 6.137444e+01 5.685443e+01
#> [311] 5.753766e+01 5.743787e+01 5.900917e+01 6.451644e+01 6.691522e+01
#> [316] 7.540169e+02 8.592929e+01 6.924214e+01 5.796221e+01 2.094116e+02
#> [321] 7.229403e+01 6.099018e+01 5.934046e+01 5.740882e+01 5.627231e+01
#> [326] 7.061667e+01 1.361975e+02 6.489729e+01 5.784827e+01 7.008946e+01
#> [331] 6.363114e+01 5.351429e+01 1.355886e+02 3.637919e+05 5.563704e+01
#> [336] 7.781343e+01 8.782831e+01 5.991536e+01 6.187865e+01 6.897172e+01
#> [341] 4.196318e+03 5.427482e+01 6.642353e+01 7.298966e+01 6.604656e+01
#> [346] 5.061137e+01 1.450340e+02 7.006884e+01 5.610701e+02 6.051897e+01
#> [351] 6.427711e+01 9.162129e+01 7.051575e+01 5.177434e+01 5.911415e+01
#> [356] 6.051273e+01 6.615777e+01 5.836406e+01 8.272472e+01 1.072617e+02
#> [361] 5.852556e+01 6.674978e+01 6.719395e+01 5.971639e+01 1.046464e+02
#> [366] 1.846505e+09 6.738949e+01 7.009948e+01 6.252153e+01 6.036827e+01
#> [371] 5.991687e+01 8.753698e+01 6.601732e+01 6.353939e+01 1.014294e+02
#> [376] 5.976336e+01 6.074780e+01 7.517013e+01 6.992282e+01 7.020192e+01
#> [381] 6.086017e+01 5.984293e+01 6.016331e+01 1.311417e+02 6.028861e+01
#> [386] 1.529286e+02 6.141148e+01 6.240933e+01 5.590302e+01 8.519963e+01
#> [391] 6.587924e+01 4.920468e+01 8.116678e+01 5.998781e+01 6.365258e+01
#> [396] 9.095644e+01 6.035755e+01 7.785688e+01 5.598986e+01 7.507263e+01
#> [401] 1.553529e+02 6.295534e+01 5.219567e+01 7.072703e+01 5.897472e+01
#> [406] 6.237662e+01 7.854131e+01 6.055703e+01 8.206223e+01 7.953141e+01
#> [411] 8.667073e+01 1.052919e+02 6.599023e+01 7.462882e+01 6.108220e+01
#> [416] 6.188316e+01 6.844449e+01 5.484290e+01 5.556917e+01 5.803917e+01
#> [421] 5.690958e+01 7.534518e+01 5.761432e+01 5.907502e+01 6.563781e+01
#> [426] 6.048329e+01 1.274183e+02 7.346700e+01 6.028778e+01 6.089598e+01
#> [431] 6.414782e+01 5.883365e+01 5.771026e+01 6.316525e+01 1.104388e+02
#> [436] 9.628509e+01 6.282969e+01 5.708022e+01 5.955346e+01 1.309671e+02
#> [441] 5.618947e+01 5.570422e+01 3.547277e+02 6.360595e+01 1.048499e+02
#> [446] 6.587115e+01 6.010183e+01 8.918503e+01 5.380364e+01 1.363249e+02
#> [451] 5.492759e+01 5.627865e+01 7.837888e+01 7.134122e+01 1.010305e+02
#> [456] 6.119157e+01 6.584267e+01 5.916221e+01 6.282890e+01 6.256774e+01
#> [461] 6.000943e+01 2.495772e+03 6.895405e+01 2.434340e+02 6.123091e+01
#> [466] 1.488765e+02 6.872754e+01 7.292922e+01 5.243667e+01 5.856285e+01
#> [471] 6.828378e+01 5.999455e+01 6.211605e+01 7.750210e+01 1.171105e+02
#> [476] 6.597556e+01 8.628999e+01 5.833815e+01 6.018628e+01 6.104023e+01
#> [481] 6.114327e+01 6.609227e+01 6.367705e+01 7.583796e+01 5.932652e+01
#> [486] 6.803865e+01 8.579033e+01 6.949777e+01 6.085018e+01 1.573856e+02
#> [491] 2.684143e+03 1.130405e+02 5.853834e+01 1.383229e+02 6.189558e+01
#> [496] 6.966657e+01 9.514532e+01 6.939555e+01 5.899652e+01 5.719124e+01
#> [501] 8.593538e+01 6.774926e+01 6.938859e+01 7.595106e+01 6.554095e+01
#> [506] 1.071605e+02 6.047097e+01 5.987037e+01 6.182948e+01 6.028691e+01
#> [511] 6.634157e+01 6.084190e+01 6.769899e+01 6.727378e+01 6.002802e+01
#> [516] 6.216582e+01 5.724203e+01 6.129069e+01 5.223527e+01 5.380042e+01
#> [521] 7.365160e+01 6.688502e+01 5.644974e+02 6.875789e+01 9.114609e+01
#> [526] 9.695951e+01 6.185235e+01 5.546736e+01 6.850133e+01 9.919688e+01
#> [531] 1.369792e+02 9.247545e+01 5.701949e+01 6.128064e+01 6.753431e+01
#> [536] 7.232118e+01 6.299107e+01 5.664075e+01 6.411214e+01 6.704872e+01
#> [541] 5.952429e+01 8.753926e+01 5.720980e+01 5.518012e+01 7.003265e+01
#> [546] 5.273923e+01 5.568058e+01 1.345462e+02 8.332930e+01 1.061374e+02
#> [551] 8.327783e+01 5.996898e+01 6.160429e+01 8.274087e+01 6.050962e+01
#> [556] 5.908346e+01 6.584626e+01 6.274532e+01 7.805507e+01 6.086920e+01
#> [561] 6.275771e+01 2.011422e+02 1.424677e+02 6.367399e+01 6.734300e+01
#> [566] 7.359689e+01 8.693518e+01 7.035107e+01 7.535776e+01 5.411613e+01
#> [571] 5.741245e+01 5.532289e+01 6.131214e+01 8.090942e+01 1.141682e+02
#> [576] 5.680270e+01 6.340729e+01 5.184213e+01 6.107060e+01 3.704064e+02
#> [581] 6.412454e+01 6.486253e+01 6.349834e+01 8.004968e+01 7.793983e+01
#> [586] 7.076962e+01 9.945847e+01 5.943827e+01 5.729985e+01 8.094039e+01
#> [591] 6.614704e+01 6.398105e+01 7.666535e+01 5.947388e+01 6.245965e+01
#> [596] 6.912653e+01 5.935517e+01 1.932466e+02 9.890019e+01 6.403732e+01
#> [601] 6.176980e+01 6.295437e+01 7.515961e+01 6.235376e+01 6.608409e+01
#> [606] 6.378938e+01 6.553304e+01 6.844220e+01 8.968366e+01 2.600136e+02
#> [611] 6.997832e+01 5.980173e+01 8.410107e+01 5.762242e+01 5.628561e+01
#> [616] 1.451250e+02 6.091050e+01 6.594856e+01 7.638499e+01 5.728114e+01
#> [621] 7.167909e+01 8.610394e+01 7.451105e+01 8.723994e+01 5.560665e+01
#> [626] 6.453940e+01 6.474605e+01 6.371185e+01 6.120827e+01 6.832128e+01
#> [631] 6.430248e+01 4.868926e+01 5.669519e+01 5.970211e+01 6.119974e+01
#> [636] 6.483770e+01 3.196268e+03 1.869365e+02 5.124253e+01 4.708875e+01
#> [641] 3.194600e+02 5.303637e+01 7.497370e+01 7.121920e+01 6.666690e+01
#> [646] 6.325285e+01 5.923121e+01 5.857163e+01 6.368858e+01 5.703676e+01
#> [651] 6.390088e+01 6.202632e+01 6.889725e+01 8.436909e+01 7.743693e+01
#> [656] 3.446410e+02 7.192367e+01 6.952040e+01 9.254813e+01 6.468145e+01
#> [661] 6.795079e+01 8.116979e+01 1.294236e+02 6.908372e+01 5.461185e+01
#> [666] 6.646353e+01 7.059175e+01 7.434668e+01 5.803133e+01 1.535518e+02
#> [671] 9.018541e+01 5.163666e+01 5.712476e+01 8.710046e+01 5.842226e+01
#> [676] 6.190998e+02 7.230136e+01 6.674638e+01 5.807606e+01 5.936340e+01
#> [681] 5.841060e+01 6.285911e+01 6.286962e+01 7.234579e+01 6.364165e+01
#> [686] 6.236875e+01 3.606478e+02 6.051006e+01 7.401579e+01 5.050852e+01
#> [691] 6.308584e+01 8.185562e+03 6.402281e+01 6.407454e+01 5.976358e+01
#> [696] 9.347711e+01 6.037843e+01 5.815949e+01 6.215544e+01 6.064624e+01
#> [701] 6.820259e+01 6.279191e+01 1.160713e+02 7.471057e+01 9.572550e+01
#> [706] 6.867091e+01 6.350893e+01 6.495325e+01 7.199821e+01 5.926126e+01
#> [711] 9.776032e+01 6.250795e+01 7.260989e+01 7.857486e+01 6.231552e+01
#> [716] 9.231278e+01 7.327080e+01 5.679954e+01 5.055086e+01 5.913801e+01
#> [721] 7.627248e+01 5.196310e+01 9.461919e+01 5.876593e+01 6.277601e+01
#> [726] 6.232119e+01 6.178064e+01 5.480892e+01 5.392763e+01 5.717233e+01
#> [731] 6.784223e+01 5.495751e+01 6.974214e+01 5.708353e+01 6.501525e+01
#> [736] 5.708952e+01 5.801089e+01 8.348956e+02 1.085690e+02 6.884230e+01
#> [741] 5.834068e+01 6.065851e+01 7.885706e+01 6.121873e+01 8.940547e+01
#> [746] 7.980381e+01 5.881113e+01 5.630149e+01 5.661802e+01 9.706992e+01
#> [751] 8.349503e+01 7.218255e+01 2.375055e+02 9.470819e+01 6.743874e+01
#> [756] 7.512688e+01 1.238303e+02 5.640774e+01 5.986344e+01 7.677157e+01
#> [761] 6.082927e+01 6.430102e+01 6.620045e+01 4.119377e+02 5.852717e+01
#> [766] 7.094951e+01 6.941064e+01 5.114587e+02 7.039502e+01 5.536690e+01
#> [771] 9.007149e+01 6.724895e+01 6.245180e+01 6.282322e+01 5.555120e+01
#> [776] 6.326747e+01 5.741322e+01 6.461535e+01 6.794875e+01 6.450646e+01
#> [781] 6.942527e+01 5.743312e+01 8.094618e+01 1.557238e+02 1.877947e+03
#> [786] 6.129941e+01 6.002053e+01 1.125558e+02 5.604225e+01 1.732240e+02
#> [791] 6.505392e+01 5.541663e+01 5.715997e+01 5.591692e+01 6.918582e+01
#> [796] 5.984030e+01 6.529056e+01 1.313749e+02 1.052211e+02 2.202994e+02
#> [801] 1.665199e+03 9.732619e+01 5.780130e+01 6.077962e+01 6.000947e+01
#> [806] 6.772966e+01 6.134597e+01 7.905640e+01 6.795008e+01 1.857893e+02
#> [811] 6.616848e+01 5.707093e+01 6.614196e+01 6.377919e+01 7.146226e+01
#> [816] 6.595652e+01 6.013756e+01 6.619093e+01 2.689056e+03 7.072695e+01
#> [821] 1.005614e+02 5.732632e+01 1.137125e+02 7.628497e+01 1.008369e+02
#> [826] 7.373274e+01 5.933978e+01 7.091237e+01 5.650056e+01 1.297887e+02
#> [831] 5.639182e+01 5.799184e+01 5.879136e+01 5.621966e+01 1.184280e+02
#> [836] 2.440315e+02 1.440273e+02 3.453120e+03 7.214102e+01 6.268630e+01
#> [841] 6.448745e+01 8.554343e+01 5.832913e+01 8.722194e+01 6.582989e+01
#> [846] 7.431132e+01 2.639199e+02 1.356705e+02 5.601449e+01 8.476130e+01
#> [851] 7.441778e+01 6.950004e+01 5.982583e+01 7.052997e+01 6.358196e+01
#> [856] 6.417141e+01 7.645163e+01 2.629087e+03 1.080227e+02 5.886514e+01
#> [861] 7.138050e+01 6.259860e+01 7.133997e+01 5.953554e+01 5.836041e+01
#> [866] 1.366696e+02 5.870734e+01 6.525156e+01 1.667772e+04 4.264704e+03
#> [871] 6.514155e+01 7.391791e+01 6.287422e+01 6.639691e+01 6.856083e+01
#> [876] 6.504994e+01 6.437646e+01 2.442219e+02 6.509416e+01 6.822472e+01
#> [881] 6.698168e+01 7.186125e+01 5.969576e+01 6.754058e+01 5.748479e+01
#> [886] 6.534351e+01 5.760685e+01 6.458894e+01 1.010981e+02 7.915152e+01
#> [891] 5.536337e+01 8.812646e+01 5.832976e+01 8.364636e+01 1.408843e+02
#> [896] 5.974992e+01 5.561478e+01 7.757402e+01 2.272371e+02 6.827818e+01
#> [901] 6.333192e+01 1.365949e+02 1.317013e+02 9.777280e+01 7.022338e+01
#> [906] 8.957526e+01 7.253355e+01 5.804759e+01 8.103131e+01 6.216290e+01
#> [911] 8.310076e+01 7.587367e+01 5.847828e+01 7.183789e+01 5.712948e+01
#> [916] 5.966182e+01 5.217882e+01 5.594073e+01 5.675262e+01 1.876435e+02
#> [921] 5.896353e+01 7.606439e+01 6.580450e+01 7.685641e+01 5.593961e+01
#> [926] 6.090338e+01 5.890497e+01 5.915237e+01 6.712425e+01 5.873618e+01
#> [931] 6.396964e+01 7.140063e+01 7.183264e+01 7.483129e+01 6.718147e+01
#> [936] 5.818393e+01 5.793089e+01 5.938519e+01 6.067106e+01 3.129139e+02
#> [941] 4.812229e+01 5.425802e+01 7.943130e+01 2.146868e+02 6.671391e+01
#> [946] 6.853714e+01 6.892495e+01 6.266153e+01 3.684980e+01 5.431476e+01
#> [951] 6.083440e+01 6.598951e+01 8.327134e+01 6.099779e+01 5.656510e+01
#> [956] 1.201095e+02 5.596546e+01 5.758087e+01 5.795682e+01 6.633953e+01
#> [961] 6.234446e+01 6.173168e+01 5.957658e+01 5.661757e+01 4.220235e+01
#> [966] 6.155284e+01 6.095779e+01 5.870534e+01 6.210972e+01 6.251193e+01
#> [971] 6.730319e+01 7.724341e+01 5.503548e+01 5.873886e+01 5.975218e+01
#> [976] 6.189738e+01 6.357470e+01 6.167089e+01 6.141181e+01 6.267319e+01
#> [981] 7.717960e+01 6.054202e+01 5.761316e+01 6.308088e+01 7.145061e+01
#> [986] 5.954127e+01 5.679304e+01 7.425641e+01 7.454707e+01 6.820912e+01
#> [991] 2.065972e+02 7.337769e+01 1.084037e+02 6.106227e+01 7.409658e+01
#> [996] 6.660890e+01 5.798989e+01 7.733150e+01 7.232025e+01 8.788145e+01