![[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): Return a dose function for aGeneralModel.doseFunction(ModelPseudo): Return a dose function for aModelPseudo.doseFunction(LogisticLogNormalOrdinal): Return a grade-specific dose function for aLogisticLogNormalOrdinalmodel.
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.319571e+01 8.855370e+02 1.741886e+02 1.168427e+02 1.108437e+02
#> [6] 9.541036e+01 8.954650e+01 7.546828e+01 8.486949e+01 7.482167e+01
#> [11] 7.087064e+01 7.382850e+01 8.569430e+02 8.929326e+01 8.754033e+01
#> [16] 8.571346e+01 2.126977e+02 7.415752e+01 1.735163e+02 1.160713e+02
#> [21] 1.180213e+02 8.842502e+01 7.122360e+01 9.458791e+01 6.890318e+01
#> [26] 1.098908e+02 8.312398e+01 9.592483e+02 3.753505e+02 3.341925e+03
#> [31] 8.653001e+01 6.619673e+01 9.226934e+01 6.962170e+01 9.247202e+01
#> [36] 8.099447e+01 7.826299e+02 8.638343e+01 2.914006e+02 7.362387e+01
#> [41] 2.021030e+02 2.003336e+02 1.005326e+03 1.110668e+02 9.435108e+01
#> [46] 9.587448e+01 3.062984e+02 6.831658e+01 1.048565e+02 1.048011e+02
#> [51] 7.021923e+01 8.683912e+01 1.230537e+02 4.462352e+02 7.456285e+01
#> [56] 7.697681e+01 8.002713e+01 2.091165e+02 1.284178e+02 6.412143e+01
#> [61] 8.755517e+01 7.289982e+01 8.385831e+01 9.553292e+01 1.045262e+03
#> [66] 1.122562e+02 6.972717e+01 7.051008e+01 9.515429e+01 9.561535e+01
#> [71] 6.632509e+01 7.320648e+01 1.346111e+02 2.761857e+02 2.706160e+02
#> [76] 6.634422e+01 8.496211e+01 7.296854e+01 7.112947e+01 1.594395e+02
#> [81] 1.456930e+02 6.870178e+01 6.773515e+01 6.758395e+01 7.145361e+01
#> [86] 6.899058e+01 9.608007e+01 6.560344e+01 1.366253e+03 6.401571e+01
#> [91] 7.157139e+01 1.094302e+02 7.730845e+01 2.211784e+02 8.703954e+01
#> [96] 5.455957e+02 1.268206e+03 6.840143e+01 1.815700e+02 1.519196e+02
#> [101] 1.269368e+02 3.617568e+02 7.564492e+01 8.949823e+01 8.425495e+01
#> [106] 6.558573e+01 2.585686e+07 8.705206e+01 9.441405e+01 1.119732e+02
#> [111] 6.533859e+01 6.582401e+01 5.419652e+02 6.395982e+01 7.698404e+01
#> [116] 1.200418e+02 7.403808e+01 9.048494e+01 8.968448e+01 7.552362e+01
#> [121] 7.876783e+01 7.562726e+01 6.502762e+01 7.430755e+01 7.885873e+01
#> [126] 1.181234e+02 1.037360e+02 8.325985e+01 1.987359e+02 9.563603e+01
#> [131] 7.134147e+01 7.150768e+01 6.169901e+02 6.651598e+01 6.655636e+01
#> [136] 4.055481e+02 7.123229e+01 7.027358e+01 6.689442e+01 9.535171e+01
#> [141] 8.697228e+01 6.732989e+01 7.315967e+01 1.008726e+02 4.942240e+03
#> [146] 1.770996e+04 7.890181e+01 1.502403e+03 6.436437e+02 6.945812e+02
#> [151] 9.908159e+01 8.221052e+01 1.354994e+02 1.306748e+02 1.238105e+02
#> [156] 8.840100e+01 7.692386e+01 7.653745e+01 6.489840e+01 1.878071e+02
#> [161] 6.464251e+01 1.573562e+02 1.342645e+02 9.240737e+01 7.583449e+01
#> [166] 6.564053e+01 1.973567e+02 4.088095e+02 7.890605e+01 8.587745e+01
#> [171] 8.446599e+01 7.799347e+01 6.966612e+01 1.702691e+02 6.520680e+01
#> [176] 1.432661e+02 2.006055e+02 7.222146e+01 2.678542e+02 7.314708e+01
#> [181] 6.818573e+01 7.343518e+01 7.666367e+01 9.042533e+01 1.052232e+02
#> [186] 9.579773e+01 7.541000e+01 1.218397e+02 6.798570e+01 9.191726e+01
#> [191] 7.022424e+01 9.865116e+01 1.080050e+02 7.240069e+01 5.986701e+02
#> [196] 1.709219e+02 7.198774e+01 6.762066e+01 1.388900e+02 2.317152e+02
#> [201] 9.542955e+01 9.963489e+02 7.420018e+01 7.849400e+01 6.843132e+01
#> [206] 1.219843e+02 1.873545e+02 8.660494e+01 1.567438e+02 8.155546e+01
#> [211] 8.481571e+01 9.746208e+01 8.271732e+01 1.303591e+02 9.033480e+01
#> [216] 1.465336e+02 3.949676e+03 9.074524e+01 2.864505e+03 7.242309e+01
#> [221] 6.435558e+01 7.089777e+01 1.445731e+03 1.037601e+02 6.098158e+02
#> [226] 8.693519e+01 1.207043e+02 7.593362e+01 2.755968e+02 6.410136e+01
#> [231] 9.241114e+01 6.988896e+01 8.076544e+01 6.802488e+01 7.410666e+01
#> [236] 7.542008e+01 7.581843e+01 6.979953e+01 1.794540e+02 8.069347e+01
#> [241] 1.719955e+02 3.412792e+02 7.169916e+01 1.540025e+02 7.741213e+01
#> [246] 2.084277e+02 1.237675e+02 7.703199e+01 4.084457e+02 2.575884e+02
#> [251] 1.212443e+02 1.102956e+02 1.688483e+02 1.253195e+02 2.925536e+02
#> [256] 8.029188e+01 8.465308e+01 6.606686e+01 6.452276e+01 8.341733e+01
#> [261] 7.164562e+01 1.334783e+02 1.013032e+02 7.249267e+03 6.637489e+01
#> [266] 7.289977e+01 7.397548e+01 1.305011e+02 8.364581e+01 1.867636e+02
#> [271] 9.371289e+01 1.012284e+02 4.331683e+03 9.933838e+01 7.565242e+01
#> [276] 2.353562e+02 6.162213e+01 6.594222e+01 1.948147e+02 8.761594e+01
#> [281] 1.343832e+02 6.151908e+02 1.143143e+02 6.134180e+03 3.865434e+04
#> [286] 1.069775e+02 1.440553e+02 9.423901e+01 6.971067e+01 7.907040e+01
#> [291] 7.438864e+01 7.773245e+01 8.153888e+01 2.495707e+04 7.654917e+01
#> [296] 1.361056e+02 6.746622e+01 1.061591e+02 8.926474e+01 7.035202e+01
#> [301] 7.566714e+01 6.720585e+01 6.335543e+01 6.779633e+01 1.012351e+02
#> [306] 7.218014e+01 8.782828e+01 7.175329e+01 8.862374e+01 9.745177e+01
#> [311] 8.989596e+01 7.577034e+01 8.145226e+01 8.641736e+01 6.629953e+01
#> [316] 1.400370e+02 6.049427e+01 9.897761e+01 9.969099e+01 9.266319e+01
#> [321] 7.404603e+01 6.450172e+01 9.249491e+01 7.830600e+01 6.741671e+01
#> [326] 8.950526e+02 2.068898e+02 7.232368e+01 7.697897e+01 7.352569e+01
#> [331] 3.713729e+02 1.365598e+02 3.258074e+03 1.115699e+02 8.337877e+01
#> [336] 7.972912e+01 6.701717e+01 6.555039e+01 2.930961e+02 1.015567e+02
#> [341] 1.794466e+02 7.679187e+01 7.370063e+01 4.576072e+02 9.658414e+01
#> [346] 7.135897e+01 7.234594e+01 7.563145e+01 3.285703e+02 6.490790e+01
#> [351] 7.768492e+01 9.409324e+01 3.811795e+02 1.134003e+02 1.041053e+02
#> [356] 7.148666e+01 7.050555e+01 7.538370e+01 6.892549e+01 4.019885e+02
#> [361] 9.089109e+01 1.291126e+02 4.238695e+02 1.049294e+02 2.728321e+02
#> [366] 9.176378e+01 9.253813e+01 9.804413e+01 9.996457e+01 6.624504e+02
#> [371] 6.708343e+01 1.649978e+02 1.785134e+02 1.127883e+02 7.098581e+01
#> [376] 1.664971e+02 6.940716e+01 2.997743e+02 1.134012e+02 1.095922e+02
#> [381] 8.994731e+01 1.549061e+02 7.849659e+02 1.005882e+02 8.888413e+01
#> [386] 9.576887e+01 3.106780e+02 9.600136e+02 6.186209e+01 9.555845e+01
#> [391] 7.390985e+01 7.774422e+01 8.332882e+01 3.136676e+03 1.864725e+02
#> [396] 9.084516e+01 8.778739e+01 1.676577e+02 9.439582e+01 6.268271e+02
#> [401] 9.549024e+01 6.528033e+01 6.397370e+01 7.285200e+01 7.803500e+01
#> [406] 8.625452e+01 6.829279e+01 6.612209e+01 7.050346e+01 1.782338e+02
#> [411] 1.306762e+02 8.705559e+01 5.209199e+02 5.785436e+03 1.200733e+02
#> [416] 7.872394e+01 7.392795e+01 1.549896e+02 1.058966e+02 6.395745e+01
#> [421] 1.809407e+02 7.346257e+01 1.163799e+02 8.169846e+02 9.609838e+01
#> [426] 8.867593e+01 6.757720e+01 8.819280e+01 8.159907e+01 1.043875e+02
#> [431] 1.735695e+02 1.220952e+02 6.554234e+02 1.810089e+02 8.826962e+01
#> [436] 1.519214e+02 7.091787e+01 1.584840e+02 6.910451e+01 1.549415e+02
#> [441] 7.415348e+01 1.231565e+02 6.458029e+01 6.459801e+01 6.588156e+01
#> [446] 1.285595e+02 1.266125e+02 1.161800e+02 8.789277e+01 6.561260e+01
#> [451] 7.664607e+01 8.568418e+01 1.055992e+02 4.649511e+02 1.720891e+02
#> [456] 5.172269e+02 2.100857e+02 1.281409e+02 2.210487e+03 1.266219e+02
#> [461] 7.692563e+01 6.619133e+01 9.722768e+01 1.183468e+02 6.481315e+01
#> [466] 7.028225e+01 1.724609e+02 3.357781e+03 1.218086e+02 6.717563e+01
#> [471] 4.445543e+03 3.049080e+02 1.107106e+02 1.064546e+02 7.489146e+01
#> [476] 6.418265e+01 6.689045e+01 7.395075e+01 6.359403e+01 8.012451e+01
#> [481] 8.633258e+01 3.623032e+02 1.258462e+02 1.405378e+02 7.660990e+01
#> [486] 6.434557e+01 7.359715e+01 7.047525e+01 7.014168e+01 7.709661e+01
#> [491] 9.012150e+01 6.708801e+01 7.499557e+01 7.572933e+01 2.040461e+03
#> [496] 6.884487e+01 8.230431e+01 1.182956e+02 6.360847e+01 8.462369e+01
#> [501] 7.490395e+01 1.022979e+02 8.957785e+01 1.703462e+02 9.653557e+01
#> [506] 1.226863e+02 8.945254e+01 1.649882e+02 7.471991e+01 1.033769e+02
#> [511] 1.024436e+02 7.351359e+01 1.547731e+02 7.444417e+01 9.144098e+01
#> [516] 7.076011e+01 1.325790e+02 7.021206e+01 7.526249e+01 7.401521e+01
#> [521] 5.195285e+02 7.580207e+01 7.264610e+01 6.381057e+01 8.649597e+01
#> [526] 7.027229e+01 8.247926e+01 9.181909e+01 8.835321e+01 1.244685e+02
#> [531] 1.659137e+02 1.552653e+02 3.598972e+02 7.257582e+01 7.853503e+01
#> [536] 9.286765e+01 6.671518e+01 7.106682e+01 7.056413e+01 8.293928e+01
#> [541] 7.058544e+01 1.225640e+02 1.000166e+02 7.628610e+01 1.003634e+02
#> [546] 6.520688e+01 6.982886e+01 3.044224e+02 1.310887e+02 6.625193e+01
#> [551] 9.328544e+02 6.803223e+01 5.440401e+02 1.319224e+04 1.753417e+02
#> [556] 6.114800e+01 7.783839e+01 8.479653e+01 7.205343e+01 1.124204e+02
#> [561] 6.124007e+01 6.392676e+01 7.747685e+01 7.941164e+01 1.025722e+02
#> [566] 5.742586e+01 8.535814e+01 9.549882e+01 6.782310e+01 7.686023e+01
#> [571] 1.203893e+02 7.096614e+01 9.038936e+01 9.089813e+01 8.394137e+01
#> [576] 7.401892e+01 7.382545e+01 8.309530e+01 9.099171e+01 7.457540e+01
#> [581] 1.329500e+02 1.646834e+02 2.369656e+02 1.171586e+02 7.016939e+01
#> [586] 6.423807e+01 4.795516e+02 2.514767e+02 8.540899e+01 1.008189e+02
#> [591] 5.723048e+01 1.118075e+02 7.572737e+01 2.108185e+02 8.296942e+02
#> [596] 2.775699e+02 6.894457e+01 4.179630e+03 7.572863e+01 8.884109e+01
#> [601] 1.129231e+02 7.096880e+01 1.161298e+02 8.132777e+01 7.651394e+01
#> [606] 5.914087e+02 7.573953e+01 6.955453e+01 1.155515e+03 1.115016e+02
#> [611] 1.278178e+02 7.592521e+01 7.401927e+01 1.039717e+02 8.486277e+01
#> [616] 8.480794e+01 9.120420e+01 7.228983e+01 3.405497e+02 1.623725e+02
#> [621] 7.804172e+01 6.887781e+01 1.209007e+02 2.095904e+02 6.763749e+01
#> [626] 7.058693e+01 3.443124e+02 1.648855e+02 9.617540e+01 7.238988e+01
#> [631] 6.488693e+01 9.479002e+01 8.815126e+01 1.908394e+02 6.443532e+01
#> [636] 1.683269e+02 7.397096e+01 1.687265e+03 3.804353e+02 9.195814e+01
#> [641] 1.183420e+02 8.334719e+03 3.021014e+02 9.324077e+01 8.780866e+01
#> [646] 7.624110e+01 1.168578e+02 7.491664e+01 1.017843e+02 6.566862e+01
#> [651] 7.204375e+01 6.838185e+01 7.715616e+01 6.636928e+01 4.568826e+02
#> [656] 8.097497e+01 8.270109e+01 6.920061e+01 6.218093e+02 8.671640e+01
#> [661] 1.541343e+02 7.699119e+01 2.287374e+02 6.469159e+01 1.031063e+02
#> [666] 1.336032e+02 6.706884e+01 7.156461e+01 7.076222e+01 1.001191e+02
#> [671] 2.453137e+02 6.692400e+01 6.625998e+01 8.754064e+01 8.024504e+01
#> [676] 1.108317e+02 1.083453e+02 1.063231e+02 6.420154e+01 6.990458e+01
#> [681] 7.574307e+01 1.039062e+02 7.003157e+01 8.081050e+01 8.379581e+01
#> [686] 8.758705e+01 8.479546e+01 1.585491e+02 1.496248e+02 7.658957e+01
#> [691] 6.281847e+01 2.389132e+02 6.945899e+01 7.887881e+01 6.973354e+01
#> [696] 8.144844e+01 6.749245e+01 8.289981e+01 7.227735e+01 1.033140e+02
#> [701] 2.212664e+02 8.911413e+02 6.551876e+01 8.819922e+01 8.309615e+01
#> [706] 2.216908e+02 1.562325e+02 6.628116e+02 8.574582e+03 7.798253e+05
#> [711] 1.194162e+02 7.505193e+01 7.078455e+01 1.140824e+02 8.145940e+01
#> [716] 1.785029e+02 7.740008e+01 7.241507e+01 8.195524e+01 9.156507e+01
#> [721] 1.693770e+02 8.243545e+01 6.689303e+01 7.540210e+01 6.569126e+01
#> [726] 1.071252e+02 8.366097e+01 6.910465e+01 7.292091e+01 8.084237e+01
#> [731] 1.013026e+02 1.460563e+02 6.808362e+01 3.412965e+03 4.364702e+02
#> [736] 8.429662e+01 1.116897e+02 6.457475e+01 1.047843e+02 4.575927e+02
#> [741] 2.285889e+05 7.921513e+01 1.121869e+02 4.864028e+02 1.051440e+02
#> [746] 1.166355e+02 9.322530e+01 1.154826e+02 6.732598e+01 6.331507e+01
#> [751] 9.441467e+01 2.719286e+02 8.254864e+01 1.534472e+02 6.611118e+01
#> [756] 6.762674e+01 1.073085e+02 1.046462e+02 7.455111e+01 6.469783e+01
#> [761] 2.616021e+03 1.242720e+02 1.192228e+02 7.880773e+01 1.877099e+03
#> [766] 7.376887e+01 1.679863e+02 1.660010e+02 1.003615e+02 6.332417e+01
#> [771] 6.962726e+01 9.771135e+01 8.051067e+01 6.896655e+01 6.403535e+01
#> [776] 1.036684e+02 9.507623e+01 1.268789e+02 1.093551e+02 8.386242e+01
#> [781] 7.377085e+01 6.022625e+01 6.298784e+01 7.103519e+01 2.077944e+02
#> [786] 6.566476e+01 3.540200e+02 3.893038e+02 7.933588e+01 1.898745e+02
#> [791] 1.052683e+02 6.345030e+01 6.194290e+01 6.562413e+01 1.376026e+02
#> [796] 8.058583e+01 1.281861e+02 8.873976e+01 6.421130e+01 8.205088e+01
#> [801] 9.486082e+01 1.404194e+02 7.792973e+01 8.727637e+01 6.587025e+01
#> [806] 7.588321e+01 8.310069e+01 7.243000e+01 3.350426e+02 1.037754e+02
#> [811] 8.505449e+01 7.445106e+01 9.331286e+01 1.294595e+02 9.558444e+01
#> [816] 1.243913e+02 1.155909e+03 6.640124e+01 1.530152e+02 1.608581e+02
#> [821] 6.854209e+01 8.072214e+01 7.407034e+01 7.136733e+01 6.999944e+01
#> [826] 1.198785e+02 8.400038e+01 1.114386e+02 1.081549e+02 2.394609e+02
#> [831] 6.532070e+01 6.178561e+01 4.719792e+02 1.134058e+02 7.427188e+01
#> [836] 8.809740e+01 7.366467e+01 9.561342e+01 1.230884e+02 1.423515e+02
#> [841] 7.008259e+01 9.563096e+01 7.555904e+01 6.897877e+01 6.428295e+01
#> [846] 1.404956e+02 8.490612e+01 9.156847e+01 8.905832e+01 9.610285e+01
#> [851] 3.881882e+02 1.760968e+02 2.690706e+02 6.548348e+01 7.624725e+01
#> [856] 7.631403e+01 1.139014e+02 6.456285e+01 1.086320e+02 7.191786e+01
#> [861] 8.355759e+01 7.106905e+01 6.864441e+01 7.201282e+01 1.132934e+02
#> [866] 1.246431e+02 1.003870e+02 1.053613e+02 1.074786e+03 8.066035e+01
#> [871] 9.810387e+01 8.989812e+02 2.120190e+03 2.925826e+02 7.563680e+01
#> [876] 1.105376e+02 7.555809e+01 6.632454e+01 7.881377e+01 1.667776e+02
#> [881] 3.194217e+02 7.063290e+01 6.905786e+01 2.581947e+02 1.196661e+02
#> [886] 6.415987e+01 7.666189e+01 1.004274e+02 9.362803e+01 7.384928e+01
#> [891] 8.026386e+01 7.059019e+01 6.959827e+01 8.521227e+01 6.400017e+01
#> [896] 7.990201e+01 7.861177e+01 1.097870e+02 9.450723e+01 9.796856e+01
#> [901] 8.693884e+01 9.068481e+01 8.641997e+01 1.015563e+02 6.735746e+01
#> [906] 1.289830e+02 1.426205e+03 1.141572e+03 7.265354e+01 1.648681e+02
#> [911] 6.226428e+01 1.313727e+02 7.288479e+01 7.193695e+01 6.947555e+01
#> [916] 7.460988e+01 6.665537e+01 8.015137e+01 1.202613e+02 7.460876e+01
#> [921] 7.537070e+01 7.132405e+01 7.576627e+01 1.758657e+02 1.035694e+02
#> [926] 1.579116e+02 6.787382e+01 7.725139e+01 7.347958e+01 7.866402e+01
#> [931] 7.353411e+01 1.701068e+02 6.665365e+01 6.738149e+01 7.776741e+01
#> [936] 1.334200e+02 7.653308e+01 1.119394e+02 3.510009e+02 8.551199e+01
#> [941] 6.430647e+01 6.410124e+01 1.183767e+02 6.398676e+01 6.433971e+01
#> [946] 8.306748e+01 6.975016e+01 1.319418e+03 2.661404e+03 1.728835e+02
#> [951] 7.692587e+01 9.837784e+02 9.344909e+01 9.753933e+01 1.712361e+02
#> [956] 8.038963e+01 8.635581e+01 1.006932e+02 8.040605e+01 7.126012e+01
#> [961] 8.486921e+01 1.040678e+02 7.043633e+01 1.262335e+02 6.592122e+01
#> [966] 1.503402e+02 7.289700e+01 7.715561e+01 8.266465e+01 6.681610e+01
#> [971] 7.672397e+01 1.039564e+03 1.031294e+05 6.392822e+01 9.643839e+01
#> [976] 8.999677e+01 6.831271e+01 8.553510e+01 6.691016e+01 6.806960e+01
#> [981] 2.288954e+02 1.090339e+02 5.963711e+01 1.594993e+02 2.424593e+02
#> [986] 8.199821e+03 6.138687e+02 8.984883e+01 9.140660e+01 1.130124e+02
#> [991] 1.942040e+02 1.384294e+02 3.655950e+03 7.309670e+01 9.089079e+01
#> [996] 7.336830e+01 7.536545e+01 8.298200e+01 8.681917e+01 7.070481e+01
doseFunction(
model,
alpha2 = samples@data$alpha2,
beta = samples@data$beta,
grade = 2L
)(x = 0.25)
#> [1] 63.54391 368.50452 75.73909 76.64513 68.84074 66.72431
#> [7] 70.61497 61.40550 62.48280 63.46367 60.59056 61.28721
#> [13] 158.35116 66.59646 66.30933 69.48178 101.76935 63.05207
#> [19] 86.93709 67.08541 68.71718 63.90152 62.03920 66.94148
#> [25] 59.57242 54.66769 47.00173 144.65068 144.99736 406.68738
#> [31] 65.49863 56.96792 61.82909 59.90974 70.26764 60.02263
#> [37] 221.93207 62.96909 132.05499 60.45612 89.20512 111.20601
#> [43] 168.25169 75.46670 63.09410 71.29598 173.91352 60.36464
#> [49] 59.85704 59.90618 60.01478 59.68193 61.77098 221.80429
#> [55] 61.50983 48.87806 62.10102 121.52627 81.09018 57.70046
#> [61] 71.22704 58.40453 60.28683 67.04028 105.32381 69.84193
#> [67] 59.20925 59.98350 75.08081 62.46476 55.64193 61.11418
#> [73] 70.92512 37.51616 91.87323 55.34344 70.50427 61.84276
#> [79] 60.02023 86.46333 80.49563 60.03287 58.89720 60.06337
#> [85] 59.34235 57.67003 74.51686 59.11058 444.80256 58.22964
#> [91] 63.35431 76.97509 66.95926 115.64015 59.38286 125.99905
#> [97] 106.05777 55.78946 89.53169 63.39087 80.06549 138.92215
#> [103] 58.76437 68.17835 62.25610 58.20093 36483.06860 65.31482
#> [109] 61.55053 61.72814 54.41417 56.96307 214.48312 57.40263
#> [115] 64.62670 81.49583 61.68015 71.83259 70.22021 60.91694
#> [121] 66.76654 64.41332 56.58676 56.08686 61.93784 74.27949
#> [127] 52.70444 55.24351 82.23215 65.36751 54.28089 53.99296
#> [133] 201.13721 58.55151 57.75040 199.46193 59.22578 59.12350
#> [139] 60.82376 68.18064 71.62970 58.89266 60.71343 72.40696
#> [145] 302.67719 2128.71351 61.17100 278.96048 214.98427 99.98525
#> [151] 75.67501 59.83783 86.45928 89.24685 85.73360 60.81430
#> [157] 62.70801 64.39001 58.96453 109.97072 56.62848 96.76775
#> [163] 64.52939 70.62560 57.56633 57.15516 82.78148 100.20222
#> [169] 56.02534 60.45021 63.05571 67.83253 61.85764 103.46718
#> [175] 55.33284 95.39860 113.36405 63.14564 132.32244 56.49724
#> [181] 52.21165 59.33454 63.39767 64.10673 55.61014 68.63642
#> [187] 60.69752 81.32460 56.08820 67.96019 61.23329 65.20826
#> [193] 78.99125 64.53220 270.05582 90.22419 56.41598 56.68091
#> [199] 82.69165 85.31412 78.10724 138.04041 57.97156 65.87653
#> [205] 53.81364 65.83511 97.62379 62.76005 82.69390 65.39563
#> [211] 65.50743 63.36447 68.25642 88.40218 68.73522 65.89664
#> [217] 906.10255 74.37854 243.94549 60.33958 57.57891 61.04660
#> [223] 371.71147 81.53081 120.14502 58.58538 78.85486 54.53393
#> [229] 89.13556 54.95814 72.76784 58.05817 66.88669 58.33097
#> [235] 66.51409 65.01806 60.01535 57.90402 95.50022 60.08558
#> [241] 107.88466 107.59362 57.98122 101.20713 62.64703 70.60320
#> [247] 75.05221 61.30593 149.26326 111.46011 82.30679 62.41808
#> [253] 92.61884 42.10239 117.04211 69.58896 67.39609 56.11135
#> [259] 57.37282 63.44387 62.37468 75.65055 74.49048 483.43435
#> [265] 56.89941 57.87289 57.39240 94.32344 65.62458 99.95849
#> [271] 55.56866 64.10501 551.39213 60.29839 57.13962 159.17194
#> [277] 56.32333 58.61230 116.55436 62.84540 87.55654 176.68873
#> [283] 85.85320 1040.51254 5904.08110 62.32587 81.64137 54.24470
#> [289] 59.81819 66.29120 65.56767 67.62667 70.00540 3503.89736
#> [295] 61.57099 86.90656 55.17506 72.41663 68.90514 56.23281
#> [301] 55.97295 59.29985 57.46276 59.48414 74.79696 63.35265
#> [307] 54.28005 56.87581 64.74544 67.68570 65.64023 59.00302
#> [313] 68.83656 59.78464 60.45285 71.78394 53.89749 71.23675
#> [319] 62.30258 50.19278 61.49108 59.45583 63.32197 68.12958
#> [325] 61.75928 222.59728 99.26716 61.71291 65.98900 59.80339
#> [331] 103.21719 70.19286 812.44321 83.80117 56.81371 68.10798
#> [337] 60.70823 56.74271 92.41249 58.75655 94.49557 63.02977
#> [343] 60.26130 138.98536 67.47095 61.41656 54.90811 62.59706
#> [349] 138.80496 52.97335 58.53523 72.66149 155.75669 66.72650
#> [355] 47.75343 64.09492 60.00507 61.40272 52.28276 136.58688
#> [361] 64.48271 78.69516 152.75645 62.29448 102.56591 63.28079
#> [367] 58.34986 63.06293 75.41320 187.95242 57.53947 105.88412
#> [373] 113.62996 56.45648 55.73723 78.75024 59.21679 110.74336
#> [379] 61.29378 73.63312 67.95056 89.86144 214.37300 48.16516
#> [385] 53.06716 64.22517 125.39297 273.84821 55.35184 67.65566
#> [391] 59.96092 64.90967 66.35357 303.74497 112.74646 71.23549
#> [397] 57.87847 108.55320 63.61598 101.90939 63.28405 60.12779
#> [403] 59.70096 63.61528 65.55598 66.51417 57.63742 53.88886
#> [409] 61.28845 113.52483 84.34788 64.98658 137.45609 216.31470
#> [415] 76.68958 60.47421 54.08959 91.14311 75.71735 58.41753
#> [421] 100.68553 64.78346 81.81448 284.71502 79.04504 65.43628
#> [427] 58.09876 66.37306 62.06892 72.76955 80.90486 76.51550
#> [433] 101.99539 89.35223 63.81739 86.04666 58.49780 94.90294
#> [439] 49.00533 80.92571 57.06311 85.16562 57.57364 56.90219
#> [445] 58.61551 91.05558 53.58351 81.11937 64.46690 58.15945
#> [451] 66.97116 60.23360 59.04352 88.43335 87.08926 90.27032
#> [457] 130.59295 92.10309 587.07039 79.82292 61.45267 59.16678
#> [463] 71.37441 74.30417 56.37480 61.33721 91.69655 408.35677
#> [469] 72.46069 56.23680 540.71604 119.07066 74.23717 70.02855
#> [475] 65.46857 57.58616 59.67516 62.34771 58.11762 57.53635
#> [481] 39.24958 196.30128 92.52800 82.63003 57.00592 58.03785
#> [487] 63.17596 62.29647 61.05522 62.16486 67.71742 58.95649
#> [493] 64.93903 64.29826 593.78797 60.85672 60.44359 43.08845
#> [499] 47.87298 64.18528 65.62638 82.93834 67.53440 97.29500
#> [505] 66.77294 78.37587 68.21701 78.37458 62.92656 82.58186
#> [511] 81.74729 63.00588 98.02423 54.00573 66.20460 54.51327
#> [517] 78.29010 51.45569 49.29737 60.44218 229.96940 60.78924
#> [523] 60.20595 56.22064 63.74312 56.44710 66.17650 65.73346
#> [529] 59.98803 78.74191 89.27767 87.03991 143.23707 58.35595
#> [535] 60.95003 67.42272 58.88189 62.55732 61.66584 65.98537
#> [541] 60.68558 83.40359 74.89846 60.78910 64.86349 59.28809
#> [547] 62.67668 157.76733 96.99390 54.92601 126.31028 58.59443
#> [553] 146.21495 139.41985 61.90790 50.97710 66.94612 64.63998
#> [559] 60.88456 73.43466 56.52146 58.12980 60.24227 63.48206
#> [565] 71.45891 53.01508 58.70818 62.95270 60.82091 64.98426
#> [571] 73.62417 58.94430 70.10213 67.77949 71.30342 61.24449
#> [577] 59.56187 68.98930 75.25602 64.34566 90.53867 93.83366
#> [583] 54.34316 80.28585 56.54716 56.22421 146.84931 103.17429
#> [589] 71.69306 79.01912 46.30480 87.93181 66.34489 98.05678
#> [595] 71.75285 117.57090 55.22993 434.94996 55.99201 61.22376
#> [601] 70.05957 59.77398 72.88863 59.32811 65.62771 234.37066
#> [607] 63.28299 61.75760 398.26291 79.89100 88.93926 44.85969
#> [613] 55.99895 83.03332 65.59845 58.10360 72.54316 55.27297
#> [619] 59.70071 93.18740 51.57728 60.60780 80.21974 115.50697
#> [625] 59.25579 61.73548 86.92699 80.42023 78.03147 61.64894
#> [631] 55.15714 69.61260 66.53630 78.77967 57.35921 105.18741
#> [637] 64.78381 354.54800 139.10597 65.90683 68.29499 272.43653
#> [643] 108.97343 60.73186 65.89386 60.85356 76.73468 61.75008
#> [649] 67.80939 58.63498 62.44710 60.42062 62.44688 57.53329
#> [655] 223.65347 60.81615 67.24969 57.90360 129.09842 64.15548
#> [661] 85.43353 60.82194 96.35664 52.99041 81.03180 56.43737
#> [667] 54.49476 58.91199 61.27251 70.07353 111.52314 55.34813
#> [673] 59.43840 71.04766 65.22754 62.68011 63.52163 66.88617
#> [679] 56.64266 60.10585 64.32727 80.45664 58.94388 66.50201
#> [685] 64.77660 66.23940 64.98305 89.12868 55.17222 58.40337
#> [691] 55.19594 107.81479 58.59145 55.20567 61.20694 69.06240
#> [697] 59.50156 69.28418 63.76155 74.52912 94.27234 298.56744
#> [703] 57.79221 74.31831 63.15203 92.76062 85.53226 120.44334
#> [709] 415.37382 2752.31391 69.93471 61.41468 60.11325 75.40929
#> [715] 58.17229 59.30405 66.10023 58.49136 62.20742 62.40194
#> [721] 72.91947 45.49546 59.26668 66.49671 55.60718 68.89735
#> [727] 62.73073 61.06370 51.99651 66.44160 68.96675 101.91424
#> [733] 49.14694 662.58670 164.87532 58.44302 61.36975 56.45583
#> [739] 77.79174 145.48051 298.14109 64.43292 72.42921 77.73543
#> [745] 78.57079 59.53665 69.49569 83.30538 59.72551 56.66215
#> [751] 73.57941 98.39935 52.02365 67.17572 58.67629 60.72470
#> [757] 84.39788 71.29286 62.00134 56.89142 234.79645 83.57952
#> [763] 79.44482 58.01828 365.52725 61.21201 84.50582 109.45586
#> [769] 75.68387 57.24023 60.50978 69.02533 62.26353 61.94859
#> [775] 58.00983 79.52591 66.14986 79.68319 68.05019 59.15210
#> [781] 56.29077 54.55493 56.65139 63.27665 80.94148 50.72626
#> [787] 157.90667 173.52921 66.53275 120.69659 64.43076 55.65141
#> [793] 55.52114 58.92944 89.17589 66.52180 65.13499 51.07720
#> [799] 57.65693 65.75856 74.22537 78.77375 64.44467 66.36361
#> [805] 59.21759 62.01656 66.89208 61.35747 128.41575 73.13410
#> [811] 63.45161 61.83009 69.06661 91.37768 67.50091 78.46533
#> [817] 256.57226 57.52230 110.82500 101.86378 58.53693 64.80621
#> [823] 63.82428 63.61650 61.89551 89.74130 65.14876 74.46827
#> [829] 78.29345 74.37315 55.06481 53.32877 140.74059 84.89692
#> [835] 61.65036 63.42738 58.58342 67.53310 68.19601 48.60582
#> [841] 60.62840 71.80310 61.47381 57.03178 56.61877 91.24950
#> [847] 62.51798 66.19086 73.68933 55.45910 173.04320 117.35445
#> [853] 128.11051 58.95766 62.64630 53.98771 70.50179 51.77061
#> [859] 78.80546 57.04740 61.38853 61.42427 62.21268 63.83533
#> [865] 77.92282 71.23109 62.28017 74.87377 379.69351 65.74705
#> [871] 69.38541 202.71318 275.37773 83.22846 56.89279 52.80945
#> [877] 59.85899 57.82899 66.00752 90.77465 107.62937 54.16490
#> [883] 57.97550 65.94674 78.28211 56.37684 64.75285 77.07404
#> [889] 59.45131 62.89652 66.09172 58.72717 62.53264 71.58606
#> [895] 59.21477 66.65127 66.07406 81.16637 71.11051 65.88938
#> [901] 70.16516 70.49907 66.26591 78.78447 58.14950 93.41018
#> [907] 474.47719 257.92662 60.45285 60.93893 52.09350 93.69644
#> [913] 64.06590 60.99684 60.41475 64.31892 60.18601 64.52942
#> [919] 75.12758 55.31311 61.36428 57.69772 56.01890 80.21764
#> [925] 72.16787 89.20888 55.46961 57.74418 59.31447 62.95632
#> [931] 58.14193 92.07945 54.47517 56.76946 63.81612 65.38471
#> [937] 54.32713 81.49105 114.53898 62.76802 58.03340 55.79954
#> [943] 83.65230 57.76700 58.29728 71.02727 61.68623 290.29357
#> [949] 407.78542 67.25245 56.56480 325.58183 62.56539 59.21602
#> [955] 97.57802 57.61837 53.85645 68.81045 62.65571 63.20877
#> [961] 67.72495 67.72347 60.26829 80.82234 56.21289 47.03960
#> [967] 62.32642 66.91910 63.06483 59.57222 60.24422 217.66059
#> [973] 1838.34993 57.04000 71.80408 64.63546 57.88427 63.97291
#> [979] 48.73827 49.49163 48.72252 84.36203 53.19090 86.79845
#> [985] 120.47690 465.32329 174.30046 62.58361 69.04008 73.77577
#> [991] 121.36551 88.54213 288.40691 63.38775 70.57615 60.66321
#> [997] 67.63785 67.26060 65.31551 58.46773