![[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] 5.512345e+04 3.512849e+03 3.952637e+02 6.627826e+01 6.951973e+01
#> [6] 2.764341e+02 1.692522e+02 3.736906e+02 9.880643e+01 1.568918e+02
#> [11] 2.894455e+04 1.340297e+03 7.857375e+01 9.168609e+01 1.016785e+02
#> [16] 8.755627e+01 9.370234e+01 7.452098e+01 7.400284e+01 9.744250e+01
#> [21] 7.715467e+01 7.715413e+01 7.684655e+01 1.933761e+03 8.505584e+06
#> [26] 1.124677e+02 7.389674e+01 8.776250e+01 9.063644e+01 7.056090e+01
#> [31] 1.215242e+02 7.949333e+01 3.085457e+03 6.445858e+01 9.009171e+01
#> [36] 1.348066e+02 1.212317e+02 1.177273e+02 6.507555e+02 3.675543e+02
#> [41] 7.648316e+01 7.747100e+01 6.834272e+01 8.421208e+01 5.850113e+02
#> [46] 8.172697e+01 1.396074e+02 6.769281e+01 6.478779e+01 7.172683e+01
#> [51] 1.647267e+02 1.090105e+02 6.795408e+01 7.436359e+01 6.747884e+01
#> [56] 7.552709e+01 1.141463e+02 1.674681e+02 1.164365e+02 1.293188e+02
#> [61] 4.310246e+02 9.920507e+01 7.671287e+01 7.180851e+01 6.673378e+01
#> [66] 6.690047e+01 4.460876e+02 7.997484e+01 6.800488e+01 6.518076e+01
#> [71] 6.267010e+01 7.022380e+01 7.172782e+01 1.529298e+02 1.263691e+02
#> [76] 7.918369e+01 2.060693e+02 1.387693e+02 6.863904e+01 1.900947e+02
#> [81] 6.054539e+01 6.381597e+01 7.550344e+01 7.140911e+01 1.311522e+02
#> [86] 1.395887e+02 7.814323e+01 7.644729e+01 5.648974e+02 2.452926e+02
#> [91] 1.240601e+02 1.737720e+02 7.152607e+01 9.611892e+01 6.646463e+01
#> [96] 1.043958e+02 6.814688e+01 2.699002e+02 6.961470e+01 8.219710e+01
#> [101] 7.768785e+01 1.571079e+02 1.075077e+02 9.824013e+01 1.037315e+02
#> [106] 3.263957e+02 6.397387e+01 9.464390e+01 7.057478e+01 6.730691e+01
#> [111] 7.963645e+01 8.154756e+01 6.316492e+01 7.710790e+01 1.791602e+02
#> [116] 1.647305e+02 7.285437e+01 7.153816e+01 7.360728e+01 7.613251e+01
#> [121] 1.348418e+03 9.403635e+01 9.852243e+01 1.444875e+02 3.329229e+02
#> [126] 8.868636e+01 7.231291e+01 7.850726e+01 9.870897e+01 7.717032e+01
#> [131] 1.605981e+02 3.128949e+02 6.492506e+01 7.374996e+01 7.721241e+01
#> [136] 2.083253e+02 1.386759e+02 8.771631e+01 6.937489e+01 7.120460e+01
#> [141] 6.489761e+01 6.376392e+01 7.520345e+01 8.201508e+01 7.091566e+01
#> [146] 6.474103e+01 2.299929e+02 2.258406e+03 5.361219e+04 5.428891e+07
#> [151] 1.629670e+02 1.063816e+02 6.578874e+01 1.723142e+02 1.417047e+02
#> [156] 6.289540e+01 4.677515e+02 3.537169e+02 4.878033e+02 1.095798e+02
#> [161] 1.159087e+02 6.948035e+01 8.085978e+01 6.578330e+01 9.256576e+01
#> [166] 1.318555e+02 1.079420e+02 1.096735e+02 8.113343e+01 2.294138e+02
#> [171] 1.043548e+02 8.310484e+01 6.729770e+01 1.856734e+02 7.395384e+01
#> [176] 1.746134e+02 8.920845e+01 7.733359e+01 4.922173e+02 1.450055e+02
#> [181] 7.316249e+01 7.393412e+01 8.860900e+01 7.534251e+01 1.763907e+02
#> [186] 9.050079e+01 7.547185e+01 1.147694e+02 1.077311e+02 6.682350e+01
#> [191] 3.771160e+02 1.982664e+03 9.738749e+01 1.610648e+02 6.905758e+01
#> [196] 2.647170e+02 8.418956e+01 8.953001e+01 8.755309e+01 7.912320e+01
#> [201] 8.744128e+01 6.553232e+01 1.059392e+02 7.486837e+01 1.304072e+02
#> [206] 7.350736e+01 1.543569e+02 1.041840e+02 1.988121e+02 2.217972e+04
#> [211] 1.244941e+02 7.103537e+01 1.177963e+02 7.842938e+01 8.058101e+01
#> [216] 6.842504e+01 1.713127e+02 6.908785e+01 7.235522e+01 7.275634e+01
#> [221] 7.347683e+01 8.574216e+01 1.266866e+02 8.914085e+01 5.995769e+01
#> [226] 6.794599e+01 2.813127e+02 9.080187e+01 1.848040e+02 9.722758e+01
#> [231] 1.437379e+02 7.694490e+01 2.217727e+02 8.601789e+01 1.095195e+02
#> [236] 1.912229e+02 1.011475e+02 7.442513e+01 7.044076e+01 6.677950e+01
#> [241] 8.869605e+01 7.035289e+01 1.225983e+02 6.933631e+01 1.004561e+02
#> [246] 7.769261e+01 9.532171e+01 6.781803e+01 1.540930e+02 3.779787e+02
#> [251] 1.404103e+02 7.188253e+01 9.645647e+01 1.003069e+02 6.282792e+01
#> [256] 1.195433e+02 7.258870e+01 7.443629e+01 1.303244e+02 6.223176e+01
#> [261] 2.242292e+02 1.127627e+02 4.080072e+02 2.921442e+02 1.172407e+02
#> [266] 6.692565e+01 8.190923e+01 6.743187e+01 7.102213e+01 1.173219e+02
#> [271] 7.228473e+01 1.108394e+02 8.522312e+01 7.207101e+01 6.698379e+01
#> [276] 7.458109e+01 3.621679e+02 6.851666e+01 2.657226e+02 9.230988e+01
#> [281] 6.843109e+01 6.470726e+01 8.232096e+01 7.625025e+01 9.203725e+01
#> [286] 7.018718e+01 7.146046e+01 6.806177e+01 7.564029e+01 6.861819e+01
#> [291] 8.211927e+01 1.030103e+02 7.809336e+01 1.082263e+02 1.695690e+02
#> [296] 1.035563e+02 2.109815e+02 8.589334e+01 6.316482e+01 7.930361e+01
#> [301] 6.990283e+01 8.310795e+01 2.157224e+02 6.457525e+01 1.781446e+02
#> [306] 1.017309e+02 1.710172e+02 1.362077e+02 3.193528e+02 7.946703e+01
#> [311] 7.974624e+01 6.687932e+01 7.535270e+01 8.445437e+01 7.189969e+01
#> [316] 1.730309e+02 8.241494e+01 1.578703e+02 1.566445e+02 7.492768e+01
#> [321] 6.618529e+01 1.022567e+02 8.232875e+01 7.351856e+01 6.124759e+01
#> [326] 6.538003e+01 8.339077e+01 6.921710e+04 8.296281e+01 8.062845e+01
#> [331] 7.301701e+01 7.759173e+01 7.348397e+01 7.478590e+01 8.093657e+01
#> [336] 7.381704e+01 9.671360e+01 6.937115e+01 1.120574e+02 9.947712e+01
#> [341] 8.322167e+01 7.971531e+01 1.145147e+03 1.045592e+02 2.089356e+02
#> [346] 8.098956e+01 2.872363e+03 1.566789e+02 1.191387e+02 8.449501e+01
#> [351] 1.335626e+02 1.025857e+03 6.600549e+01 6.871603e+01 2.345023e+02
#> [356] 3.034574e+02 1.540573e+02 1.791934e+03 1.067231e+02 1.235014e+02
#> [361] 1.306135e+02 9.635400e+01 9.286615e+01 6.779096e+01 6.427938e+01
#> [366] 9.343377e+01 2.652276e+02 9.827316e+02 6.862107e+01 6.596627e+01
#> [371] 8.637854e+01 1.366543e+02 9.224221e+01 1.047574e+02 6.038550e+02
#> [376] 1.337200e+03 3.861870e+02 4.266622e+02 7.138438e+01 8.087956e+01
#> [381] 1.858984e+02 2.754760e+02 1.959297e+02 1.076774e+02 6.878239e+01
#> [386] 8.069925e+01 1.269786e+02 9.337103e+01 1.612414e+02 8.557757e+01
#> [391] 9.804417e+01 7.648965e+01 1.071763e+02 9.663175e+01 1.509253e+02
#> [396] 7.638251e+01 7.872881e+01 7.141864e+01 8.420817e+01 6.580254e+01
#> [401] 1.463466e+02 1.463735e+02 7.035239e+01 1.119347e+02 9.244133e+01
#> [406] 7.601117e+01 2.176637e+02 9.637956e+01 8.938637e+02 1.863645e+03
#> [411] 8.129157e+02 7.672034e+01 9.653159e+01 9.506452e+01 7.986473e+01
#> [416] 2.189302e+05 6.804423e+01 1.031572e+02 8.358739e+01 1.587462e+02
#> [421] 7.855935e+01 7.889375e+01 1.074076e+02 1.021951e+02 2.167600e+02
#> [426] 7.650119e+01 1.107676e+02 7.890302e+01 8.081799e+01 4.387337e+02
#> [431] 1.389266e+02 6.648138e+01 6.502715e+01 3.060520e+02 1.037200e+02
#> [436] 1.096604e+02 8.521574e+01 2.246003e+02 6.670343e+01 6.427704e+01
#> [441] 6.489824e+01 6.570931e+01 6.247008e+01 6.733550e+01 7.120498e+01
#> [446] 6.901477e+01 7.846169e+01 1.261559e+02 8.585509e+01 1.535468e+02
#> [451] 8.507192e+01 1.099777e+02 1.584618e+02 8.139243e+01 2.375463e+02
#> [456] 8.654881e+01 6.502208e+01 1.134848e+02 1.953207e+02 1.039650e+02
#> [461] 7.965322e+01 1.105211e+02 1.151640e+02 1.190741e+02 1.011444e+02
#> [466] 1.011704e+02 8.772401e+01 6.071165e+01 8.284578e+01 1.063453e+02
#> [471] 6.598019e+01 1.192638e+02 1.330565e+02 4.246528e+02 2.115862e+02
#> [476] 1.822519e+02 9.091514e+01 1.138796e+02 7.107884e+01 8.545814e+01
#> [481] 7.335162e+01 1.269803e+02 7.022577e+01 1.372159e+02 6.586312e+02
#> [486] 7.716553e+01 6.779051e+01 9.900585e+01 8.101322e+03 1.144484e+02
#> [491] 2.797032e+03 1.861084e+02 2.051081e+02 6.779865e+01 8.101812e+01
#> [496] 6.996077e+01 1.151794e+02 6.078604e+01 7.236822e+01 4.608243e+03
#> [501] 7.172387e+01 8.125925e+01 1.074487e+02 8.117716e+01 7.245984e+01
#> [506] 1.053673e+02 1.048663e+02 1.170144e+02 9.850815e+01 1.008498e+02
#> [511] 1.411348e+03 9.429167e+01 1.512125e+02 8.198586e+01 9.009311e+01
#> [516] 1.027453e+02 1.020108e+02 7.862052e+01 9.079603e+01 7.408996e+01
#> [521] 8.395788e+01 1.443052e+04 1.836818e+04 1.849736e+02 1.140041e+02
#> [526] 1.734589e+02 2.331429e+02 8.209915e+01 6.530947e+01 8.105743e+01
#> [531] 1.003117e+02 8.503011e+01 1.537605e+02 8.184752e+01 2.892898e+02
#> [536] 3.567195e+02 8.549183e+01 1.358452e+02 7.689391e+01 6.783788e+01
#> [541] 8.525550e+01 8.265740e+01 1.347327e+02 3.993115e+03 1.919146e+02
#> [546] 1.494484e+02 1.291528e+02 1.027780e+02 1.458452e+02 8.177772e+01
#> [551] 2.133015e+02 6.176099e+01 2.849614e+02 3.026538e+02 6.716028e+01
#> [556] 1.022864e+02 8.333115e+01 1.129874e+02 6.484905e+01 7.303268e+01
#> [561] 7.112062e+01 6.883592e+01 1.369978e+02 2.851550e+02 7.327313e+01
#> [566] 1.409012e+02 8.600531e+01 6.498397e+01 8.929645e+01 8.445436e+01
#> [571] 7.616219e+01 1.375828e+02 7.815377e+01 2.612880e+02 7.330341e+02
#> [576] 1.034848e+02 1.019831e+02 8.677413e+01 1.012755e+02 7.794022e+01
#> [581] 9.486984e+01 2.021886e+02 1.176736e+02 1.032712e+02 7.860790e+01
#> [586] 7.719548e+01 1.453394e+02 1.411379e+02 1.405576e+02 1.476535e+02
#> [591] 8.279628e+01 8.061006e+01 1.195255e+02 8.463532e+01 1.188608e+02
#> [596] 6.917699e+01 7.166406e+01 9.913183e+01 1.108962e+02 7.082456e+02
#> [601] 5.801688e+01 6.893054e+01 1.476441e+02 1.473253e+02 7.309982e+01
#> [606] 6.150077e+01 9.380960e+01 6.585166e+02 2.508466e+02 4.312931e+02
#> [611] 1.863480e+02 1.630457e+02 1.061546e+02 1.110315e+02 8.163273e+01
#> [616] 9.888768e+01 1.020626e+02 8.182882e+01 1.241293e+02 8.902868e+01
#> [621] 3.631130e+03 2.405010e+02 1.671730e+03 3.470351e+02 1.222081e+02
#> [626] 7.040240e+01 1.227469e+02 7.019671e+01 7.328672e+01 1.308774e+03
#> [631] 5.041903e+05 9.364703e+04 9.153832e+01 1.034652e+02 2.296139e+02
#> [636] 8.677488e+02 8.148609e+01 6.665630e+01 6.378775e+01 9.974424e+01
#> [641] 1.606759e+02 3.124489e+03 1.036965e+02 1.391432e+02 6.746037e+01
#> [646] 1.537767e+02 1.017745e+02 1.089039e+02 1.562899e+02 8.661379e+01
#> [651] 8.495098e+01 1.070457e+02 1.208422e+02 8.047104e+01 1.548815e+02
#> [656] 8.057000e+01 9.163342e+01 2.019348e+04 7.855750e+02 6.511618e+01
#> [661] 8.054399e+01 6.750213e+01 6.549183e+01 7.093745e+01 1.456652e+02
#> [666] 1.174072e+02 1.627195e+02 1.785320e+02 7.067338e+01 8.268033e+01
#> [671] 7.233787e+02 1.197035e+02 4.567151e+02 7.272004e+01 8.389210e+01
#> [676] 1.948613e+02 1.392547e+02 1.141931e+02 1.291212e+02 9.944362e+01
#> [681] 7.165256e+01 6.594707e+01 6.879386e+01 7.852713e+01 6.439065e+01
#> [686] 8.723112e+01 5.618329e+02 8.114601e+06 7.767709e+01 7.457762e+01
#> [691] 9.049034e+01 5.666220e+02 7.120638e+01 8.876408e+01 6.704394e+01
#> [696] 7.813742e+01 8.995230e+01 8.952970e+01 9.043287e+01 8.931460e+01
#> [701] 6.244139e+01 9.414957e+01 1.495621e+02 1.721067e+02 6.396689e+01
#> [706] 4.795772e+02 1.176705e+02 9.551763e+01 2.700347e+02 4.886275e+02
#> [711] 1.008736e+02 6.931101e+01 7.380569e+01 7.399255e+01 9.103689e+01
#> [716] 1.221938e+02 7.319877e+01 7.774630e+01 8.649183e+02 1.641204e+02
#> [721] 7.221853e+01 1.002123e+02 5.823748e+03 1.307218e+02 1.081701e+02
#> [726] 5.899106e+02 1.076045e+03 1.029755e+02 6.899023e+01 7.069789e+01
#> [731] 8.694544e+01 6.549117e+01 9.288956e+01 2.402150e+02 1.438245e+02
#> [736] 6.741450e+01 2.155109e+05 1.060708e+03 2.863881e+02 7.256597e+01
#> [741] 1.014757e+02 1.862585e+02 6.628917e+01 1.831870e+02 7.098413e+01
#> [746] 7.565679e+01 6.536785e+01 7.858598e+01 7.192076e+02 6.444357e+01
#> [751] 1.134397e+02 5.167525e+02 2.254305e+02 8.500375e+02 7.144755e+01
#> [756] 6.726448e+01 1.046389e+02 8.676248e+01 7.201778e+01 2.200983e+03
#> [761] 6.566536e+01 1.155108e+02 8.674594e+01 8.007906e+01 9.802401e+01
#> [766] 7.375789e+01 2.806488e+03 1.016904e+03 6.190212e+01 8.010508e+01
#> [771] 9.476719e+01 6.435227e+01 6.789965e+01 6.659051e+01 7.547143e+01
#> [776] 8.043557e+01 8.907352e+01 1.082823e+02 4.242705e+02 1.624278e+02
#> [781] 1.398525e+02 2.770266e+02 7.408064e+01 2.413557e+02 3.078963e+02
#> [786] 2.294266e+02 9.677482e+01 1.173062e+02 1.063638e+02 7.792608e+01
#> [791] 9.268428e+01 1.075807e+02 1.096890e+02 1.054151e+02 7.826496e+01
#> [796] 4.218465e+02 8.079390e+01 8.436101e+01 1.088796e+02 1.894037e+02
#> [801] 1.080453e+02 8.436989e+01 7.502816e+01 1.078861e+02 7.324554e+01
#> [806] 1.485379e+02 1.136401e+02 1.039689e+02 1.403179e+02 1.197276e+02
#> [811] 7.373836e+01 6.959761e+01 7.263090e+01 1.444986e+02 3.788744e+02
#> [816] 9.643586e+01 8.265082e+01 3.642359e+02 6.822007e+01 8.053214e+01
#> [821] 3.930730e+02 1.600347e+03 5.424447e+04 2.327108e+03 8.021796e+01
#> [826] 1.552812e+02 9.053785e+01 2.048851e+02 9.759615e+01 2.292208e+02
#> [831] 1.965087e+02 6.560729e+01 6.892496e+01 6.383694e+01 7.849955e+01
#> [836] 8.221888e+01 9.311223e+01 6.818259e+01 6.294908e+01 6.937712e+01
#> [841] 7.512027e+01 2.438663e+02 1.103333e+02 7.248892e+01 8.340210e+01
#> [846] 7.542616e+01 1.407912e+02 4.683372e+02 7.112723e+01 7.072387e+01
#> [851] 8.195207e+01 9.480624e+01 9.134818e+01 7.569272e+01 6.362658e+01
#> [856] 6.545446e+01 7.221938e+01 6.698968e+01 7.354910e+01 7.476950e+01
#> [861] 9.415774e+01 6.993441e+01 2.454750e+02 7.521238e+01 8.687786e+01
#> [866] 1.518948e+02 1.177453e+02 9.048133e+01 8.571775e+01 7.237960e+01
#> [871] 6.477090e+01 6.523130e+01 7.449526e+01 1.003895e+02 1.020623e+02
#> [876] 3.805822e+02 2.134236e+03 1.146587e+02 1.278200e+02 7.325034e+01
#> [881] 8.742722e+01 9.677988e+01 6.351291e+01 6.316306e+01 1.130115e+02
#> [886] 1.018917e+02 2.620430e+02 8.708681e+01 3.645875e+02 2.612691e+02
#> [891] 6.151136e+01 7.268526e+01 6.210152e+01 7.110729e+01 1.268769e+02
#> [896] 3.345002e+02 8.757786e+01 8.429146e+01 1.585384e+02 8.007854e+01
#> [901] 6.894731e+01 9.211997e+01 7.408080e+01 8.889145e+01 7.786517e+01
#> [906] 9.538429e+01 6.309363e+01 6.865261e+01 8.973019e+01 1.704520e+02
#> [911] 6.287529e+01 1.101670e+02 7.753127e+01 6.891450e+01 1.977652e+02
#> [916] 8.632271e+01 7.719842e+01 9.109733e+01 7.106141e+01 7.264582e+01
#> [921] 5.896380e+02 6.291733e+02 6.311263e+01 8.488404e+01 6.653940e+01
#> [926] 8.805296e+01 6.913135e+01 6.330285e+01 7.393683e+01 9.735937e+01
#> [931] 7.301202e+02 7.906620e+01 7.168955e+01 8.745009e+01 6.475598e+01
#> [936] 6.802162e+01 3.077951e+02 6.668939e+01 7.876824e+01 6.004178e+01
#> [941] 8.315621e+01 8.816048e+01 9.560387e+01 6.335918e+01 1.051511e+02
#> [946] 1.087411e+02 8.034630e+01 8.804840e+01 7.124100e+01 1.766118e+02
#> [951] 1.577217e+02 7.332385e+01 1.320019e+02 1.238724e+02 1.143710e+02
#> [956] 1.403680e+02 1.016137e+02 7.310679e+01 6.752307e+01 1.013149e+02
#> [961] 8.814089e+01 2.951475e+02 8.080363e+01 1.001647e+02 7.543145e+01
#> [966] 1.720751e+02 9.793800e+01 6.967744e+01 7.835543e+01 7.352316e+01
#> [971] 1.315117e+02 1.405251e+03 8.695294e+02 9.026669e+01 8.947349e+01
#> [976] 8.081011e+01 9.680537e+01 6.898962e+01 1.356802e+02 7.336701e+01
#> [981] 6.724626e+01 7.059191e+01 9.776969e+01 7.261770e+01 1.991323e+02
#> [986] 8.759126e+01 7.181677e+01 2.722298e+02 7.628483e+01 8.327146e+01
#> [991] 9.745430e+02 7.419334e+01 6.987643e+01 7.020212e+01 1.487281e+02
#> [996] 6.141143e+01 7.278439e+01 1.199585e+02 1.480431e+02 7.088834e+01
doseFunction(
model,
alpha2 = samples@data$alpha2,
beta = samples@data$beta,
grade = 2L
)(x = 0.25)
#> [1] 1896.72755 155.28485 60.00073 55.36907 61.72185
#> [6] 164.65382 119.34676 167.13303 76.18897 92.87348
#> [11] 946.35797 173.82467 65.92123 68.52815 73.25081
#> [16] 72.87386 75.53381 59.33714 64.91135 55.74278
#> [21] 63.32251 61.31299 66.73278 771.76666 274988.72900
#> [26] 66.60272 61.42112 63.86143 71.07511 61.39649
#> [31] 84.69537 66.28474 293.71735 55.45661 64.65798
#> [36] 90.93827 70.36765 74.21012 275.53494 125.70573
#> [41] 64.64954 59.49895 60.90803 63.40653 171.02957
#> [46] 61.43471 61.81572 59.77033 59.03011 62.05560
#> [51] 115.14450 85.07439 59.43212 63.12599 59.44026
#> [56] 64.75186 87.34891 63.16875 73.46534 77.03802
#> [61] 205.28365 74.92257 66.41532 58.64489 57.53235
#> [66] 56.38190 214.36393 66.76321 59.77451 58.09040
#> [71] 56.82173 62.36850 62.21494 77.31329 83.48457
#> [76] 61.51542 104.36643 69.73556 57.61479 123.98847
#> [81] 56.00765 56.59733 64.81286 59.90466 86.94533
#> [86] 82.21153 61.69201 63.29951 177.58431 69.20755
#> [91] 69.39350 103.42175 54.15506 61.77005 58.62849
#> [96] 68.76402 61.03207 153.64869 60.14996 68.44474
#> [101] 65.31933 101.83776 58.95800 70.93793 60.07360
#> [106] 153.70394 56.85167 65.73720 59.55723 56.31335
#> [111] 68.22073 68.96759 56.95467 63.58715 107.97975
#> [116] 53.78130 62.39766 59.84566 60.95623 61.51284
#> [121] 344.13252 75.26440 62.50586 73.60487 162.57160
#> [126] 71.37435 59.28019 63.83522 72.92608 64.62596
#> [131] 104.57663 131.26862 57.84853 64.79544 66.13305
#> [136] 85.85919 58.14645 56.01313 59.36877 59.32337
#> [141] 57.03114 57.36438 65.76410 66.35167 61.86400
#> [146] 56.66514 112.57333 340.06227 2268.61172 264.96848
#> [151] 90.88152 80.18448 54.85169 93.10808 78.71280
#> [156] 54.07147 108.13343 94.56385 110.56286 61.96707
#> [161] 78.12453 59.75553 60.38977 58.09862 71.72951
#> [166] 65.99612 75.03442 76.21432 57.73157 98.16090
#> [171] 76.95749 42.82498 60.55462 99.93115 60.92040
#> [176] 92.63742 70.44196 58.34613 184.98338 70.94277
#> [181] 61.35841 62.27272 69.65224 60.79893 105.05487
#> [186] 54.92452 59.98741 73.15338 76.22429 54.63272
#> [191] 116.29439 474.86028 64.06346 104.21411 58.93552
#> [196] 114.72789 64.62475 60.85270 58.29988 57.20105
#> [201] 55.36209 57.22424 78.80773 55.99727 77.27439
#> [206] 39.32999 93.15557 68.99880 90.14985 616.11503
#> [211] 78.13457 60.91073 70.39802 62.37555 66.47625
#> [216] 59.24640 111.44685 60.37449 61.88809 61.48824
#> [221] 60.53276 66.45920 76.16609 53.87231 55.67567
#> [226] 63.25364 188.03786 77.75678 115.73201 56.41356
#> [231] 94.06962 56.21644 71.92439 65.15069 83.98674
#> [236] 106.71383 65.86404 63.82462 60.96036 59.88048
#> [241] 68.40506 57.97472 67.95107 58.46591 61.82197
#> [246] 61.63371 71.69116 58.89745 68.27540 143.83182
#> [251] 98.03561 58.48221 64.74782 63.17530 56.25008
#> [256] 70.40352 60.01913 63.63940 77.49081 53.41824
#> [261] 148.36243 51.35845 132.24567 74.64336 85.26913
#> [266] 56.09375 63.04965 60.88775 59.10176 81.41130
#> [271] 58.63125 77.13181 62.29382 58.04523 54.19144
#> [276] 58.06386 107.29267 59.78378 129.01254 63.48924
#> [281] 59.19662 59.08204 70.69799 64.10386 72.27985
#> [286] 58.03345 58.54243 58.27913 58.84793 56.41198
#> [291] 64.90904 66.68757 55.27337 80.97400 85.03516
#> [296] 79.04792 91.40785 66.14272 56.85715 61.32114
#> [301] 54.97623 57.35663 92.91361 52.62373 107.22146
#> [306] 66.83823 75.02942 83.83680 123.03387 63.16307
#> [311] 63.80902 55.20783 59.96188 68.97293 58.36657
#> [316] 77.07202 63.53676 98.52587 60.01137 52.80600
#> [321] 57.57612 70.95160 66.57684 59.68201 55.65837
#> [326] 60.60246 73.59763 13388.07170 72.49136 66.91710
#> [331] 59.17904 62.09870 60.42231 62.94978 66.55285
#> [336] 60.67304 71.81678 53.72081 57.52651 69.24063
#> [341] 51.34343 51.69484 338.72848 59.20829 82.13245
#> [346] 60.03621 335.00326 85.99633 81.48387 53.93142
#> [351] 75.13255 371.05617 58.99549 60.48245 136.08262
#> [356] 88.98494 86.91121 347.47291 73.95949 89.07584
#> [361] 75.81110 72.07498 75.35552 61.12529 59.06691
#> [366] 73.69431 167.89896 192.20224 58.54286 60.88683
#> [371] 69.39087 91.81385 54.16204 72.56961 135.26595
#> [376] 166.56010 112.61842 184.97945 61.35695 59.89901
#> [381] 109.23593 75.35014 69.91729 70.66405 59.74942
#> [386] 68.21808 87.67034 56.59288 59.76772 55.68104
#> [391] 78.11833 58.99822 71.83613 71.33403 83.85986
#> [396] 57.75055 63.35302 55.64132 67.96814 56.51117
#> [401] 76.68693 86.97976 58.36077 66.69940 65.85841
#> [406] 61.37297 91.36331 75.16142 312.87837 376.94905
#> [411] 98.59141 62.61347 61.01696 54.39212 56.49729
#> [416] 5007.16470 57.87672 75.90644 60.29212 104.58079
#> [421] 63.43002 63.07861 75.35284 80.81398 55.56892
#> [426] 63.37197 79.35910 56.55481 62.84521 169.01874
#> [431] 99.50291 60.52836 57.65760 126.90440 58.07663
#> [436] 67.98298 63.12877 125.38671 54.57874 57.16751
#> [441] 55.69415 59.42719 57.28302 61.87055 63.49061
#> [446] 61.18913 63.38967 73.04837 63.84889 91.51820
#> [451] 55.97337 65.10199 75.33921 63.65785 86.20454
#> [456] 50.60834 54.39879 81.24128 95.09616 76.04900
#> [461] 65.12345 76.12011 67.50093 73.26619 71.77971
#> [466] 73.52571 66.79985 55.58082 70.24303 82.91930
#> [471] 54.99330 79.40031 84.49571 103.93273 102.24065
#> [476] 105.01940 73.00448 86.48943 59.52781 68.86481
#> [481] 60.83577 62.69608 59.77935 81.46686 113.23837
#> [486] 62.92814 56.39462 63.16420 87.01306 70.86558
#> [491] 353.14800 108.69461 93.25400 55.90868 61.10670
#> [496] 58.21379 71.38906 55.28086 60.44225 919.52604
#> [501] 59.05836 49.39028 83.10176 68.98011 59.98615
#> [506] 57.80082 73.72874 67.37413 72.24800 73.09134
#> [511] 474.38579 74.14741 95.80088 63.47105 71.99568
#> [516] 77.64883 63.00023 60.04986 57.70114 62.69179
#> [521] 65.75189 485.02645 481.23213 98.27131 70.27263
#> [526] 95.69735 72.62709 56.30804 58.79789 66.50598
#> [531] 71.91945 64.11476 90.18014 66.88676 124.64024
#> [536] 152.02950 62.71436 87.23756 64.65859 58.80338
#> [541] 59.66451 59.75515 52.07186 1012.64705 55.10091
#> [546] 42.84118 73.56905 61.40813 95.73782 62.54464
#> [551] 116.33517 55.45540 97.85859 124.69403 55.67522
#> [556] 51.22574 66.10902 64.31375 59.08389 60.40214
#> [561] 59.14953 63.22506 105.11868 167.06206 59.69747
#> [566] 85.93998 66.01190 58.88430 72.93693 68.41326
#> [571] 53.51454 75.00747 60.00049 122.01619 236.11450
#> [576] 70.42472 76.99765 65.47583 77.22742 62.01347
#> [581] 70.29751 122.01516 85.17818 72.70300 61.16683
#> [586] 50.66832 87.09813 81.61315 80.28504 88.64913
#> [591] 55.32055 63.12622 90.10328 53.45472 65.14984
#> [596] 60.04376 62.66495 72.99764 71.65011 216.11983
#> [601] 50.73347 57.20711 82.87916 90.44845 58.50395
#> [606] 56.16454 68.23962 118.75091 89.54606 154.23287
#> [611] 98.03793 61.81570 64.86340 71.23884 59.16679
#> [616] 67.37279 61.62917 58.63368 71.72972 62.25224
#> [621] 557.46060 58.76888 88.12448 125.35209 85.57482
#> [626] 60.22509 85.09772 59.63206 65.14056 306.29608
#> [631] 3988.81295 577.86455 52.82403 69.98743 72.81757
#> [636] 201.91611 62.42747 59.41085 56.10150 60.21844
#> [641] 94.39140 345.04964 61.17740 86.17795 59.13871
#> [646] 87.53912 69.61812 63.43252 100.26109 68.49852
#> [651] 63.45818 66.11837 75.88846 52.46402 78.84945
#> [656] 60.50476 71.00419 2494.22002 296.81745 55.67839
#> [661] 66.07282 57.43038 56.16843 58.39215 81.18187
#> [666] 71.76019 93.42870 96.45137 58.25021 66.64881
#> [671] 146.04387 75.96197 151.72021 58.21764 66.51093
#> [676] 135.27869 99.25639 60.73242 68.38959 68.63745
#> [681] 59.59093 56.88182 61.85973 63.74032 57.99620
#> [686] 70.62388 176.94504 22449.56810 51.07408 57.18007
#> [691] 66.40702 240.60095 54.67973 67.77032 56.89974
#> [696] 63.30936 73.11641 67.63591 69.75228 65.29869
#> [701] 57.45063 73.08164 82.22096 81.13634 52.84967
#> [706] 241.21513 71.19815 73.68788 134.02667 271.51944
#> [711] 77.92650 61.34407 65.43816 65.40342 71.79216
#> [716] 83.79447 58.57767 65.03331 256.79367 104.21391
#> [721] 52.16680 74.13885 401.47239 76.22496 70.55441
#> [726] 176.02542 215.75821 78.00490 59.99312 59.77506
#> [731] 68.04645 58.35266 72.16512 124.31168 95.26785
#> [736] 56.79833 6138.49970 106.82764 101.87062 58.02995
#> [741] 61.45933 89.81436 56.97606 88.99937 59.84762
#> [746] 63.14530 57.08105 64.82358 151.87748 55.73025
#> [751] 72.40609 160.49352 114.18566 253.81432 56.50540
#> [756] 56.13423 78.04407 70.32156 60.02719 466.79358
#> [761] 59.64680 88.13822 69.27628 67.39204 60.17812
#> [766] 59.61549 407.28545 198.76069 56.83236 62.44532
#> [771] 65.80643 57.72985 59.84919 57.47396 59.77187
#> [776] 65.50890 61.61136 74.92566 221.80272 103.42545
#> [781] 74.57890 77.48976 56.57439 96.67975 129.39940
#> [786] 93.73342 64.13690 77.18328 74.76123 61.37145
#> [791] 69.69624 70.88017 68.43703 52.87480 61.10768
#> [796] 101.79501 54.99109 60.69504 80.43853 103.17510
#> [801] 65.04215 63.83652 62.59251 79.77167 63.63638
#> [806] 83.58644 65.86579 64.11761 84.96374 67.46139
#> [811] 60.09570 58.95017 60.64176 92.34200 132.50036
#> [816] 65.46325 67.54240 74.92019 55.23259 62.40157
#> [821] 102.55495 424.64485 1755.19817 418.00817 62.68065
#> [826] 57.65911 62.42592 106.58802 71.31732 98.17852
#> [831] 79.37806 58.36678 60.06225 59.64552 56.82109
#> [836] 63.97346 67.20826 61.80624 57.85310 63.34913
#> [841] 63.40264 105.12102 74.96255 58.16360 65.62397
#> [846] 60.81876 73.79168 130.43921 60.56647 58.72711
#> [851] 68.99248 71.05657 71.77903 61.77005 56.50951
#> [856] 59.73638 58.71988 53.15657 61.35347 59.56142
#> [861] 70.33570 57.72637 123.72550 60.07877 64.51881
#> [866] 75.09352 77.30385 68.15095 66.48312 62.84837
#> [871] 57.62085 57.84568 61.25417 73.31391 59.51875
#> [876] 141.53266 447.44645 34.16494 97.67266 64.82945
#> [881] 72.96112 62.68435 49.13980 50.31106 76.78335
#> [886] 68.22227 127.40836 61.71966 106.52428 101.98363
#> [891] 53.73109 60.47766 57.38188 60.12693 45.73199
#> [896] 154.30591 63.56421 64.06589 89.14230 65.42448
#> [901] 61.69080 74.29688 59.36245 62.44583 58.89711
#> [906] 63.97176 55.02885 58.34407 65.12192 83.66481
#> [911] 56.10509 81.32820 66.60610 62.43406 92.34754
#> [916] 65.53390 64.00037 60.42932 61.36218 58.31374
#> [921] 179.29069 169.15449 55.77313 67.42859 57.67916
#> [926] 72.83524 59.39749 56.65427 63.12807 57.99314
#> [931] 225.04144 59.92336 59.45851 67.96346 58.87652
#> [936] 62.00431 172.67829 60.77258 62.76196 53.29119
#> [941] 63.18726 62.61606 75.21155 57.50503 72.36847
#> [946] 66.16597 65.86367 68.64277 55.55584 93.25625
#> [951] 88.21595 61.26636 81.72861 80.09875 68.31337
#> [956] 74.72576 60.47282 58.54651 59.45552 78.50831
#> [961] 66.89630 129.29980 52.61949 74.28963 64.40480
#> [966] 101.66832 61.75401 57.92528 67.29697 64.58432
#> [971] 62.43225 144.15365 115.00451 70.80825 68.00358
#> [976] 64.96202 72.21824 59.43283 81.51669 56.48314
#> [981] 56.60173 59.72828 71.98865 57.34299 101.83166
#> [986] 69.98901 60.69317 114.80969 63.95260 63.71609
#> [991] 260.59410 58.47553 58.82786 58.50512 84.09734
#> [996] 55.35348 59.99475 61.44179 87.02367 60.13394