![[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 GeneralModel
doseFunction(model, ...)
# S4 method for ModelPseudo
doseFunction(model, ...)
# S4 method for LogisticLogNormalOrdinal
doseFunction(model, grade, ...)Arguments
- model
(
GeneralModelorModelPseudo)
the model.- ...
model specific parameters.
- grade
(
integer)
the toxicity grade for which the dose function is required
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.293732e+01 7.472328e+01 4.063285e+02 8.822810e+01 9.295479e+01
#> [6] 8.113170e+01 1.214556e+02 1.800308e+02 2.717524e+02 8.129628e+01
#> [11] 7.516753e+01 2.014847e+02 1.575024e+02 7.895441e+01 8.081570e+01
#> [16] 6.159598e+01 1.803469e+02 1.657403e+02 9.122808e+01 8.631712e+01
#> [21] 1.168641e+02 1.225334e+02 9.842611e+01 9.014967e+01 7.340275e+01
#> [26] 7.755605e+01 1.267725e+02 3.485873e+02 3.267553e+02 1.997349e+03
#> [31] 7.993418e+01 6.336397e+01 7.118573e+01 8.842151e+01 2.044683e+02
#> [36] 4.127848e+02 1.080945e+02 3.662112e+02 2.724302e+02 4.833247e+02
#> [41] 3.516097e+02 1.262178e+02 1.091739e+02 7.935587e+01 8.420038e+01
#> [46] 7.023228e+01 1.035512e+02 2.921311e+02 2.471110e+05 3.606052e+02
#> [51] 7.703772e+01 7.607436e+01 1.104535e+02 7.896314e+01 8.675551e+01
#> [56] 9.235756e+01 7.387744e+01 3.535280e+02 6.830161e+01 2.704661e+02
#> [61] 9.487532e+01 7.905392e+01 6.967285e+01 1.647693e+02 6.914029e+01
#> [66] 8.899217e+01 8.126441e+01 1.729898e+03 2.113755e+02 5.107217e+02
#> [71] 7.596176e+01 7.179265e+01 2.585496e+02 7.182347e+01 9.081209e+01
#> [76] 6.786687e+01 1.098155e+03 2.085101e+03 1.564800e+03 1.248346e+02
#> [81] 6.864314e+01 7.094992e+01 8.112140e+01 8.535036e+01 2.045517e+02
#> [86] 2.416547e+02 1.263768e+02 1.500285e+02 1.181752e+02 6.652044e+01
#> [91] 6.828949e+01 1.014080e+02 1.520325e+02 9.713380e+01 3.643880e+02
#> [96] 3.104328e+02 6.289760e+01 6.485401e+03 7.701826e+01 1.408405e+02
#> [101] 6.250235e+01 1.374756e+02 6.439078e+01 7.256889e+01 9.531237e+01
#> [106] 9.762754e+01 1.116826e+02 6.255115e+01 9.209966e+01 6.872343e+01
#> [111] 8.346596e+01 1.290638e+02 2.293784e+03 1.225527e+02 1.164454e+02
#> [116] 6.997944e+01 7.470897e+01 3.324348e+02 6.550206e+01 6.971675e+01
#> [121] 1.213902e+02 6.967839e+01 6.576415e+01 1.047544e+02 7.864652e+01
#> [126] 7.005541e+01 6.474083e+01 1.035920e+02 3.276756e+02 8.380610e+01
#> [131] 7.533122e+01 6.634472e+01 7.764499e+01 8.878790e+02 1.624034e+02
#> [136] 6.225442e+01 3.031866e+02 1.033441e+02 8.881061e+01 1.615895e+02
#> [141] 6.689411e+01 7.522861e+01 7.736045e+01 6.867980e+01 7.535080e+01
#> [146] 6.236076e+01 8.977880e+01 1.469256e+02 1.243142e+02 1.255014e+02
#> [151] 1.112047e+02 6.266245e+01 8.928944e+01 1.151437e+02 8.807068e+01
#> [156] 8.822658e+01 7.890879e+01 6.296141e+01 6.046143e+01 5.466926e+02
#> [161] 8.349074e+01 4.250878e+02 1.434689e+04 1.220656e+02 9.533173e+01
#> [166] 7.973502e+01 7.005828e+01 1.711766e+02 1.076149e+02 1.707305e+02
#> [171] 2.807670e+03 6.603191e+01 8.183633e+01 2.635461e+02 6.556748e+01
#> [176] 1.819158e+02 1.195195e+02 9.598960e+01 8.340769e+01 9.672410e+01
#> [181] 7.345276e+01 8.654069e+01 8.277042e+01 2.201745e+02 1.295534e+02
#> [186] 3.035358e+02 1.689634e+02 6.179814e+01 1.577603e+02 7.390344e+01
#> [191] 7.503493e+01 8.948786e+01 7.396980e+01 1.089257e+02 1.188691e+02
#> [196] 6.987607e+01 8.786079e+01 7.284396e+01 1.532608e+02 1.313739e+02
#> [201] 1.206010e+02 1.615653e+02 6.880084e+01 1.825900e+02 3.802076e+02
#> [206] 1.026807e+02 7.536190e+01 7.400961e+01 1.074252e+02 7.320284e+01
#> [211] 7.115590e+01 9.542177e+01 2.638690e+02 6.887344e+01 8.069066e+01
#> [216] 7.823809e+01 8.899935e+01 8.992997e+01 1.052868e+02 8.114569e+01
#> [221] 1.359148e+02 2.038658e+02 1.123833e+02 9.351801e+01 1.605926e+02
#> [226] 2.363509e+02 7.509383e+01 7.287074e+01 8.138933e+01 3.078663e+02
#> [231] 6.797300e+01 7.068773e+01 6.352865e+01 1.604111e+02 1.036796e+02
#> [236] 6.266028e+01 8.617142e+01 1.159396e+02 6.703695e+01 9.661960e+01
#> [241] 1.475566e+02 6.762664e+01 1.319390e+02 7.015366e+01 1.070314e+02
#> [246] 1.149326e+02 1.523585e+02 1.085236e+02 2.209774e+02 6.719343e+01
#> [251] 6.408886e+01 8.871572e+01 7.084525e+01 9.518029e+02 8.100511e+01
#> [256] 9.782468e+01 8.838508e+01 7.140246e+01 7.913321e+01 1.129113e+02
#> [261] 1.254790e+03 2.067932e+02 8.007655e+01 6.809994e+01 1.358389e+02
#> [266] 4.575745e+02 5.350791e+02 8.836542e+01 1.878045e+02 9.520286e+01
#> [271] 7.153175e+01 8.546887e+01 6.684558e+01 8.448223e+01 3.317389e+02
#> [276] 1.466778e+02 7.829854e+01 7.774109e+01 1.095682e+02 6.479566e+01
#> [281] 7.782471e+01 7.397105e+01 9.056114e+01 6.511593e+01 8.086521e+01
#> [286] 9.142032e+01 6.027380e+01 1.058506e+02 3.419389e+02 1.080248e+02
#> [291] 7.397134e+01 1.051854e+02 7.201871e+01 7.387708e+01 7.293290e+01
#> [296] 7.679247e+01 6.808511e+01 1.315168e+02 7.977762e+01 7.022998e+01
#> [301] 9.837098e+01 7.305776e+01 6.637706e+01 7.613288e+01 6.460624e+01
#> [306] 6.802044e+01 7.121856e+01 6.676397e+01 7.487707e+01 9.388451e+01
#> [311] 6.555087e+01 7.358283e+01 7.414757e+01 1.133459e+03 1.006123e+02
#> [316] 7.421131e+01 7.791229e+01 2.327466e+02 4.014782e+02 7.711119e+01
#> [321] 7.154084e+01 7.228584e+01 7.262521e+01 8.255918e+01 8.593481e+01
#> [326] 3.847682e+02 7.888796e+01 6.242444e+01 1.201480e+02 7.787583e+01
#> [331] 1.732531e+02 6.635833e+01 9.086261e+01 7.592969e+01 1.343950e+02
#> [336] 1.364261e+02 7.400731e+01 7.308213e+01 1.204399e+02 7.161057e+01
#> [341] 6.922727e+01 6.830513e+01 6.990612e+01 7.134277e+01 2.045965e+02
#> [346] 7.564540e+01 1.438484e+02 6.499126e+02 2.592839e+03 7.788985e+01
#> [351] 6.960700e+01 7.667338e+01 1.433576e+02 1.516527e+02 6.994603e+01
#> [356] 2.618628e+02 4.080491e+02 1.054654e+05 8.822194e+01 6.566745e+01
#> [361] 8.457592e+01 1.113693e+02 1.011755e+02 7.732232e+01 8.108537e+01
#> [366] 1.022381e+02 1.193248e+02 1.774424e+04 8.436342e+01 8.440142e+01
#> [371] 1.123842e+02 7.108896e+01 1.212300e+02 7.801519e+01 9.295792e+01
#> [376] 1.879760e+02 1.042920e+02 8.443276e+01 8.224034e+01 1.102746e+02
#> [381] 7.388798e+01 1.083988e+02 8.990031e+01 6.454799e+01 7.133607e+01
#> [386] 7.188710e+01 1.101318e+02 6.662599e+01 8.379382e+01 2.840726e+03
#> [391] 1.433940e+02 1.048502e+02 1.131610e+02 7.460512e+01 1.176093e+02
#> [396] 7.985495e+01 2.291690e+03 6.616168e+01 9.602745e+01 8.049483e+01
#> [401] 6.253066e+01 6.734393e+01 7.188287e+01 1.103603e+02 1.505116e+02
#> [406] 8.390228e+01 1.190148e+02 6.996597e+01 1.100051e+03 1.304612e+02
#> [411] 6.776542e+01 7.581152e+01 2.821316e+02 1.650218e+02 6.416961e+01
#> [416] 1.053333e+02 9.663479e+01 1.258943e+02 8.505179e+01 7.800192e+01
#> [421] 8.690995e+01 1.077594e+02 8.283617e+01 7.494307e+01 7.426037e+01
#> [426] 7.728844e+01 6.351371e+01 7.708099e+01 1.775608e+02 2.220976e+02
#> [431] 6.293360e+01 1.499129e+03 3.367619e+02 4.518726e+02 7.466714e+01
#> [436] 9.446724e+01 1.054320e+02 3.416661e+03 9.561688e+01 1.362424e+02
#> [441] 6.669596e+01 7.788350e+01 8.725265e+01 6.474408e+01 8.657256e+01
#> [446] 6.555413e+01 6.475331e+01 7.424069e+01 1.143760e+02 1.896664e+02
#> [451] 8.607080e+01 7.536749e+01 9.606224e+01 9.205982e+01 6.556985e+01
#> [456] 8.040009e+01 6.942089e+01 3.675378e+02 7.305838e+02 1.402227e+02
#> [461] 1.411695e+02 1.114370e+02 8.048905e+01 1.108211e+02 9.704158e+01
#> [466] 8.251619e+01 8.235584e+01 8.735988e+01 6.900207e+01 1.219325e+02
#> [471] 9.286778e+01 6.892364e+01 7.246959e+01 8.113300e+01 6.958964e+01
#> [476] 3.620179e+02 6.927144e+01 8.131980e+01 7.508606e+01 7.073189e+01
#> [481] 7.522109e+01 6.674642e+01 6.249979e+02 7.258507e+01 1.060988e+02
#> [486] 2.595025e+02 4.154148e+03 6.321571e+02 6.926014e+01 9.954787e+01
#> [491] 5.773356e+02 8.741730e+01 6.869435e+01 8.715090e+01 7.064994e+01
#> [496] 6.278748e+01 7.242870e+01 7.625401e+01 1.405795e+02 1.967149e+03
#> [501] 9.343045e+01 1.076609e+02 7.506453e+01 6.579748e+01 6.920599e+01
#> [506] 6.799117e+01 8.117431e+01 1.888350e+02 1.034722e+02 1.120430e+02
#> [511] 1.089332e+02 1.498021e+02 7.139917e+01 1.066702e+02 6.967630e+01
#> [516] 8.149248e+01 6.462814e+01 6.947027e+01 2.250090e+02 7.884526e+01
#> [521] 8.744269e+01 1.145779e+02 9.139947e+01 1.289848e+02 9.982267e+01
#> [526] 8.415115e+01 7.487630e+01 9.424168e+01 2.005116e+02 8.020250e+01
#> [531] 1.074972e+02 7.067224e+01 1.290164e+02 7.462489e+01 8.639301e+01
#> [536] 7.798241e+01 9.569597e+01 2.536812e+02 7.285010e+01 6.448697e+01
#> [541] 6.672836e+01 8.461922e+01 8.863292e+01 6.848172e+01 2.273363e+02
#> [546] 6.716525e+01 6.781580e+01 8.097244e+01 8.307359e+01 1.128464e+02
#> [551] 1.278724e+02 7.076454e+01 1.305211e+02 7.683655e+01 7.322982e+01
#> [556] 6.798174e+01 9.611305e+01 7.293789e+01 2.165447e+03 8.941377e+02
#> [561] 2.777000e+03 1.283842e+02 1.716588e+02 8.393449e+01 1.249844e+02
#> [566] 6.498136e+01 1.632894e+02 7.706011e+01 8.487474e+01 9.602467e+01
#> [571] 7.083707e+01 1.246977e+03 2.732967e+02 1.091006e+02 6.720708e+01
#> [576] 2.030226e+02 1.282787e+02 6.717127e+01 6.362020e+01 8.108666e+01
#> [581] 6.911048e+01 7.622640e+01 1.148653e+02 9.369395e+01 1.534292e+02
#> [586] 7.372382e+01 2.940213e+02 8.923829e+01 6.602821e+01 6.623422e+01
#> [591] 8.380563e+03 7.991405e+01 1.514691e+02 8.592883e+01 6.283505e+02
#> [596] 2.471753e+02 8.364931e+01 9.171681e+01 9.246885e+01 8.015742e+01
#> [601] 7.338597e+01 1.146066e+02 8.743374e+01 1.428435e+02 1.001516e+02
#> [606] 9.483453e+01 1.044764e+02 6.777740e+01 6.172868e+01 7.893793e+01
#> [611] 7.750782e+01 7.765846e+01 3.226235e+03 9.257734e+01 8.714518e+01
#> [616] 8.921799e+01 7.547914e+01 1.138798e+02 1.013884e+02 2.506718e+02
#> [621] 2.497773e+02 7.848158e+01 6.870878e+01 7.391417e+01 7.348673e+01
#> [626] 8.700160e+01 6.926840e+01 6.922758e+01 6.523706e+01 8.791755e+01
#> [631] 2.729081e+06 1.008091e+03 5.439911e+02 6.416607e+01 1.017094e+02
#> [636] 6.839453e+01 1.475062e+02 6.815294e+01 6.465213e+01 6.472405e+01
#> [641] 1.308088e+02 1.099207e+02 6.332721e+01 7.333247e+01 7.915080e+01
#> [646] 6.813777e+01 6.340112e+01 2.176725e+02 3.510876e+02 2.978491e+12
#> [651] 2.642961e+02 7.164995e+01 7.422670e+01 1.028426e+02 1.401024e+02
#> [656] 9.211802e+01 7.908395e+01 7.692078e+01 7.838944e+01 7.549900e+01
#> [661] 1.461714e+02 1.549190e+02 1.207507e+02 6.398182e+01 7.410251e+01
#> [666] 1.188413e+02 9.159508e+01 6.252377e+03 1.151852e+03 9.242667e+01
#> [671] 8.766879e+01 1.163243e+02 3.485784e+02 7.903065e+01 1.208668e+02
#> [676] 6.833143e+01 6.584139e+01 7.888006e+01 3.829602e+02 7.881048e+02
#> [681] 1.888891e+02 7.214727e+01 7.334496e+01 8.430787e+01 7.748631e+01
#> [686] 2.216145e+02 7.236554e+01 7.744825e+01 4.547633e+02 6.216735e+03
#> [691] 7.034626e+01 1.126281e+02 1.050901e+02 2.421625e+03 9.271546e+01
#> [696] 4.842936e+03 2.623487e+02 9.108641e+02 2.961860e+03 1.177811e+02
#> [701] 1.181567e+02 1.219189e+02 7.370645e+01 6.570537e+01 6.476789e+01
#> [706] 1.182706e+02 6.576910e+01 6.577622e+01 9.835541e+09 7.239986e+01
#> [711] 7.227295e+01 1.024030e+02 1.266336e+02 4.552174e+02 1.087566e+02
#> [716] 8.161572e+01 9.576671e+01 7.788116e+01 6.880183e+01 7.395470e+01
#> [721] 8.252934e+01 9.754692e+01 9.326956e+01 1.295212e+02 7.664937e+01
#> [726] 6.961899e+01 6.490467e+02 8.446385e+01 6.894852e+01 7.290828e+01
#> [731] 6.331154e+01 1.670322e+02 6.584420e+01 6.136170e+01 8.249502e+01
#> [736] 1.030396e+02 8.737005e+01 6.746305e+01 1.275639e+02 7.098541e+01
#> [741] 8.250211e+01 7.051530e+01 1.228509e+02 6.317155e+01 8.453323e+01
#> [746] 8.390735e+01 1.511098e+02 6.413954e+01 1.660734e+02 1.374750e+02
#> [751] 8.284884e+01 2.210906e+02 1.110813e+02 9.013381e+01 1.060801e+02
#> [756] 7.584127e+01 4.228039e+02 9.940178e+01 6.901853e+01 9.487175e+01
#> [761] 7.599069e+01 3.447642e+02 9.525486e+01 9.064056e+01 8.410833e+01
#> [766] 8.177176e+01 2.658371e+02 3.790149e+02 2.211468e+02 4.008660e+02
#> [771] 7.715451e+01 6.225836e+01 8.398705e+01 7.652535e+01 1.025963e+02
#> [776] 7.636880e+01 8.066236e+01 7.127757e+01 7.811726e+01 6.920160e+01
#> [781] 3.340310e+02 6.662128e+01 5.097135e+03 9.057630e+01 6.994082e+01
#> [786] 8.710203e+01 7.943194e+02 7.958523e+01 9.627050e+01 8.277627e+01
#> [791] 6.813077e+01 6.482821e+01 7.255848e+01 2.075917e+02 6.747610e+01
#> [796] 9.112598e+01 7.126157e+01 6.345110e+01 7.094969e+01 6.501310e+01
#> [801] 6.913367e+01 6.296444e+01 4.145744e+02 7.294537e+01 9.099910e+01
#> [806] 8.940953e+01 8.856775e+01 1.085144e+02 6.753051e+01 1.286534e+02
#> [811] 6.928353e+01 1.058522e+02 1.038054e+02 7.525829e+01 1.122147e+02
#> [816] 1.293124e+02 8.328995e+01 8.633172e+01 8.555077e+01 7.763049e+01
#> [821] 1.265667e+02 1.369555e+02 2.193849e+02 1.040503e+02 1.530282e+02
#> [826] 9.656102e+01 9.129639e+01 7.830115e+01 2.644794e+02 2.807066e+02
#> [831] 1.428392e+02 2.272879e+02 7.526338e+01 8.764798e+01 1.015035e+02
#> [836] 1.006310e+02 6.884980e+01 1.164135e+02 6.893201e+01 7.375770e+01
#> [841] 1.047085e+02 1.882856e+02 7.217964e+01 6.574756e+01 1.119520e+02
#> [846] 1.351267e+02 7.739825e+01 7.081999e+01 8.936292e+01 1.757157e+02
#> [851] 3.255632e+02 1.275380e+02 9.698007e+01 6.894761e+01 1.310723e+02
#> [856] 8.041579e+01 6.839267e+01 2.569734e+02 1.346939e+02 1.226509e+03
#> [861] 7.123514e+01 4.451517e+02 1.018107e+02 4.903065e+02 1.716978e+02
#> [866] 7.814327e+01 8.546549e+01 9.626218e+01 6.681861e+01 6.682943e+01
#> [871] 6.709475e+01 8.526134e+01 2.534309e+02 3.066141e+02 1.371490e+02
#> [876] 7.062599e+01 8.577968e+01 9.708876e+01 8.521467e+01 8.952073e+01
#> [881] 8.643881e+01 2.012446e+02 8.880657e+01 9.173687e+01 1.041512e+02
#> [886] 6.567107e+01 9.886862e+01 1.036631e+02 4.428624e+02 1.323945e+02
#> [891] 7.212259e+01 2.739452e+02 3.200065e+02 2.272460e+02 1.075773e+02
#> [896] 1.415299e+02 7.319792e+01 1.294229e+02 7.076012e+01 6.806606e+01
#> [901] 1.348042e+02 1.176085e+02 1.043181e+02 4.537493e+03 1.579722e+02
#> [906] 1.419754e+02 8.054663e+01 9.344211e+01 9.569949e+01 1.278546e+02
#> [911] 6.142432e+01 8.875322e+01 7.157323e+01 6.266166e+01 8.244128e+01
#> [916] 9.912450e+01 8.478800e+01 1.592937e+02 2.242221e+02 7.007857e+01
#> [921] 7.626251e+01 7.045964e+01 9.389911e+01 1.169687e+02 7.320945e+01
#> [926] 6.541601e+01 8.798191e+01 7.763227e+01 8.752979e+01 8.322523e+01
#> [931] 7.326929e+01 2.926207e+02 6.149776e+01 7.289355e+01 3.740953e+02
#> [936] 9.881166e+02 1.268597e+02 3.799072e+02 1.554213e+02 1.468826e+02
#> [941] 4.469570e+02 6.771639e+01 4.299143e+02 1.020299e+02 8.221919e+01
#> [946] 7.062705e+01 1.292297e+02 6.455160e+02 1.548960e+02 1.361687e+02
#> [951] 7.000707e+01 6.885831e+01 7.052760e+01 6.830676e+01 9.873094e+01
#> [956] 1.308835e+02 6.905023e+01 7.886537e+01 2.564397e+02 1.453932e+02
#> [961] 7.504213e+01 7.636090e+01 1.044476e+02 8.155772e+01 7.918273e+01
#> [966] 9.266461e+01 9.289715e+01 1.718613e+02 7.324029e+01 8.413825e+01
#> [971] 1.116264e+02 7.543429e+01 1.098378e+02 2.469716e+02 1.589115e+02
#> [976] 8.155026e+01 8.390755e+01 1.043726e+02 6.896569e+01 2.991655e+03
#> [981] 1.118382e+02 6.273980e+01 6.842387e+01 1.173037e+02 1.008692e+02
#> [986] 6.777895e+01 6.222932e+01 6.795762e+01 9.722765e+01 8.971749e+01
#> [991] 1.124500e+02 8.474998e+01 8.883785e+01 8.578052e+01 7.429611e+01
#> [996] 7.783596e+01 9.682662e+01 2.306855e+02 6.898650e+01 6.954046e+01
doseFunction(model, alpha2 = samples@data$alpha2, beta = samples@data$beta, grade = 2L)(x = 0.25)
#> [1] 6.166616e+01 6.272767e+01 1.613896e+02 6.274464e+01 6.064978e+01
#> [6] 6.371446e+01 7.547644e+01 8.855888e+01 1.177778e+02 6.610311e+01
#> [11] 5.229946e+01 9.496709e+01 1.025823e+02 5.819889e+01 6.100643e+01
#> [16] 5.451126e+01 7.911931e+01 8.148291e+01 5.877700e+01 6.961519e+01
#> [21] 6.484467e+01 5.886881e+01 7.509199e+01 7.029466e+01 6.116145e+01
#> [26] 5.849555e+01 8.503840e+01 1.397334e+02 1.501034e+02 2.856198e+02
#> [31] 6.579348e+01 5.844474e+01 6.465693e+01 6.719642e+01 1.170925e+02
#> [36] 1.423861e+02 7.358947e+01 1.432218e+02 1.098561e+02 1.387641e+02
#> [41] 1.352401e+02 7.444393e+01 6.642682e+01 5.857453e+01 6.817279e+01
#> [46] 6.053618e+01 7.142997e+01 1.047276e+02 2.627165e+03 1.249939e+02
#> [51] 6.041209e+01 5.733894e+01 7.526119e+01 5.329991e+01 6.081406e+01
#> [56] 7.163310e+01 5.862304e+01 1.475288e+02 6.055224e+01 1.340716e+02
#> [61] 6.524800e+01 6.406265e+01 5.896397e+01 5.891437e+01 4.472066e+01
#> [66] 7.148165e+01 6.269812e+01 1.580953e+02 8.668503e+01 1.590298e+02
#> [71] 5.557021e+01 6.076428e+01 1.337439e+02 6.555589e+01 7.414437e+01
#> [76] 5.852825e+01 1.263348e+02 2.682222e+02 2.367450e+02 8.689120e+01
#> [81] 5.713283e+01 5.790182e+01 6.331365e+01 6.544584e+01 9.034015e+01
#> [86] 6.774956e+01 7.126238e+01 1.080320e+02 7.592893e+01 5.543653e+01
#> [91] 5.758122e+01 6.178648e+01 9.636360e+01 7.505029e+01 8.274844e+01
#> [96] 1.035646e+02 5.485044e+01 5.997775e+02 6.396347e+01 7.340632e+01
#> [101] 5.591224e+01 8.312463e+01 5.795288e+01 6.324081e+01 6.858362e+01
#> [106] 6.361478e+01 7.162657e+01 5.759415e+01 6.952222e+01 5.726891e+01
#> [111] 6.170763e+01 9.108114e+01 4.503156e+02 8.398238e+01 8.014330e+01
#> [116] 6.198962e+01 6.482962e+01 1.232193e+02 5.706479e+01 5.887995e+01
#> [121] 7.661230e+01 5.544025e+01 6.008812e+01 7.561046e+01 6.287680e+01
#> [126] 5.440557e+01 5.612850e+01 8.209198e+01 1.133400e+02 6.681937e+01
#> [131] 6.201213e+01 5.713220e+01 6.582742e+01 2.190840e+02 6.046856e+01
#> [136] 5.209193e+01 1.425249e+02 8.219242e+01 6.936108e+01 7.908026e+01
#> [141] 5.300213e+01 6.223218e+01 6.685108e+01 6.095172e+01 6.250957e+01
#> [146] 5.767137e+01 7.425429e+01 1.004070e+02 7.961573e+01 5.148972e+01
#> [151] 6.393819e+01 5.607614e+01 6.989968e+01 6.242515e+01 6.807494e+01
#> [156] 6.403502e+01 6.820588e+01 5.625397e+01 5.557489e+01 2.056985e+02
#> [161] 6.639867e+01 7.231204e+01 4.746986e+02 6.632455e+01 6.871039e+01
#> [166] 6.436448e+01 6.119836e+01 1.085158e+02 7.722412e+01 1.199448e+02
#> [171] 8.186471e+02 5.984021e+01 6.161465e+01 1.237997e+02 5.595658e+01
#> [176] 9.555766e+01 7.728002e+01 7.108792e+01 6.291335e+01 6.793169e+01
#> [181] 6.217068e+01 6.653884e+01 6.455480e+01 8.723752e+01 7.333592e+01
#> [186] 9.596596e+01 4.080783e+01 4.998082e+01 8.413519e+01 5.931954e+01
#> [191] 6.262887e+01 6.218482e+01 6.062026e+01 6.406595e+01 8.861746e+01
#> [196] 6.026015e+01 6.140239e+01 6.216841e+01 8.839810e+01 7.818113e+01
#> [201] 5.362891e+01 8.106817e+01 5.739188e+01 8.432729e+01 1.039742e+02
#> [206] 7.150312e+01 6.063551e+01 5.633735e+01 6.592514e+01 5.897385e+01
#> [211] 6.260193e+01 7.293045e+01 9.636090e+01 5.629522e+01 6.810541e+01
#> [216] 5.777316e+01 6.187731e+01 6.022020e+01 7.644430e+01 6.519512e+01
#> [221] 8.583308e+01 9.081244e+01 7.376730e+01 7.262770e+01 8.757835e+01
#> [226] 8.577796e+01 6.546202e+01 6.053814e+01 6.705829e+01 1.672305e+02
#> [231] 6.173829e+01 5.772660e+01 5.644785e+01 7.197474e+01 6.750975e+01
#> [236] 5.660854e+01 6.981632e+01 7.927338e+01 5.907278e+01 7.032952e+01
#> [241] 8.203590e+01 5.585481e+01 6.779257e+01 5.995531e+01 5.930973e+01
#> [246] 8.286547e+01 9.959730e+01 7.233039e+01 1.204272e+02 5.304811e+01
#> [251] 5.716907e+01 7.457997e+01 6.347968e+01 1.751616e+02 7.022511e+01
#> [256] 6.647446e+01 6.949053e+01 5.811249e+01 5.634743e+01 6.803837e+01
#> [261] 1.992248e+02 9.215445e+01 5.180519e+01 5.656654e+01 6.787560e+01
#> [266] 1.291089e+02 1.535614e+02 7.407191e+01 9.646543e+01 6.661531e+01
#> [271] 5.805893e+01 5.380150e+01 5.957791e+01 5.973608e+01 1.612463e+02
#> [276] 8.104052e+01 5.889584e+01 6.268806e+01 7.210890e+01 5.496860e+01
#> [281] 6.227626e+01 6.376101e+01 6.806352e+01 5.938705e+01 5.467708e+01
#> [286] 5.894196e+01 5.355967e+01 7.296932e+01 1.738319e+02 8.447933e+01
#> [291] 6.350069e+01 8.061942e+01 5.910874e+01 6.164721e+01 6.003044e+01
#> [296] 6.227806e+01 5.196782e+01 8.881304e+01 6.586120e+01 6.144075e+01
#> [301] 6.647694e+01 5.708748e+01 5.666620e+01 6.299938e+01 5.947797e+01
#> [306] 5.763102e+01 6.321103e+01 5.902518e+01 6.445691e+01 6.839782e+01
#> [311] 5.786556e+01 6.103424e+01 4.229993e+01 3.255071e+01 7.273030e+01
#> [316] 6.035626e+01 6.792004e+01 1.194177e+02 1.140072e+02 6.357347e+01
#> [321] 6.209316e+01 6.104655e+01 5.900157e+01 6.838148e+01 7.082931e+01
#> [326] 1.476851e+02 6.259507e+01 5.814792e+01 8.073844e+01 5.825930e+01
#> [331] 9.234718e+01 5.358229e+01 5.750041e+01 6.194647e+01 9.202476e+01
#> [336] 8.965512e+01 6.243741e+01 5.694734e+01 7.172048e+01 6.090907e+01
#> [341] 6.018765e+01 6.007797e+01 6.082542e+01 5.970080e+01 8.925274e+01
#> [346] 6.105822e+01 9.721788e+01 7.641267e+01 1.904001e+02 6.497970e+01
#> [351] 6.043127e+01 6.720183e+01 1.008501e+02 1.041269e+02 5.587288e+01
#> [356] 1.205868e+02 8.601598e+01 1.123422e+04 5.626622e+01 5.920960e+01
#> [361] 6.318747e+01 5.761293e+01 6.758339e+01 6.181922e+01 5.895373e+01
#> [366] 6.895156e+01 6.685052e+01 8.441761e+02 7.171148e+01 6.651114e+01
#> [371] 7.941139e+01 5.372410e+01 7.697293e+01 6.362318e+01 5.408895e+01
#> [376] 1.065313e+02 5.696449e+01 6.186718e+01 5.876557e+01 6.731921e+01
#> [381] 6.192001e+01 8.499530e+01 6.997736e+01 5.725769e+01 5.746358e+01
#> [386] 5.960724e+01 8.006551e+01 5.937864e+01 6.880978e+01 4.245284e+02
#> [391] 6.942636e+01 7.255438e+01 7.223202e+01 6.040835e+01 7.445031e+01
#> [396] 6.183005e+01 2.157062e+02 5.642127e+01 7.144864e+01 6.671841e+01
#> [401] 5.777082e+01 6.099962e+01 5.995780e+01 7.307325e+01 8.040093e+01
#> [406] 6.833588e+01 8.287097e+01 6.044927e+01 2.797953e+02 7.572017e+01
#> [411] 5.640833e+01 6.027076e+01 9.365280e+01 8.547549e+01 5.048618e+01
#> [416] 6.916311e+01 6.420808e+01 9.188574e+01 6.333089e+01 6.308308e+01
#> [421] 6.373639e+01 7.621757e+01 5.973693e+01 6.316186e+01 6.314009e+01
#> [426] 6.228715e+01 5.390317e+01 6.033982e+01 9.869039e+01 1.154293e+02
#> [431] 5.558795e+01 3.828609e+02 6.714010e+01 1.132671e+02 5.746862e+01
#> [436] 5.895620e+01 6.454907e+01 4.775241e+02 7.704117e+01 1.008607e+02
#> [441] 6.166734e+01 6.986372e+01 7.631474e+01 5.901634e+01 7.349846e+01
#> [446] 5.952449e+01 5.817861e+01 6.354043e+01 8.416470e+01 1.036289e+02
#> [451] 6.882235e+01 6.469199e+01 7.192701e+01 6.730192e+01 5.835671e+01
#> [456] 6.547880e+01 5.554929e+01 1.694899e+02 2.683121e+02 9.247396e+01
#> [461] 8.992679e+01 7.500381e+01 6.672893e+01 8.636290e+01 5.557918e+01
#> [466] 6.693734e+01 6.293500e+01 5.829946e+01 5.426766e+01 8.174653e+01
#> [471] 5.254897e+01 6.036758e+01 6.162072e+01 6.713681e+01 5.377536e+01
#> [476] 8.289863e+01 5.798106e+01 5.658337e+01 6.414093e+01 6.243319e+01
#> [481] 6.368792e+01 5.180247e+01 7.918583e+01 5.481537e+01 5.821372e+01
#> [486] 9.975828e+01 3.839229e+02 2.434104e+02 5.357953e+01 7.221722e+01
#> [491] 1.431999e+02 5.743567e+01 5.916072e+01 7.063249e+01 5.795681e+01
#> [496] 5.509036e+01 6.367895e+01 6.572263e+01 1.031608e+02 1.371715e+02
#> [501] 6.999149e+01 7.096497e+01 6.083051e+01 5.692351e+01 5.903657e+01
#> [506] 5.703782e+01 6.441448e+01 9.951252e+01 5.554696e+01 5.739146e+01
#> [511] 7.473701e+01 6.614448e+01 5.494061e+01 6.921177e+01 5.783673e+01
#> [516] 7.074689e+01 5.902063e+01 6.155010e+01 9.606980e+01 6.369188e+01
#> [521] 6.548931e+01 8.085709e+01 6.413451e+01 9.010860e+01 7.742464e+01
#> [526] 6.566837e+01 6.327471e+01 7.041533e+01 1.006520e+02 5.653891e+01
#> [531] 7.271557e+01 5.565499e+01 7.108218e+01 6.007324e+01 7.148691e+01
#> [536] 6.746875e+01 7.777631e+01 1.273016e+02 6.097654e+01 5.769791e+01
#> [541] 5.891102e+01 7.327023e+01 6.522213e+01 5.659102e+01 9.760098e+01
#> [546] 5.017656e+01 5.371771e+01 5.355676e+01 6.658442e+01 5.394595e+01
#> [551] 6.638257e+01 6.204495e+01 9.525936e+01 6.179709e+01 5.699409e+01
#> [556] 5.561765e+01 7.062710e+01 5.641408e+01 4.207005e+02 2.544768e+02
#> [561] 6.525494e+02 8.126782e+01 7.906937e+01 5.956151e+01 7.826121e+01
#> [566] 5.418018e+01 1.001225e+02 5.837996e+01 6.579558e+01 7.093549e+01
#> [571] 5.733985e+01 5.388320e+01 1.124081e+02 6.427119e+01 5.809660e+01
#> [576] 1.169354e+02 8.532910e+01 5.891835e+01 5.563898e+01 6.228147e+01
#> [581] 5.895379e+01 6.147209e+01 6.086103e+01 7.520456e+01 8.365279e+01
#> [586] 6.480301e+01 9.367726e+01 6.979230e+01 6.084854e+01 5.957724e+01
#> [591] 1.266569e+03 6.863073e+01 8.959754e+01 7.045489e+01 1.017540e+02
#> [596] 1.125036e+02 6.968225e+01 6.168025e+01 6.487279e+01 6.009201e+01
#> [601] 6.063579e+01 8.372881e+01 6.866689e+01 7.228660e+01 5.480008e+01
#> [606] 7.679530e+01 7.246533e+01 5.960100e+01 5.652453e+01 6.386044e+01
#> [611] 6.533114e+01 6.463657e+01 6.572224e+02 6.220472e+01 6.379685e+01
#> [616] 6.895683e+01 6.510352e+01 7.752414e+01 6.999003e+01 9.897779e+01
#> [621] 1.137789e+02 5.930440e+01 5.450903e+01 6.409724e+01 6.398002e+01
#> [626] 6.551080e+01 6.016553e+01 6.188793e+01 5.832260e+01 7.240693e+01
#> [631] 1.201604e+04 1.610792e+02 1.620423e+02 5.296821e+01 8.180165e+01
#> [636] 6.167494e+01 9.889412e+01 5.912478e+01 5.849946e+01 5.957495e+01
#> [641] 9.346523e+01 7.208188e+01 5.511886e+01 5.911976e+01 6.334112e+01
#> [646] 5.544636e+01 5.300568e+01 1.233031e+02 1.540608e+02 8.289451e+08
#> [651] 8.365745e+01 5.613948e+01 5.822077e+01 6.680585e+01 8.294866e+01
#> [656] 6.296911e+01 6.112962e+01 5.940152e+01 6.041496e+01 5.756941e+01
#> [661] 6.692678e+01 5.697033e+01 8.101846e+01 5.799533e+01 6.397009e+01
#> [666] 7.872461e+01 7.407537e+01 1.475225e+03 3.007868e+02 7.237789e+01
#> [671] 7.137037e+01 7.341092e+01 1.314785e+02 6.534589e+01 4.768929e+01
#> [676] 5.768140e+01 5.876114e+01 6.337168e+01 1.734678e+02 2.484033e+02
#> [681] 1.132767e+02 6.213857e+01 6.431794e+01 6.947791e+01 6.706629e+01
#> [686] 1.338523e+02 5.831983e+01 5.957812e+01 1.991926e+02 7.765346e+02
#> [691] 5.944916e+01 8.150611e+01 7.887247e+01 1.132249e+02 6.537295e+01
#> [696] 2.395678e+02 8.786979e+01 3.109627e+02 7.215439e+02 8.164989e+01
#> [701] 9.033068e+01 9.348889e+01 6.354194e+01 5.707011e+01 5.698042e+01
#> [706] 8.760070e+01 6.101097e+01 6.140747e+01 1.605537e+07 5.977149e+01
#> [711] 6.291423e+01 6.165308e+01 6.804916e+01 6.103382e+01 6.098060e+01
#> [716] 5.614140e+01 5.894804e+01 6.525630e+01 6.141123e+01 6.238697e+01
#> [721] 6.935394e+01 7.610714e+01 6.113189e+01 7.116612e+01 5.909293e+01
#> [726] 5.797699e+01 2.055478e+02 5.985790e+01 6.010471e+01 6.185347e+01
#> [731] 5.837336e+01 1.142872e+02 5.909010e+01 5.663438e+01 6.771030e+01
#> [736] 7.478370e+01 6.624914e+01 5.940013e+01 7.564317e+01 5.934687e+01
#> [741] 6.136237e+01 6.222524e+01 8.895361e+01 5.548197e+01 6.622706e+01
#> [746] 4.833397e+01 7.279589e+01 5.637378e+01 1.026841e+02 8.470041e+01
#> [751] 6.339452e+01 1.152502e+02 7.981606e+01 6.491899e+01 7.693418e+01
#> [756] 6.158231e+01 2.312940e+02 6.308693e+01 5.999018e+01 6.675783e+01
#> [761] 5.872352e+01 1.339011e+02 6.960878e+01 6.288175e+01 6.743221e+01
#> [766] 6.748353e+01 1.428269e+02 1.312745e+02 1.117181e+02 1.352126e+02
#> [771] 6.427375e+01 5.693402e+01 6.450869e+01 6.592759e+01 8.156739e+01
#> [776] 6.387847e+01 6.599210e+01 6.010777e+01 5.877795e+01 5.902681e+01
#> [781] 9.290823e+01 5.336359e+01 5.813762e+02 6.979208e+01 5.995837e+01
#> [786] 6.647898e+01 1.001754e+02 6.060986e+01 6.288131e+01 6.346774e+01
#> [791] 5.534771e+01 5.663170e+01 6.375231e+01 8.743889e+01 5.869745e+01
#> [796] 7.569488e+01 5.934050e+01 5.743031e+01 5.905897e+01 5.909812e+01
#> [801] 5.954635e+01 5.587964e+01 1.792342e+02 5.690333e+01 7.067454e+01
#> [806] 7.337377e+01 6.695918e+01 6.802922e+01 5.932456e+01 7.959170e+01
#> [811] 5.536273e+01 6.412158e+01 5.739717e+01 6.379292e+01 7.424977e+01
#> [816] 8.033798e+01 6.113901e+01 7.066499e+01 6.136227e+01 5.613695e+01
#> [821] 6.917953e+01 8.104699e+01 6.679911e+01 6.732307e+01 9.352440e+01
#> [826] 6.813021e+01 6.988472e+01 6.217479e+01 1.074968e+02 1.059395e+02
#> [831] 8.853631e+01 1.340082e+02 6.506569e+01 7.303032e+01 6.669160e+01
#> [836] 6.783282e+01 5.286214e+01 6.728503e+01 6.036036e+01 5.746459e+01
#> [841] 7.608554e+01 8.090402e+01 6.316479e+01 5.726953e+01 7.697150e+01
#> [846] 7.776877e+01 5.979900e+01 5.999429e+01 6.812370e+01 8.336964e+01
#> [851] 1.194217e+02 7.239341e+01 5.838041e+01 5.803320e+01 8.286245e+01
#> [856] 5.926719e+01 5.702614e+01 1.353698e+02 5.842404e+01 1.734409e+02
#> [861] 6.242330e+01 1.175358e+02 6.949616e+01 2.647901e+02 7.884167e+01
#> [866] 5.500634e+01 5.967960e+01 6.907920e+01 5.845302e+01 5.942173e+01
#> [871] 5.897678e+01 7.010664e+01 1.513249e+02 1.338265e+02 8.214985e+01
#> [876] 5.630540e+01 6.479634e+01 7.675074e+01 6.763463e+01 6.665962e+01
#> [881] 5.179981e+01 8.526876e+01 5.769923e+01 6.643518e+01 7.236619e+01
#> [886] 5.628737e+01 6.743709e+01 6.940042e+01 1.140780e+02 5.812955e+01
#> [891] 6.002293e+01 1.281598e+02 5.015892e+01 1.031094e+02 6.602952e+01
#> [896] 8.661929e+01 6.179487e+01 9.189142e+01 5.742447e+01 5.881808e+01
#> [901] 8.892244e+01 8.441233e+01 6.822366e+01 6.014261e+02 8.080438e+01
#> [906] 7.581181e+01 6.519013e+01 6.955521e+01 7.815142e+01 8.635785e+01
#> [911] 5.655157e+01 6.839699e+01 6.235592e+01 5.615372e+01 6.643563e+01
#> [916] 7.444370e+01 6.304225e+01 6.027408e+01 9.942767e+01 6.074523e+01
#> [921] 6.291793e+01 6.055734e+01 6.765971e+01 8.177601e+01 6.016827e+01
#> [926] 5.760209e+01 7.389319e+01 6.797234e+01 6.473447e+01 6.152691e+01
#> [931] 6.388699e+01 1.491231e+02 5.267645e+01 5.935580e+01 1.643341e+02
#> [936] 2.038893e+02 8.555184e+01 1.686822e+02 6.543638e+01 8.434140e+01
#> [941] 1.092483e+02 5.944164e+01 1.616521e+02 6.850311e+01 5.940367e+01
#> [946] 5.875259e+01 6.607917e+01 2.076980e+02 1.005667e+02 5.943843e+01
#> [951] 5.910691e+01 5.724742e+01 5.976775e+01 5.294117e+01 5.915243e+01
#> [956] 8.049441e+01 5.237390e+01 5.794765e+01 9.388163e+01 7.665081e+01
#> [961] 5.981217e+01 6.210536e+01 7.069693e+01 6.383472e+01 6.564845e+01
#> [966] 6.729330e+01 6.336205e+01 1.024854e+02 6.044407e+01 6.496134e+01
#> [971] 8.191239e+01 6.338331e+01 7.241423e+01 9.218072e+01 9.475968e+01
#> [976] 6.517921e+01 6.763390e+01 7.333051e+01 5.877759e+01 2.956068e+02
#> [981] 7.333630e+01 5.320524e+01 6.217501e+01 8.005741e+01 7.273110e+01
#> [986] 6.119346e+01 5.771636e+01 6.084738e+01 6.880534e+01 6.275819e+01
#> [991] 7.523009e+01 5.958197e+01 6.556391e+01 6.241506e+01 6.142553e+01
#> [996] 6.471703e+01 7.153491e+01 1.196430e+02 6.013673e+01 5.890647e+01