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
(
GeneralModel
orModelPseudo
)
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.445180e+02 1.322162e+02 1.339096e+02 1.723122e+02 9.651866e+01
#> [6] 7.400532e+02 9.734743e+01 1.690455e+02 8.130419e+01 8.041935e+01
#> [11] 6.610001e+01 9.187052e+01 9.191528e+01 1.297635e+02 9.019012e+01
#> [16] 1.864280e+02 1.512489e+02 8.921574e+01 1.553881e+02 9.469229e+01
#> [21] 1.361454e+02 1.193003e+02 1.794959e+02 7.764095e+01 6.501217e+01
#> [26] 7.651283e+01 7.975129e+01 7.691988e+01 6.444792e+01 7.474535e+01
#> [31] 6.074717e+01 9.189681e+01 1.202842e+02 1.114340e+02 8.945025e+01
#> [36] 1.377218e+02 6.311969e+01 4.373525e+03 8.436218e+01 7.057643e+01
#> [41] 8.239678e+01 7.651307e+01 7.206465e+01 1.027939e+02 7.760176e+01
#> [46] 7.881910e+01 6.518850e+01 8.667574e+01 9.055783e+01 9.061699e+01
#> [51] 8.412353e+01 1.036849e+02 2.408859e+04 1.845982e+02 6.889109e+01
#> [56] 1.109446e+02 1.094271e+02 1.743556e+02 8.168450e+01 6.265998e+01
#> [61] 7.828203e+01 7.633054e+01 7.409357e+01 8.933097e+01 9.382607e+01
#> [66] 1.065747e+02 8.295922e+01 7.639031e+01 6.746773e+01 8.612754e+01
#> [71] 1.213064e+02 1.718720e+02 8.113671e+01 4.885124e+02 1.367041e+02
#> [76] 6.581957e+01 2.938535e+02 8.225225e+01 8.278012e+01 7.770174e+01
#> [81] 8.155909e+01 1.037727e+02 1.004179e+02 5.199129e+02 1.377487e+02
#> [86] 6.449028e+01 7.237359e+01 8.663363e+01 1.895074e+02 7.444990e+01
#> [91] 8.846258e+01 1.132332e+04 1.197173e+02 6.274647e+01 1.068803e+03
#> [96] 2.144085e+02 6.541421e+01 6.201634e+01 2.219412e+02 1.240176e+02
#> [101] 8.297499e+01 1.863713e+02 9.913548e+01 7.899058e+01 7.832842e+01
#> [106] 1.280972e+02 6.766797e+02 5.769860e+02 1.367287e+02 9.103786e+01
#> [111] 9.660107e+01 7.144744e+01 8.181445e+01 1.106086e+02 1.111699e+02
#> [116] 8.360774e+01 7.607234e+01 9.403165e+01 7.309523e+01 1.609369e+02
#> [121] 2.775650e+02 3.931628e+02 7.424184e+01 9.388809e+01 4.327575e+02
#> [126] 1.474900e+02 2.068206e+02 6.152035e+01 9.133130e+01 8.223233e+01
#> [131] 8.769164e+01 7.649757e+01 1.348076e+02 2.395477e+02 7.764552e+01
#> [136] 7.918407e+01 8.896788e+01 6.736526e+01 1.816681e+02 1.220760e+02
#> [141] 3.544916e+02 7.426578e+01 2.307264e+03 7.864027e+01 1.894437e+02
#> [146] 8.117553e+01 7.795980e+01 8.085816e+01 7.390503e+01 6.855546e+01
#> [151] 7.362052e+01 9.168349e+01 8.207773e+01 7.586167e+01 8.397670e+01
#> [156] 2.346227e+02 1.084754e+02 1.850142e+02 2.207476e+02 7.236167e+01
#> [161] 1.222373e+02 7.766725e+01 6.397812e+01 8.794011e+01 8.148554e+01
#> [166] 1.118388e+02 7.749884e+01 8.951964e+01 9.004225e+01 4.025684e+02
#> [171] 8.474765e+02 8.695172e+01 1.930002e+02 8.918524e+01 8.498154e+01
#> [176] 8.429673e+01 7.799966e+01 7.181959e+01 1.272417e+02 7.509511e+01
#> [181] 7.620533e+01 1.524398e+02 8.533800e+01 6.818560e+01 7.626969e+01
#> [186] 9.989723e+01 8.747778e+01 1.038203e+02 7.095278e+01 6.337852e+01
#> [191] 8.953438e+01 1.141471e+02 1.063945e+02 6.574085e+01 6.411226e+01
#> [196] 1.819414e+02 1.271099e+02 8.416499e+01 6.545649e+02 8.536475e+01
#> [201] 7.168938e+01 1.019516e+02 3.333981e+02 7.781354e+01 7.998697e+01
#> [206] 7.551017e+01 9.673733e+01 8.522521e+01 1.327299e+02 6.582310e+01
#> [211] 1.013469e+02 9.241205e+01 7.763522e+01 1.000721e+02 6.073401e+01
#> [216] 2.674642e+03 1.080461e+02 9.514770e+01 8.366710e+01 8.991243e+01
#> [221] 7.405701e+01 2.133880e+02 6.319647e+01 5.358738e+02 1.281074e+02
#> [226] 7.888456e+01 2.804828e+02 5.494701e+02 1.053584e+03 1.500269e+04
#> [231] 2.231671e+02 6.725259e+02 6.958821e+02 1.919647e+02 4.981651e+02
#> [236] 6.134287e+01 7.321128e+01 9.062958e+01 9.994984e+01 7.751697e+01
#> [241] 1.302061e+02 9.582316e+01 2.221320e+03 1.322969e+02 1.179794e+02
#> [246] 1.111769e+02 8.588027e+01 2.222594e+02 1.782539e+02 1.130825e+02
#> [251] 9.374893e+01 1.183625e+02 7.929055e+01 1.104899e+02 8.423604e+01
#> [256] 9.417530e+01 3.269263e+02 7.648655e+01 7.103860e+01 7.382061e+01
#> [261] 6.688389e+01 6.323264e+01 6.859591e+01 6.097774e+01 9.425516e+01
#> [266] 7.835773e+01 1.099707e+02 6.475271e+01 7.437602e+01 7.611694e+01
#> [271] 7.528788e+01 6.574246e+01 8.314190e+01 1.517744e+02 1.069607e+03
#> [276] 1.897717e+04 1.139813e+02 9.091508e+01 1.955644e+02 1.184012e+02
#> [281] 1.006542e+02 7.971854e+01 7.333653e+01 6.174992e+01 7.877589e+01
#> [286] 1.138732e+02 8.145192e+01 9.828532e+01 7.744543e+01 1.661016e+02
#> [291] 2.090349e+02 7.534226e+01 1.105010e+02 7.288376e+01 8.059513e+01
#> [296] 6.718881e+01 1.145717e+02 1.210318e+02 2.403392e+02 2.857365e+02
#> [301] 1.139020e+02 8.243739e+01 1.098740e+02 8.119777e+01 3.602400e+02
#> [306] 1.032133e+03 6.791726e+01 1.353131e+02 3.674924e+02 8.684867e+01
#> [311] 1.762141e+02 7.564600e+01 6.466501e+01 6.783229e+01 1.009082e+02
#> [316] 2.688530e+02 1.671224e+03 2.295926e+04 6.147969e+01 1.218279e+02
#> [321] 1.455399e+02 2.963080e+02 1.627023e+02 8.405865e+01 8.635836e+01
#> [326] 8.852390e+01 8.951051e+01 6.922514e+01 6.645404e+01 7.229236e+03
#> [331] 8.173037e+01 6.337354e+02 1.381243e+02 7.443833e+01 4.959756e+02
#> [336] 7.002008e+01 6.160190e+01 6.472334e+01 4.435723e+02 1.200342e+02
#> [341] 7.173985e+01 1.989115e+02 7.805490e+01 1.220730e+02 7.554781e+01
#> [346] 8.996809e+01 9.600648e+01 2.306039e+02 7.129557e+01 6.648480e+01
#> [351] 6.096163e+01 1.936244e+02 1.691569e+02 7.408455e+01 7.358427e+01
#> [356] 9.295133e+01 2.838085e+02 1.989292e+02 7.121369e+01 9.342579e+01
#> [361] 7.342902e+01 6.343180e+01 1.192647e+02 2.151348e+02 8.096312e+01
#> [366] 7.410189e+01 6.633888e+01 7.029197e+01 9.614390e+01 7.260177e+01
#> [371] 8.622289e+01 6.974797e+01 9.810058e+01 6.252921e+01 7.189592e+01
#> [376] 8.579073e+01 7.068698e+01 8.867667e+01 7.287081e+01 1.141546e+02
#> [381] 1.405751e+02 1.667878e+02 6.514783e+01 9.660924e+01 6.924543e+01
#> [386] 8.295018e+01 1.891342e+02 4.139004e+02 1.222244e+02 5.540952e+03
#> [391] 9.137079e+01 8.409199e+01 9.090832e+02 7.113120e+01 6.471931e+01
#> [396] 1.575024e+02 6.213849e+01 6.778271e+01 6.849377e+01 8.343947e+01
#> [401] 8.760392e+01 1.305523e+02 7.218500e+01 3.380503e+02 6.856884e+01
#> [406] 9.397479e+01 7.821512e+01 6.743324e+01 1.015532e+02 7.473384e+01
#> [411] 9.881561e+01 1.155722e+02 1.080701e+02 1.110174e+02 9.013301e+01
#> [416] 1.424997e+02 7.062049e+01 3.577570e+02 7.776286e+01 6.808146e+01
#> [421] 1.414678e+02 8.197652e+01 6.118849e+01 6.396353e+01 6.120050e+01
#> [426] 7.839613e+01 6.693003e+01 5.110252e+05 1.405826e+02 2.418423e+02
#> [431] 8.550268e+01 6.790963e+01 7.739154e+01 3.256600e+02 2.152761e+02
#> [436] 2.167500e+02 1.315260e+02 7.352665e+01 1.674532e+02 7.108183e+01
#> [441] 6.958678e+01 7.673394e+01 1.047379e+03 2.046528e+02 6.730911e+01
#> [446] 9.093018e+01 6.942513e+01 6.617782e+01 8.137791e+01 7.046747e+01
#> [451] 1.075565e+02 8.979266e+01 1.752648e+02 6.999356e+01 9.278142e+01
#> [456] 7.537419e+01 8.453987e+01 7.242167e+01 1.650871e+02 6.455079e+01
#> [461] 1.452776e+02 6.871726e+01 1.306991e+02 6.897587e+01 6.726045e+01
#> [466] 7.265431e+01 8.030916e+01 6.638953e+01 7.698561e+01 7.030061e+01
#> [471] 2.080913e+02 8.953680e+01 1.959462e+02 8.810912e+01 8.904719e+01
#> [476] 7.068419e+01 1.050573e+02 1.137272e+02 1.495695e+02 7.513493e+01
#> [481] 8.837192e+01 7.131750e+01 1.158094e+02 6.161128e+01 7.922060e+01
#> [486] 7.581021e+01 8.190160e+01 6.879798e+01 6.486017e+01 2.090750e+02
#> [491] 6.226960e+01 7.186024e+01 8.428019e+01 7.438688e+01 6.020422e+01
#> [496] 6.756395e+01 6.982990e+01 7.732411e+01 6.229580e+01 1.880501e+02
#> [501] 6.796904e+01 1.161699e+02 1.208918e+02 9.419265e+02 1.185031e+02
#> [506] 7.086954e+01 1.008569e+02 8.979660e+01 9.465602e+01 6.709939e+01
#> [511] 6.493021e+01 6.453958e+01 1.146231e+02 6.022856e+01 6.356410e+01
#> [516] 6.881069e+01 6.587036e+01 9.164368e+01 8.882913e+01 7.555138e+01
#> [521] 6.738491e+01 7.505573e+01 4.470554e+03 8.034245e+01 5.938715e+02
#> [526] 2.362877e+02 1.015162e+02 9.717427e+01 1.047782e+02 7.264140e+01
#> [531] 9.565446e+01 5.380080e+02 5.660703e+03 7.654595e+01 7.992901e+01
#> [536] 9.926119e+01 7.276836e+01 6.708635e+01 7.376456e+01 1.443000e+02
#> [541] 1.337784e+02 1.919070e+02 8.181696e+01 8.056279e+02 1.257724e+02
#> [546] 9.055170e+01 8.831768e+01 1.280129e+02 1.227136e+02 1.583916e+02
#> [551] 7.533237e+01 8.638621e+01 1.354703e+02 1.588981e+02 7.978153e+01
#> [556] 4.985115e+02 1.132225e+02 6.913535e+01 3.297624e+02 1.755687e+02
#> [561] 1.190583e+02 6.814409e+01 7.575536e+01 1.077858e+02 7.975640e+01
#> [566] 2.115534e+02 7.747865e+01 1.056179e+02 1.381232e+02 1.777683e+02
#> [571] 2.504126e+02 7.848117e+01 8.350991e+01 6.840401e+01 7.569935e+01
#> [576] 6.422764e+01 6.641595e+01 2.232626e+02 1.043386e+02 7.201145e+01
#> [581] 9.209173e+01 7.543569e+01 8.187405e+01 1.574686e+04 6.575877e+01
#> [586] 7.939951e+01 6.820039e+01 8.639475e+01 5.597544e+02 7.842689e+01
#> [591] 8.372084e+01 1.982896e+02 6.274801e+01 9.942507e+01 8.906499e+01
#> [596] 7.016559e+01 5.973971e+01 6.873554e+01 6.282311e+01 8.219100e+01
#> [601] 8.216195e+01 6.502518e+01 6.780539e+01 1.024224e+02 1.040644e+02
#> [606] 6.936723e+01 7.670966e+01 7.359815e+01 4.286484e+02 1.319777e+02
#> [611] 9.224334e+01 1.642704e+02 1.602811e+02 8.969170e+01 1.041115e+02
#> [616] 9.962604e+01 8.390440e+01 6.991526e+01 8.739419e+01 1.902587e+02
#> [621] 1.122033e+02 7.959589e+03 6.496562e+01 1.008525e+02 6.718791e+01
#> [626] 6.355190e+01 7.321977e+01 6.295591e+01 4.992561e+02 7.891591e+01
#> [631] 7.308628e+01 1.266285e+02 6.840231e+01 1.742463e+04 5.506767e+02
#> [636] 7.462424e+01 1.649372e+02 9.838200e+01 7.634589e+01 1.137771e+02
#> [641] 6.368727e+01 6.300091e+01 6.794239e+01 8.457307e+01 1.155144e+02
#> [646] 1.304843e+02 7.203424e+01 8.564167e+01 6.853984e+01 8.050863e+01
#> [651] 6.739812e+01 1.606498e+02 1.079515e+02 6.872714e+01 6.796674e+01
#> [656] 7.008907e+01 1.036131e+02 3.026201e+02 6.725227e+01 6.978607e+01
#> [661] 1.016243e+02 1.186339e+02 4.051233e+02 1.018569e+02 2.019300e+02
#> [666] 2.221115e+02 8.358633e+01 7.026835e+02 6.188212e+01 9.089988e+01
#> [671] 7.483591e+01 1.235172e+02 7.789139e+01 1.572906e+02 6.672748e+02
#> [676] 1.058006e+02 1.203621e+02 1.872436e+02 6.468983e+01 7.512011e+01
#> [681] 6.647291e+01 9.223351e+01 3.904449e+02 2.983836e+02 1.582551e+02
#> [686] 1.385867e+02 2.097798e+02 7.103623e+01 8.630478e+02 1.353197e+02
#> [691] 8.301085e+02 5.373223e+02 1.041516e+02 6.673326e+01 1.343685e+02
#> [696] 5.586872e+03 1.555997e+02 6.737058e+01 7.055456e+01 1.871247e+02
#> [701] 7.434800e+01 6.936980e+01 1.007101e+02 1.239291e+02 7.411310e+01
#> [706] 1.465045e+02 8.732640e+01 1.709527e+02 6.364129e+01 6.849882e+01
#> [711] 7.092288e+01 1.079399e+02 6.994559e+01 1.690061e+02 1.082065e+02
#> [716] 7.144079e+01 6.605763e+01 8.008986e+01 7.200933e+01 6.774310e+01
#> [721] 7.922452e+01 1.033994e+02 1.214764e+02 6.452459e+01 6.733882e+01
#> [726] 1.139184e+02 8.224328e+01 6.622711e+01 9.614114e+01 2.806863e+02
#> [731] 8.895260e+01 8.471828e+01 1.352582e+02 1.109150e+02 1.280254e+02
#> [736] 8.615450e+01 1.123023e+02 9.141361e+04 7.295377e+01 7.440399e+01
#> [741] 1.181801e+02 2.619327e+03 6.492213e+01 1.131437e+02 7.796791e+01
#> [746] 1.011402e+02 6.594227e+01 8.766476e+01 6.834433e+01 1.464531e+02
#> [751] 7.663853e+01 1.042957e+02 3.422264e+02 8.435080e+01 4.778787e+02
#> [756] 2.647439e+03 4.168433e+02 7.278268e+01 6.249580e+01 7.239274e+01
#> [761] 6.698512e+01 3.938782e+02 6.263894e+01 6.380713e+01 8.116752e+01
#> [766] 8.946713e+01 2.036803e+02 1.139016e+02 7.594399e+01 6.994401e+01
#> [771] 1.163848e+02 9.692350e+01 8.899104e+01 8.124947e+01 7.810457e+01
#> [776] 8.224505e+01 1.664651e+02 8.939775e+01 7.643458e+01 1.303933e+02
#> [781] 1.763180e+02 1.872332e+02 1.068254e+02 8.850163e+01 7.045455e+01
#> [786] 1.218115e+02 1.121976e+02 9.070110e+01 7.528669e+01 3.675962e+02
#> [791] 4.010441e+02 6.853554e+01 4.641323e+02 1.217362e+02 7.604334e+01
#> [796] 1.313732e+02 3.841233e+02 2.137318e+02 2.058403e+02 1.828106e+02
#> [801] 8.224121e+01 7.771789e+01 6.694912e+01 9.233833e+01 7.166842e+01
#> [806] 1.177872e+02 6.753134e+01 1.340324e+02 7.139789e+01 1.927966e+02
#> [811] 1.002546e+02 1.165795e+02 1.395710e+02 9.369397e+01 9.429940e+01
#> [816] 7.522409e+01 1.318666e+02 6.389838e+01 9.475487e+01 7.607062e+01
#> [821] 7.848718e+01 9.514664e+01 2.421159e+02 8.508166e+01 2.696681e+02
#> [826] 7.586968e+03 7.763695e+01 8.044823e+01 8.835257e+01 1.012958e+02
#> [831] 8.699420e+01 8.108165e+01 7.962797e+01 3.803531e+02 7.293559e+01
#> [836] 6.892427e+01 9.033595e+01 6.289781e+01 6.000284e+01 2.027951e+02
#> [841] 8.194778e+01 6.433167e+01 7.264121e+01 7.492673e+01 9.952306e+01
#> [846] 7.482153e+01 7.036155e+01 1.136883e+02 1.575012e+02 1.503684e+02
#> [851] 7.931387e+01 6.681333e+01 6.925821e+01 1.145160e+02 7.914495e+01
#> [856] 6.098727e+01 6.234427e+01 8.157287e+01 6.541819e+01 1.300380e+02
#> [861] 1.180920e+02 6.718740e+01 7.482390e+01 8.203852e+01 8.134664e+01
#> [866] 9.966378e+01 1.250388e+02 1.162430e+02 8.033296e+01 8.072443e+01
#> [871] 1.307418e+02 8.430948e+01 8.402531e+01 6.195878e+02 3.499104e+02
#> [876] 2.219712e+02 9.728031e+01 2.512503e+02 8.139398e+01 8.463314e+01
#> [881] 6.999033e+01 6.266723e+01 1.284515e+02 1.332188e+02 8.532073e+01
#> [886] 8.036978e+01 6.403519e+01 8.424994e+01 6.528122e+01 6.841069e+02
#> [891] 1.069715e+02 2.110523e+02 7.069631e+01 1.317424e+02 1.510415e+09
#> [896] 1.528229e+02 8.463265e+01 2.316636e+02 1.354618e+02 1.105332e+02
#> [901] 7.636239e+01 6.830247e+01 6.590757e+01 6.920198e+01 2.643052e+02
#> [906] 4.481689e+02 4.767598e+02 9.001419e+01 7.065491e+01 7.527391e+02
#> [911] 1.888639e+02 9.284198e+01 2.557901e+02 5.459913e+02 1.568705e+02
#> [916] 1.012532e+02 6.697952e+01 7.043165e+01 1.998968e+03 7.220093e+01
#> [921] 8.586554e+01 1.706202e+02 2.781022e+04 7.269842e+02 6.508055e+01
#> [926] 1.649614e+02 9.367832e+01 1.332676e+02 1.162470e+02 9.348431e+01
#> [931] 7.194679e+01 9.015832e+01 6.566238e+01 2.417118e+02 8.700938e+01
#> [936] 1.040789e+02 1.672312e+02 1.209230e+02 8.326404e+01 1.221055e+02
#> [941] 9.967239e+01 1.094427e+02 1.399351e+02 1.039416e+02 7.613433e+01
#> [946] 9.736438e+01 1.743620e+02 9.143383e+02 1.250758e+02 1.627898e+02
#> [951] 1.248082e+02 7.506660e+01 8.574092e+01 9.342270e+01 7.605888e+01
#> [956] 7.435600e+01 6.784481e+01 7.099382e+01 6.777541e+01 8.103795e+01
#> [961] 6.183805e+01 2.195096e+02 1.023290e+02 9.160571e+01 7.057991e+01
#> [966] 6.575897e+01 8.140294e+01 2.046752e+02 7.201108e+01 8.674704e+01
#> [971] 1.904977e+03 1.124767e+02 1.100037e+02 1.121464e+02 8.408601e+03
#> [976] 1.801361e+02 7.499410e+01 7.356427e+01 7.897452e+01 1.935241e+02
#> [981] 6.724304e+01 8.336447e+01 7.178331e+01 9.571655e+01 6.964949e+01
#> [986] 6.316340e+01 9.185011e+01 2.509429e+02 1.817731e+02 9.513923e+01
#> [991] 1.525598e+02 1.820891e+02 1.027475e+02 6.467627e+01 1.255053e+02
#> [996] 7.026786e+01 7.448383e+01 7.938968e+01 7.038913e+01 7.110957e+01
doseFunction(model, alpha2 = samples@data$alpha2, beta = samples@data$beta, grade = 2L)(x = 0.25)
#> [1] 195.09822 81.35926 78.39750 81.44608 68.93855
#> [6] 162.63669 69.04521 81.34231 63.95363 68.37721
#> [11] 53.81473 65.12864 54.31416 67.79721 60.32905
#> [16] 98.07203 64.22577 69.28725 72.41814 68.51948
#> [21] 73.05225 68.47864 102.97638 62.27034 58.24754
#> [26] 64.26832 61.53277 60.94975 54.60279 60.70782
#> [31] 53.60791 68.48563 77.25165 82.05295 66.48677
#> [36] 72.20224 55.82553 291.09081 61.24806 58.46750
#> [41] 68.30035 63.77295 59.44158 73.81259 60.23703
#> [46] 62.71348 59.10920 71.61619 70.53832 63.79891
#> [51] 64.64780 78.51206 1979.68523 120.66563 59.28485
#> [56] 77.58349 68.44731 95.50322 64.06216 56.56979
#> [61] 60.70523 64.25099 60.66848 67.78512 65.56572
#> [66] 71.49784 65.18275 65.43645 60.54549 71.27396
#> [71] 71.80616 118.69343 70.44959 230.15042 98.32975
#> [76] 56.45410 154.20532 69.43134 65.65335 60.79941
#> [81] 59.32414 66.40238 66.41336 164.82325 86.34462
#> [86] 56.12344 62.22918 62.66062 83.31135 55.54442
#> [91] 67.93910 743.64529 74.48804 51.75642 426.43246
#> [96] 89.07419 55.48867 53.82756 57.04417 65.39266
#> [101] 50.92038 108.71317 58.27425 64.69253 53.99822
#> [106] 69.73224 171.02275 188.29795 91.79656 70.61068
#> [111] 69.26589 60.83556 67.81910 69.79844 64.62605
#> [116] 62.75724 56.82709 72.18397 59.84202 69.16977
#> [121] 126.68817 110.30336 57.72803 61.64118 87.68286
#> [126] 64.10886 62.58381 50.93645 66.79452 56.06319
#> [131] 57.05465 58.73764 54.08904 118.28609 62.47993
#> [136] 64.87631 73.38893 57.23035 113.18318 50.61843
#> [141] 103.05003 57.04836 470.89914 66.80742 100.56790
#> [146] 63.79778 59.88316 64.45505 62.65453 61.26891
#> [151] 63.65785 65.28891 61.26339 59.77439 65.23968
#> [156] 132.01777 82.91693 109.90889 103.12045 58.71372
#> [161] 75.95353 55.42189 55.04196 65.90402 63.53184
#> [166] 77.08601 60.85743 72.52468 65.33598 121.83096
#> [171] 155.17250 64.64455 105.55817 69.20718 62.33315
#> [176] 58.45554 64.41438 63.04386 70.54090 60.17751
#> [181] 60.49310 92.64854 67.67116 58.81511 60.51088
#> [186] 76.73256 67.10166 72.29449 64.73425 58.12949
#> [191] 73.26656 88.82555 80.82882 57.56652 58.39024
#> [196] 109.00627 87.88494 64.29290 81.97713 53.56448
#> [201] 52.02911 60.22842 150.93239 61.13064 61.90907
#> [206] 58.13507 63.29029 56.77190 70.94586 57.49607
#> [211] 66.84172 66.41600 66.26775 73.84313 53.56607
#> [216] 209.31198 80.09723 69.38911 67.09530 69.42672
#> [221] 61.18311 76.09846 53.02713 211.36361 88.78953
#> [226] 60.71570 88.90149 106.13761 117.28045 677.97249
#> [231] 56.08779 78.85594 138.07315 52.11260 156.85519
#> [236] 55.49411 62.25953 65.00230 79.61982 58.30206
#> [241] 74.42401 68.42003 304.26868 72.02201 69.54267
#> [246] 61.37749 66.92676 91.03095 104.10297 74.49584
#> [251] 71.01330 84.13059 63.12147 73.96928 65.90901
#> [256] 63.63875 165.06914 65.03680 61.30117 64.63748
#> [261] 56.43101 58.10901 60.80893 55.97045 72.69177
#> [266] 63.52876 83.78822 58.39758 65.69093 67.71173
#> [271] 63.35023 58.97385 66.98088 83.85597 300.08177
#> [276] 2436.94152 70.54051 60.55341 88.26043 73.16669
#> [281] 63.06323 54.47935 64.37076 56.70649 65.53001
#> [286] 89.54751 63.50863 75.90352 65.85263 88.07625
#> [291] 86.49169 57.42385 73.93461 59.31826 61.98649
#> [296] 55.61814 72.15350 79.38025 100.00855 103.12515
#> [301] 82.18355 63.23430 59.17677 57.61208 93.87723
#> [306] 209.12169 59.30669 77.35849 137.51008 66.75049
#> [311] 78.98109 64.16955 57.95572 60.63492 80.12382
#> [316] 120.27846 112.13877 1810.14886 57.09717 86.88680
#> [321] 70.57662 146.77077 77.12898 60.94412 67.42624
#> [326] 72.86134 71.07549 59.73188 58.48254 277.62141
#> [331] 58.55315 130.35781 63.23255 59.27723 159.67256
#> [336] 60.66034 57.21045 59.37616 163.39161 62.69787
#> [341] 58.52594 66.03031 52.60121 55.95435 54.33978
#> [346] 67.61063 73.43528 128.47291 60.63627 61.32416
#> [351] 56.18514 122.04849 105.73667 61.53471 56.68938
#> [356] 58.18490 85.80646 86.06128 63.26885 73.36952
#> [361] 61.04625 56.52845 74.04190 95.21574 60.51788
#> [366] 60.06959 57.21510 62.32609 73.15763 60.45511
#> [371] 66.81284 52.26166 62.98902 55.66292 56.86870
#> [376] 65.54788 56.58226 70.78875 55.06033 71.83492
#> [381] 86.80669 92.81820 53.82155 73.12286 60.49066
#> [386] 59.78789 61.22890 69.77130 70.02499 176.83611
#> [391] 64.41269 59.17164 278.55525 60.93152 56.60753
#> [396] 91.21114 55.97618 59.26787 58.95062 68.77975
#> [401] 66.50543 85.86759 58.20726 104.97786 57.96097
#> [406] 70.98720 66.22681 62.01010 83.56663 67.09027
#> [411] 76.01346 74.61075 80.10886 73.97252 59.77706
#> [416] 76.03570 59.19263 167.16120 53.55509 60.22613
#> [421] 88.22557 58.61882 54.72616 56.71535 54.22518
#> [426] 64.53501 57.22555 14231.60253 76.47334 124.68698
#> [431] 69.61049 59.94233 65.00487 87.41925 94.37983
#> [436] 108.19308 92.18681 56.11078 77.75360 61.63102
#> [441] 60.16424 65.11666 206.43733 91.12752 56.41089
#> [446] 68.22656 56.09009 58.98072 65.00456 56.81297
#> [451] 73.97366 71.51640 112.80924 58.54675 66.25161
#> [456] 63.72367 65.98974 57.57970 94.01633 55.41612
#> [461] 88.84091 56.65234 71.18290 57.49921 57.49857
#> [466] 60.51227 67.07232 59.59843 65.64527 56.23253
#> [471] 109.99399 59.12906 96.44403 69.41530 62.00031
#> [476] 56.49512 79.70381 60.58573 80.60783 59.03535
#> [481] 68.38755 54.47063 59.78636 55.07267 61.61582
#> [486] 62.13170 63.92312 59.26483 55.89417 119.53903
#> [491] 56.43829 62.32164 68.54247 58.29817 52.83279
#> [496] 59.34988 62.55569 65.90173 57.96582 113.41704
#> [501] 57.89252 78.66721 71.25626 231.45669 72.66283
#> [506] 59.25231 70.98879 69.38904 73.84338 56.77521
#> [511] 58.50473 59.27136 81.01079 56.13995 60.49991
#> [516] 63.66717 61.17323 76.07587 75.25013 64.68930
#> [521] 56.07549 56.87883 960.32146 63.16259 141.49440
#> [526] 122.31231 69.62361 70.50747 73.71226 58.14421
#> [531] 76.03561 232.65479 1182.83662 63.67394 59.02901
#> [536] 76.43077 64.07941 55.61563 60.51046 76.87283
#> [541] 82.51621 101.20744 63.40050 101.24690 74.38864
#> [546] 60.51659 55.33128 66.81880 83.86279 83.90472
#> [551] 60.86377 71.88927 80.50998 78.69693 64.60185
#> [556] 192.91344 63.46010 59.89248 150.80744 94.75917
#> [561] 73.17832 60.24308 62.37683 78.15439 63.79412
#> [566] 98.26090 61.95382 67.21795 55.13129 110.85625
#> [571] 105.81815 63.14093 53.76929 55.65365 60.66876
#> [576] 54.87233 57.69044 125.70536 69.54020 61.39229
#> [581] 70.46594 63.24227 64.93641 1061.20443 58.78562
#> [586] 65.16817 57.21091 66.80628 141.19692 60.05701
#> [591] 63.82521 90.61092 52.99937 74.11336 59.32514
#> [596] 53.42600 55.66318 62.15855 58.23526 60.77076
#> [601] 66.21266 58.01485 58.92375 73.23244 56.35659
#> [606] 58.57567 66.70915 65.20805 153.84219 73.81294
#> [611] 61.84375 77.73366 94.30987 63.48550 67.40926
#> [616] 63.71531 64.06842 58.69811 72.99747 86.21251
#> [621] 74.09454 257.67369 58.28674 76.39237 59.18319
#> [626] 59.26876 65.03671 56.42952 235.41503 69.84420
#> [631] 63.41976 74.15544 57.84014 394.00038 117.52682
#> [636] 61.55604 102.88523 69.52933 63.30031 79.56868
#> [641] 58.07558 58.56307 61.28318 71.21759 80.01978
#> [646] 78.66360 52.17297 56.01253 59.27456 66.08341
#> [651] 61.27245 93.17939 64.03840 58.87088 59.78958
#> [656] 58.76324 71.09852 107.68692 60.02492 57.91883
#> [661] 67.23121 70.45092 57.93069 73.35500 104.08226
#> [666] 125.12988 65.88259 204.34619 43.67089 52.14125
#> [671] 67.13395 89.25826 54.66943 84.07140 277.35157
#> [676] 61.62626 79.34754 106.20716 56.54964 61.88715
#> [681] 57.45167 73.34313 64.31967 71.58433 80.90106
#> [686] 76.04757 91.55081 59.41739 218.67396 90.57811
#> [691] 181.79854 213.07065 79.91703 60.84688 67.56314
#> [696] 390.44395 105.07062 59.37199 57.89140 117.93182
#> [701] 63.90381 61.05714 69.70410 74.87817 61.18954
#> [706] 95.27812 56.13704 109.84679 54.81351 62.35534
#> [711] 61.32043 84.33876 61.08413 107.09763 61.81796
#> [716] 56.48257 56.82119 63.00220 58.74338 60.67276
#> [721] 68.25570 65.14718 72.40174 56.61752 62.09902
#> [726] 92.01829 66.21701 57.47584 71.10143 113.12639
#> [731] 64.38605 62.61999 90.84195 82.20618 84.21728
#> [736] 68.34291 68.29615 740.48720 59.29537 62.75208
#> [741] 68.18327 370.85533 57.71184 82.98287 66.02253
#> [746] 75.45538 57.39342 65.98219 58.78701 87.07417
#> [751] 58.74029 69.34169 156.19026 51.81600 51.69988
#> [756] 189.70927 105.14485 53.18082 56.59670 57.85074
#> [761] 57.59211 100.63091 54.24104 53.45701 64.07642
#> [766] 66.80832 93.57795 71.36456 60.16585 56.29928
#> [771] 86.56424 73.21071 66.84851 64.44249 62.80108
#> [776] 63.55346 87.16925 59.10137 59.89482 81.12978
#> [781] 82.20605 92.26910 80.96036 61.01097 55.00328
#> [786] 87.20115 62.82244 61.41057 60.83613 111.34724
#> [791] 135.49939 58.67326 186.60229 75.11957 65.63173
#> [796] 79.54436 149.44927 86.97759 104.71615 86.00923
#> [801] 64.03534 59.51455 57.25501 66.37361 55.00336
#> [806] 67.38526 58.96211 75.47121 59.80246 94.56001
#> [811] 72.22432 76.28489 90.34312 66.51562 64.82612
#> [816] 63.59723 81.76529 57.61879 76.61912 59.32883
#> [821] 63.76431 68.38508 92.17941 64.26513 147.82466
#> [826] 786.07629 55.78110 63.58670 70.56771 77.47537
#> [831] 65.57320 63.04253 63.65850 97.99106 58.81883
#> [836] 61.22887 67.55614 55.52891 53.08191 129.62325
#> [841] 65.00599 56.39161 64.83813 64.29180 60.47464
#> [846] 59.32928 59.70234 77.72723 90.42910 93.16591
#> [851] 65.41731 59.13471 61.75204 69.75720 60.51741
#> [856] 56.15272 57.21375 65.13669 59.68296 88.08576
#> [861] 79.16182 60.56731 62.05625 64.11266 62.39983
#> [866] 69.80490 74.90367 76.55930 66.19747 60.66483
#> [871] 92.05372 70.06091 69.56326 161.11876 108.38367
#> [876] 101.98973 55.66152 127.96703 68.73925 71.68209
#> [881] 64.50457 58.32711 88.44792 77.97254 68.37280
#> [886] 57.16098 57.76021 61.93737 56.99147 292.54675
#> [891] 72.65400 88.72725 49.43718 71.43325 157284.71098
#> [896] 79.65870 65.38231 100.93748 80.07907 69.18098
#> [901] 66.01727 61.46116 60.07485 59.73077 118.13878
#> [906] 86.74157 157.09370 59.15542 57.72315 209.68340
#> [911] 97.24818 61.97365 109.17770 252.21423 69.88084
#> [916] 79.74560 55.86622 57.88140 365.57732 56.49837
#> [921] 67.77161 87.74472 2366.76139 114.26590 55.69276
#> [926] 88.47877 73.20794 86.54557 92.11920 78.88806
#> [931] 61.21771 65.76634 57.98666 130.79452 71.80927
#> [936] 76.37574 87.23902 70.70204 60.94044 75.62161
#> [941] 62.05995 78.76080 92.23052 75.47903 66.33953
#> [946] 74.49468 64.88001 213.97757 70.68856 95.61366
#> [951] 76.75441 63.85556 68.05571 65.72890 63.27557
#> [956] 64.83713 58.02808 62.76788 58.29143 68.30739
#> [961] 54.65406 116.50344 75.69229 69.44459 52.57203
#> [966] 58.98019 66.38464 114.17888 60.58450 70.21126
#> [971] 222.80347 73.02650 80.07901 84.29708 1803.11676
#> [976] 102.18884 53.78736 56.07949 55.43875 119.20371
#> [981] 59.24160 67.54102 57.47851 66.31847 61.34338
#> [986] 56.79188 72.21293 132.20734 91.39971 76.77323
#> [991] 82.99018 116.12655 67.22471 55.37711 84.07354
#> [996] 57.77875 59.80113 61.68643 60.81856 59.74095