![[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] 3.159029e+02 1.758503e+02 2.433804e+02 6.507361e+01 6.556283e+01
#> [6] 5.981034e+01 1.385670e+03 8.801113e+02 2.654958e+02 8.917769e+01
#> [11] 1.878383e+02 6.327540e+01 1.364458e+02 1.091114e+03 7.194881e+01
#> [16] 1.852051e+02 7.260585e+01 7.484790e+05 8.912171e+01 6.790513e+01
#> [21] 8.038794e+01 1.051371e+02 6.816110e+01 8.425552e+01 7.157370e+01
#> [26] 1.035973e+02 2.226990e+02 2.691375e+02 7.001572e+01 6.245211e+01
#> [31] 6.312846e+01 1.322834e+02 7.546320e+01 7.220080e+01 1.948866e+02
#> [36] 1.577088e+02 6.664593e+01 1.141888e+02 1.269675e+02 1.889124e+02
#> [41] 9.264081e+01 2.547111e+02 6.845880e+01 7.347089e+01 7.687758e+01
#> [46] 6.955615e+01 6.874093e+01 1.508662e+02 7.892005e+01 1.295606e+02
#> [51] 1.452272e+02 6.805623e+01 6.397053e+01 7.138421e+01 7.875534e+01
#> [56] 5.900010e+01 7.289789e+01 6.305500e+01 6.714884e+01 1.367880e+02
#> [61] 8.138374e+01 6.818673e+01 8.664477e+01 8.307194e+01 1.961702e+02
#> [66] 8.590197e+01 8.354413e+01 1.801372e+02 9.847526e+01 2.733357e+02
#> [71] 1.352479e+03 2.766433e+02 9.987436e+01 6.627684e+01 7.911778e+01
#> [76] 8.120167e+01 7.398994e+01 7.407307e+01 7.540040e+01 7.338235e+01
#> [81] 8.153939e+01 1.601976e+02 3.109438e+02 5.097632e+04 2.547398e+04
#> [86] 2.823141e+03 4.960943e+02 2.209601e+03 8.756956e+01 7.024901e+01
#> [91] 7.363569e+01 7.548897e+01 8.204585e+01 7.829883e+01 9.173278e+01
#> [96] 1.452844e+02 8.470013e+01 8.802289e+01 2.176142e+03 1.044155e+03
#> [101] 8.347081e+01 1.014707e+02 1.125378e+02 1.897731e+02 6.930331e+01
#> [106] 7.797970e+01 6.226513e+01 2.159381e+02 6.189248e+01 8.127682e+01
#> [111] 8.163920e+01 9.118886e+01 6.816294e+01 3.427124e+02 6.541786e+01
#> [116] 1.346910e+02 6.690947e+01 1.013463e+04 8.903017e+01 8.459338e+01
#> [121] 1.331637e+02 2.532727e+02 1.083921e+02 7.045206e+01 6.272024e+01
#> [126] 8.859706e+01 2.373233e+02 1.361159e+02 1.802203e+02 1.125671e+02
#> [131] 1.047607e+03 8.216319e+01 6.770642e+01 5.976153e+02 1.759763e+02
#> [136] 1.003588e+02 6.704821e+01 7.072015e+01 1.211003e+02 9.011208e+01
#> [141] 1.108452e+02 7.177340e+01 1.171588e+02 8.068649e+01 2.227478e+02
#> [146] 1.497793e+02 7.051855e+01 1.017048e+02 7.271037e+01 8.565784e+01
#> [151] 1.325962e+02 7.976470e+01 7.857427e+01 7.770870e+01 6.948088e+01
#> [156] 1.181537e+02 1.502510e+02 8.690050e+01 1.223197e+02 2.853309e+02
#> [161] 5.923560e+01 7.769287e+01 6.712569e+01 5.581582e+02 6.917967e+01
#> [166] 7.405985e+01 6.304982e+01 6.490477e+01 2.050738e+02 2.636145e+02
#> [171] 1.809376e+02 2.274038e+02 1.841618e+02 9.619812e+01 9.997606e+01
#> [176] 7.289345e+01 6.324945e+01 7.903968e+01 3.609686e+02 7.160034e+01
#> [181] 2.484845e+02 5.719881e+05 2.659878e+02 1.107857e+02 7.898689e+01
#> [186] 6.943270e+01 1.500528e+02 9.392081e+01 8.758465e+01 8.961333e+01
#> [191] 9.134379e+01 7.550421e+01 6.883387e+01 6.296474e+01 6.468697e+01
#> [196] 6.821672e+01 7.787550e+01 7.410851e+01 6.563831e+02 6.556572e+01
#> [201] 9.411901e+01 8.309513e+01 6.600243e+01 8.310446e+01 9.434452e+01
#> [206] 9.872694e+01 1.029953e+02 1.114067e+02 3.665041e+03 8.906651e+01
#> [211] 1.170493e+02 7.020583e+01 7.098486e+01 1.303096e+02 2.208698e+02
#> [216] 3.569672e+02 6.723921e+01 1.900647e+02 8.203511e+01 8.334184e+01
#> [221] 8.828997e+01 1.096508e+02 7.872990e+01 1.088579e+02 1.280542e+02
#> [226] 7.598972e+01 1.988354e+02 9.228563e+01 7.975393e+01 6.379529e+01
#> [231] 1.090663e+02 6.795759e+01 7.439867e+01 9.471234e+01 1.241583e+02
#> [236] 6.448984e+01 7.957753e+01 9.054315e+01 9.592316e+01 9.352042e+01
#> [241] 1.136576e+02 9.960568e+01 1.870799e+02 7.369099e+02 8.284073e+01
#> [246] 2.548431e+02 6.891848e+01 8.000453e+01 8.974660e+01 9.488924e+01
#> [251] 8.546283e+01 5.564112e+02 2.030427e+02 6.755822e+01 6.990569e+01
#> [256] 1.283985e+02 2.169039e+02 2.722552e+02 8.049393e+01 7.195769e+01
#> [261] 6.299733e+01 6.459323e+01 6.365644e+01 7.323502e+01 4.785386e+03
#> [266] 7.818688e+01 9.461832e+01 6.966401e+02 3.426209e+02 2.523446e+02
#> [271] 1.353791e+02 7.102539e+01 8.356193e+01 1.018145e+02 8.236292e+01
#> [276] 1.330638e+04 8.505381e+01 6.928006e+01 1.677598e+02 6.940430e+01
#> [281] 2.501538e+02 1.194744e+02 7.526196e+01 6.706747e+01 1.310093e+02
#> [286] 1.012375e+02 7.965530e+01 6.828846e+01 6.547642e+01 1.325634e+02
#> [291] 1.235223e+02 5.708486e+02 1.347498e+02 8.396540e+01 7.495590e+01
#> [296] 1.109118e+02 6.682475e+01 8.526268e+01 9.049041e+01 3.476737e+02
#> [301] 7.472678e+01 7.663184e+01 4.518041e+02 1.075648e+02 7.842367e+01
#> [306] 1.002846e+02 6.284607e+01 7.513502e+01 1.748335e+02 9.826352e+01
#> [311] 9.219074e+01 1.036214e+02 1.570057e+02 7.081989e+01 1.375906e+02
#> [316] 8.528287e+01 1.251022e+02 7.045556e+01 8.471929e+01 1.155670e+02
#> [321] 6.156616e+01 2.196208e+02 8.137681e+01 7.172465e+01 1.133614e+02
#> [326] 6.996817e+01 9.023569e+01 1.227175e+02 2.771150e+02 7.331507e+01
#> [331] 7.358499e+01 8.283156e+01 1.150644e+02 6.343340e+01 6.977150e+01
#> [336] 9.249104e+01 1.215128e+02 1.717979e+02 1.010520e+02 7.316239e+01
#> [341] 5.444979e+02 1.048738e+02 8.192753e+01 1.420224e+02 2.375794e+02
#> [346] 1.088485e+02 1.050083e+02 3.152628e+02 8.452469e+02 1.596059e+02
#> [351] 1.964770e+02 6.971343e+01 8.993250e+01 1.078381e+02 1.293644e+02
#> [356] 1.656428e+02 6.507471e+01 8.814906e+01 1.437836e+02 7.329239e+01
#> [361] 9.104635e+01 6.228683e+01 7.103365e+01 7.442044e+01 6.479052e+01
#> [366] 7.007139e+01 7.669077e+01 6.919302e+01 9.278334e+01 9.459692e+01
#> [371] 9.268594e+01 1.471891e+02 1.112091e+02 1.082301e+02 9.867399e+01
#> [376] 7.106003e+01 7.452734e+01 6.671701e+01 1.511570e+02 6.435848e+01
#> [381] 1.112946e+02 6.184598e+01 1.751601e+02 7.206911e+01 6.778435e+01
#> [386] 8.511427e+01 1.466694e+02 7.394954e+01 7.841041e+01 7.062897e+01
#> [391] 6.591991e+01 6.995680e+01 3.132634e+02 1.045107e+02 1.805130e+02
#> [396] 4.403605e+02 5.825744e+02 6.677703e+01 6.769409e+01 7.895642e+01
#> [401] 7.742099e+01 9.788086e+01 6.488240e+01 7.963322e+01 7.466903e+01
#> [406] 7.159631e+01 7.841972e+01 9.198058e+01 5.876001e+01 6.154896e+01
#> [411] 8.730751e+01 1.339219e+02 3.393317e+02 8.699514e+01 8.731614e+01
#> [416] 8.723432e+01 8.704059e+01 6.597424e+01 1.071190e+02 1.322154e+02
#> [421] 7.406209e+01 8.354542e+01 1.926984e+02 6.663340e+01 9.135459e+01
#> [426] 8.516959e+01 9.421696e+01 6.878514e+01 2.342498e+04 1.873148e+02
#> [431] 1.940444e+03 6.309014e+02 1.670823e+03 7.307707e+01 7.500121e+01
#> [436] 7.516589e+01 1.251398e+02 1.272122e+02 6.768343e+01 1.082653e+02
#> [441] 7.671607e+01 8.502265e+02 1.563777e+03 1.280241e+02 8.882869e+01
#> [446] 1.046363e+02 6.926020e+01 1.049316e+02 1.838507e+02 6.755975e+01
#> [451] 1.132954e+02 9.835337e+01 9.813732e+01 6.735909e+01 6.235383e+01
#> [456] 9.903995e+01 1.122145e+02 1.427916e+02 1.663352e+02 1.409747e+02
#> [461] 7.392940e+01 1.034321e+03 7.402404e+01 7.097266e+01 1.075105e+02
#> [466] 1.564894e+02 6.344301e+01 6.483286e+01 7.749067e+01 1.503880e+02
#> [471] 7.029682e+01 2.328584e+02 2.035267e+02 6.861087e+01 7.605767e+01
#> [476] 9.938666e+01 6.535452e+01 8.782397e+01 6.852706e+01 6.543135e+01
#> [481] 7.638873e+01 1.200413e+02 6.472481e+01 7.650039e+01 7.256502e+01
#> [486] 9.066626e+01 2.131047e+02 9.273343e+01 7.990846e+01 6.300703e+01
#> [491] 9.160641e+01 9.106449e+01 8.845756e+01 1.334541e+02 8.181038e+01
#> [496] 8.673056e+01 1.079937e+02 7.466063e+01 1.992180e+02 1.211622e+02
#> [501] 7.763869e+01 7.257178e+01 8.432347e+01 9.176566e+01 1.438565e+03
#> [506] 7.774697e+01 9.047485e+01 6.974713e+01 6.610054e+01 1.111373e+02
#> [511] 7.804196e+01 9.461591e+01 1.157075e+02 1.080116e+02 7.794771e+01
#> [516] 9.002519e+01 1.762425e+02 1.229893e+02 7.101463e+01 6.838450e+01
#> [521] 7.776505e+01 1.402869e+02 9.112101e+01 1.127314e+02 7.637365e+01
#> [526] 1.075083e+02 1.045113e+02 9.204035e+01 1.927598e+02 7.471521e+01
#> [531] 1.093365e+02 7.099533e+01 5.788640e+01 6.611934e+01 1.058138e+02
#> [536] 1.107478e+02 7.040229e+01 8.806157e+01 2.113127e+02 6.970482e+01
#> [541] 8.063933e+01 8.552717e+01 6.266340e+01 1.867042e+02 3.485121e+02
#> [546] 2.659178e+02 2.454782e+04 4.716539e+03 2.114114e+02 7.396625e+01
#> [551] 1.008627e+02 7.412649e+01 2.154659e+02 1.538950e+02 1.062885e+02
#> [556] 3.133211e+02 2.034119e+02 1.500876e+03 1.150404e+02 7.035605e+01
#> [561] 1.047890e+02 9.298118e+01 7.560613e+01 1.385122e+02 1.667888e+02
#> [566] 1.717988e+02 6.725442e+01 7.739740e+01 7.291666e+01 9.529225e+01
#> [571] 7.848816e+01 6.376323e+01 2.127350e+02 3.328024e+04 2.668815e+03
#> [576] 7.690464e+01 1.054644e+02 1.285059e+02 8.346070e+01 5.646198e+02
#> [581] 2.220452e+05 1.371377e+05 6.230330e+01 8.307629e+02 1.114870e+02
#> [586] 1.240015e+02 6.919463e+01 8.808962e+01 6.780987e+01 8.340838e+01
#> [591] 9.955135e+01 8.206505e+01 7.983099e+01 1.110665e+02 2.278716e+02
#> [596] 7.018846e+01 1.550209e+02 7.216082e+01 7.233480e+01 1.557941e+02
#> [601] 7.063715e+01 9.072143e+01 8.871776e+01 1.583391e+02 7.764327e+01
#> [606] 7.515724e+01 3.263850e+02 1.570959e+03 6.325706e+03 7.394628e+01
#> [611] 7.506974e+01 8.358759e+01 2.730507e+02 7.797440e+01 8.263560e+01
#> [616] 7.144961e+01 9.607582e+01 9.077355e+01 1.327923e+02 3.646310e+02
#> [621] 9.958970e+01 1.569009e+02 1.777400e+02 6.886561e+01 2.268608e+02
#> [626] 9.968693e+01 1.590980e+02 1.158541e+02 9.655932e+01 9.776692e+01
#> [631] 1.443960e+02 1.127246e+02 8.941751e+01 7.011054e+01 7.587698e+01
#> [636] 2.027223e+02 7.387310e+01 2.417604e+02 1.768463e+02 1.614680e+03
#> [641] 6.820326e+01 6.948056e+01 7.435368e+01 1.231870e+02 3.195673e+02
#> [646] 7.153694e+01 8.406109e+01 1.434924e+02 7.038022e+01 1.939568e+02
#> [651] 7.270780e+01 7.994738e+01 1.168520e+02 6.780324e+01 6.161507e+01
#> [656] 5.510664e+02 1.171325e+03 3.446588e+02 3.379177e+06 9.702813e+01
#> [661] 7.649243e+01 5.474208e+02 3.064182e+02 3.282072e+02 6.819608e+01
#> [666] 7.538903e+01 7.273217e+01 6.223584e+01 1.002040e+02 1.776891e+02
#> [671] 7.499877e+01 1.411177e+02 6.587219e+01 1.382474e+02 1.419149e+03
#> [676] 7.008234e+01 7.897749e+01 1.250712e+02 1.353347e+02 9.344896e+01
#> [681] 1.650971e+02 1.284265e+02 6.165336e+01 7.213666e+01 1.264840e+02
#> [686] 1.543943e+02 1.695210e+02 1.182509e+02 6.619217e+01 9.054212e+01
#> [691] 9.208798e+01 9.131328e+01 6.530324e+01 8.305433e+01 6.617451e+01
#> [696] 7.018190e+01 8.291075e+01 8.595137e+01 1.377859e+02 1.412180e+02
#> [701] 9.431272e+01 1.152128e+02 2.203294e+02 6.269649e+01 7.191104e+01
#> [706] 9.396810e+01 1.059261e+02 2.550478e+02 6.452341e+02 7.391433e+01
#> [711] 7.606877e+01 8.483708e+01 7.635737e+01 8.014363e+01 8.674863e+01
#> [716] 8.441708e+01 1.345190e+02 1.308285e+03 9.045520e+01 7.224311e+01
#> [721] 1.989251e+02 2.577996e+03 7.096108e+01 7.918693e+01 1.037239e+02
#> [726] 7.413791e+01 8.470605e+01 8.031796e+01 7.068118e+01 6.430509e+01
#> [731] 2.872908e+02 1.409052e+02 8.307649e+01 7.811899e+01 7.167800e+01
#> [736] 1.073917e+02 1.132337e+02 1.017709e+02 7.103610e+01 7.959046e+01
#> [741] 3.208300e+02 1.734167e+03 4.562655e+02 6.604821e+01 8.859456e+01
#> [746] 6.468539e+01 7.324505e+01 1.882297e+02 1.204125e+02 1.454431e+02
#> [751] 7.267990e+01 7.040866e+01 6.925556e+01 9.500590e+01 1.560094e+03
#> [756] 1.671479e+03 1.216529e+02 6.836796e+01 6.163489e+01 2.256236e+02
#> [761] 8.732863e+01 7.048979e+01 6.658179e+01 7.506155e+01 1.046680e+02
#> [766] 8.300056e+01 7.457980e+01 7.708832e+01 8.731139e+01 8.232153e+01
#> [771] 8.038660e+01 1.415364e+02 8.368062e+01 7.886711e+01 6.178410e+01
#> [776] 6.441493e+01 1.852018e+02 7.416580e+01 2.501780e+02 2.253904e+02
#> [781] 1.137073e+02 9.149106e+01 6.296443e+01 7.095329e+01 1.070684e+02
#> [786] 9.771503e+01 7.009359e+01 7.305948e+01 7.108254e+01 9.010918e+01
#> [791] 7.903290e+01 4.419062e+03 8.476174e+01 6.263028e+01 9.618620e+01
#> [796] 9.182219e+01 1.199990e+02 6.501965e+01 1.292996e+02 7.217719e+01
#> [801] 1.122371e+02 6.254541e+01 7.850032e+01 6.915633e+01 6.628717e+01
#> [806] 8.056113e+01 1.251129e+03 6.905944e+01 2.910290e+02 2.812032e+02
#> [811] 7.049043e+01 1.237221e+02 1.710750e+02 2.397049e+03 9.593838e+01
#> [816] 9.398546e+01 7.044073e+01 6.213229e+01 6.708913e+01 6.830723e+01
#> [821] 1.221940e+02 1.726166e+02 1.447457e+04 6.581106e+01 3.505360e+02
#> [826] 8.862698e+01 7.872403e+01 1.635294e+03 1.724186e+02 2.918847e+02
#> [831] 7.050081e+01 7.998863e+01 1.042729e+02 3.108549e+02 1.976249e+02
#> [836] 4.736249e+02 1.726942e+02 1.100090e+02 9.571341e+01 1.133483e+02
#> [841] 1.079499e+02 1.117701e+02 6.332580e+02 7.993148e+01 6.024365e+01
#> [846] 8.220792e+01 7.323742e+01 6.753089e+01 8.874759e+01 7.645786e+01
#> [851] 5.731689e+08 3.705345e+02 6.487899e+01 7.492455e+01 5.549980e+02
#> [856] 1.686572e+02 7.214260e+01 1.069366e+02 8.470673e+01 1.004477e+02
#> [861] 6.995808e+01 1.341742e+02 1.734260e+02 7.782591e+01 9.074676e+01
#> [866] 7.603757e+01 7.720845e+01 6.836855e+01 1.420338e+02 7.045056e+01
#> [871] 7.512305e+01 6.818474e+01 9.174299e+01 3.308209e+02 1.391805e+02
#> [876] 1.217644e+02 7.147514e+01 8.834882e+01 8.278347e+01 5.352918e+02
#> [881] 6.663436e+01 9.531434e+01 7.068487e+01 8.671516e+01 7.516888e+01
#> [886] 8.494232e+01 7.011631e+01 9.633253e+01 9.838351e+01 4.467892e+02
#> [891] 6.640662e+01 8.771972e+01 1.377270e+02 8.887944e+01 6.245714e+01
#> [896] 1.272442e+02 1.132335e+02 1.005148e+02 2.472082e+02 8.097069e+01
#> [901] 9.309179e+01 7.112725e+01 6.460869e+01 8.873448e+01 2.116725e+02
#> [906] 9.392442e+01 7.046930e+01 9.260856e+01 6.653661e+01 6.012245e+02
#> [911] 6.996503e+01 9.611694e+01 7.278605e+01 2.159909e+02 1.826590e+02
#> [916] 2.585413e+02 6.666582e+01 1.329104e+02 6.848168e+01 6.354093e+01
#> [921] 3.681885e+02 1.328828e+02 6.845278e+01 7.207133e+01 1.107555e+02
#> [926] 7.532117e+01 8.408383e+01 2.054883e+02 8.521985e+01 1.097261e+02
#> [931] 6.774474e+01 8.274175e+01 1.791662e+02 6.916728e+01 9.303365e+01
#> [936] 1.355918e+02 1.015926e+02 1.449952e+02 2.696426e+02 6.378451e+01
#> [941] 6.757723e+01 9.605345e+01 1.827128e+02 2.096918e+02 7.757962e+01
#> [946] 8.504229e+01 6.537158e+01 1.222229e+02 1.264327e+02 1.954307e+02
#> [951] 6.605243e+01 7.727445e+01 7.786819e+02 1.479183e+02 8.802236e+02
#> [956] 3.632250e+02 1.590239e+02 9.250870e+01 8.270223e+01 1.363853e+02
#> [961] 6.035389e+02 9.298260e+01 1.772933e+02 6.723720e+01 6.965714e+01
#> [966] 1.086270e+02 2.163233e+02 1.083037e+02 1.530397e+02 1.439211e+02
#> [971] 7.589914e+01 1.730795e+02 8.015887e+01 1.304670e+02 5.233108e+08
#> [976] 9.734344e+01 7.195047e+01 1.112445e+02 9.537395e+01 1.575316e+02
#> [981] 1.024070e+02 1.004506e+02 1.395924e+02 1.651101e+02 7.527467e+01
#> [986] 8.000290e+01 1.474897e+02 8.261646e+01 6.477160e+01 1.412584e+02
#> [991] 9.878064e+01 1.068740e+04 1.125657e+06 2.343191e+05 1.850430e+02
#> [996] 8.678238e+01 3.448430e+03 2.513123e+02 8.672458e+01 6.838547e+01
doseFunction(
model,
alpha2 = samples@data$alpha2,
beta = samples@data$beta,
grade = 2L
)(x = 0.25)
#> [1] 78.56608 78.73103 142.89075 57.98986 58.57536
#> [6] 54.75326 410.19021 228.38086 122.07827 58.00360
#> [11] 87.26821 54.27537 38.21157 174.12513 56.60885
#> [16] 96.15559 58.95598 3301.07320 71.52505 58.09754
#> [21] 66.19632 75.63268 55.39225 70.97361 63.97652
#> [26] 70.36329 70.85237 106.69814 56.43263 56.60021
#> [31] 58.79106 99.23489 64.17955 58.52191 118.69827
#> [36] 80.40486 57.55206 74.34044 74.98398 100.50335
#> [41] 67.80736 114.08505 53.25237 60.94349 64.01678
#> [46] 59.07456 59.14364 97.71037 54.57535 68.83321
#> [51] 91.58263 56.46463 57.86359 63.58067 65.98891
#> [56] 51.14511 61.99835 55.91316 58.54643 88.79645
#> [61] 55.23151 56.58822 60.40066 64.81171 117.19586
#> [66] 67.03815 65.22844 108.67556 68.33810 112.30739
#> [71] 215.39076 110.22919 67.10681 56.52201 52.39259
#> [76] 61.00811 58.81213 61.83954 62.51931 61.02161
#> [81] 68.12213 105.62277 178.77622 750.85314 1277.89665
#> [86] 123.44249 155.29971 369.86927 66.69295 54.98320
#> [91] 59.75087 62.92350 66.36840 69.11421 78.03833
#> [96] 64.34211 68.67436 66.69408 665.00138 217.31210
#> [101] 64.80278 71.00363 74.49942 74.49535 59.92911
#> [106] 58.06223 57.01310 105.23544 54.94123 60.67778
#> [111] 66.02729 73.97212 59.27748 106.52889 54.78283
#> [116] 58.14423 54.29602 702.72740 70.37137 67.15494
#> [121] 83.94809 133.09934 59.60314 59.60875 56.17805
#> [126] 64.53537 104.82489 94.61215 86.50108 82.55765
#> [131] 256.92917 65.01797 57.71670 180.30165 84.58017
#> [136] 71.72277 59.11364 63.23186 95.75101 73.21970
#> [141] 69.70482 57.23529 71.67333 63.57369 55.60709
#> [146] 71.53436 59.64922 73.21487 62.67152 58.65419
#> [151] 93.39332 62.42267 57.48508 68.69827 62.70875
#> [156] 83.95238 83.85693 61.39019 76.84912 154.00006
#> [161] 52.65553 58.03138 58.53338 165.50844 59.31989
#> [166] 62.57707 57.02387 58.09972 117.20245 134.41550
#> [171] 109.10811 86.10106 91.44867 70.40135 69.33305
#> [176] 64.02532 57.63799 60.23760 114.24924 55.67836
#> [181] 64.81668 2321.88104 133.50784 84.49367 63.42347
#> [186] 59.42398 54.35534 69.30908 66.37947 71.51198
#> [191] 74.14426 62.84397 57.26223 55.53282 56.18934
#> [196] 59.79296 64.13891 57.01797 188.82425 54.90800
#> [201] 70.22797 61.07744 55.07558 67.21118 69.18081
#> [206] 72.82513 75.37008 63.17646 957.46888 70.05627
#> [211] 78.93842 61.84953 62.72886 82.34112 91.60581
#> [216] 81.74356 56.67833 103.17766 66.43280 64.79287
#> [221] 65.36557 78.13734 62.29077 84.74404 68.73389
#> [226] 57.10804 118.29743 69.48721 65.79963 57.02631
#> [231] 61.25413 58.74835 61.46745 73.47505 87.56684
#> [236] 56.43255 63.77518 62.86836 59.92079 62.49316
#> [241] 71.94875 62.04385 85.54241 224.73156 65.20385
#> [246] 130.85036 59.79764 68.63078 69.29993 68.44381
#> [251] 53.39350 227.86638 99.25318 57.86431 59.35024
#> [256] 62.65771 114.70271 120.45145 64.28712 55.15667
#> [261] 57.80303 59.14754 56.82547 64.49521 227.84421
#> [266] 57.21623 49.85948 199.20573 87.16033 68.82185
#> [271] 92.47466 55.71381 62.95925 62.80742 66.71406
#> [276] 626.83835 64.75034 57.82769 72.10570 57.20407
#> [281] 129.79901 63.50497 64.75489 59.42475 79.63196
#> [286] 75.26533 58.34138 59.33086 58.16763 79.26343
#> [291] 85.51665 220.70411 82.39653 62.83991 61.68623
#> [296] 84.04553 57.03782 66.32215 75.33568 76.94952
#> [301] 62.11117 59.33269 110.23054 79.69330 60.61629
#> [306] 72.54461 51.63916 56.32989 105.19910 69.40575
#> [311] 62.82062 75.99003 77.61907 58.64688 82.88905
#> [316] 71.46614 85.12735 61.52478 66.60254 74.81605
#> [321] 56.02975 123.29625 67.00255 58.16789 75.12806
#> [326] 61.08515 75.14550 78.25721 128.67011 58.32456
#> [331] 63.68991 63.15965 74.26246 56.06106 59.63311
#> [336] 69.18757 67.61812 86.37117 68.43513 61.06992
#> [341] 215.59149 75.45422 56.88701 78.21884 131.00006
#> [346] 59.57501 63.78233 132.97958 283.16112 78.40889
#> [351] 105.62175 57.43861 64.90389 75.41785 89.84202
#> [356] 76.18180 59.06022 67.72164 66.84292 61.47646
#> [361] 62.48256 52.76762 59.15345 50.60679 58.38412
#> [366] 59.69018 58.89812 60.56973 72.73996 66.98948
#> [371] 69.77326 86.35761 73.05829 75.57577 69.08952
#> [376] 60.14811 64.07477 56.64240 79.18146 49.85451
#> [381] 72.22404 56.69933 113.83598 60.35950 57.79267
#> [386] 63.41526 82.28484 64.08114 67.35820 59.13348
#> [391] 58.01260 54.13817 114.22466 63.39824 87.25460
#> [396] 159.22208 126.96419 56.34213 57.75195 63.21648
#> [401] 63.99826 76.92238 58.60077 67.62703 64.16046
#> [406] 62.12033 64.14715 72.23849 50.44145 55.37282
#> [411] 71.82356 96.71692 155.58057 73.36449 68.81493
#> [416] 69.20413 68.83555 57.90160 76.28523 74.54664
#> [421] 59.57211 71.64020 94.24470 57.94436 76.57631
#> [426] 74.39489 81.93507 64.20303 2744.63488 116.02271
#> [431] 413.39432 99.93624 336.77729 57.25019 55.57394
#> [436] 61.28241 87.14104 79.32681 56.85695 76.29392
#> [441] 62.54281 72.84080 234.53041 85.02224 65.63421
#> [446] 62.93827 53.22572 80.17761 87.83503 58.18026
#> [451] 73.09235 61.93735 74.07385 56.58881 56.02041
#> [456] 72.21286 76.86408 73.19574 78.43949 87.51242
#> [461] 58.50528 127.87601 62.65595 61.88774 78.84566
#> [466] 99.37057 56.86726 58.80660 64.88491 78.27800
#> [471] 61.05187 123.46120 116.36547 56.39395 60.00588
#> [476] 71.85407 51.07944 73.42712 61.27402 56.06392
#> [481] 57.12620 72.30216 57.77729 58.20558 55.16102
#> [486] 57.63687 87.98002 63.09531 64.78637 54.87390
#> [491] 70.04489 71.57388 57.47470 80.25493 55.24307
#> [496] 55.21468 64.34611 61.14475 89.25154 79.68658
#> [501] 65.44406 60.45258 68.61417 75.83316 250.13640
#> [506] 60.05546 58.51124 59.89111 58.79765 75.57280
#> [511] 62.65528 74.14397 78.95426 72.61212 60.05384
#> [516] 66.19289 69.97282 84.32533 61.65423 57.47149
#> [521] 58.62134 89.34588 58.47217 69.09605 63.73962
#> [526] 79.66187 74.37610 73.27140 121.09154 61.98271
#> [531] 83.48030 64.23332 52.54816 60.12705 83.92036
#> [536] 81.05974 59.19635 70.14119 110.11877 54.12634
#> [541] 63.86890 66.63646 55.14390 104.36561 165.33832
#> [546] 139.92732 815.38634 601.75978 82.71665 59.52127
#> [551] 67.19991 59.31341 111.20269 91.63697 68.19364
#> [556] 96.92225 66.99541 347.74171 76.39026 59.24789
#> [561] 52.29771 59.10800 60.92973 87.75284 95.30386
#> [566] 101.71896 55.76364 66.13956 64.05415 67.92655
#> [571] 60.10332 51.17858 90.50452 1893.56507 399.18554
#> [576] 59.63311 76.81030 77.74444 70.02697 133.79091
#> [581] 3341.93579 986.65015 57.70737 323.60937 77.81039
#> [586] 80.08415 58.47923 57.21595 59.26792 67.70045
#> [591] 67.40703 62.10700 61.75787 77.06984 82.22865
#> [596] 60.28443 92.57688 62.46637 62.57824 89.57792
#> [601] 57.89605 60.81552 60.57505 81.17497 60.07413
#> [606] 60.14429 71.08176 257.65939 475.48088 59.90264
#> [611] 58.71314 57.34573 110.17761 57.89188 69.04194
#> [616] 61.76245 64.77516 70.56028 57.10811 120.37531
#> [621] 51.76859 51.25816 74.07435 55.27252 129.06231
#> [626] 68.91855 73.29325 56.33618 58.67599 65.13315
#> [631] 95.94178 76.20287 60.86169 61.86212 66.28459
#> [636] 98.68353 54.44343 80.61684 98.70927 243.66547
#> [641] 57.77145 59.04167 52.45390 79.98136 119.82938
#> [646] 58.79453 64.39208 89.27624 55.12114 98.54683
#> [651] 61.21068 65.84576 71.08391 59.72049 54.81805
#> [656] 204.09060 350.44744 125.98795 4147.27684 63.87324
#> [661] 59.23942 213.33628 131.64015 123.45192 62.31392
#> [666] 63.47065 61.42152 55.61556 73.14468 95.64731
#> [671] 63.48781 96.04173 58.39039 70.15947 168.57930
#> [676] 61.69756 69.07479 87.37491 90.65986 67.30834
#> [681] 71.60648 80.44006 49.45524 59.76539 65.35054
#> [686] 107.25458 86.06261 82.68732 58.73007 68.37157
#> [691] 58.83184 60.74703 56.46877 71.00261 56.33933
#> [696] 60.24832 66.28596 60.18244 85.86029 94.76058
#> [701] 74.04342 67.95297 105.41617 55.43712 61.66119
#> [706] 70.65620 65.71093 129.49122 171.91433 60.59582
#> [711] 63.56468 70.34392 64.49171 66.61263 70.42831
#> [716] 72.32793 76.10193 502.02566 70.98718 60.72160
#> [721] 90.36813 260.99255 64.52742 64.25340 77.87597
#> [726] 65.17405 68.04309 68.04538 61.49828 58.33588
#> [731] 162.66684 71.12056 67.46346 59.10764 52.20220
#> [736] 80.70800 72.05071 68.82917 60.23462 60.77843
#> [741] 147.99618 495.36548 177.41477 57.48217 58.45182
#> [746] 56.77082 63.19748 115.60305 90.93194 103.96318
#> [751] 59.60639 60.95199 59.61951 70.72384 212.57571
#> [756] 141.53998 74.18329 60.62464 57.37070 139.83483
#> [761] 73.24228 59.73790 56.68325 60.18039 74.05484
#> [766] 61.60294 55.00323 58.53261 69.83376 66.95512
#> [771] 50.19238 94.56210 64.73982 63.30759 56.05224
#> [776] 56.21866 108.20868 64.84207 68.38372 78.10467
#> [781] 78.43881 73.27892 55.92307 56.15502 71.26666
#> [786] 71.45005 60.13642 59.45648 58.23555 64.36731
#> [791] 64.26804 693.69097 74.41381 54.82171 60.89342
#> [796] 71.67139 72.27635 54.16051 83.84377 56.75193
#> [801] 87.04627 57.84957 65.19856 58.96269 59.10073
#> [806] 62.78386 192.96761 59.95376 131.36959 136.84141
#> [811] 58.83548 89.70195 82.40877 307.72749 65.97041
#> [816] 73.53220 61.29662 56.27354 60.99397 61.67943
#> [821] 91.90204 95.73899 307.76108 54.25998 154.58149
#> [826] 66.30210 67.63340 542.98107 64.41201 126.81114
#> [831] 58.04815 59.51742 75.49524 97.50542 90.26412
#> [836] 126.56918 107.13740 76.73623 68.13029 63.66367
#> [841] 53.06166 75.98541 145.24184 68.35245 54.93406
#> [846] 61.75341 63.65842 57.77246 64.46247 61.23151
#> [851] 11759.23541 72.08983 54.16623 56.21207 168.45159
#> [856] 83.29418 62.77231 73.44044 60.47532 64.18168
#> [861] 55.31935 80.09056 108.54724 51.21057 60.79631
#> [866] 56.90275 59.82890 56.24699 87.29836 58.52266
#> [871] 62.69740 60.06647 66.43517 93.97310 83.88516
#> [876] 70.90782 57.94291 64.57490 58.65931 161.37447
#> [881] 58.25578 70.79826 60.08660 68.44888 62.42405
#> [886] 64.11954 62.10187 67.69399 66.88886 174.17454
#> [891] 57.55215 73.06869 94.71979 65.22120 53.69505
#> [896] 63.73962 74.72405 72.63910 98.67982 58.38586
#> [901] 73.72951 60.03864 58.96082 65.13716 122.07987
#> [906] 69.74198 56.90725 68.57703 57.07134 155.79712
#> [911] 60.12192 69.16012 59.34521 104.75812 110.81608
#> [916] 130.98377 59.37216 81.73337 58.04952 55.36329
#> [921] 146.39127 72.06932 55.08826 60.13069 69.49412
#> [926] 61.30376 68.35286 93.17964 68.32836 76.19179
#> [931] 56.64084 58.19276 99.96979 62.14001 63.14007
#> [936] 71.22079 67.34382 84.19186 117.12297 52.14092
#> [941] 59.05930 72.12002 66.36513 109.67208 64.83743
#> [946] 62.11450 53.26161 74.41590 78.49514 98.47208
#> [951] 49.05772 52.75193 174.19764 74.67304 142.85028
#> [956] 123.74091 101.73396 55.50526 67.41884 85.30606
#> [961] 139.75303 60.51737 104.51466 59.04707 62.03797
#> [966] 71.80082 90.22595 47.43704 80.98746 85.23247
#> [971] 61.47697 110.82960 64.15447 84.29105 214242.84037
#> [976] 73.07838 58.80663 82.61791 70.39338 86.59442
#> [981] 71.27020 62.91601 61.82828 95.83929 61.19795
#> [986] 67.73997 80.39113 59.95589 57.24045 79.67941
#> [991] 72.15522 677.94054 6643.39594 4420.54623 78.99679
#> [996] 70.83425 358.16904 100.41492 63.00446 57.61863