![[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] 1.423457e+05 1.835495e+02 8.162403e+01 7.337644e+01 7.118981e+01
#> [6] 9.945067e+01 6.692337e+01 7.159096e+01 7.329037e+01 7.459835e+01
#> [11] 6.488263e+01 6.410400e+01 7.486712e+01 6.173720e+01 1.802438e+04
#> [16] 7.229418e+01 4.267514e+02 1.752347e+02 1.163487e+02 7.515198e+01
#> [21] 1.640667e+02 7.159544e+01 8.453224e+01 1.070889e+02 2.541382e+02
#> [26] 1.329870e+02 2.136092e+02 7.816437e+01 4.170615e+02 8.989573e+01
#> [31] 6.398752e+01 6.060075e+02 2.265689e+03 1.924708e+02 1.110233e+02
#> [36] 9.592029e+01 1.058829e+02 6.580108e+01 1.056868e+02 8.979572e+01
#> [41] 1.829590e+02 9.263000e+01 6.589516e+01 1.146246e+02 4.820375e+04
#> [46] 1.352398e+03 3.317249e+11 3.197594e+03 9.163185e+02 1.368836e+02
#> [51] 1.594115e+03 7.040432e+01 7.969901e+01 3.940009e+02 7.550860e+01
#> [56] 9.064809e+01 1.006343e+02 6.672303e+01 1.132836e+02 9.331714e+01
#> [61] 2.897642e+02 1.470642e+02 7.886501e+01 1.277141e+02 6.307545e+01
#> [66] 6.030007e+01 8.383764e+01 8.610886e+01 7.246375e+01 8.908044e+01
#> [71] 7.357020e+01 1.243909e+02 7.327769e+01 6.368555e+01 1.423172e+02
#> [76] 2.325108e+03 9.615272e+01 1.184034e+02 8.968470e+01 6.792433e+01
#> [81] 1.180424e+02 1.307932e+02 1.291011e+02 9.374879e+01 8.755295e+01
#> [86] 1.178481e+02 3.372366e+02 2.710769e+02 9.313598e+02 8.818035e+01
#> [91] 1.441263e+02 6.526176e+01 7.200859e+01 6.404929e+01 6.637330e+01
#> [96] 8.583599e+01 7.913978e+01 7.045036e+01 8.108795e+01 1.085093e+02
#> [101] 7.594455e+01 8.614321e+01 7.818382e+01 3.627758e+02 1.607405e+02
#> [106] 2.903762e+03 1.254284e+02 6.930794e+01 9.741448e+01 7.138673e+01
#> [111] 1.106909e+02 6.998952e+01 6.668146e+01 1.676708e+02 8.045478e+01
#> [116] 1.186248e+02 5.151051e+02 2.511858e+12 1.510678e+02 1.298347e+03
#> [121] 9.943293e+02 5.241912e+02 1.547439e+02 6.284167e+01 7.764019e+01
#> [126] 2.612067e+02 7.324045e+01 8.661346e+01 1.177638e+02 5.979187e+01
#> [131] 1.165185e+02 6.325233e+03 2.779068e+02 1.338893e+02 2.473047e+02
#> [136] 9.550337e+01 7.909660e+01 8.284204e+01 7.616152e+01 8.142573e+01
#> [141] 7.012396e+01 3.220422e+02 1.482398e+04 1.302612e+04 2.705962e+03
#> [146] 4.498257e+02 2.400582e+02 7.883186e+01 9.415735e+01 8.476993e+01
#> [151] 9.586769e+01 1.342671e+02 2.428195e+02 1.043032e+02 9.191663e+01
#> [156] 1.618028e+02 9.407531e+01 7.033120e+01 6.275770e+01 7.092376e+01
#> [161] 1.077026e+02 9.992572e+01 7.195176e+01 8.009655e+01 6.528124e+01
#> [166] 9.344058e+01 3.323202e+02 1.420167e+02 9.281441e+01 2.098485e+02
#> [171] 1.002633e+02 6.756606e+01 2.809604e+02 7.228959e+01 6.968629e+01
#> [176] 9.101395e+01 1.395833e+02 9.721237e+01 7.501189e+01 7.158226e+01
#> [181] 8.586032e+01 7.695245e+01 3.077413e+02 6.788226e+01 7.096221e+01
#> [186] 6.367785e+01 1.113472e+02 1.225126e+02 9.411350e+01 8.190180e+01
#> [191] 1.020243e+02 1.635554e+02 2.144411e+02 3.450576e+02 8.116054e+01
#> [196] 1.115072e+02 1.165407e+02 8.498052e+01 7.290496e+01 1.472063e+02
#> [201] 1.387891e+02 2.007359e+03 6.469024e+02 1.058006e+02 1.254654e+02
#> [206] 2.190912e+02 5.722031e+02 1.483527e+02 8.678884e+01 6.989619e+01
#> [211] 9.062710e+01 2.614686e+02 9.287234e+01 7.496456e+01 1.161884e+02
#> [216] 7.423064e+01 1.165951e+02 6.371086e+01 7.307784e+01 9.840517e+01
#> [221] 2.030378e+02 7.099121e+01 1.031698e+02 7.607130e+01 1.167750e+02
#> [226] 5.877829e+01 9.125737e+02 9.747197e+01 9.764530e+01 8.876434e+02
#> [231] 7.572428e+01 1.088331e+02 8.135553e+01 8.405619e+01 1.253903e+02
#> [236] 7.557025e+01 1.643650e+02 6.774280e+01 7.001739e+01 1.237564e+02
#> [241] 9.879883e+01 1.021201e+02 6.770850e+01 7.407286e+01 1.030683e+02
#> [246] 4.726377e+02 1.130083e+03 3.352299e+02 8.248158e+01 1.051860e+02
#> [251] 7.529236e+01 7.331700e+01 1.101868e+02 8.463458e+01 7.997818e+01
#> [256] 7.647482e+01 9.877741e+01 1.073933e+02 7.328267e+01 8.947516e+01
#> [261] 7.608102e+01 1.004780e+02 7.771878e+01 6.537970e+02 6.777790e+01
#> [266] 2.829989e+02 2.748873e+02 7.643109e+01 3.329455e+02 1.572274e+02
#> [271] 8.249761e+01 6.980857e+01 7.508492e+01 2.356970e+02 8.144363e+02
#> [276] 1.830636e+06 6.446863e+01 8.397226e+01 7.999670e+01 8.153226e+01
#> [281] 8.779979e+01 6.440607e+01 8.127815e+01 2.429653e+02 9.586368e+01
#> [286] 2.188885e+02 3.723690e+02 1.042997e+02 9.474818e+01 9.828620e+01
#> [291] 1.377086e+02 3.826779e+02 1.306416e+02 9.511597e+01 8.650451e+01
#> [296] 6.944430e+01 7.200725e+01 9.500559e+01 6.273949e+01 9.239607e+01
#> [301] 1.158652e+02 1.878246e+02 2.381370e+02 6.755691e+01 1.647752e+07
#> [306] 5.920343e+02 7.062549e+01 9.914575e+01 1.081339e+02 1.154179e+02
#> [311] 1.227293e+02 1.092395e+02 1.036553e+02 6.291665e+01 6.359109e+01
#> [316] 6.388081e+01 6.875621e+01 6.697893e+01 7.006084e+01 7.084952e+01
#> [321] 8.539177e+01 8.279920e+01 1.624677e+02 1.290479e+02 6.031483e+02
#> [326] 1.059807e+02 7.714391e+01 3.787980e+02 7.765923e+01 8.765855e+01
#> [331] 5.894790e+01 1.128577e+02 6.831263e+01 8.601111e+01 6.599475e+01
#> [336] 6.465211e+01 7.330603e+01 7.422697e+01 1.354445e+02 7.472874e+01
#> [341] 7.316252e+01 1.011210e+02 1.227006e+02 7.824218e+01 9.709316e+01
#> [346] 1.127772e+02 7.344941e+01 8.188705e+01 6.498615e+01 7.564058e+01
#> [351] 1.471231e+02 1.481020e+02 1.606154e+02 1.147554e+03 6.759858e+01
#> [356] 6.832623e+01 6.289997e+01 7.076672e+01 6.899247e+01 2.280385e+08
#> [361] 6.467602e+01 7.740988e+01 1.097820e+02 7.521149e+01 7.810817e+01
#> [366] 6.340488e+01 1.078457e+02 6.626721e+01 8.789202e+01 6.345872e+01
#> [371] 6.427737e+01 7.027655e+01 6.865506e+01 6.681016e+01 6.351972e+01
#> [376] 7.538738e+01 6.739541e+01 8.697141e+05 9.614196e+01 1.304254e+03
#> [381] 1.442735e+02 8.228819e+01 7.027436e+01 7.107649e+01 7.937943e+01
#> [386] 1.275111e+02 7.648712e+01 1.024568e+02 8.084786e+01 9.376636e+01
#> [391] 1.139291e+02 6.125777e+02 7.842591e+01 6.433062e+01 2.065951e+03
#> [396] 3.830488e+02 1.912654e+02 8.675184e+01 1.231133e+02 1.557378e+02
#> [401] 6.498353e+01 3.686842e+03 1.226176e+02 1.181014e+02 6.506921e+01
#> [406] 3.028460e+02 2.936731e+02 6.611902e+01 8.490061e+01 1.262455e+02
#> [411] 3.130016e+02 3.072423e+02 1.242236e+02 1.461067e+02 7.334969e+01
#> [416] 6.686789e+01 6.495824e+01 2.320670e+02 8.192262e+01 8.947807e+01
#> [421] 6.553923e+01 8.471659e+01 1.616762e+02 2.013976e+02 1.053494e+02
#> [426] 5.700088e+02 1.626263e+02 6.569344e+01 9.262651e+01 9.530898e+01
#> [431] 7.832403e+01 3.415184e+02 6.740739e+01 1.665372e+02 6.560833e+01
#> [436] 1.069004e+02 7.476624e+01 6.307194e+01 8.902951e+01 6.404579e+01
#> [441] 6.834968e+01 6.311129e+01 8.091682e+01 6.938984e+01 8.718388e+01
#> [446] 7.435309e+01 7.385399e+01 1.082591e+02 9.890265e+01 6.734576e+01
#> [451] 1.485000e+02 6.938974e+01 8.942225e+01 1.142122e+02 1.094485e+02
#> [456] 1.671633e+02 8.882374e+01 9.511879e+01 7.305982e+01 7.981336e+01
#> [461] 1.453291e+02 2.870237e+04 2.050652e+02 7.103932e+01 8.766903e+01
#> [466] 9.422777e+01 7.299921e+01 5.985045e+02 7.215169e+01 7.512364e+01
#> [471] 7.753989e+01 6.608052e+01 1.244225e+02 1.439766e+02 7.970548e+01
#> [476] 7.974252e+01 6.522669e+01 8.841486e+01 6.863561e+01 9.771016e+01
#> [481] 1.540526e+02 7.415581e+01 1.045636e+02 1.220220e+02 6.620502e+01
#> [486] 7.790322e+01 1.123000e+02 9.602169e+01 7.898949e+01 1.329220e+02
#> [491] 1.458263e+02 6.631644e+01 1.717663e+02 7.936887e+01 7.563687e+01
#> [496] 1.317247e+02 9.685664e+01 6.395396e+01 6.729328e+01 8.226933e+01
#> [501] 9.583051e+01 8.123179e+01 6.603096e+01 8.166833e+01 9.534416e+01
#> [506] 7.479518e+01 8.941666e+01 2.652630e+03 6.853687e+02 7.603995e+01
#> [511] 7.182443e+01 7.351562e+01 6.420767e+01 6.540612e+01 6.814163e+01
#> [516] 7.118934e+01 6.705539e+01 6.783003e+01 1.020870e+02 8.049179e+01
#> [521] 6.336036e+01 7.552180e+01 7.524355e+01 1.071391e+03 7.692908e+02
#> [526] 7.582135e+01 6.215677e+01 6.627066e+01 6.369027e+01 6.366029e+01
#> [531] 6.841149e+01 6.395978e+01 1.055324e+02 8.762721e+01 7.436197e+01
#> [536] 6.892574e+01 1.714359e+02 7.375749e+01 1.734593e+02 1.085592e+02
#> [541] 9.332053e+01 8.455645e+01 7.162417e+01 6.845966e+01 3.398956e+02
#> [546] 6.689121e+01 1.462649e+02 6.888376e+01 7.828973e+01 6.394159e+01
#> [551] 7.210700e+01 7.442499e+01 2.131563e+02 6.249770e+01 5.221885e+02
#> [556] 3.211450e+02 1.085235e+02 9.683881e+01 2.058931e+02 3.214072e+02
#> [561] 8.239223e+01 8.857543e+01 7.193743e+01 8.851649e+01 7.818673e+01
#> [566] 7.469181e+01 7.534703e+01 7.694348e+01 7.815842e+01 1.349556e+02
#> [571] 8.204498e+01 7.012355e+01 6.616020e+01 9.112752e+01 1.223094e+02
#> [576] 8.405050e+01 8.960626e+01 3.397330e+02 7.830416e+01 8.931703e+01
#> [581] 7.487886e+01 1.085070e+02 7.870136e+01 9.454342e+01 7.102361e+01
#> [586] 1.113638e+02 7.547809e+01 7.952469e+01 1.394289e+02 7.002473e+01
#> [591] 7.631599e+01 7.679162e+01 6.515385e+01 6.063544e+02 2.718746e+02
#> [596] 9.536368e+01 1.721589e+02 1.505850e+03 1.365445e+02 1.217510e+02
#> [601] 6.543588e+01 9.109663e+01 9.083918e+01 7.536294e+01 8.046477e+01
#> [606] 1.281005e+02 1.079334e+02 3.352045e+02 4.533984e+02 7.091470e+01
#> [611] 7.236715e+01 7.468380e+01 3.062332e+02 1.184928e+02 9.918626e+01
#> [616] 1.093950e+02 1.041965e+02 1.961782e+02 1.042989e+02 8.969589e+04
#> [621] 7.078192e+05 7.097854e+01 7.751499e+01 6.892130e+01 1.108707e+02
#> [626] 2.232476e+02 2.006620e+02 7.932013e+01 1.159676e+03 1.883037e+02
#> [631] 6.726310e+01 9.542041e+01 2.443368e+02 1.043137e+03 2.983515e+02
#> [636] 7.777657e+01 8.070879e+01 6.403229e+01 8.888162e+01 8.538830e+01
#> [641] 1.247818e+02 1.399037e+02 1.089952e+02 6.629702e+01 1.080189e+02
#> [646] 1.132571e+02 1.005687e+02 1.860959e+02 1.208713e+02 9.258811e+01
#> [651] 6.594061e+01 9.379116e+02 7.697684e+01 6.894358e+01 6.191327e+01
#> [656] 8.313893e+01 1.092612e+02 6.995375e+01 7.133640e+01 7.197575e+01
#> [661] 7.672857e+01 8.499236e+01 8.685402e+01 6.048245e+01 6.357940e+01
#> [666] 6.202909e+01 8.940386e+01 7.570720e+01 6.631564e+01 8.376930e+01
#> [671] 6.180111e+01 6.758016e+01 1.077578e+02 7.898430e+01 6.296639e+01
#> [676] 7.440182e+01 7.834852e+01 1.481117e+02 2.367684e+02 7.299526e+01
#> [681] 9.040394e+01 1.058387e+02 1.191788e+02 6.549132e+01 7.424083e+01
#> [686] 7.172023e+01 6.408927e+01 6.243297e+01 1.329574e+02 8.056273e+01
#> [691] 1.197647e+02 1.530409e+02 1.038433e+02 8.166770e+01 1.253467e+02
#> [696] 1.220505e+02 3.236457e+03 2.599411e+02 6.745421e+01 3.414145e+04
#> [701] 9.280191e+01 9.388428e+01 7.385027e+01 1.040451e+02 6.772765e+01
#> [706] 7.069772e+01 2.458752e+02 1.037016e+02 8.885122e+01 2.029264e+02
#> [711] 6.428868e+01 1.641644e+02 2.253427e+02 8.709608e+01 5.984575e+01
#> [716] 6.374235e+01 6.808142e+01 1.657057e+02 6.482579e+01 7.508153e+01
#> [721] 7.215192e+01 6.393913e+01 1.023408e+02 1.012799e+02 1.246391e+02
#> [726] 1.358544e+02 7.305769e+01 1.510989e+02 7.086984e+01 9.640582e+01
#> [731] 1.114223e+02 1.151677e+02 5.424175e+02 8.139288e+01 6.498051e+01
#> [736] 1.853839e+02 8.735235e+01 6.151357e+01 1.184437e+02 7.846750e+01
#> [741] 1.888292e+02 6.669805e+01 1.671514e+02 7.130596e+01 7.488092e+01
#> [746] 1.228597e+02 1.786718e+02 7.530245e+01 6.905881e+01 7.231150e+01
#> [751] 7.652166e+01 2.609376e+02 3.756595e+02 1.111810e+02 9.367800e+01
#> [756] 7.708475e+01 1.191537e+02 7.292673e+01 7.098604e+01 8.090175e+01
#> [761] 7.226130e+01 9.768445e+01 1.920278e+02 9.518280e+01 8.027478e+01
#> [766] 6.360850e+01 7.142285e+01 6.445069e+01 6.277446e+01 6.845923e+01
#> [771] 8.710731e+01 6.911309e+01 8.537515e+01 6.524240e+01 6.303661e+01
#> [776] 4.255921e+02 1.956847e+02 1.425612e+02 9.839049e+01 4.267303e+02
#> [781] 6.413809e+03 3.770298e+03 7.420378e+01 7.215922e+01 7.769070e+01
#> [786] 8.641536e+01 8.939648e+01 6.336229e+01 7.469905e+01 7.917920e+01
#> [791] 1.113769e+02 1.018322e+03 7.373018e+01 9.776016e+01 1.276342e+02
#> [796] 8.638706e+01 6.503516e+01 6.795793e+01 3.053523e+02 6.835251e+01
#> [801] 9.122353e+01 8.763866e+03 3.447374e+02 1.125528e+02 1.388075e+02
#> [806] 7.587302e+01 7.117863e+01 1.978684e+02 5.393260e+02 1.834707e+02
#> [811] 1.642752e+03 2.662979e+02 7.812013e+01 7.228698e+01 1.010694e+02
#> [816] 1.026655e+02 3.499768e+02 1.263577e+02 7.415225e+01 5.918865e+02
#> [821] 6.309299e+01 7.211715e+01 6.198020e+01 7.342231e+01 3.566684e+02
#> [826] 1.122463e+02 6.955022e+01 1.216275e+02 1.308598e+03 6.969307e+01
#> [831] 6.436507e+01 1.532354e+02 7.243903e+01 6.625231e+01 7.756719e+01
#> [836] 8.507288e+01 7.887272e+01 9.122139e+01 7.294194e+01 6.874402e+01
#> [841] 6.930561e+01 8.889144e+01 6.864367e+01 3.611316e+03 8.019284e+01
#> [846] 6.901827e+01 1.072320e+02 7.209109e+01 7.543683e+01 8.090216e+01
#> [851] 6.740235e+01 6.618130e+01 7.228172e+01 7.250018e+01 6.306170e+01
#> [856] 6.276427e+01 8.831652e+01 7.107250e+01 1.897550e+02 9.905587e+01
#> [861] 6.910906e+01 1.393682e+02 4.213533e+02 1.146287e+02 8.047136e+01
#> [866] 6.667827e+01 7.364462e+01 7.154574e+01 1.010718e+02 2.111903e+02
#> [871] 8.719962e+01 8.130098e+01 6.143689e+01 7.053530e+01 1.015365e+02
#> [876] 1.096776e+02 7.505575e+01 2.010695e+03 4.764197e+02 1.230764e+02
#> [881] 7.924228e+01 6.964553e+01 3.469246e+02 1.129306e+04 5.283912e+02
#> [886] 8.470742e+01 7.963693e+01 7.432975e+01 8.174614e+01 3.558594e+02
#> [891] 2.026117e+02 2.860567e+02 6.240347e+01 8.912034e+01 6.912312e+01
#> [896] 6.516002e+01 9.790650e+01 1.524927e+02 1.454723e+02 7.285340e+01
#> [901] 1.072139e+02 9.152749e+01 6.425220e+01 7.396011e+01 7.038702e+01
#> [906] 1.378688e+02 1.117388e+02 2.376573e+02 5.516574e+03 6.303718e+03
#> [911] 7.306213e+01 7.725080e+01 7.292899e+01 9.077480e+01 8.946691e+01
#> [916] 1.154797e+02 1.344876e+02 6.695959e+01 8.432211e+01 2.546537e+02
#> [921] 1.924867e+02 1.124726e+02 1.056092e+02 4.507556e+02 1.248037e+02
#> [926] 7.455828e+01 9.101557e+01 7.360711e+01 6.428359e+01 7.415501e+01
#> [931] 2.525721e+02 7.930411e+01 2.088275e+02 3.209373e+02 9.224608e+01
#> [936] 6.118172e+01 7.567622e+01 9.554154e+01 1.540784e+03 9.768333e+01
#> [941] 8.712167e+01 1.396945e+02 8.074237e+01 1.809050e+02 3.740310e+02
#> [946] 3.272693e+02 6.846996e+01 7.114949e+01 6.817119e+01 7.515374e+01
#> [951] 6.831056e+01 6.729589e+01 6.528424e+01 6.527961e+01 6.721524e+01
#> [956] 7.301140e+01 9.578640e+01 1.097611e+02 7.088212e+01 7.289435e+01
#> [961] 1.466408e+02 1.821299e+02 7.578862e+01 7.266793e+01 7.700324e+01
#> [966] 1.005007e+02 8.330161e+01 1.594924e+02 9.464745e+01 7.622715e+01
#> [971] 6.245023e+01 2.535291e+02 7.485090e+01 6.303593e+01 7.058788e+01
#> [976] 7.260208e+01 7.962407e+02 6.249805e+01 6.432167e+01 7.391316e+01
#> [981] 8.779484e+01 1.300962e+02 8.051008e+01 9.094027e+01 7.560507e+01
#> [986] 7.713010e+01 4.096403e+03 1.143959e+02 8.726464e+01 9.503593e+01
#> [991] 1.601539e+02 8.468683e+01 7.088685e+01 1.783087e+02 6.450481e+01
#> [996] 1.464373e+02 4.917125e+02 9.098289e+01 7.431634e+01 9.711856e+01
doseFunction(
model,
alpha2 = samples@data$alpha2,
beta = samples@data$beta,
grade = 2L
)(x = 0.25)
#> [1] 3.326068e+03 8.689574e+01 5.738979e+01 5.745355e+01 6.133776e+01
#> [6] 7.665158e+01 5.963980e+01 6.122455e+01 5.357155e+01 5.781265e+01
#> [11] 5.754635e+01 5.783688e+01 6.592526e+01 5.517784e+01 2.572881e+03
#> [16] 6.145467e+01 1.245207e+02 7.544382e+01 8.536159e+01 6.102780e+01
#> [21] 8.333319e+01 4.903585e+01 4.939718e+01 7.093240e+01 1.109160e+02
#> [26] 7.750067e+01 1.092889e+02 6.286564e+01 1.358358e+02 6.923076e+01
#> [31] 5.704840e+01 1.812908e+02 2.106545e+02 9.631713e+01 7.238544e+01
#> [36] 7.161385e+01 7.093204e+01 5.684382e+01 6.967699e+01 6.256409e+01
#> [41] 8.330844e+01 7.224826e+01 5.836196e+01 8.060180e+01 1.752463e+03
#> [46] 2.306543e+02 2.239749e+06 2.938212e+02 2.006710e+02 8.991320e+01
#> [51] 4.330831e+02 6.227873e+01 6.442332e+01 1.267931e+02 6.323095e+01
#> [56] 6.969568e+01 7.546278e+01 5.826588e+01 7.984885e+01 6.995754e+01
#> [61] 1.414124e+02 9.870623e+01 5.821526e+01 8.577280e+01 5.688042e+01
#> [66] 5.516091e+01 6.751767e+01 6.817521e+01 5.821865e+01 6.677147e+01
#> [71] 5.826930e+01 8.508404e+01 6.261571e+01 5.608686e+01 8.392292e+01
#> [76] 5.469850e+02 7.217584e+01 7.140792e+01 7.119037e+01 5.810246e+01
#> [81] 7.757801e+01 8.847147e+01 5.982239e+01 6.758833e+01 6.253872e+01
#> [86] 8.056877e+01 1.456553e+02 1.455999e+02 2.615592e+02 5.839864e+01
#> [91] 9.771553e+01 5.485931e+01 6.102589e+01 5.754851e+01 5.753442e+01
#> [96] 6.313779e+01 6.305522e+01 5.982269e+01 6.639310e+01 8.115673e+01
#> [101] 6.182993e+01 6.679523e+01 6.474195e+01 1.328874e+02 6.821791e+01
#> [106] 8.987904e+02 8.909541e+01 5.705202e+01 6.564566e+01 5.963010e+01
#> [111] 7.444350e+01 6.026915e+01 6.110965e+01 9.333618e+01 6.430912e+01
#> [116] 7.430710e+01 1.492377e+02 3.544201e+07 8.018933e+01 9.853304e+01
#> [121] 2.913133e+02 1.403177e+02 8.469539e+01 5.135918e+01 6.588585e+01
#> [126] 1.029623e+02 5.924486e+01 6.702149e+01 6.802721e+01 5.483772e+01
#> [131] 8.131391e+01 3.399250e+02 1.298091e+02 3.790825e+01 7.141046e+01
#> [136] 5.706397e+01 6.160783e+01 5.966245e+01 6.405967e+01 6.570378e+01
#> [141] 5.760955e+01 1.142786e+02 5.054755e+02 2.471068e+02 2.506299e+02
#> [146] 1.442829e+02 8.239291e+01 5.713013e+01 6.939088e+01 6.825895e+01
#> [151] 6.131906e+01 8.099906e+01 1.047359e+02 7.146355e+01 7.233008e+01
#> [156] 8.173221e+01 6.764028e+01 6.062751e+01 5.555863e+01 6.198217e+01
#> [161] 7.297516e+01 7.812162e+01 5.720035e+01 5.961841e+01 5.565993e+01
#> [166] 7.233721e+01 1.772452e+02 1.006419e+02 5.864602e+01 8.747927e+01
#> [171] 7.187308e+01 5.862663e+01 1.076023e+02 6.173395e+01 5.741562e+01
#> [176] 6.591952e+01 8.699453e+01 7.217618e+01 5.141689e+01 6.385798e+01
#> [181] 7.032578e+01 6.572853e+01 1.013966e+02 5.577007e+01 5.995933e+01
#> [186] 5.599402e+01 7.019028e+01 7.997825e+01 7.323318e+01 6.041394e+01
#> [191] 6.544836e+01 7.782087e+01 9.501940e+01 1.287018e+02 5.677770e+01
#> [196] 7.787891e+01 8.093025e+01 6.055683e+01 5.808126e+01 8.280956e+01
#> [201] 9.384905e+01 4.697290e+02 1.836870e+02 6.480092e+01 7.996989e+01
#> [206] 9.360591e+01 1.596381e+02 8.662077e+01 6.512850e+01 6.146271e+01
#> [211] 6.571255e+01 1.179396e+02 6.534958e+01 5.576694e+01 8.184481e+01
#> [216] 6.315499e+01 8.733513e+01 5.630238e+01 6.128278e+01 7.814631e+01
#> [221] 1.248898e+02 5.635317e+01 6.658296e+01 6.218628e+01 9.483958e+01
#> [226] 5.409081e+01 3.219369e+02 7.181792e+01 6.066047e+01 1.411903e+02
#> [231] 5.677942e+01 7.938855e+01 7.078819e+01 7.117180e+01 8.453161e+01
#> [236] 5.604901e+01 7.145111e+01 5.834623e+01 6.027975e+01 6.530039e+01
#> [241] 6.105682e+01 7.865382e+01 5.727155e+01 6.031480e+01 7.327273e+01
#> [246] 9.311674e+01 1.331279e+02 9.933112e+01 5.981795e+01 7.520236e+01
#> [251] 6.054898e+01 5.741720e+01 7.028683e+01 6.595634e+01 6.136064e+01
#> [256] 6.310972e+01 7.659808e+01 8.198580e+01 6.388752e+01 6.206344e+01
#> [261] 6.255768e+01 6.503342e+01 6.701371e+01 5.008721e+01 5.775356e+01
#> [266] 9.713564e+01 1.116067e+02 5.880344e+01 1.185391e+02 1.036320e+02
#> [271] 6.222184e+01 6.203844e+01 6.494445e+01 8.681241e+01 1.285260e+02
#> [276] 1.611569e+04 5.834000e+01 6.654066e+01 6.662860e+01 6.170055e+01
#> [281] 5.591627e+01 5.780996e+01 6.335565e+01 1.143059e+02 6.213440e+01
#> [286] 1.010717e+02 1.320433e+02 6.487739e+01 6.262107e+01 7.779153e+01
#> [291] 7.307012e+01 1.287274e+02 8.639943e+01 7.112479e+01 6.860058e+01
#> [296] 5.629886e+01 6.297740e+01 7.449239e+01 5.628328e+01 7.308740e+01
#> [301] 8.299837e+01 6.979097e+01 8.450201e+01 5.566373e+01 1.282063e+04
#> [306] 1.416082e+02 6.239232e+01 6.183617e+01 7.175234e+01 8.011504e+01
#> [311] 8.588807e+01 6.431623e+01 6.906836e+01 5.350579e+01 5.966706e+01
#> [316] 5.983891e+01 6.249558e+01 6.016789e+01 6.386456e+01 6.285876e+01
#> [321] 6.607580e+01 6.381709e+01 7.389719e+01 7.677530e+01 1.296524e+02
#> [326] 7.316221e+01 6.284770e+01 1.338876e+02 6.315836e+01 6.233732e+01
#> [331] 5.389151e+01 8.300722e+01 5.973550e+01 6.340218e+01 5.816915e+01
#> [336] 5.541716e+01 6.095329e+01 6.210350e+01 6.926455e+01 5.797151e+01
#> [341] 6.090434e+01 7.405780e+01 6.832635e+01 6.094469e+01 7.384540e+01
#> [346] 7.881391e+01 6.299980e+01 6.129910e+01 5.342289e+01 6.397800e+01
#> [351] 6.593126e+01 8.669863e+01 7.577220e+01 2.158903e+02 5.924292e+01
#> [356] 5.935412e+01 5.588769e+01 6.134346e+01 5.938591e+01 8.095554e+05
#> [361] 5.173936e+01 6.298899e+01 7.222080e+01 6.136841e+01 5.165428e+01
#> [366] 4.976463e+01 8.032376e+01 5.751724e+01 6.566880e+01 5.563012e+01
#> [371] 5.628820e+01 5.918566e+01 6.124220e+01 5.889698e+01 5.773260e+01
#> [376] 6.610159e+01 5.860148e+01 1.462183e+03 5.319977e+01 3.105920e+02
#> [381] 7.080844e+01 6.015031e+01 5.840938e+01 6.045268e+01 6.360880e+01
#> [386] 8.364453e+01 6.292388e+01 6.967311e+01 6.232672e+01 6.952278e+01
#> [391] 7.826374e+01 1.054964e+02 6.173600e+01 5.519718e+01 9.595903e+01
#> [396] 1.298624e+02 9.140448e+01 5.896212e+01 7.854950e+01 8.703302e+01
#> [401] 5.584812e+01 6.114817e+02 7.232430e+01 6.716577e+01 5.511042e+01
#> [406] 9.521400e+01 9.657829e+01 5.775760e+01 6.979765e+01 7.853240e+01
#> [411] 8.901112e+01 9.783245e+01 8.969895e+01 8.743662e+01 5.595633e+01
#> [416] 5.732655e+01 5.815792e+01 1.151365e+02 5.794843e+01 6.928243e+01
#> [421] 5.868118e+01 7.053839e+01 1.028565e+02 1.136351e+02 6.694370e+01
#> [426] 2.241597e+02 6.422838e+01 5.562033e+01 6.500243e+01 6.504918e+01
#> [431] 6.769583e+01 1.110836e+02 5.726780e+01 8.296974e+01 5.786153e+01
#> [436] 7.915928e+01 6.372838e+01 5.650291e+01 7.364965e+01 5.663012e+01
#> [441] 5.973355e+01 5.617873e+01 6.427448e+01 6.022021e+01 6.804277e+01
#> [446] 6.205809e+01 6.407502e+01 8.292507e+01 7.249879e+01 5.649771e+01
#> [451] 9.207264e+01 5.676682e+01 7.127809e+01 7.638262e+01 5.372200e+01
#> [456] 7.792071e+01 5.877612e+01 7.333861e+01 6.196786e+01 5.851068e+01
#> [461] 7.585811e+01 3.520852e+02 9.528817e+01 5.730095e+01 6.787374e+01
#> [466] 6.972242e+01 5.947596e+01 2.347543e+02 6.016511e+01 6.126380e+01
#> [471] 5.939919e+01 5.836965e+01 6.732076e+01 8.268968e+01 6.392372e+01
#> [476] 6.517217e+01 5.810783e+01 7.343178e+01 6.186279e+01 7.725541e+01
#> [481] 8.781358e+01 5.651984e+01 5.991199e+01 6.627216e+01 5.699839e+01
#> [486] 6.486247e+01 8.440707e+01 7.427769e+01 6.390685e+01 8.577185e+01
#> [491] 8.054236e+01 5.597120e+01 9.866445e+01 6.802435e+01 6.305985e+01
#> [496] 8.609351e+01 7.718342e+01 5.712246e+01 5.912258e+01 6.589637e+01
#> [501] 7.408863e+01 6.410904e+01 5.504450e+01 5.905312e+01 6.900084e+01
#> [506] 6.242377e+01 6.443977e+01 1.569908e+02 1.059974e+02 5.982093e+01
#> [511] 5.128571e+01 6.466630e+01 5.576014e+01 5.750840e+01 6.072055e+01
#> [516] 6.371798e+01 6.080483e+01 6.047417e+01 7.252623e+01 6.381078e+01
#> [521] 5.548130e+01 6.420363e+01 5.976615e+01 1.670750e+02 1.048721e+02
#> [526] 6.193241e+01 5.555494e+01 5.820047e+01 5.672451e+01 5.765174e+01
#> [531] 6.162111e+01 5.896178e+01 7.567773e+01 7.087674e+01 6.273300e+01
#> [536] 5.868401e+01 1.088536e+02 6.071672e+01 8.999997e+01 7.221984e+01
#> [541] 6.111027e+01 6.060747e+01 5.869948e+01 5.760989e+01 1.650412e+02
#> [546] 5.629464e+01 8.625241e+01 4.828178e+01 6.319842e+01 5.751873e+01
#> [551] 6.072938e+01 6.529165e+01 1.157024e+02 5.741447e+01 1.103868e+02
#> [556] 1.170334e+02 7.113333e+01 7.942601e+01 1.083817e+02 1.788806e+02
#> [561] 6.991320e+01 6.449929e+01 6.345065e+01 7.215893e+01 6.470696e+01
#> [566] 5.676961e+01 6.373398e+01 6.344087e+01 6.761562e+01 5.492745e+01
#> [571] 6.199420e+01 6.170279e+01 5.790703e+01 7.140964e+01 7.699444e+01
#> [576] 6.921227e+01 6.894288e+01 1.354757e+02 5.719327e+01 5.713640e+01
#> [581] 5.301407e+01 7.942515e+01 5.999562e+01 7.143871e+01 5.890882e+01
#> [586] 8.098626e+01 6.495603e+01 5.219688e+01 9.776994e+01 5.592681e+01
#> [591] 5.582470e+01 6.134853e+01 5.828395e+01 1.925614e+02 1.719362e+02
#> [596] 8.107999e+01 1.265000e+02 5.689855e+02 9.865074e+01 8.883959e+01
#> [601] 5.789532e+01 6.870132e+01 5.869721e+01 6.067833e+01 5.961409e+01
#> [606] 7.529863e+01 7.102736e+01 1.026593e+02 9.328064e+01 6.318849e+01
#> [611] 5.718823e+01 5.853751e+01 1.073772e+02 7.804495e+01 6.038157e+01
#> [616] 6.164477e+01 7.207352e+01 8.986288e+01 5.694010e+01 6.081123e+03
#> [621] 2.672484e+04 5.434573e+01 6.173388e+01 5.777153e+01 7.343878e+01
#> [626] 9.573433e+01 7.202090e+01 6.246688e+01 2.377217e+02 9.125608e+01
#> [631] 5.294113e+01 7.616766e+01 1.250071e+02 3.275706e+02 1.461289e+02
#> [636] 6.532337e+01 6.507898e+01 5.509298e+01 6.800455e+01 5.315737e+01
#> [641] 8.221035e+01 8.165948e+01 8.319403e+01 5.980723e+01 6.146091e+01
#> [646] 5.919537e+01 7.361225e+01 9.905614e+01 8.948780e+01 7.521572e+01
#> [651] 5.916435e+01 2.734806e+02 6.410985e+01 5.914460e+01 5.546756e+01
#> [656] 6.922141e+01 8.228981e+01 6.052121e+01 6.163465e+01 5.569918e+01
#> [661] 6.417501e+01 7.050764e+01 7.037056e+01 5.602408e+01 5.765065e+01
#> [666] 5.557140e+01 7.204310e+01 6.479864e+01 5.779530e+01 6.871418e+01
#> [671] 5.606956e+01 5.780809e+01 7.267804e+01 6.589843e+01 5.792434e+01
#> [676] 6.325835e+01 6.731484e+01 7.700888e+01 9.814127e+01 5.593448e+01
#> [681] 5.427161e+01 7.482230e+01 8.077568e+01 5.703643e+01 6.461621e+01
#> [686] 6.119038e+01 5.797026e+01 5.668396e+01 8.157438e+01 6.338748e+01
#> [691] 6.560771e+01 8.135137e+01 6.951573e+01 5.822465e+01 6.483246e+01
#> [696] 4.899744e+01 4.990645e+02 1.368719e+02 5.752302e+01 7.229764e+02
#> [701] 5.354077e+01 6.231494e+01 5.409423e+01 7.516283e+01 5.286434e+01
#> [706] 6.180574e+01 1.098597e+02 6.581498e+01 6.397875e+01 1.050614e+02
#> [711] 5.671290e+01 1.164550e+02 9.836061e+01 5.834369e+01 5.372669e+01
#> [716] 5.797957e+01 6.235852e+01 8.884896e+01 5.515446e+01 6.407031e+01
#> [721] 5.831878e+01 5.717913e+01 7.953518e+01 6.249503e+01 7.311136e+01
#> [726] 9.721996e+01 6.436719e+01 9.811971e+01 5.795483e+01 6.588695e+01
#> [731] 6.774384e+01 6.358932e+01 9.777595e+01 6.575215e+01 5.784146e+01
#> [736] 1.073274e+02 6.161947e+01 5.702895e+01 8.390822e+01 5.587246e+01
#> [741] 4.936154e+01 5.692004e+01 9.831820e+01 5.492021e+01 5.983924e+01
#> [746] 7.717665e+01 7.528805e+01 6.113892e+01 6.364355e+01 6.065931e+01
#> [751] 5.936294e+01 1.419307e+02 1.267910e+02 7.998564e+01 6.017358e+01
#> [756] 6.694718e+01 7.533295e+01 5.740655e+01 5.984968e+01 6.661903e+01
#> [761] 6.284542e+01 7.703571e+01 9.303535e+01 6.616285e+01 5.932576e+01
#> [766] 5.429731e+01 6.468890e+01 6.052068e+01 5.921090e+01 6.129513e+01
#> [771] 6.856310e+01 5.545237e+01 6.857595e+01 5.751209e+01 5.790405e+01
#> [776] 1.866840e+02 1.151287e+02 8.062475e+01 5.901535e+01 1.158154e+02
#> [781] 1.348507e+03 9.345415e+02 5.919588e+01 6.290790e+01 6.241786e+01
#> [786] 6.515739e+01 7.041450e+01 5.704245e+01 6.041994e+01 6.383386e+01
#> [791] 6.783052e+01 1.767941e+02 6.105619e+01 6.157850e+01 8.572432e+01
#> [796] 6.610467e+01 5.838076e+01 6.100275e+01 1.356351e+02 5.918639e+01
#> [801] 6.324240e+01 7.152503e+02 2.078644e+02 8.241754e+01 6.987751e+01
#> [806] 5.684536e+01 5.924408e+01 9.884786e+01 1.601229e+02 1.140859e+02
#> [811] 2.972598e+02 1.030804e+02 6.020393e+01 6.117296e+01 6.615539e+01
#> [816] 7.956736e+01 1.452727e+02 8.571953e+01 5.296427e+01 1.100919e+02
#> [821] 5.304164e+01 6.507674e+01 5.637719e+01 6.461755e+01 1.102533e+02
#> [826] 7.044239e+01 6.040601e+01 9.193635e+01 4.014692e+02 6.161170e+01
#> [831] 5.825315e+01 9.833424e+01 6.259565e+01 5.799878e+01 6.420102e+01
#> [836] 6.339223e+01 5.954208e+01 6.508400e+01 6.023263e+01 5.940993e+01
#> [841] 5.793901e+01 7.278850e+01 5.509427e+01 8.203677e+02 6.571138e+01
#> [846] 6.116840e+01 7.869053e+01 6.318044e+01 6.292358e+01 6.140627e+01
#> [851] 5.878804e+01 6.086314e+01 6.517163e+01 6.517257e+01 5.876082e+01
#> [856] 5.700151e+01 6.467404e+01 5.991833e+01 1.000759e+02 4.499068e+01
#> [861] 5.093554e+01 8.589470e+01 7.790802e+01 6.432596e+01 6.361459e+01
#> [866] 5.847177e+01 6.119116e+01 6.223772e+01 6.602223e+01 9.110682e+01
#> [871] 6.074605e+01 6.306553e+01 5.548498e+01 5.722982e+01 6.988671e+01
#> [876] 7.376517e+01 5.927068e+01 1.986438e+02 1.498318e+02 8.785360e+01
#> [881] 4.840251e+01 5.091576e+01 4.217802e+01 1.432927e+03 1.472376e+02
#> [886] 6.843964e+01 6.114480e+01 6.599398e+01 6.964106e+01 1.197920e+02
#> [891] 7.817386e+01 1.686437e+02 5.757097e+01 7.546698e+01 6.032653e+01
#> [896] 5.825968e+01 7.794619e+01 9.891694e+01 8.981153e+01 6.078933e+01
#> [901] 7.786282e+01 6.428728e+01 5.701945e+01 6.007559e+01 6.210923e+01
#> [906] 8.186686e+01 6.282994e+01 6.482669e+01 4.466501e+02 7.010108e+02
#> [911] 5.507889e+01 6.036548e+01 6.116921e+01 6.139383e+01 6.035606e+01
#> [916] 7.955232e+01 7.798700e+01 5.522490e+01 6.455534e+01 1.195475e+02
#> [921] 1.026484e+02 7.277187e+01 8.101304e+01 2.266926e+02 8.876736e+01
#> [926] 5.506969e+01 7.014742e+01 6.018653e+01 5.828254e+01 6.206229e+01
#> [931] 1.095573e+02 6.204358e+01 9.956227e+01 1.587706e+02 7.215842e+01
#> [936] 5.562925e+01 6.360373e+01 7.392358e+01 4.127161e+02 6.823474e+01
#> [941] 6.596732e+01 7.265823e+01 5.146095e+01 8.420060e+01 1.539283e+02
#> [946] 1.815373e+02 5.461980e+01 6.048157e+01 5.730622e+01 5.951630e+01
#> [951] 5.743253e+01 5.677938e+01 5.723026e+01 5.927441e+01 5.796557e+01
#> [956] 6.346922e+01 6.835851e+01 6.678093e+01 5.304158e+01 5.898211e+01
#> [961] 7.093879e+01 9.817586e+01 5.934210e+01 5.922897e+01 6.005375e+01
#> [966] 6.157553e+01 6.873406e+01 9.049453e+01 6.899595e+01 5.851136e+01
#> [971] 5.413877e+01 1.079205e+02 5.445946e+01 5.556303e+01 6.183820e+01
#> [976] 6.483781e+01 2.198868e+02 5.567938e+01 5.672187e+01 5.729337e+01
#> [981] 7.012645e+01 8.201787e+01 5.359570e+01 6.632177e+01 6.335573e+01
#> [986] 6.449737e+01 4.747765e+02 8.050749e+01 6.893047e+01 7.046261e+01
#> [991] 8.125179e+01 5.114394e+01 5.731434e+01 9.999049e+01 5.788842e+01
#> [996] 8.501921e+01 1.758998e+02 6.932987e+01 5.602188e+01 7.072449e+01