![[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] 7.938356e+01 9.280500e+01 9.575981e+01 7.120207e+01 6.828746e+01
#> [6] 1.662719e+02 1.174134e+02 9.616218e+01 1.111584e+02 7.222707e+01
#> [11] 1.584912e+02 5.011943e+02 9.534363e+01 2.711630e+02 8.647937e+01
#> [16] 1.514253e+02 1.548075e+02 6.585510e+01 6.647620e+01 6.695945e+01
#> [21] 9.018683e+01 7.270942e+01 6.757715e+01 1.076267e+02 1.172573e+02
#> [26] 6.757067e+01 7.541249e+01 7.418129e+01 6.812183e+01 7.711120e+01
#> [31] 7.008104e+01 5.486061e+02 8.182787e+01 1.188729e+02 7.917431e+01
#> [36] 7.736121e+01 7.069200e+01 1.231558e+02 9.538364e+01 2.231258e+02
#> [41] 7.697597e+01 9.392205e+01 1.008399e+02 6.345903e+01 7.738378e+01
#> [46] 7.565894e+01 7.756619e+01 1.210196e+02 1.908142e+02 1.367465e+02
#> [51] 9.244816e+01 1.325759e+02 8.014804e+03 1.233411e+02 1.453260e+02
#> [56] 2.073387e+02 1.179911e+02 6.062267e+01 3.644117e+03 2.472105e+12
#> [61] 8.484446e+01 1.198999e+02 6.197902e+02 1.860897e+03 8.284705e+02
#> [66] 1.834415e+02 9.961159e+01 1.480386e+02 7.146083e+01 3.418432e+08
#> [71] 8.321075e+03 9.016177e+02 1.005722e+04 2.614859e+04 6.707012e+01
#> [76] 7.769243e+01 1.309789e+02 2.190886e+02 1.859172e+03 7.926414e+01
#> [81] 7.146857e+01 8.611062e+01 6.996398e+02 1.039822e+02 6.487624e+01
#> [86] 8.897777e+01 1.127727e+02 2.781751e+02 1.014825e+03 3.716138e+02
#> [91] 8.837325e+01 6.356739e+01 7.929188e+02 9.904055e+01 1.073354e+02
#> [96] 1.073413e+02 6.654278e+01 6.393244e+01 1.366063e+02 6.428773e+01
#> [101] 7.634977e+01 9.901327e+01 9.277344e+01 9.561475e+02 1.161823e+02
#> [106] 6.611331e+01 1.674585e+02 2.345559e+02 7.464157e+01 6.435542e+01
#> [111] 6.832181e+01 6.934556e+01 1.354139e+02 6.565337e+01 6.844597e+01
#> [116] 5.101300e+02 1.354586e+02 3.437455e+02 8.062240e+01 7.434363e+01
#> [121] 1.020001e+02 1.532563e+02 7.000262e+01 7.273181e+01 8.696890e+01
#> [126] 7.844611e+01 6.715286e+01 8.628481e+01 4.155919e+02 1.582264e+02
#> [131] 9.474528e+01 6.448843e+01 6.931013e+01 7.279706e+01 1.103836e+02
#> [136] 1.161549e+02 7.353809e+01 6.577921e+01 7.394978e+01 6.475459e+01
#> [141] 8.361631e+01 1.268783e+02 1.913751e+02 9.806621e+01 1.367899e+02
#> [146] 1.783074e+02 9.806182e+02 8.475957e+01 7.619342e+01 9.060826e+01
#> [151] 9.748383e+01 1.038391e+02 8.121305e+01 1.568408e+02 9.845318e+01
#> [156] 8.885988e+01 2.768995e+02 1.870343e+02 3.799618e+02 2.024878e+02
#> [161] 9.625297e+01 7.197084e+01 9.309669e+01 9.381134e+01 1.239915e+02
#> [166] 3.869183e+02 6.285305e+01 2.271680e+03 1.319946e+02 1.147073e+03
#> [171] 9.039303e+01 6.771348e+01 6.278941e+01 1.008181e+02 1.121452e+02
#> [176] 7.967181e+01 6.535633e+01 7.803584e+01 2.255843e+02 1.459928e+02
#> [181] 2.636986e+02 6.411207e+01 1.312930e+02 1.641781e+02 6.841938e+02
#> [186] 6.627664e+01 6.315861e+01 6.725812e+01 1.046016e+02 6.993920e+01
#> [191] 6.340268e+01 2.086251e+02 6.476927e+01 8.320569e+01 8.238909e+01
#> [196] 7.155053e+01 8.917607e+01 6.206305e+01 6.652574e+01 7.119370e+01
#> [201] 7.201309e+01 9.810868e+01 1.153135e+02 1.481779e+02 8.448667e+03
#> [206] 1.579148e+02 4.476169e+09 2.133185e+02 1.084949e+03 2.391532e+02
#> [211] 7.538916e+01 1.839576e+03 5.800960e+02 2.736086e+02 6.810668e+01
#> [216] 7.920573e+01 7.232357e+01 8.640053e+01 1.108484e+02 1.437399e+02
#> [221] 9.937977e+01 7.158076e+01 8.439039e+01 9.842934e+01 1.375785e+02
#> [226] 7.215102e+01 7.283381e+01 8.101594e+02 1.689646e+03 6.126201e+02
#> [231] 7.331015e+01 1.013621e+02 6.822293e+01 6.769485e+01 1.401457e+02
#> [236] 1.177196e+02 3.145271e+02 6.994557e+01 1.812776e+02 1.011232e+02
#> [241] 1.303855e+02 9.777868e+01 8.263740e+01 7.641638e+01 8.919421e+01
#> [246] 2.109451e+02 2.997530e+02 7.572122e+01 1.220231e+02 8.556468e+01
#> [251] 1.006859e+02 4.144079e+03 6.993956e+01 8.423861e+01 2.848982e+02
#> [256] 1.392407e+02 1.062946e+02 6.357939e+01 7.058059e+01 6.259852e+01
#> [261] 7.286834e+01 8.515936e+01 1.097182e+02 6.996925e+01 6.192267e+01
#> [266] 6.959548e+01 7.881238e+01 1.996038e+02 1.667545e+02 7.046819e+01
#> [271] 1.095017e+02 2.185821e+02 1.145330e+02 1.934218e+02 1.872387e+02
#> [276] 6.113846e+01 7.337956e+01 8.614386e+01 1.189560e+02 6.961309e+01
#> [281] 1.006936e+02 9.285344e+01 9.054758e+01 6.897012e+01 7.768786e+02
#> [286] 9.030803e+03 7.478300e+01 8.382799e+01 8.385771e+01 6.631432e+01
#> [291] 2.012342e+02 6.045359e+01 8.417418e+01 1.048996e+03 6.560971e+01
#> [296] 8.663805e+01 7.040224e+01 9.926535e+01 8.700190e+01 7.213526e+01
#> [301] 3.181655e+02 7.102286e+01 7.437387e+01 1.748685e+02 6.039052e+01
#> [306] 1.360214e+02 6.232619e+01 9.400373e+03 8.797463e+01 7.806854e+01
#> [311] 1.275430e+02 7.364544e+01 1.030646e+02 8.593872e+01 7.055373e+01
#> [316] 1.829354e+02 7.762133e+01 2.397443e+02 6.968630e+01 8.234007e+01
#> [321] 6.465610e+01 8.351943e+01 2.828372e+03 2.036477e+02 7.358991e+01
#> [326] 2.342614e+02 8.701192e+01 6.381894e+01 1.062612e+02 1.441124e+08
#> [331] 1.528788e+02 7.700240e+01 2.495412e+02 7.533459e+01 7.185995e+01
#> [336] 9.720295e+01 7.157253e+01 1.040266e+02 6.149046e+01 7.745765e+01
#> [341] 7.797970e+01 1.595783e+02 1.472828e+02 7.126596e+01 9.737070e+01
#> [346] 8.858664e+01 1.393746e+02 7.596843e+01 7.269070e+01 9.529689e+01
#> [351] 6.653735e+01 6.022279e+01 7.012528e+01 1.114026e+02 6.282169e+01
#> [356] 7.930911e+01 3.487969e+02 1.157677e+02 6.906644e+01 7.182453e+01
#> [361] 8.866728e+01 7.623240e+01 1.241924e+02 2.805389e+02 1.899479e+02
#> [366] 1.785837e+02 6.787123e+01 1.625325e+02 1.131651e+02 7.430483e+01
#> [371] 9.075020e+01 6.228683e+01 8.492206e+01 7.950428e+01 9.546058e+01
#> [376] 7.188246e+01 4.411151e+02 5.675448e+02 1.008730e+03 2.206935e+02
#> [381] 6.644583e+02 1.074591e+02 9.075959e+01 6.479304e+01 6.896726e+01
#> [386] 6.381211e+01 5.053721e+02 1.250329e+02 7.592254e+01 8.010857e+01
#> [391] 8.879941e+01 2.593163e+02 7.319344e+01 9.873828e+01 8.181642e+01
#> [396] 8.369438e+01 6.872355e+01 7.680162e+01 7.869762e+01 1.746163e+02
#> [401] 8.697626e+01 9.064164e+01 9.512344e+01 1.137213e+02 1.145106e+02
#> [406] 7.259231e+01 6.328858e+01 6.416550e+01 8.226283e+01 1.400078e+02
#> [411] 6.688591e+01 6.764526e+01 1.725731e+02 6.423693e+01 4.006822e+02
#> [416] 9.108078e+01 8.957918e+01 7.781035e+01 1.830402e+02 8.224552e+01
#> [421] 2.217669e+02 8.364128e+01 1.115176e+02 8.010760e+01 1.682898e+02
#> [426] 1.232463e+02 7.939076e+01 8.500414e+01 6.850237e+01 6.870452e+01
#> [431] 3.674266e+02 2.301760e+02 6.410473e+01 6.521290e+01 9.342455e+01
#> [436] 7.991755e+01 9.798023e+01 2.034085e+02 1.295988e+02 8.994013e+01
#> [441] 7.427812e+01 8.372559e+01 6.650725e+01 8.894196e+01 7.214471e+01
#> [446] 2.413889e+02 9.245542e+01 1.246026e+02 7.696368e+01 1.172738e+02
#> [451] 6.518962e+01 8.768995e+01 1.091758e+02 8.045738e+01 1.163151e+02
#> [456] 9.133900e+01 8.325074e+01 6.422731e+01 7.597912e+01 7.389868e+01
#> [461] 7.991449e+01 1.103641e+02 2.139346e+02 1.903914e+02 9.914173e+01
#> [466] 6.147677e+02 3.355098e+02 1.128730e+02 6.551926e+01 7.945097e+01
#> [471] 8.606111e+02 2.878272e+02 5.137956e+05 1.312263e+02 8.980968e+01
#> [476] 1.228814e+03 1.080733e+02 6.750885e+01 8.216451e+01 8.722150e+01
#> [481] 8.823342e+01 6.957774e+01 6.531950e+01 7.280979e+01 9.118471e+02
#> [486] 7.173806e+01 8.556296e+01 1.119770e+02 1.300783e+02 1.155870e+02
#> [491] 2.555082e+02 8.933226e+01 6.906827e+01 8.069714e+01 1.326950e+02
#> [496] 7.207910e+01 6.904078e+01 4.552569e+02 1.469509e+02 1.086182e+02
#> [501] 7.398000e+01 6.840537e+01 6.744847e+01 6.370424e+01 7.848055e+01
#> [506] 7.664171e+01 1.057561e+02 7.735076e+01 1.651787e+02 7.474445e+01
#> [511] 1.128658e+02 1.340915e+02 6.507446e+01 6.794967e+01 1.196263e+02
#> [516] 7.571347e+01 6.521099e+01 1.135963e+02 7.339712e+01 1.680237e+02
#> [521] 3.367368e+02 6.191011e+01 1.216486e+02 8.709575e+01 8.608195e+01
#> [526] 7.915276e+01 7.882521e+01 8.542705e+01 1.337681e+02 1.635037e+02
#> [531] 1.250305e+02 1.008941e+02 7.319906e+01 1.088903e+02 1.146834e+02
#> [536] 2.479170e+02 7.204001e+01 1.169926e+02 1.547451e+02 8.423102e+01
#> [541] 8.045672e+01 2.324351e+03 1.496840e+02 8.178565e+01 7.836285e+01
#> [546] 1.003076e+02 1.200288e+02 7.797614e+01 7.175609e+01 1.098945e+02
#> [551] 9.200185e+01 6.156500e+01 6.546448e+01 7.050973e+01 6.839012e+01
#> [556] 8.943080e+01 1.128977e+02 1.492892e+02 9.244873e+01 4.323006e+02
#> [561] 8.224909e+01 6.977957e+01 8.176617e+01 7.926616e+01 7.003673e+01
#> [566] 8.089727e+01 8.258013e+01 1.032218e+02 3.116068e+02 1.001046e+02
#> [571] 8.965753e+01 7.933766e+01 7.552484e+01 6.888800e+01 7.658775e+01
#> [576] 2.114994e+02 2.513838e+02 1.276666e+02 1.294036e+02 1.856967e+02
#> [581] 2.127898e+02 6.697397e+01 1.783574e+02 3.040218e+02 4.203457e+02
#> [586] 9.232637e+01 8.480293e+01 7.839322e+01 1.016246e+02 7.350609e+01
#> [591] 2.745148e+02 7.139467e+01 7.859528e+02 1.245012e+02 8.265515e+01
#> [596] 7.234251e+01 6.735657e+01 6.659282e+01 9.066793e+01 9.598171e+01
#> [601] 1.035886e+02 7.394050e+01 1.137762e+02 1.028304e+03 4.315287e+08
#> [606] 1.226164e+02 2.118073e+02 1.119697e+02 2.088377e+02 6.677747e+01
#> [611] 7.417991e+01 1.034188e+02 8.250986e+01 7.491627e+01 3.652912e+03
#> [616] 2.080570e+03 6.636838e+01 6.999457e+01 6.694797e+01 6.759779e+01
#> [621] 8.284906e+01 6.550123e+01 1.170914e+02 1.318637e+02 9.157992e+01
#> [626] 7.262152e+01 1.107027e+02 8.712652e+01 9.417236e+01 1.208764e+02
#> [631] 1.559598e+02 7.323778e+01 6.918686e+01 1.241329e+02 8.232395e+03
#> [636] 1.233266e+02 9.027258e+01 1.123178e+02 8.803861e+01 1.023977e+02
#> [641] 6.626420e+01 6.708583e+01 7.349089e+01 8.829444e+01 7.074863e+01
#> [646] 9.884451e+01 7.048659e+01 9.427703e+01 1.576280e+02 2.523960e+02
#> [651] 1.310035e+02 8.329790e+01 9.944101e+01 6.765154e+02 2.968303e+02
#> [656] 9.397508e+01 7.081119e+01 3.040338e+02 1.374247e+02 8.934466e+01
#> [661] 6.578030e+01 1.503893e+02 7.018881e+01 6.298257e+01 6.801468e+01
#> [666] 7.923224e+01 1.135022e+03 1.221907e+04 1.567866e+02 1.863307e+02
#> [671] 1.152022e+02 1.004493e+02 9.465060e+01 9.242222e+01 4.412235e+02
#> [676] 8.000603e+01 7.450822e+01 8.280635e+01 2.028318e+02 1.534604e+02
#> [681] 8.427760e+01 6.861058e+01 6.945803e+01 6.735805e+01 6.845450e+01
#> [686] 1.512843e+02 7.892489e+01 6.525835e+01 5.998194e+01 6.733367e+01
#> [691] 6.109488e+01 8.157150e+01 7.563220e+01 8.615367e+01 7.243937e+01
#> [696] 6.945854e+01 8.995774e+01 6.185692e+01 8.871263e+01 8.066857e+01
#> [701] 1.037086e+02 7.920481e+01 1.636445e+02 9.378139e+01 7.802700e+01
#> [706] 1.307628e+02 7.103752e+01 1.387570e+02 1.138933e+02 3.299471e+03
#> [711] 1.362702e+02 1.882969e+02 6.984689e+01 1.015989e+02 2.528170e+02
#> [716] 7.800093e+01 9.203287e+01 1.090185e+02 1.078488e+02 1.154565e+02
#> [721] 7.986354e+01 7.880205e+01 1.146025e+02 5.338355e+02 6.590486e+01
#> [726] 8.290876e+01 3.303507e+03 1.037123e+04 3.407871e+03 1.032388e+04
#> [731] 1.340971e+02 1.041026e+02 1.610622e+02 1.570023e+02 6.173166e+01
#> [736] 6.929398e+01 6.866049e+01 6.828349e+01 9.917242e+01 1.458089e+02
#> [741] 9.192229e+01 8.218465e+01 8.446976e+01 9.002381e+01 1.458635e+02
#> [746] 6.355646e+01 9.153432e+01 7.781004e+01 6.845284e+01 6.759327e+01
#> [751] 7.970791e+01 7.337773e+01 8.312079e+01 7.330687e+01 7.386999e+01
#> [756] 7.558951e+01 7.510568e+01 1.282230e+02 7.488247e+01 1.525074e+02
#> [761] 1.266261e+02 2.461214e+02 6.982210e+01 6.918671e+01 1.708496e+02
#> [766] 6.611151e+01 8.299337e+01 1.243924e+02 7.717521e+02 1.475260e+02
#> [771] 8.099069e+01 9.930509e+01 7.084569e+01 8.072445e+01 7.077363e+01
#> [776] 1.701382e+02 1.205369e+02 2.044775e+02 1.125106e+02 6.349363e+01
#> [781] 6.943145e+01 1.458828e+02 6.901527e+01 7.389880e+01 6.811370e+01
#> [786] 1.336375e+02 7.595606e+01 1.108487e+02 1.145272e+02 6.644359e+01
#> [791] 1.135604e+02 9.942345e+01 9.157605e+01 7.917644e+01 3.692200e+02
#> [796] 1.395901e+02 8.806359e+01 1.783515e+02 1.745154e+02 1.074247e+02
#> [801] 1.158036e+02 6.544501e+01 6.319507e+01 6.452533e+01 7.558181e+01
#> [806] 8.946509e+01 2.148486e+02 1.547113e+02 9.642982e+01 1.614226e+02
#> [811] 1.146669e+02 8.241789e+01 6.326824e+01 6.492677e+01 1.873229e+02
#> [816] 7.997568e+01 9.560067e+01 1.227075e+02 1.622421e+02 1.546364e+02
#> [821] 1.081916e+02 1.108237e+02 8.370952e+01 1.545682e+02 1.058007e+02
#> [826] 4.334423e+02 1.497373e+03 8.692040e+01 8.549387e+01 1.166170e+02
#> [831] 7.997433e+01 8.483729e+01 9.482232e+01 2.355204e+03 3.298158e+02
#> [836] 1.080230e+02 6.101497e+01 6.995670e+01 7.788334e+01 7.637274e+01
#> [841] 6.444030e+01 9.947283e+01 8.488273e+01 9.536892e+01 8.029387e+01
#> [846] 6.445195e+01 7.615709e+01 1.324741e+02 7.197997e+01 8.371326e+01
#> [851] 1.012360e+02 7.921667e+02 7.096226e+01 9.254614e+01 8.169771e+01
#> [856] 9.429936e+01 6.987731e+01 7.231678e+01 4.420223e+02 7.247378e+01
#> [861] 7.712583e+01 7.599752e+01 6.435761e+01 2.236039e+02 6.798402e+01
#> [866] 6.286344e+01 6.569294e+01 7.355191e+01 1.919717e+02 2.967789e+02
#> [871] 1.363270e+02 6.593043e+01 6.703813e+01 9.162821e+01 1.632998e+02
#> [876] 1.094411e+02 8.793011e+01 6.929628e+02 6.364229e+01 1.378675e+02
#> [881] 2.941522e+02 7.107045e+01 1.814234e+02 7.997805e+01 7.116687e+01
#> [886] 6.566635e+01 5.909226e+02 1.462069e+02 1.059556e+02 7.515524e+01
#> [891] 6.769710e+01 1.190685e+02 7.600393e+01 9.487184e+01 6.086540e+01
#> [896] 5.770034e+02 1.242801e+02 1.046489e+02 7.270100e+01 8.781356e+01
#> [901] 8.583326e+01 6.944616e+01 8.137234e+01 7.138196e+01 6.747048e+01
#> [906] 8.725883e+01 6.248545e+01 1.004478e+02 1.487221e+02 3.554384e+02
#> [911] 1.161525e+03 8.149860e+07 7.643433e+01 2.302238e+02 1.158812e+02
#> [916] 7.628526e+01 7.794296e+01 6.580619e+01 1.203351e+02 7.214499e+01
#> [921] 2.642645e+02 3.446787e+02 3.042935e+03 6.483072e+01 6.311306e+01
#> [926] 6.851241e+01 7.792757e+01 3.999608e+02 6.323909e+01 2.865310e+02
#> [931] 1.687215e+03 8.319281e+02 2.868000e+02 8.622963e+01 6.195736e+01
#> [936] 6.820448e+01 6.619180e+01 6.959584e+01 1.131543e+02 3.535602e+02
#> [941] 1.347163e+02 7.042548e+01 9.655453e+01 3.785938e+02 2.693509e+02
#> [946] 9.034774e+01 8.872926e+01 8.946538e+01 6.120311e+01 8.155471e+01
#> [951] 9.343884e+01 9.210575e+01 8.350345e+01 8.734335e+01 3.019534e+02
#> [956] 1.737670e+02 1.905819e+02 2.559929e+02 8.210756e+01 1.024322e+02
#> [961] 1.846653e+02 1.636825e+02 1.084228e+02 8.783175e+01 6.826636e+01
#> [966] 1.335906e+02 1.118169e+02 6.124274e+01 8.821347e+01 8.516886e+01
#> [971] 4.433300e+03 9.956494e+01 1.214237e+02 5.858519e+05 1.013826e+02
#> [976] 6.936044e+01 5.874184e+02 7.321977e+01 1.317136e+02 9.139549e+03
#> [981] 7.915550e+01 7.696809e+01 1.377439e+02 6.958423e+01 7.839881e+01
#> [986] 1.455461e+02 7.006874e+01 6.529496e+01 9.098148e+01 1.238506e+02
#> [991] 9.115327e+02 7.417864e+01 8.389176e+01 1.195032e+02 7.524849e+01
#> [996] 1.164377e+02 7.547232e+01 8.089417e+01 8.881327e+01 7.168877e+01
doseFunction(
model,
alpha2 = samples@data$alpha2,
beta = samples@data$beta,
grade = 2L
)(x = 0.25)
#> [1] 6.081849e+01 6.703126e+01 7.396868e+01 6.307079e+01 6.106598e+01
#> [6] 1.134020e+02 6.484842e+01 7.069962e+01 7.474332e+01 6.219438e+01
#> [11] 7.866325e+01 1.203591e+02 6.483552e+01 1.267282e+02 6.425162e+01
#> [16] 1.066852e+02 9.518114e+01 5.739829e+01 5.845476e+01 5.896932e+01
#> [21] 6.250804e+01 6.414805e+01 5.996269e+01 6.321501e+01 8.266910e+01
#> [26] 5.878078e+01 5.590521e+01 5.783672e+01 5.922261e+01 6.603773e+01
#> [31] 6.215115e+01 2.156017e+02 6.280971e+01 7.434861e+01 5.705582e+01
#> [36] 6.291528e+01 6.008121e+01 7.224300e+01 6.759389e+01 8.239930e+01
#> [41] 6.437367e+01 7.162315e+01 7.070634e+01 5.939444e+01 6.653785e+01
#> [46] 6.270563e+01 6.232332e+01 7.840734e+01 1.038061e+02 7.661622e+01
#> [51] 6.293267e+01 9.746986e+01 4.839904e+02 7.778618e+01 8.114279e+01
#> [56] 1.155014e+02 6.501898e+01 5.462510e+01 5.798701e+02 3.231983e+08
#> [61] 6.828586e+01 7.722148e+01 1.339813e+02 2.699512e+02 7.921213e+01
#> [66] 8.993150e+01 6.947812e+01 9.507581e+01 5.676004e+01 1.061725e+05
#> [71] 6.695056e+02 1.464322e+02 4.300062e+02 1.141489e+03 5.431923e+01
#> [76] 5.556266e+01 8.839973e+01 6.259750e+01 4.284049e+02 4.080457e+01
#> [81] 6.178342e+01 6.935686e+01 1.894328e+02 6.178598e+01 5.786219e+01
#> [86] 6.627669e+01 5.941324e+01 8.588200e+01 4.321603e+02 8.781638e+01
#> [91] 5.796261e+01 5.670715e+01 1.246259e+02 5.942598e+01 8.011698e+01
#> [96] 6.935071e+01 5.648289e+01 5.664606e+01 8.307323e+01 5.345224e+01
#> [101] 5.798612e+01 6.223365e+01 7.437285e+01 2.410099e+02 8.711155e+01
#> [106] 5.414443e+01 8.483200e+01 1.011139e+02 5.399649e+01 5.769426e+01
#> [111] 6.046720e+01 5.985601e+01 9.011879e+01 5.641713e+01 5.760148e+01
#> [116] 1.317635e+02 9.366234e+01 1.565151e+02 6.631242e+01 6.185871e+01
#> [121] 7.837401e+01 1.039454e+02 6.062137e+01 6.341757e+01 7.420606e+01
#> [126] 6.948676e+01 5.954600e+01 6.104543e+01 1.163494e+02 1.026878e+02
#> [131] 7.037259e+01 5.759748e+01 5.470433e+01 6.079543e+01 7.247766e+01
#> [136] 7.260169e+01 6.290623e+01 5.736700e+01 5.951800e+01 5.452900e+01
#> [141] 6.292816e+01 7.477147e+01 1.112679e+02 6.595780e+01 7.659696e+01
#> [146] 1.111703e+02 2.164392e+02 6.480119e+01 6.559964e+01 7.192225e+01
#> [151] 7.043105e+01 8.129693e+01 6.687758e+01 9.436431e+01 8.062419e+01
#> [156] 4.925292e+01 1.279783e+02 8.872405e+01 8.214863e+01 6.087648e+01
#> [161] 6.227133e+01 6.349794e+01 7.252500e+01 7.143690e+01 6.698528e+01
#> [166] 1.330479e+02 5.749928e+01 5.041408e+02 6.742334e+01 2.032701e+02
#> [171] 6.275912e+01 5.956098e+01 5.638496e+01 7.536276e+01 8.512980e+01
#> [176] 6.467683e+01 5.934677e+01 6.705857e+01 1.156805e+02 8.003528e+01
#> [181] 8.540802e+01 5.537680e+01 8.762618e+01 9.279451e+01 2.974972e+02
#> [186] 6.050236e+01 5.355016e+01 5.943615e+01 7.614982e+01 5.949364e+01
#> [191] 5.659401e+01 1.057106e+02 5.739500e+01 6.947884e+01 6.629104e+01
#> [196] 5.941812e+01 6.964858e+01 5.536880e+01 5.825067e+01 5.971138e+01
#> [201] 6.161348e+01 6.778483e+01 6.870751e+01 7.734248e+01 1.890799e+02
#> [206] 6.953150e+01 6.084718e+03 7.118160e+01 2.628996e+02 8.732632e+01
#> [211] 5.660196e+01 4.056055e+02 6.345305e+01 9.998142e+01 5.975802e+01
#> [216] 6.411316e+01 6.093090e+01 7.100737e+01 8.163652e+01 8.410090e+01
#> [221] 6.725020e+01 5.746578e+01 6.009951e+01 7.424103e+01 9.287487e+01
#> [226] 6.426590e+01 5.216809e+01 9.036654e+01 4.927196e+02 9.610693e+01
#> [231] 5.610321e+01 7.073842e+01 5.844317e+01 5.807221e+01 7.697789e+01
#> [236] 6.982202e+01 1.149524e+02 5.792601e+01 1.057468e+02 6.063530e+01
#> [241] 7.079064e+01 6.202216e+01 6.956304e+01 6.327332e+01 7.435606e+01
#> [246] 1.334456e+02 1.612915e+02 4.864962e+01 7.401176e+01 5.254538e+01
#> [251] 5.968362e+01 2.923144e+01 6.180268e+01 6.620697e+01 7.856265e+01
#> [256] 9.339450e+01 6.511658e+01 5.698867e+01 6.337364e+01 5.844923e+01
#> [261] 6.618038e+01 7.193185e+01 6.597290e+01 6.062017e+01 5.705427e+01
#> [266] 6.148770e+01 6.152204e+01 7.225818e+01 9.192817e+01 5.902030e+01
#> [271] 6.941037e+01 8.521854e+01 7.970376e+01 1.200731e+02 1.180372e+02
#> [276] 5.529287e+01 6.768019e+01 7.546758e+01 7.639191e+01 5.847399e+01
#> [281] 5.795384e+01 6.751487e+01 6.866518e+01 5.366878e+01 1.981144e+02
#> [286] 2.469278e+02 6.203115e+01 6.435812e+01 6.498041e+01 5.678416e+01
#> [291] 6.414558e+01 5.510787e+01 6.910139e+01 3.703657e+02 5.904232e+01
#> [296] 7.270777e+01 6.205706e+01 6.434994e+01 6.662449e+01 6.101364e+01
#> [301] 1.135985e+02 4.145272e+01 5.798770e+01 8.825481e+01 5.295552e+01
#> [306] 8.900370e+01 5.745418e+01 2.946293e+03 6.756476e+01 5.806519e+01
#> [311] 7.004554e+01 6.209995e+01 7.523346e+01 6.417698e+01 5.592056e+01
#> [316] 1.162457e+02 5.691218e+01 1.352869e+02 6.254240e+01 6.978238e+01
#> [321] 5.605309e+01 6.459904e+01 2.652112e+02 1.091962e+02 5.799187e+01
#> [326] 9.676868e+01 6.256892e+01 5.439841e+01 7.330506e+01 1.005134e+06
#> [331] 6.823998e+01 5.014207e+01 1.390655e+02 6.555037e+01 6.351877e+01
#> [336] 7.209180e+01 6.169184e+01 8.006456e+01 5.481925e+01 6.499206e+01
#> [341] 6.042610e+01 8.422164e+01 8.247796e+01 5.442187e+01 6.045364e+01
#> [346] 5.474825e+01 7.025818e+01 6.237587e+01 5.924804e+01 6.910092e+01
#> [351] 5.973307e+01 5.636187e+01 6.193906e+01 8.732078e+01 5.751152e+01
#> [356] 6.696221e+01 1.079472e+02 6.776702e+01 5.769111e+01 5.967415e+01
#> [361] 5.753995e+01 5.841107e+01 7.727733e+01 1.272504e+02 9.297035e+01
#> [366] 5.647894e+01 5.614392e+01 1.057620e+02 7.997978e+01 6.011082e+01
#> [371] 6.724296e+01 5.675779e+01 6.528735e+01 6.249318e+01 6.458251e+01
#> [376] 6.007815e+01 2.027305e+02 3.177661e+02 1.902677e+02 8.409077e+01
#> [381] 2.768490e+02 7.761780e+01 7.035068e+01 6.046391e+01 6.103512e+01
#> [386] 5.480323e+01 1.598545e+02 8.130785e+01 6.033764e+01 5.600768e+01
#> [391] 6.748885e+01 8.690592e+01 5.666369e+01 6.882136e+01 5.936702e+01
#> [396] 5.670977e+01 5.662084e+01 5.799803e+01 6.433808e+01 7.882526e+01
#> [401] 6.205142e+01 6.764427e+01 7.342060e+01 5.912278e+01 7.301394e+01
#> [406] 6.150453e+01 5.347452e+01 5.744296e+01 6.982845e+01 8.343008e+01
#> [411] 5.722522e+01 5.682821e+01 8.587654e+01 5.432633e+01 1.563573e+02
#> [416] 6.814130e+01 7.082714e+01 6.548185e+01 1.074736e+02 6.028043e+01
#> [421] 9.248256e+01 6.295656e+01 7.567898e+01 6.414651e+01 9.834632e+01
#> [426] 7.946655e+01 5.525253e+01 6.713862e+01 5.884357e+01 6.135265e+01
#> [431] 1.558142e+02 9.927286e+01 5.571308e+01 5.695896e+01 6.217080e+01
#> [436] 6.472928e+01 4.956508e+01 1.337150e+02 8.796826e+01 7.141026e+01
#> [441] 4.979280e+01 7.087165e+01 6.069729e+01 7.215311e+01 6.355417e+01
#> [446] 1.486978e+02 5.983824e+01 7.893327e+01 5.766005e+01 9.124899e+01
#> [451] 6.029031e+01 5.761745e+01 7.211919e+01 6.343812e+01 8.102133e+01
#> [456] 7.305671e+01 6.251207e+01 5.760241e+01 6.721179e+01 6.327914e+01
#> [461] 6.228831e+01 7.462994e+01 7.941144e+01 9.570791e+01 6.619047e+01
#> [466] 1.438420e+02 1.280329e+02 7.663962e+01 5.460685e+01 5.642236e+01
#> [471] 2.869870e+02 9.099107e+01 1.852641e+04 8.940088e+01 6.922552e+01
#> [476] 3.236668e+02 6.915189e+01 5.440966e+01 6.713429e+01 6.826660e+01
#> [481] 5.879644e+01 6.094280e+01 5.783137e+01 6.002195e+01 3.124887e+02
#> [486] 5.754929e+01 6.602592e+01 7.267604e+01 9.313561e+01 5.852569e+01
#> [491] 1.418477e+02 5.976184e+01 6.087720e+01 6.921068e+01 9.248533e+01
#> [496] 6.244720e+01 6.005778e+01 1.262838e+02 8.209956e+01 6.024103e+01
#> [501] 6.243239e+01 6.044960e+01 5.890966e+01 5.791600e+01 6.925095e+01
#> [506] 6.778723e+01 7.866986e+01 5.640309e+01 9.460555e+01 6.275665e+01
#> [511] 6.910719e+01 8.024626e+01 5.797663e+01 5.673878e+01 8.439439e+01
#> [516] 6.321479e+01 5.430249e+01 7.585915e+01 5.766067e+01 7.802607e+01
#> [521] 1.263879e+02 5.809498e+01 8.979134e+01 7.027212e+01 6.891226e+01
#> [526] 6.132022e+01 6.110326e+01 6.279486e+01 8.447344e+01 9.563183e+01
#> [531] 2.415404e+01 5.851744e+01 5.644787e+01 6.518132e+01 7.923319e+01
#> [536] 1.165246e+02 5.874782e+01 8.836582e+01 9.521180e+01 6.387137e+01
#> [541] 6.809799e+01 2.579606e+02 9.898788e+01 6.174607e+01 6.168396e+01
#> [546] 6.261695e+01 8.637257e+01 6.659228e+01 6.309850e+01 7.206905e+01
#> [551] 7.072031e+01 5.585158e+01 5.839825e+01 5.871376e+01 5.824667e+01
#> [556] 6.203457e+01 8.452309e+01 8.313432e+01 7.459225e+01 1.101199e+02
#> [561] 6.232256e+01 5.826223e+01 6.734158e+01 6.087587e+01 5.919812e+01
#> [566] 6.086868e+01 6.035031e+01 6.285283e+01 1.240015e+02 7.577757e+01
#> [571] 7.350915e+01 6.487866e+01 6.298476e+01 5.520314e+01 5.756180e+01
#> [576] 1.102874e+02 1.282750e+02 8.605468e+01 8.141086e+01 9.428146e+01
#> [581] 8.330977e+01 5.737757e+01 8.044113e+01 1.161113e+02 1.998714e+02
#> [586] 6.103663e+01 6.003571e+01 6.516027e+01 7.668447e+01 5.631315e+01
#> [591] 9.078524e+01 5.747010e+01 1.482824e+02 9.124455e+01 6.957683e+01
#> [596] 6.327823e+01 5.917027e+01 5.879641e+01 6.864089e+01 6.616051e+01
#> [601] 6.394360e+01 6.266867e+01 8.091874e+01 2.074288e+02 2.390188e+05
#> [606] 7.743203e+01 1.149503e+02 8.438859e+01 1.216673e+02 5.616844e+01
#> [611] 6.299647e+01 7.231980e+01 6.708028e+01 6.112677e+01 6.195883e+02
#> [616] 3.548503e+02 5.823663e+01 5.813660e+01 5.820518e+01 5.785983e+01
#> [621] 5.833587e+01 5.670905e+01 7.322799e+01 5.748657e+01 6.637027e+01
#> [626] 6.213467e+01 7.127117e+01 6.485768e+01 6.534735e+01 8.909062e+01
#> [631] 9.940669e+01 5.352717e+01 5.694030e+01 7.598929e+01 6.425413e+02
#> [636] 7.523594e+01 5.912156e+01 7.133758e+01 6.518324e+01 6.846030e+01
#> [641] 5.802058e+01 5.791316e+01 6.236844e+01 5.981052e+01 6.043273e+01
#> [646] 6.031376e+01 5.531848e+01 5.935553e+01 7.088252e+01 1.282715e+02
#> [651] 5.948458e+01 6.638462e+01 7.577987e+01 1.986679e+02 1.082486e+02
#> [656] 7.150779e+01 5.542972e+01 1.199001e+02 8.547551e+01 6.610362e+01
#> [661] 5.718496e+01 8.962610e+01 5.855406e+01 5.520378e+01 6.101362e+01
#> [666] 6.200074e+01 1.407862e+02 5.090035e+02 7.241117e+01 9.590012e+01
#> [671] 7.929257e+01 6.332784e+01 7.040286e+01 5.707080e+01 1.234455e+02
#> [676] 6.084646e+01 5.990444e+01 6.193089e+01 8.666969e+01 8.491322e+01
#> [681] 7.117280e+01 6.283382e+01 6.046126e+01 5.858942e+01 5.699123e+01
#> [686] 8.661850e+01 6.667410e+01 5.766620e+01 5.505559e+01 6.017975e+01
#> [691] 5.523804e+01 6.860701e+01 5.908755e+01 6.388519e+01 6.057634e+01
#> [696] 5.895241e+01 7.266349e+01 5.308011e+01 6.567535e+01 5.913225e+01
#> [701] 6.320294e+01 6.163036e+01 9.603652e+01 6.142340e+01 5.802572e+01
#> [706] 7.074244e+01 5.930207e+01 6.712482e+01 7.294651e+01 5.150636e+02
#> [711] 7.947739e+01 7.929383e+01 5.742613e+01 7.050025e+01 1.440390e+02
#> [716] 5.210899e+01 6.810085e+01 6.917113e+01 6.770867e+01 8.003432e+01
#> [721] 6.160358e+01 6.041156e+01 6.410767e+01 2.401939e+02 5.201572e+01
#> [726] 6.216956e+01 1.307031e+02 2.740170e+02 9.544737e+01 3.268001e+02
#> [731] 8.048767e+01 6.785020e+01 7.756505e+01 7.085377e+01 5.772790e+01
#> [736] 6.137691e+01 6.073948e+01 5.872361e+01 6.822887e+01 7.422642e+01
#> [741] 6.874407e+01 6.064674e+01 6.261531e+01 6.265863e+01 8.366355e+01
#> [746] 5.765734e+01 7.429466e+01 5.945714e+01 5.869159e+01 6.072131e+01
#> [751] 6.425931e+01 6.189828e+01 6.896146e+01 6.329176e+01 5.636156e+01
#> [756] 6.200044e+01 5.821870e+01 8.353407e+01 5.445148e+01 8.963334e+01
#> [761] 7.602408e+01 1.251497e+02 6.063736e+01 5.776728e+01 7.865290e+01
#> [766] 5.920155e+01 6.408420e+01 5.939942e+01 1.040276e+02 5.540256e+01
#> [771] 5.493961e+01 7.556892e+01 6.197499e+01 6.343823e+01 5.593636e+01
#> [776] 9.693804e+01 8.297070e+01 8.071850e+01 7.544130e+01 5.802264e+01
#> [781] 6.124757e+01 9.548976e+01 5.654379e+01 6.299798e+01 5.873281e+01
#> [786] 9.073658e+01 6.187716e+01 7.330303e+01 8.122947e+01 5.980051e+01
#> [791] 7.688607e+01 7.692708e+01 6.006606e+01 5.399370e+01 9.569308e+01
#> [796] 7.150264e+01 6.129142e+01 1.093976e+02 9.596031e+01 5.733385e+01
#> [801] 8.809390e+01 5.781708e+01 5.872599e+01 5.994580e+01 6.148602e+01
#> [806] 6.287969e+01 1.132889e+02 9.298410e+01 6.280828e+01 8.055917e+01
#> [811] 7.318072e+01 6.441625e+01 5.706922e+01 5.775243e+01 1.102400e+02
#> [816] 6.636330e+01 7.605338e+01 7.696935e+01 7.721495e+01 8.654275e+01
#> [821] 7.613998e+01 7.815203e+01 7.087850e+01 9.866307e+01 6.936200e+01
#> [826] 1.867545e+02 3.187427e+02 6.701437e+01 6.673461e+01 8.735972e+01
#> [831] 6.057823e+01 6.507306e+01 7.339518e+01 4.675528e+02 9.531812e+01
#> [836] 7.245101e+01 5.534734e+01 5.939982e+01 6.012918e+01 6.621723e+01
#> [841] 5.730977e+01 7.091370e+01 6.341758e+01 6.912805e+01 5.837659e+01
#> [846] 5.738657e+01 6.698910e+01 7.227672e+01 5.135548e+01 6.106852e+01
#> [851] 7.628790e+01 2.810464e+02 5.388985e+01 7.121770e+01 6.362477e+01
#> [856] 6.677215e+01 5.969526e+01 5.909933e+01 1.707473e+02 6.106141e+01
#> [861] 6.581667e+01 6.343501e+01 5.653408e+01 1.023070e+02 5.963018e+01
#> [866] 5.419764e+01 5.487705e+01 5.741543e+01 6.095320e+01 1.236359e+02
#> [871] 8.203340e+01 5.874071e+01 5.852409e+01 7.028116e+01 1.047467e+02
#> [876] 5.426765e+01 6.538778e+01 2.743184e+02 5.639182e+01 8.181985e+01
#> [881] 1.038173e+02 5.997094e+01 9.757172e+01 4.927210e+01 5.719085e+01
#> [886] 5.731540e+01 1.640227e+02 7.101526e+01 6.536442e+01 6.443390e+01
#> [891] 5.611825e+01 7.750368e+01 5.979424e+01 7.706807e+01 5.399199e+01
#> [896] 1.787258e+02 7.751496e+01 6.688427e+01 5.944458e+01 6.379455e+01
#> [901] 6.539275e+01 6.103986e+01 6.905559e+01 6.244155e+01 6.027004e+01
#> [906] 7.021485e+01 5.497540e+01 6.753754e+01 8.870550e+01 9.785344e+01
#> [911] 2.224355e+02 3.243479e+04 6.426553e+01 7.306104e+01 7.963043e+01
#> [916] 6.506338e+01 6.663967e+01 5.987556e+01 9.054815e+01 5.939179e+01
#> [921] 1.140199e+02 1.135018e+02 3.513165e+02 5.976155e+01 5.758537e+01
#> [926] 5.812480e+01 6.302911e+01 1.177548e+02 5.783218e+01 1.702908e+02
#> [931] 4.373453e+02 2.666452e+02 1.266857e+02 6.213143e+01 5.537726e+01
#> [936] 6.142211e+01 6.066340e+01 6.213891e+01 9.307516e+01 1.936126e+02
#> [941] 9.754972e+01 5.981553e+01 7.445156e+01 1.624627e+02 9.717378e+01
#> [946] 7.139592e+01 7.267399e+01 6.089640e+01 5.426398e+01 6.463731e+01
#> [951] 7.094833e+01 7.213250e+01 7.072809e+01 6.628473e+01 1.544755e+02
#> [956] 7.782077e+01 7.439772e+01 9.991529e+01 6.537224e+01 7.742958e+01
#> [961] 1.049651e+02 7.291247e+01 6.042193e+01 5.504434e+01 5.900651e+01
#> [966] 7.720475e+01 6.590764e+01 5.427878e+01 6.746266e+01 5.582711e+01
#> [971] 2.187786e+02 7.129174e+01 7.511970e+01 3.968267e+03 7.815874e+01
#> [976] 6.297766e+01 2.277017e+02 6.320530e+01 9.394167e+01 1.028890e+03
#> [981] 6.462315e+01 6.095550e+01 6.894586e+01 5.349056e+01 6.544230e+01
#> [986] 7.105805e+01 6.185790e+01 5.955956e+01 6.876143e+01 9.334425e+01
#> [991] 4.138910e+02 5.854251e+01 6.681189e+01 5.483406e+01 6.188273e+01
#> [996] 6.118029e+01 5.510007e+01 5.807351e+01 6.326679e+01 5.740576e+01