![[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] 6.476251e+01 7.020646e+01 1.026954e+02 6.527192e+02 7.468853e+01
#> [6] 7.763991e+01 8.182519e+01 2.977043e+02 1.776080e+02 7.027963e+01
#> [11] 8.286808e+01 8.423098e+01 6.637516e+01 7.079726e+01 7.108300e+01
#> [16] 6.667132e+01 1.671552e+02 7.629784e+01 1.232707e+02 1.598916e+02
#> [21] 8.242496e+01 9.029389e+01 6.891669e+01 1.370842e+02 5.876431e+01
#> [26] 9.477470e+01 7.225244e+02 7.769562e+01 1.191198e+02 2.823788e+02
#> [31] 4.564827e+02 7.327728e+01 6.796814e+01 8.255846e+01 6.646266e+01
#> [36] 2.920236e+02 7.227287e+01 7.519243e+01 2.085610e+02 2.580755e+02
#> [41] 8.480305e+01 7.265777e+01 7.556544e+01 8.031796e+01 1.803595e+03
#> [46] 3.474450e+02 8.863897e+01 7.068657e+01 6.725077e+01 8.075429e+01
#> [51] 8.343060e+01 6.706104e+01 6.247578e+01 6.712644e+01 7.927549e+01
#> [56] 6.881728e+01 6.800607e+02 2.141552e+02 2.144643e+02 3.056750e+02
#> [61] 1.225896e+02 8.743284e+01 6.356930e+01 2.868397e+03 8.398313e+02
#> [66] 1.355752e+02 7.800288e+01 7.711572e+01 3.291242e+02 2.394430e+05
#> [71] 6.287547e+01 7.573049e+01 1.260875e+02 6.280173e+01 6.265153e+01
#> [76] 6.609855e+01 1.533493e+02 6.842277e+01 7.971971e+01 1.275430e+02
#> [81] 1.358853e+02 9.152592e+01 8.033473e+01 8.467619e+01 1.389135e+02
#> [86] 7.178960e+01 8.431796e+01 6.573025e+01 7.335199e+01 8.974734e+01
#> [91] 1.926628e+02 1.628679e+02 1.966695e+02 7.129323e+01 1.027634e+02
#> [96] 8.427243e+01 9.486081e+01 8.353327e+02 7.604163e+01 7.124741e+01
#> [101] 6.043877e+01 6.000616e+01 5.753750e+03 9.863653e+01 7.418418e+01
#> [106] 7.024584e+01 6.832579e+01 8.648597e+01 8.021685e+02 3.539852e+03
#> [111] 8.372947e+03 2.642727e+03 1.447266e+02 6.599536e+01 9.996779e+01
#> [116] 8.693677e+01 9.440446e+01 7.833483e+01 9.219222e+01 6.911644e+01
#> [121] 9.358932e+01 1.494444e+02 1.585323e+02 8.341322e+01 7.463349e+01
#> [126] 6.362146e+01 9.592502e+01 7.093337e+02 9.192672e+01 8.946778e+01
#> [131] 5.980252e+01 1.000228e+03 9.166786e+01 7.021912e+01 5.890885e+01
#> [136] 2.968271e+03 6.425116e+01 6.648387e+01 7.750626e+01 3.246747e+02
#> [141] 1.184969e+02 7.165461e+01 1.037573e+02 7.835013e+01 7.264219e+01
#> [146] 7.257275e+01 8.438229e+01 2.052609e+02 5.198852e+02 1.335608e+02
#> [151] 1.207047e+02 1.199424e+02 7.294916e+01 1.317885e+02 4.452180e+02
#> [156] 7.034524e+01 2.554794e+03 7.225849e+01 7.172388e+01 1.574827e+02
#> [161] 1.245372e+02 1.092224e+02 1.031200e+02 1.132205e+02 1.598628e+02
#> [166] 9.658427e+01 7.003287e+01 2.987661e+03 7.860633e+01 2.230881e+02
#> [171] 7.030230e+01 9.676165e+01 6.298325e+01 6.734468e+01 6.374984e+01
#> [176] 9.311227e+01 1.587036e+02 1.256471e+02 6.531835e+01 1.112954e+02
#> [181] 2.775570e+02 6.848930e+01 7.771765e+01 1.043859e+02 7.030245e+01
#> [186] 6.794825e+01 1.256519e+02 3.205317e+02 6.319923e+02 5.634501e+02
#> [191] 1.904699e+02 9.318418e+01 1.122644e+02 7.552015e+01 1.914717e+02
#> [196] 7.033098e+01 6.760343e+01 9.704743e+01 1.407161e+02 2.326310e+02
#> [201] 7.339075e+01 7.234706e+01 6.300598e+01 4.000158e+02 6.880332e+01
#> [206] 7.703466e+01 8.309903e+01 9.120664e+01 1.113324e+02 7.426031e+01
#> [211] 1.153639e+02 8.900335e+01 9.080817e+01 7.519074e+01 9.706121e+01
#> [216] 7.354564e+01 1.197236e+02 6.045756e+02 1.530463e+04 2.631650e+02
#> [221] 2.161309e+02 8.744225e+01 7.202865e+01 7.757442e+01 6.964913e+01
#> [226] 1.706712e+02 1.029810e+02 7.895322e+01 7.850384e+01 1.232780e+02
#> [231] 7.254858e+01 1.322970e+02 1.031159e+02 8.323600e+01 3.854374e+02
#> [236] 6.576895e+01 1.081852e+02 1.688635e+02 4.621466e+02 8.590151e+01
#> [241] 7.509940e+01 9.640959e+01 2.395409e+02 1.106646e+02 4.066443e+02
#> [246] 8.982108e+01 9.981636e+01 9.955784e+01 7.255498e+01 1.073206e+02
#> [251] 1.776600e+02 7.007138e+01 8.540737e+01 1.085202e+02 1.378787e+02
#> [256] 1.198315e+02 6.840600e+01 1.311628e+02 7.438064e+01 9.752026e+01
#> [261] 1.151052e+02 1.600134e+02 1.342287e+02 1.179601e+02 7.542873e+01
#> [266] 6.988952e+01 7.594995e+01 6.219385e+02 1.292598e+02 8.656759e+01
#> [271] 7.287268e+01 7.210576e+01 9.624833e+01 8.264845e+01 1.230917e+02
#> [276] 2.404407e+02 8.730283e+01 1.697776e+02 7.726292e+01 7.966270e+01
#> [281] 8.832609e+01 6.984392e+01 6.597544e+01 9.737927e+01 6.671812e+01
#> [286] 8.472270e+01 6.593578e+01 7.603272e+01 7.869095e+01 1.033121e+02
#> [291] 1.412480e+02 5.952452e+02 7.138917e+01 2.387871e+02 1.946971e+02
#> [296] 6.630375e+01 3.241381e+02 1.966470e+02 8.274373e+01 1.071227e+03
#> [301] 7.621197e+01 1.078055e+02 1.367346e+02 7.553395e+01 7.366797e+01
#> [306] 8.088116e+01 6.940618e+01 9.221128e+01 7.235684e+01 7.508614e+01
#> [311] 7.458566e+01 8.217835e+01 7.315025e+01 1.120193e+02 9.046045e+01
#> [316] 2.527451e+02 8.696823e+01 1.616356e+03 1.180767e+02 7.093718e+01
#> [321] 6.807500e+01 1.116030e+02 6.684084e+01 9.273868e+02 7.556772e+01
#> [326] 6.904416e+01 1.995757e+02 7.787755e+02 7.927591e+01 1.988528e+02
#> [331] 3.526387e+03 7.136034e+01 7.289621e+01 6.706962e+01 7.056716e+01
#> [336] 9.975627e+01 7.221224e+01 7.582036e+01 8.113972e+01 9.283163e+01
#> [341] 6.664558e+01 5.266005e+02 1.922375e+03 6.509812e+02 6.342684e+01
#> [346] 8.764742e+01 7.335957e+01 8.159780e+01 6.914691e+02 7.977497e+01
#> [351] 6.957800e+01 7.198690e+01 2.061605e+02 8.487623e+01 6.674621e+01
#> [356] 1.190199e+02 9.192097e+01 8.791275e+01 1.053327e+03 1.150932e+02
#> [361] 8.269952e+01 3.622734e+02 1.006814e+02 1.716417e+02 9.814982e+01
#> [366] 1.010119e+02 1.096099e+02 6.315563e+02 6.050156e+02 9.156623e+01
#> [371] 7.045469e+01 8.188485e+01 1.363645e+02 1.300010e+02 2.745445e+02
#> [376] 8.967484e+01 2.138200e+02 1.356541e+02 7.207300e+01 6.744537e+01
#> [381] 2.902520e+02 7.095738e+01 6.624846e+01 1.340992e+02 1.110007e+02
#> [386] 7.362015e+01 7.366446e+01 7.716376e+01 8.425694e+01 6.365690e+01
#> [391] 6.826704e+01 7.893412e+01 7.258273e+01 1.211212e+02 1.476798e+02
#> [396] 8.580473e+01 3.108651e+02 2.018485e+03 9.459428e+01 9.232246e+01
#> [401] 1.218563e+03 7.008947e+01 6.825397e+01 9.123023e+01 1.300888e+02
#> [406] 1.245295e+02 9.939609e+01 1.113804e+02 9.107012e+01 1.151779e+02
#> [411] 7.191102e+01 1.278534e+02 7.692362e+01 9.218521e+01 7.857749e+01
#> [416] 7.448155e+01 7.869074e+01 8.269622e+01 6.961199e+01 8.926039e+01
#> [421] 1.891462e+02 8.547447e+02 1.839421e+07 9.123907e+01 8.637468e+01
#> [426] 1.186730e+02 9.222487e+01 8.023388e+01 7.708123e+01 1.933891e+02
#> [431] 2.646217e+02 9.362213e+01 1.136706e+02 1.017685e+02 9.680191e+02
#> [436] 5.181860e+02 7.611615e+01 6.709059e+01 1.102885e+02 1.664242e+02
#> [441] 8.401461e+01 4.096219e+04 6.599995e+01 1.955914e+02 7.310118e+01
#> [446] 7.954811e+01 6.414742e+01 7.550895e+01 8.265965e+01 7.419311e+01
#> [451] 1.479239e+02 2.684554e+02 1.051750e+02 4.814836e+02 3.889347e+02
#> [456] 1.158477e+02 8.275736e+01 1.050049e+02 2.552087e+12 2.232629e+05
#> [461] 2.821861e+02 5.527371e+02 1.144077e+02 1.100589e+02 8.261165e+01
#> [466] 1.520600e+02 9.971257e+01 9.714708e+01 8.410557e+01 1.275026e+02
#> [471] 7.264661e+01 8.615116e+01 1.260144e+02 7.663157e+01 7.678731e+01
#> [476] 1.299681e+02 9.791286e+01 1.090046e+02 8.103903e+01 6.765774e+01
#> [481] 6.621661e+01 7.307031e+01 7.342505e+01 9.237110e+01 2.286888e+02
#> [486] 7.117028e+01 6.885373e+01 7.773038e+01 2.248738e+04 6.616060e+01
#> [491] 7.701058e+01 1.460720e+02 1.083373e+03 1.497727e+02 6.677760e+01
#> [496] 6.521872e+01 6.407083e+01 6.712234e+01 7.212623e+01 2.367360e+02
#> [501] 7.893887e+01 2.327485e+02 6.439337e+01 6.909377e+01 7.349074e+01
#> [506] 8.627779e+01 1.081731e+02 7.539935e+01 6.656907e+01 8.024938e+01
#> [511] 7.670209e+01 6.678080e+01 7.446643e+01 2.866849e+03 2.716408e+03
#> [516] 7.977677e+01 1.127490e+02 1.685218e+02 1.108350e+02 9.048950e+01
#> [521] 6.802498e+01 8.454243e+01 7.979917e+01 9.133722e+01 1.563184e+02
#> [526] 9.541816e+01 9.950284e+01 7.733543e+01 7.471863e+01 1.003600e+02
#> [531] 3.060566e+02 7.303199e+01 8.917069e+01 7.877892e+01 3.755507e+02
#> [536] 8.237973e+01 9.165244e+01 6.843459e+01 6.711198e+01 6.670480e+01
#> [541] 7.866440e+01 7.896379e+01 1.470705e+02 7.404153e+01 6.264231e+01
#> [546] 8.555847e+01 1.934829e+02 8.838226e+01 1.477464e+02 1.184625e+02
#> [551] 4.967784e+02 8.352806e+02 1.772184e+02 1.950931e+02 4.330350e+02
#> [556] 7.486616e+01 8.606744e+01 7.917713e+01 6.272284e+01 6.489527e+01
#> [561] 1.482257e+02 1.841087e+02 9.428785e+01 6.384638e+01 1.303425e+02
#> [566] 7.390472e+01 5.894070e+01 7.099472e+02 2.938529e+03 8.295572e+01
#> [571] 1.494001e+02 1.408643e+02 7.975813e+01 7.330429e+01 1.420639e+02
#> [576] 6.971246e+01 7.326847e+01 1.144354e+02 7.464008e+01 7.397261e+01
#> [581] 7.236846e+01 1.132664e+02 7.560027e+01 2.051720e+02 9.058927e+01
#> [586] 6.693492e+01 1.260070e+02 1.291954e+02 8.546920e+01 1.955033e+02
#> [591] 7.011071e+01 6.808239e+01 6.451530e+01 6.846130e+01 8.503864e+01
#> [596] 8.523224e+01 8.468107e+01 3.825429e+02 2.860469e+02 6.874036e+01
#> [601] 8.728728e+01 1.412427e+02 2.398708e+02 8.197599e+01 6.557159e+01
#> [606] 7.019079e+01 9.121373e+01 8.719222e+01 6.589947e+01 5.867673e+01
#> [611] 1.239231e+02 2.673266e+02 9.388798e+01 6.836665e+04 6.690831e+01
#> [616] 7.512522e+01 9.659924e+01 8.296116e+01 1.401504e+02 1.128624e+02
#> [621] 9.190680e+01 9.094205e+01 6.928032e+01 7.278008e+01 1.920225e+02
#> [626] 9.649438e+01 1.161009e+02 1.306185e+02 6.521609e+01 1.058154e+02
#> [631] 8.515874e+01 8.756251e+01 8.391631e+01 6.867602e+01 6.120606e+01
#> [636] 7.055924e+01 7.005975e+01 1.405778e+02 7.798313e+01 7.582124e+01
#> [641] 6.812230e+01 7.020362e+01 1.052281e+02 6.752251e+01 7.180126e+01
#> [646] 1.579241e+02 6.748426e+01 7.697277e+01 1.439411e+02 4.446080e+02
#> [651] 7.524346e+01 3.365191e+02 2.225744e+02 2.424803e+02 6.146741e+01
#> [656] 6.805589e+01 8.514257e+01 7.405848e+01 7.501347e+01 1.057884e+02
#> [661] 2.535432e+02 6.942509e+02 8.561715e+01 1.020867e+02 7.896359e+01
#> [666] 1.579661e+02 1.447813e+02 1.070500e+02 1.362495e+02 6.626812e+01
#> [671] 1.019840e+02 7.171061e+01 8.583847e+01 7.056693e+01 7.098070e+01
#> [676] 6.731104e+01 6.817918e+01 6.870085e+01 7.575269e+01 7.041339e+01
#> [681] 7.590515e+01 9.125332e+01 6.759940e+01 1.286329e+02 6.981580e+01
#> [686] 8.888443e+01 7.549215e+01 6.718845e+01 1.507989e+02 1.597189e+04
#> [691] 7.823211e+01 1.780799e+02 1.148220e+02 8.038298e+01 1.176202e+02
#> [696] 7.837969e+01 1.105966e+02 1.481993e+02 1.071765e+02 4.168827e+03
#> [701] 6.156080e+01 7.836077e+01 7.243672e+01 1.095720e+02 7.285830e+01
#> [706] 6.351042e+03 2.263596e+03 1.174887e+02 6.286835e+01 7.212207e+01
#> [711] 7.816874e+01 7.210322e+01 1.392357e+02 6.470951e+01 1.376082e+02
#> [716] 1.106330e+02 6.317829e+01 8.396330e+01 7.804915e+01 7.633964e+01
#> [721] 6.798748e+01 4.568015e+02 3.867052e+02 7.395202e+01 7.025596e+01
#> [726] 2.641158e+02 4.336532e+02 9.656947e+01 1.017080e+02 1.231990e+02
#> [731] 7.023373e+01 1.725995e+02 8.199261e+01 6.694853e+01 9.153835e+01
#> [736] 1.122834e+02 9.457242e+01 7.243026e+01 2.051793e+02 6.718328e+01
#> [741] 7.338936e+01 7.760182e+01 6.696705e+01 6.358401e+01 7.334841e+01
#> [746] 2.963830e+02 7.655457e+01 7.047785e+01 7.393036e+01 3.848600e+03
#> [751] 7.095550e+01 6.644483e+01 7.316487e+01 7.305330e+01 7.482229e+01
#> [756] 7.247179e+01 8.094455e+01 6.901456e+01 7.205877e+01 1.473187e+02
#> [761] 8.692556e+01 7.167042e+01 6.593484e+01 2.296907e+02 8.298217e+01
#> [766] 1.093840e+02 1.491154e+02 1.959047e+02 1.137396e+02 1.001381e+02
#> [771] 8.562824e+01 6.568650e+01 6.298056e+01 8.975045e+01 1.536704e+02
#> [776] 9.301727e+01 9.482071e+01 9.573606e+01 7.321411e+01 7.322994e+01
#> [781] 8.239558e+01 1.616931e+02 6.968275e+01 8.254443e+01 7.547460e+01
#> [786] 1.314900e+02 1.203879e+02 8.333395e+01 1.058254e+02 1.766565e+02
#> [791] 7.257177e+01 7.176990e+01 9.811580e+01 9.076037e+01 5.173111e+02
#> [796] 3.919850e+03 6.880501e+01 7.295209e+01 9.669641e+01 7.529535e+01
#> [801] 6.349947e+01 8.006598e+03 6.772538e+01 1.747961e+02 3.237314e+02
#> [806] 1.516888e+02 7.104155e+01 6.525103e+01 6.789917e+01 1.112154e+02
#> [811] 7.486833e+01 2.973333e+02 1.004305e+02 7.865261e+01 7.367113e+01
#> [816] 7.952102e+01 1.274164e+02 1.794451e+02 7.525776e+01 6.558396e+01
#> [821] 7.304270e+01 1.103166e+02 8.464695e+01 2.472006e+02 1.476314e+02
#> [826] 1.023231e+02 1.341314e+02 6.503005e+01 6.413659e+01 8.533802e+01
#> [831] 6.360389e+01 7.593925e+01 6.727113e+01 8.470445e+01 2.830300e+04
#> [836] 7.132341e+01 7.010999e+01 7.882742e+01 6.872165e+01 7.162650e+01
#> [841] 2.011198e+02 7.250900e+01 7.389692e+01 2.151218e+02 1.427923e+02
#> [846] 1.164695e+02 6.361953e+01 7.621872e+01 1.191559e+02 1.180951e+02
#> [851] 1.652520e+03 3.584034e+02 9.864879e+01 7.495652e+01 6.361912e+01
#> [856] 7.096141e+01 6.901178e+01 1.187852e+02 1.024054e+02 7.427953e+01
#> [861] 1.113455e+02 9.822803e+01 1.196447e+02 6.386018e+01 1.080505e+02
#> [866] 1.119142e+02 1.391341e+02 7.153981e+01 9.021990e+01 7.771148e+01
#> [871] 1.239016e+02 1.217254e+02 8.187280e+01 7.665072e+01 1.913104e+02
#> [876] 3.822255e+02 1.695695e+02 7.047012e+01 6.724265e+01 7.091045e+01
#> [881] 8.244235e+01 1.992220e+02 6.994172e+01 7.512953e+01 8.096522e+01
#> [886] 1.697935e+02 7.906813e+01 2.792316e+02 6.312285e+01 1.185751e+02
#> [891] 8.965397e+01 1.153090e+02 9.978653e+01 1.647554e+02 8.112190e+01
#> [896] 8.580645e+01 1.668913e+02 1.634801e+02 8.723752e+01 1.908545e+02
#> [901] 1.475755e+02 1.581648e+02 7.490527e+01 8.094834e+01 1.225124e+02
#> [906] 1.450097e+02 8.187567e+01 2.176894e+02 6.588697e+01 8.334761e+01
#> [911] 7.436465e+01 1.002861e+02 7.635719e+02 3.944005e+02 7.943397e+01
#> [916] 8.279255e+01 6.803038e+01 7.418373e+01 4.228045e+02 1.401139e+02
#> [921] 6.909322e+01 8.830791e+02 8.105783e+01 7.841475e+01 7.630375e+01
#> [926] 1.454726e+02 8.836129e+01 1.295526e+02 2.540202e+02 3.791047e+07
#> [931] 1.854005e+03 6.825947e+01 7.376837e+01 6.920301e+01 6.658237e+01
#> [936] 7.033689e+01 6.620245e+01 6.954614e+01 6.679907e+01 2.536059e+02
#> [941] 6.048020e+01 7.770336e+01 1.819938e+02 6.390682e+01 6.908883e+01
#> [946] 6.357661e+01 7.048418e+01 1.566138e+02 6.939159e+01 6.767975e+01
#> [951] 7.220129e+01 1.867074e+02 5.993023e+01 9.248633e+01 9.393210e+02
#> [956] 1.438969e+02 7.056514e+01 5.285100e+02 1.200269e+02 2.087922e+02
#> [961] 7.907742e+01 4.539235e+02 6.609407e+01 1.125178e+02 7.621259e+01
#> [966] 2.138795e+02 7.868205e+01 7.550897e+01 6.625928e+01 2.412229e+02
#> [971] 1.473391e+02 2.979393e+03 2.249013e+04 5.862795e+02 3.344835e+02
#> [976] 6.760064e+01 2.199278e+02 1.094825e+02 1.281854e+02 8.658067e+01
#> [981] 1.639206e+02 8.335935e+01 1.063132e+02 7.313996e+01 3.661010e+02
#> [986] 6.809954e+03 6.797540e+01 2.441035e+02 1.470205e+02 9.888653e+01
#> [991] 1.506150e+02 7.447405e+01 7.215397e+02 1.473473e+02 9.153377e+01
#> [996] 6.748758e+01 9.634901e+01 6.897188e+01 1.176767e+02 7.756156e+01
doseFunction(
model,
alpha2 = samples@data$alpha2,
beta = samples@data$beta,
grade = 2L
)(x = 0.25)
#> [1] 5.857635e+01 5.950914e+01 6.804319e+01 1.329049e+02 5.868427e+01
#> [6] 6.668461e+01 6.304226e+01 1.475808e+02 7.465004e+01 5.496307e+01
#> [11] 6.139697e+01 6.604199e+01 6.067577e+01 6.238878e+01 6.001156e+01
#> [16] 5.671890e+01 9.500395e+01 6.413186e+01 7.641492e+01 9.194001e+01
#> [21] 6.769784e+01 5.729458e+01 5.815707e+01 8.060630e+01 4.950793e+01
#> [26] 7.173976e+01 1.330001e+02 6.347664e+01 7.119729e+01 8.992277e+01
#> [31] 1.815744e+02 6.390770e+01 6.028012e+01 6.557779e+01 5.651742e+01
#> [36] 1.107701e+02 6.207333e+01 6.235126e+01 1.338936e+02 1.347786e+02
#> [41] 6.353866e+01 5.924441e+01 6.301367e+01 6.633344e+01 2.009661e+02
#> [46] 9.207738e+01 6.513705e+01 5.534972e+01 5.734456e+01 6.520738e+01
#> [51] 6.033151e+01 6.115720e+01 5.849110e+01 6.094881e+01 6.123348e+01
#> [56] 5.872289e+01 1.228998e+02 1.084612e+02 1.188495e+02 1.611437e+02
#> [61] 9.096187e+01 7.430817e+01 5.711553e+01 1.209288e+02 1.766504e+02
#> [66] 9.713013e+01 6.213325e+01 6.464404e+01 1.450240e+02 1.092411e+04
#> [71] 5.524444e+01 6.525118e+01 9.551948e+01 5.709017e+01 5.505618e+01
#> [76] 5.923471e+01 9.217307e+01 6.081928e+01 6.736872e+01 6.451640e+01
#> [81] 7.260870e+01 6.334997e+01 5.743027e+01 6.038257e+01 8.342683e+01
#> [86] 5.201926e+01 6.092956e+01 5.441539e+01 6.144519e+01 7.659968e+01
#> [91] 8.750334e+01 7.151385e+01 1.066772e+02 5.359997e+01 7.722627e+01
#> [96] 5.615865e+01 6.766214e+01 1.982847e+02 6.044887e+01 6.271206e+01
#> [101] 5.620680e+01 5.627450e+01 1.719668e+03 7.329397e+01 6.273301e+01
#> [106] 6.126750e+01 5.878408e+01 6.502715e+01 1.258915e+02 4.783583e+02
#> [111] 2.790607e+02 3.711949e+02 9.563700e+01 5.854954e+01 7.770115e+01
#> [116] 7.167065e+01 7.138936e+01 5.395151e+01 6.226772e+01 4.998944e+01
#> [121] 7.365179e+01 7.527417e+01 8.886228e+01 6.963574e+01 5.692371e+01
#> [126] 5.704322e+01 6.462417e+01 1.869993e+02 6.617380e+01 5.907355e+01
#> [131] 5.419411e+01 1.601964e+02 6.738646e+01 6.323709e+01 5.497291e+01
#> [136] 6.680922e+02 5.754773e+01 5.917471e+01 6.363825e+01 9.748716e+01
#> [141] 6.306420e+01 5.579311e+01 7.299703e+01 6.473669e+01 6.051524e+01
#> [146] 6.330457e+01 6.564767e+01 8.121008e+01 1.972454e+02 8.926431e+01
#> [151] 7.768637e+01 8.239524e+01 5.321813e+01 8.118085e+01 1.592877e+02
#> [156] 5.569545e+01 3.460614e+02 6.021584e+01 6.161404e+01 9.082777e+01
#> [161] 8.526470e+01 8.268133e+01 6.350279e+01 7.719627e+01 7.042913e+01
#> [166] 7.332361e+01 5.291495e+01 3.358850e+02 6.214227e+01 1.022138e+02
#> [171] 5.117692e+01 7.097493e+01 5.687378e+01 5.841731e+01 5.752356e+01
#> [176] 7.250331e+01 9.040186e+01 8.672969e+01 5.709283e+01 7.132944e+01
#> [181] 1.155268e+02 5.369641e+01 6.272180e+01 5.874752e+01 5.278985e+01
#> [186] 6.155075e+01 9.253738e+01 1.560543e+02 2.662106e+02 2.282465e+02
#> [191] 1.049245e+02 6.258168e+01 7.567591e+01 6.192265e+01 9.352828e+01
#> [196] 5.618774e+01 5.335920e+01 7.221118e+01 8.043034e+01 1.344026e+02
#> [201] 5.956845e+01 6.121771e+01 5.655313e+01 1.052990e+02 5.805081e+01
#> [206] 6.127505e+01 6.688361e+01 6.993708e+01 8.065011e+01 6.204912e+01
#> [211] 6.414234e+01 6.367195e+01 6.814429e+01 6.207743e+01 7.646545e+01
#> [216] 6.221026e+01 8.178566e+01 1.250229e+02 2.582012e+03 8.783972e+01
#> [221] 5.802705e+01 6.138652e+01 6.304802e+01 6.123323e+01 5.858296e+01
#> [226] 1.067249e+02 7.318027e+01 6.018296e+01 6.257935e+01 7.218348e+01
#> [231] 5.923639e+01 8.209970e+01 7.359394e+01 7.015650e+01 1.029663e+02
#> [236] 4.933399e+01 5.844666e+01 8.924583e+01 2.001971e+02 7.198084e+01
#> [241] 6.305125e+01 6.410042e+01 8.194567e+01 7.183335e+01 1.881462e+02
#> [246] 5.871849e+01 7.342659e+01 7.030569e+01 6.117076e+01 5.680746e+01
#> [251] 1.039343e+02 6.374904e+01 6.207846e+01 8.099096e+01 9.643438e+01
#> [256] 7.273778e+01 5.556628e+01 8.105872e+01 5.853567e+01 6.458148e+01
#> [261] 7.615344e+01 7.474980e+01 7.536544e+01 7.596971e+01 5.680298e+01
#> [266] 5.681623e+01 6.115851e+01 1.738702e+02 7.530079e+01 6.547979e+01
#> [271] 6.114313e+01 5.814686e+01 6.374601e+01 6.153830e+01 7.780389e+01
#> [276] 1.015457e+02 6.550951e+01 9.091306e+01 5.540245e+01 6.462411e+01
#> [281] 6.383766e+01 6.136519e+01 5.932671e+01 7.434665e+01 5.874399e+01
#> [286] 7.221323e+01 6.014369e+01 6.576054e+01 6.228406e+01 5.712380e+01
#> [291] 8.877560e+01 1.534928e+02 5.367515e+01 9.050614e+01 7.425412e+01
#> [296] 5.500456e+01 1.521116e+02 7.928967e+01 6.188083e+01 3.389607e+02
#> [301] 6.127564e+01 7.405402e+01 9.923676e+01 6.125190e+01 6.164271e+01
#> [306] 6.414712e+01 6.219274e+01 6.970312e+01 6.025990e+01 6.116108e+01
#> [311] 5.981686e+01 6.234394e+01 5.735951e+01 7.252570e+01 5.477932e+01
#> [316] 1.019097e+02 6.615449e+01 3.148117e+02 7.382670e+01 6.168218e+01
#> [321] 6.066755e+01 8.041522e+01 5.692844e+01 3.168179e+02 5.568502e+01
#> [326] 5.212485e+01 1.066123e+02 2.015438e+02 6.213376e+01 1.199582e+02
#> [331] 6.005567e+02 5.980110e+01 6.146607e+01 5.803326e+01 5.646105e+01
#> [336] 7.256754e+01 5.818711e+01 6.482535e+01 6.102569e+01 7.168198e+01
#> [341] 5.886716e+01 1.459154e+02 1.674834e+02 1.744400e+02 5.717890e+01
#> [346] 6.466582e+01 6.259770e+01 5.656907e+01 2.018099e+02 5.054999e+01
#> [351] 5.083259e+01 5.492892e+01 9.960297e+01 6.598939e+01 5.886720e+01
#> [356] 8.069788e+01 7.068792e+01 7.014869e+01 2.624202e+02 7.522239e+01
#> [361] 6.656468e+01 1.591900e+02 5.823212e+01 8.073720e+01 6.103184e+01
#> [366] 6.601832e+01 7.632493e+01 1.127694e+02 1.825644e+02 7.179232e+01
#> [371] 5.894842e+01 6.513098e+01 6.785892e+01 7.208026e+01 1.089295e+02
#> [376] 7.185292e+01 1.086734e+02 8.233903e+01 5.568119e+01 5.685721e+01
#> [381] 1.177762e+02 6.261852e+01 5.905322e+01 7.626688e+01 7.496666e+01
#> [386] 5.978359e+01 6.194311e+01 6.386572e+01 6.199893e+01 5.715085e+01
#> [391] 5.959714e+01 6.263950e+01 5.635345e+01 8.674245e+01 6.891411e+01
#> [396] 6.792147e+01 1.294620e+02 2.612369e+02 6.886056e+01 7.164965e+01
#> [401] 4.035259e+02 6.276005e+01 6.079195e+01 6.809216e+01 8.571916e+01
#> [406] 8.940930e+01 7.586592e+01 7.270447e+01 6.166911e+01 7.237672e+01
#> [411] 6.098144e+01 7.756944e+01 6.096968e+01 7.449199e+01 6.382403e+01
#> [416] 5.494872e+01 5.662449e+01 5.692857e+01 5.563508e+01 6.903034e+01
#> [421] 1.205357e+02 2.787622e+02 3.695218e+05 7.467402e+01 4.390497e+01
#> [426] 7.752463e+01 7.260743e+01 6.303123e+01 6.429685e+01 8.351991e+01
#> [431] 8.334671e+01 6.930407e+01 7.146322e+01 6.981601e+01 1.378435e+02
#> [436] 1.720216e+02 5.803188e+01 5.849316e+01 7.937810e+01 1.028790e+02
#> [441] 5.654648e+01 9.249319e+02 5.707049e+01 1.001083e+02 5.493964e+01
#> [446] 5.533566e+01 5.380498e+01 6.457421e+01 6.327480e+01 6.317958e+01
#> [451] 6.706066e+01 1.024632e+02 7.152435e+01 1.137936e+02 1.043452e+02
#> [456] 6.634275e+01 6.618085e+01 6.956460e+01 1.466970e+08 1.876964e+03
#> [461] 1.215178e+02 1.397900e+02 7.502657e+01 7.661552e+01 7.134397e+01
#> [466] 1.045001e+02 7.206405e+01 7.321884e+01 6.326301e+01 8.923133e+01
#> [471] 6.074627e+01 6.631372e+01 8.122382e+01 5.787990e+01 5.970300e+01
#> [476] 7.979207e+01 7.568003e+01 7.034411e+01 6.355587e+01 5.905573e+01
#> [481] 5.990939e+01 6.123711e+01 5.914689e+01 7.163122e+01 9.353779e+01
#> [486] 6.007979e+01 5.889453e+01 6.498205e+01 1.377839e+03 5.881165e+01
#> [491] 6.350282e+01 9.452259e+01 2.221021e+02 7.569185e+01 5.835575e+01
#> [496] 5.598046e+01 5.841122e+01 5.681541e+01 5.943432e+01 8.514976e+01
#> [501] 6.729625e+01 1.393104e+02 5.809176e+01 6.243699e+01 6.508079e+01
#> [506] 6.236240e+01 6.742186e+01 5.936370e+01 5.646944e+01 6.983797e+01
#> [511] 6.522697e+01 5.827362e+01 6.480997e+01 4.500734e+02 7.374324e+02
#> [516] 6.438787e+01 6.361896e+01 8.335973e+01 7.292868e+01 6.303252e+01
#> [521] 5.625444e+01 6.459398e+01 6.400002e+01 6.833484e+01 9.252694e+01
#> [526] 6.762219e+01 7.206915e+01 6.084034e+01 5.556505e+01 4.874147e+01
#> [531] 1.231858e+02 4.723802e+01 7.116542e+01 5.951125e+01 1.075604e+02
#> [536] 6.595569e+01 6.947112e+01 5.843692e+01 6.145166e+01 5.883891e+01
#> [541] 6.612987e+01 6.600226e+01 9.036358e+01 6.041129e+01 5.816679e+01
#> [546] 6.980110e+01 5.421384e+01 6.647152e+01 8.434604e+01 7.118191e+01
#> [551] 5.591001e+01 2.874229e+02 1.028201e+02 9.848949e+01 2.000531e+02
#> [556] 6.316398e+01 6.322897e+01 5.749676e+01 5.676732e+01 5.856448e+01
#> [561] 1.000319e+02 8.376900e+01 6.438183e+01 5.873171e+01 9.693745e+01
#> [566] 6.261031e+01 5.431709e+01 1.812846e+02 1.278338e+02 6.688233e+01
#> [571] 8.228016e+01 8.848916e+01 6.836154e+01 5.907595e+01 8.642414e+01
#> [576] 5.823467e+01 5.939695e+01 7.855000e+01 6.048925e+01 6.661460e+01
#> [581] 6.616986e+01 8.456325e+01 6.079361e+01 1.137518e+02 6.794780e+01
#> [586] 5.538515e+01 8.008050e+01 7.743742e+01 5.822299e+01 9.653154e+01
#> [591] 5.893391e+01 5.957162e+01 5.659056e+01 5.801067e+01 6.020868e+01
#> [596] 6.388025e+01 6.750406e+01 8.969379e+01 9.963059e+01 5.275919e+01
#> [601] 7.106189e+01 9.364310e+01 1.103451e+02 6.338482e+01 5.620706e+01
#> [606] 5.967691e+01 7.070007e+01 6.765860e+01 5.667941e+01 5.169442e+01
#> [611] 9.302274e+01 9.479713e+01 6.838982e+01 2.593333e+03 5.934419e+01
#> [616] 6.380237e+01 7.443398e+01 6.695472e+01 1.008313e+02 8.633462e+01
#> [621] 7.350636e+01 6.954046e+01 6.090687e+01 6.016896e+01 1.116502e+02
#> [626] 5.899141e+01 6.009424e+01 6.752386e+01 5.900920e+01 8.220116e+01
#> [631] 5.809506e+01 6.775254e+01 5.478979e+01 5.798875e+01 5.539574e+01
#> [636] 6.136380e+01 6.166679e+01 1.004126e+02 6.265158e+01 6.339032e+01
#> [641] 5.574499e+01 5.831825e+01 7.658425e+01 6.044959e+01 5.699418e+01
#> [646] 8.533982e+01 5.696727e+01 6.643362e+01 9.136857e+01 7.430449e+01
#> [651] 6.132141e+01 1.415628e+02 1.271391e+02 9.395374e+01 5.505526e+01
#> [656] 5.686317e+01 6.575101e+01 6.392414e+01 6.078567e+01 7.629952e+01
#> [661] 7.471471e+01 2.670674e+02 6.602497e+01 6.048226e+01 5.634526e+01
#> [666] 6.821838e+01 5.943489e+01 6.340464e+01 7.426510e+01 5.493243e+01
#> [671] 6.370188e+01 6.186470e+01 6.349606e+01 6.083230e+01 6.146209e+01
#> [676] 5.726064e+01 6.045611e+01 6.097132e+01 6.319975e+01 5.478788e+01
#> [681] 6.249822e+01 6.921303e+01 5.460653e+01 7.085409e+01 5.980472e+01
#> [686] 7.126733e+01 6.462795e+01 5.900564e+01 8.545101e+01 6.743906e+02
#> [691] 6.704846e+01 1.041532e+02 7.266695e+01 5.691018e+01 7.056193e+01
#> [696] 6.206409e+01 7.822329e+01 8.640239e+01 6.943773e+01 6.006969e+02
#> [701] 5.248130e+01 6.390549e+01 5.885487e+01 8.181340e+01 5.714237e+01
#> [706] 9.486896e+02 2.835441e+02 7.908795e+01 5.618050e+01 6.313284e+01
#> [711] 6.610706e+01 6.125895e+01 8.112120e+01 5.966421e+01 8.909662e+01
#> [716] 6.951880e+01 5.635488e+01 6.623152e+01 6.142021e+01 6.565514e+01
#> [721] 6.120578e+01 1.810153e+02 1.191534e+02 6.032278e+01 5.718793e+01
#> [726] 1.072620e+02 1.554299e+02 6.765676e+01 6.916696e+01 7.431965e+01
#> [731] 5.815677e+01 8.449213e+01 6.583867e+01 5.633072e+01 6.000183e+01
#> [736] 6.739196e+01 6.946396e+01 6.342801e+01 1.174546e+02 5.687567e+01
#> [741] 6.402537e+01 6.382660e+01 5.899282e+01 5.793715e+01 6.425953e+01
#> [746] 1.404873e+02 6.090626e+01 5.513034e+01 5.835051e+01 1.327126e+02
#> [751] 5.457017e+01 6.018373e+01 6.570882e+01 6.578825e+01 6.515394e+01
#> [756] 6.141776e+01 6.577384e+01 5.991868e+01 6.208698e+01 9.359091e+01
#> [761] 6.948753e+01 6.085532e+01 5.922835e+01 9.736101e+01 6.201384e+01
#> [766] 7.901833e+01 7.427654e+01 9.465716e+01 7.403869e+01 7.293401e+01
#> [771] 6.735415e+01 5.499749e+01 5.533445e+01 7.250015e+01 8.810280e+01
#> [776] 6.928060e+01 7.553136e+01 6.805299e+01 5.995823e+01 6.273295e+01
#> [781] 6.511074e+01 1.011719e+02 5.455906e+01 6.032432e+01 6.226765e+01
#> [786] 8.694680e+01 8.373858e+01 6.181420e+01 7.249491e+01 8.824544e+01
#> [791] 5.940338e+01 5.727227e+01 6.443032e+01 6.473497e+01 1.180696e+02
#> [796] 4.459050e+02 5.730293e+01 6.014879e+01 7.066813e+01 6.540435e+01
#> [801] 5.581966e+01 9.772392e+02 5.334218e+01 1.017487e+02 5.901460e+01
#> [806] 6.320803e+01 5.557057e+01 5.723956e+01 6.203894e+01 8.564227e+01
#> [811] 6.574739e+01 9.285773e+01 5.906092e+01 5.305247e+01 6.104994e+01
#> [816] 6.408076e+01 8.651251e+01 1.027969e+02 5.750630e+01 5.916314e+01
#> [821] 6.449011e+01 8.046035e+01 6.694585e+01 1.307476e+02 5.664498e+01
#> [826] 5.407988e+01 8.108785e+01 5.875920e+01 5.782515e+01 6.454834e+01
#> [831] 5.601574e+01 6.670093e+01 5.910370e+01 6.944143e+01 1.467007e+03
#> [836] 5.994143e+01 5.539006e+01 6.467108e+01 5.553459e+01 6.246597e+01
#> [841] 9.291348e+01 5.522811e+01 6.124507e+01 1.164734e+02 7.715267e+01
#> [846] 6.453918e+01 5.782424e+01 6.416391e+01 6.925780e+01 7.027869e+01
#> [851] 2.559199e+02 1.015668e+02 7.522204e+01 6.259794e+01 5.707412e+01
#> [856] 5.943795e+01 5.752351e+01 6.304175e+01 6.645541e+01 5.587849e+01
#> [861] 7.681644e+01 6.722510e+01 6.041482e+01 5.329760e+01 7.569825e+01
#> [866] 5.327479e+01 6.978427e+01 5.977258e+01 7.058903e+01 6.074703e+01
#> [871] 7.730173e+01 7.760300e+01 7.050734e+01 6.407673e+01 9.137905e+01
#> [876] 1.544926e+02 5.637498e+01 5.693777e+01 5.627048e+01 6.150023e+01
#> [881] 6.454188e+01 1.010579e+02 5.525939e+01 6.345360e+01 6.504600e+01
#> [886] 8.336739e+01 6.650094e+01 1.454224e+02 5.745640e+01 8.940629e+01
#> [891] 7.115993e+01 6.837838e+01 6.764201e+01 9.748848e+01 6.407113e+01
#> [896] 6.101948e+01 1.089872e+02 8.152832e+01 5.999231e+01 9.194786e+01
#> [901] 7.058685e+01 7.100387e+01 5.839228e+01 5.068353e+01 8.839945e+01
#> [906] 9.251719e+01 6.756522e+01 1.191730e+02 5.825830e+01 6.800003e+01
#> [911] 6.071320e+01 6.989539e+01 1.370388e+02 1.387538e+02 6.318876e+01
#> [916] 6.827892e+01 5.953897e+01 6.115249e+01 1.416959e+02 7.413196e+01
#> [921] 5.885407e+01 1.840567e+02 6.499691e+01 6.175786e+01 6.512872e+01
#> [926] 6.748893e+01 5.713393e+01 7.908572e+01 1.121321e+02 1.083917e+05
#> [931] 3.667664e+02 6.130812e+01 6.602890e+01 6.103968e+01 5.982414e+01
#> [936] 6.172818e+01 5.822724e+01 5.967893e+01 5.943757e+01 1.329857e+02
#> [941] 5.595950e+01 6.417712e+01 8.497285e+01 5.560513e+01 5.737362e+01
#> [946] 5.768673e+01 6.122626e+01 1.135061e+02 6.087020e+01 5.853122e+01
#> [951] 6.284406e+01 9.519691e+01 5.649028e+01 7.400450e+01 2.106500e+02
#> [956] 8.252403e+01 6.167314e+01 1.674999e+02 7.324549e+01 9.479260e+01
#> [961] 6.680035e+01 1.232072e+02 5.365973e+01 8.435316e+01 6.085888e+01
#> [966] 1.035138e+02 6.646295e+01 6.487825e+01 5.859286e+01 1.035052e+02
#> [971] 8.008271e+01 2.053470e+02 1.665933e+03 1.673108e+02 1.695359e+02
#> [976] 5.844863e+01 9.552561e+01 5.903820e+01 8.610566e+01 6.733044e+01
#> [981] 9.144368e+01 5.964209e+01 6.829870e+01 6.013796e+01 1.213110e+02
#> [986] 1.048029e+03 5.766267e+01 9.801841e+01 8.557152e+01 4.957760e+01
#> [991] 6.170138e+01 5.559639e+01 2.954021e+02 8.369758e+01 6.280356e+01
#> [996] 5.859026e+01 7.077654e+01 5.724957e+01 6.735217e+01 4.940820e+01