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] 6.503084e+01 8.498407e+01 1.722265e+05 1.084608e+02 9.078398e+01
#> [6] 1.181502e+02 7.027493e+01 7.434530e+01 1.324969e+02 6.394657e+01
#> [11] 1.124658e+02 1.091526e+02 7.120897e+01 6.327834e+01 6.402245e+01
#> [16] 7.769486e+01 6.891442e+01 6.453608e+01 7.005443e+01 6.153396e+01
#> [21] 8.923055e+01 6.450479e+01 7.760224e+01 6.985643e+01 8.563263e+01
#> [26] 6.384704e+01 6.466711e+01 6.653645e+01 8.308099e+01 7.840831e+01
#> [31] 8.379316e+01 1.238230e+02 6.395280e+01 7.695206e+01 5.856745e+01
#> [36] 6.771935e+01 8.929541e+01 7.751109e+01 7.365970e+01 1.578591e+02
#> [41] 1.031466e+02 7.131603e+01 8.747489e+01 2.472047e+02 1.084568e+02
#> [46] 6.090196e+01 8.630686e+01 7.371421e+01 7.697519e+01 1.165463e+02
#> [51] 8.345645e+01 5.639928e+02 2.101983e+02 9.786847e+01 8.185470e+01
#> [56] 7.669657e+01 1.791204e+02 2.796758e+02 6.848033e+01 6.687694e+02
#> [61] 1.917765e+03 7.379869e+01 9.362901e+01 7.127939e+01 1.003004e+02
#> [66] 9.207894e+01 1.053104e+02 1.137647e+02 2.473878e+02 6.325475e+01
#> [71] 8.353038e+01 5.700847e+02 1.056529e+02 2.220344e+02 7.537199e+01
#> [76] 6.536868e+01 8.322058e+01 1.316995e+02 6.258951e+01 7.727648e+01
#> [81] 5.580148e+02 1.222888e+02 7.684023e+01 6.940840e+01 1.829506e+02
#> [86] 9.330323e+01 1.017972e+02 1.009770e+02 1.184292e+02 1.535989e+02
#> [91] 9.198828e+01 7.276377e+01 1.093656e+02 6.975872e+01 8.332459e+01
#> [96] 9.252075e+01 7.698194e+01 1.291220e+02 6.209793e+01 1.150519e+04
#> [101] 6.729705e+01 1.327842e+02 9.013334e+01 5.972891e+02 7.441872e+01
#> [106] 8.070157e+01 8.907183e+01 2.618728e+02 4.245492e+02 1.149803e+03
#> [111] 1.848492e+02 6.718782e+01 6.955089e+01 8.917288e+01 1.532461e+02
#> [116] 8.127902e+01 1.556842e+02 6.747357e+01 8.603749e+02 7.582119e+01
#> [121] 7.830263e+01 7.855383e+01 8.011230e+01 1.113070e+02 1.057506e+02
#> [126] 1.427723e+02 7.027705e+01 1.285299e+02 1.557013e+02 1.104287e+02
#> [131] 6.629406e+01 7.818412e+01 7.646418e+01 2.477438e+02 2.919791e+02
#> [136] 3.789521e+02 3.758356e+02 7.422668e+01 1.030223e+02 1.231875e+02
#> [141] 1.136044e+02 5.720861e+02 6.529358e+01 1.152568e+02 1.436768e+02
#> [146] 1.109442e+02 7.159707e+01 1.027552e+02 7.826992e+01 8.248965e+01
#> [151] 1.136497e+02 1.415924e+02 3.731579e+02 2.335800e+02 2.285431e+02
#> [156] 1.272584e+02 7.363271e+01 6.942581e+01 3.482417e+02 1.532114e+02
#> [161] 1.271255e+02 6.524294e+01 8.122912e+01 8.158074e+01 1.222241e+02
#> [166] 1.599526e+02 2.856903e+02 7.599746e+01 1.332819e+02 8.176940e+01
#> [171] 9.259351e+01 7.028315e+01 9.820258e+01 8.488094e+01 6.883887e+01
#> [176] 9.495413e+01 2.122723e+02 1.207067e+03 7.432917e+01 1.060834e+02
#> [181] 7.410315e+01 9.485360e+01 2.862606e+02 1.277007e+02 6.184486e+01
#> [186] 2.220644e+02 6.278293e+01 3.106339e+02 3.008494e+03 6.720637e+01
#> [191] 1.999328e+02 1.885044e+02 9.685263e+01 7.010346e+01 6.258795e+03
#> [196] 7.459464e+01 6.245622e+01 8.148639e+01 6.918182e+01 8.059190e+01
#> [201] 6.728309e+01 7.093447e+01 6.739035e+01 1.347430e+02 9.137718e+01
#> [206] 8.057651e+01 7.666685e+01 8.992839e+01 8.200693e+01 9.541110e+01
#> [211] 5.816691e+03 1.228327e+02 6.181771e+01 6.586417e+01 1.161272e+02
#> [216] 7.909117e+01 8.846124e+01 9.592866e+01 2.188370e+02 2.346454e+02
#> [221] 7.509577e+01 6.835090e+01 7.621722e+01 9.874818e+01 7.094198e+01
#> [226] 9.709109e+01 1.258394e+02 1.019503e+02 7.501154e+01 1.556508e+02
#> [231] 8.438213e+01 2.819699e+02 1.418858e+02 7.023149e+01 1.952700e+02
#> [236] 1.191799e+02 2.219140e+02 1.232989e+02 7.852397e+01 7.187521e+01
#> [241] 7.363604e+01 7.369231e+01 2.086001e+02 6.744804e+02 1.400254e+02
#> [246] 1.063645e+02 8.228697e+01 1.139842e+02 9.281450e+01 1.366660e+02
#> [251] 2.320574e+02 2.307426e+02 1.089490e+02 9.650974e+01 7.888544e+01
#> [256] 1.011646e+02 2.109782e+02 1.181977e+02 8.340524e+01 6.861892e+01
#> [261] 6.539948e+01 1.338659e+02 6.326899e+01 2.574105e+02 8.496280e+01
#> [266] 6.952776e+01 1.156693e+02 1.435487e+02 6.809930e+01 1.039647e+02
#> [271] 1.134985e+02 1.028105e+02 2.173632e+02 8.080296e+01 7.245273e+01
#> [276] 8.069801e+01 8.118066e+01 7.216292e+01 9.182271e+01 1.381751e+02
#> [281] 1.149461e+02 8.718183e+01 1.804098e+02 6.727142e+01 1.048407e+02
#> [286] 7.171278e+01 6.947606e+01 1.045978e+02 1.239224e+02 5.015247e+02
#> [291] 6.855526e+01 1.352003e+02 1.021910e+02 1.375786e+02 8.951363e+01
#> [296] 6.640493e+01 3.112651e+02 4.127177e+03 3.041731e+02 7.035787e+01
#> [301] 1.015686e+02 7.364512e+01 6.796977e+01 5.943621e+01 8.923011e+01
#> [306] 1.577525e+02 1.629621e+02 8.218035e+01 9.135446e+01 7.986531e+01
#> [311] 9.670040e+01 4.853841e+02 8.192245e+01 6.720938e+01 1.068993e+02
#> [316] 2.597964e+02 7.510108e+01 7.207820e+01 1.980493e+02 1.286661e+02
#> [321] 2.136060e+02 6.841335e+01 1.177406e+02 1.791463e+02 9.074243e+01
#> [326] 2.702403e+02 9.715915e+01 8.290010e+01 1.259079e+02 7.857410e+01
#> [331] 6.045756e+01 7.301783e+01 6.358667e+01 8.649604e+01 6.980643e+01
#> [336] 6.388502e+01 3.132787e+02 6.698287e+01 1.966805e+02 6.854826e+01
#> [341] 9.058445e+01 6.831299e+01 1.327763e+02 7.936325e+01 2.039416e+02
#> [346] 1.247671e+02 6.931938e+01 4.062341e+03 8.208494e+01 8.075341e+01
#> [351] 1.516907e+02 9.224136e+01 1.251799e+02 6.336984e+01 1.611782e+02
#> [356] 7.360378e+01 1.576472e+02 1.112493e+02 4.703930e+03 3.626039e+02
#> [361] 9.220015e+01 8.836313e+01 8.002635e+01 1.367394e+02 3.389248e+02
#> [366] 1.308288e+02 7.963813e+01 9.076880e+01 1.744151e+02 6.520522e+01
#> [371] 6.929842e+01 1.457986e+02 6.658277e+01 2.697688e+02 1.663014e+02
#> [376] 7.842516e+01 9.521462e+01 7.574806e+01 2.253561e+02 3.693896e+03
#> [381] 1.020538e+02 1.230892e+02 1.490889e+02 1.063988e+02 1.374652e+02
#> [386] 5.152618e+02 5.754561e+02 5.066471e+02 8.476937e+01 7.325551e+01
#> [391] 2.257513e+03 6.965076e+01 8.062518e+01 8.458528e+01 1.338855e+02
#> [396] 7.285312e+01 2.082680e+02 9.952664e+01 6.566586e+01 7.791174e+01
#> [401] 7.779504e+01 6.864016e+01 7.248550e+01 6.423324e+01 7.424535e+01
#> [406] 1.603833e+02 1.062326e+02 1.072339e+02 6.616223e+01 7.243314e+01
#> [411] 6.690499e+01 8.163824e+01 7.277554e+01 8.352850e+01 1.179314e+02
#> [416] 6.805935e+01 8.401525e+01 1.101491e+02 2.335976e+02 7.657881e+01
#> [421] 7.606132e+01 6.527077e+01 6.875819e+01 1.021488e+02 6.264788e+01
#> [426] 1.038294e+02 1.433259e+02 6.639518e+01 2.413614e+02 1.261039e+02
#> [431] 6.956612e+01 7.645859e+01 6.974067e+01 6.546398e+01 8.177875e+01
#> [436] 1.047113e+02 6.124354e+01 6.790119e+01 6.616305e+01 7.224999e+01
#> [441] 1.157478e+02 1.415222e+02 1.811120e+02 6.160772e+03 1.752306e+02
#> [446] 3.917627e+02 1.260762e+02 7.492817e+01 9.686728e+01 1.179086e+02
#> [451] 7.978477e+01 7.517379e+01 7.211045e+01 1.784269e+02 1.611152e+02
#> [456] 2.556675e+02 1.062776e+02 7.380160e+01 6.896856e+01 7.064532e+01
#> [461] 6.270236e+01 7.097065e+01 6.343107e+01 9.220503e+01 1.650248e+02
#> [466] 4.627322e+02 1.111316e+02 8.922099e+01 6.821423e+01 8.164998e+01
#> [471] 6.391668e+01 6.365634e+01 8.325735e+01 7.684420e+01 9.749441e+01
#> [476] 1.390125e+02 1.239318e+02 5.634085e+02 6.152209e+01 9.083780e+01
#> [481] 6.802851e+01 7.634628e+01 8.627444e+01 7.977474e+01 8.794724e+01
#> [486] 6.526277e+01 7.312176e+01 6.939773e+01 7.674708e+01 7.291890e+01
#> [491] 6.434157e+01 6.454216e+01 7.392111e+01 2.493895e+02 6.819520e+01
#> [496] 8.919941e+01 1.903033e+02 8.139815e+01 7.570364e+01 1.042115e+02
#> [501] 1.372302e+02 9.100375e+01 7.689756e+01 8.163560e+01 7.762703e+01
#> [506] 6.085078e+01 7.070630e+01 6.177928e+01 6.235916e+01 7.049648e+01
#> [511] 7.107035e+01 8.209730e+01 1.268741e+02 9.644858e+01 1.073805e+02
#> [516] 8.612022e+01 9.308050e+01 1.518829e+02 9.127490e+01 7.056902e+01
#> [521] 8.097765e+01 6.473212e+01 6.442351e+01 7.042769e+01 1.689398e+02
#> [526] 8.189825e+01 1.031744e+02 1.815059e+02 9.062112e+01 6.669275e+01
#> [531] 8.708299e+01 6.476419e+01 7.790505e+01 8.322525e+01 1.976810e+02
#> [536] 9.557449e+01 7.557802e+01 8.054657e+01 1.021680e+02 1.032896e+02
#> [541] 8.950749e+01 9.795868e+01 6.694175e+01 1.458403e+02 1.139801e+02
#> [546] 5.693004e+04 5.022792e+02 8.388899e+01 1.781849e+02 8.330739e+01
#> [551] 1.085215e+02 9.429007e+01 7.018223e+01 1.291206e+02 6.226269e+01
#> [556] 7.755749e+01 1.437325e+02 7.281637e+01 6.699461e+01 6.988115e+01
#> [561] 1.528451e+02 7.107760e+01 9.871291e+01 4.452729e+02 6.532630e+01
#> [566] 6.419435e+01 4.342461e+03 8.373562e+01 6.807697e+01 2.172637e+02
#> [571] 3.746384e+02 1.161669e+02 8.790881e+01 6.777657e+01 6.517945e+01
#> [576] 7.682197e+01 9.404811e+01 1.181210e+02 1.126046e+02 1.067075e+02
#> [581] 1.148417e+02 8.749214e+01 9.033284e+01 6.845147e+01 9.704365e+01
#> [586] 7.858482e+01 9.497891e+01 7.916923e+01 1.308223e+02 9.974410e+01
#> [591] 7.656605e+01 8.610852e+01 1.247100e+02 6.745917e+01 1.639858e+02
#> [596] 7.006606e+01 8.470136e+01 8.361961e+01 1.193119e+02 4.923638e+02
#> [601] 7.217552e+01 6.662284e+01 9.513087e+01 9.631167e+01 7.387117e+01
#> [606] 4.173352e+02 9.479093e+01 1.100325e+02 7.291247e+01 6.589898e+01
#> [611] 2.615167e+04 9.960697e+01 6.984198e+01 6.914031e+01 9.320362e+01
#> [616] 7.256220e+01 7.809848e+01 9.959519e+01 3.622340e+02 8.230438e+01
#> [621] 1.642883e+02 8.762355e+01 1.107786e+02 7.106525e+01 1.104771e+02
#> [626] 7.630980e+01 7.194540e+01 7.399865e+01 8.127346e+01 6.320093e+01
#> [631] 1.497209e+02 6.258932e+01 7.940668e+01 1.022388e+02 2.286406e+02
#> [636] 1.033732e+02 2.481919e+02 3.264940e+02 2.347604e+03 7.284574e+01
#> [641] 3.242767e+02 8.211219e+01 2.429147e+02 7.488144e+01 1.025128e+02
#> [646] 6.871923e+01 9.680284e+01 8.008450e+01 1.146113e+02 8.699601e+01
#> [651] 2.137526e+02 2.440099e+02 7.529815e+01 7.097522e+01 7.353278e+01
#> [656] 7.533854e+01 7.203241e+01 6.779574e+01 7.730931e+01 7.127327e+01
#> [661] 9.319975e+01 8.695814e+01 2.291523e+02 8.296521e+01 7.513481e+01
#> [666] 7.614569e+01 8.207810e+01 7.122129e+01 7.820029e+01 7.032229e+01
#> [671] 6.639852e+01 6.953890e+01 6.190981e+01 9.036169e+01 6.142469e+01
#> [676] 1.449840e+02 9.800269e+01 7.122956e+01 1.098650e+02 2.545892e+02
#> [681] 7.618136e+01 3.679906e+02 6.299330e+01 5.085663e+02 1.069025e+03
#> [686] 3.847088e+02 1.362140e+02 1.884074e+02 8.899823e+01 7.013242e+01
#> [691] 7.645829e+01 6.967685e+01 6.502035e+01 6.837303e+01 7.270277e+01
#> [696] 1.603585e+02 4.558941e+02 1.231059e+02 2.622426e+02 6.818516e+01
#> [701] 7.356042e+01 2.960673e+02 1.266247e+02 8.165447e+01 7.789208e+01
#> [706] 7.978521e+01 8.645060e+01 9.217071e+01 6.830740e+01 7.019757e+01
#> [711] 8.311072e+01 1.355279e+02 6.405546e+01 6.552920e+01 6.497825e+01
#> [716] 1.000158e+02 1.043084e+02 1.461556e+03 6.311240e+01 1.383562e+02
#> [721] 1.415537e+02 1.042629e+02 8.413007e+01 9.601770e+01 6.108690e+01
#> [726] 1.896029e+02 1.080678e+02 7.539670e+01 1.148878e+02 6.244713e+01
#> [731] 6.541991e+01 6.357502e+01 7.176918e+01 1.857017e+02 2.408103e+02
#> [736] 7.888159e+01 6.730535e+01 6.871393e+01 6.755808e+01 9.540366e+01
#> [741] 6.644645e+01 7.545907e+01 1.729810e+02 6.873647e+01 7.037858e+01
#> [746] 7.324532e+02 5.382044e+03 5.398474e+02 7.796392e+01 6.479670e+01
#> [751] 6.670543e+01 1.660615e+02 1.002886e+02 1.712647e+02 7.962938e+01
#> [756] 1.397672e+04 2.498252e+02 1.158853e+02 6.994027e+01 5.278646e+02
#> [761] 2.647985e+02 7.701714e+01 8.533231e+01 1.093354e+02 8.072386e+01
#> [766] 8.976489e+01 2.477205e+02 8.787242e+01 9.298816e+01 1.201584e+02
#> [771] 8.330275e+01 9.374226e+01 2.487042e+02 7.353362e+01 6.603288e+01
#> [776] 7.319286e+01 7.118204e+01 6.401546e+01 5.982924e+02 3.180768e+03
#> [781] 8.138040e+01 6.466809e+01 7.979642e+01 8.437622e+01 1.179181e+02
#> [786] 9.138436e+01 8.357337e+01 1.822999e+02 8.735049e+01 7.206720e+01
#> [791] 1.110641e+02 4.418442e+02 7.807747e+01 1.312790e+02 8.393741e+01
#> [796] 7.392201e+01 1.036455e+02 2.548565e+02 9.643718e+01 9.197963e+01
#> [801] 9.447140e+01 8.025034e+01 6.710468e+01 1.407838e+02 3.274664e+02
#> [806] 8.000429e+01 1.413720e+02 9.985223e+02 1.379333e+02 6.835176e+01
#> [811] 7.719082e+01 8.473138e+01 2.676810e+02 8.757422e+01 9.038539e+01
#> [816] 7.451720e+01 2.798526e+02 9.543949e+01 1.245684e+02 1.176589e+02
#> [821] 1.856986e+02 8.984961e+01 6.986043e+01 8.460862e+01 8.768672e+01
#> [826] 2.529359e+02 2.991581e+02 9.471875e+01 8.506682e+01 1.076897e+02
#> [831] 1.279349e+02 1.291407e+02 1.229130e+02 6.983130e+01 7.905521e+01
#> [836] 3.638210e+08 1.045398e+03 7.995731e+01 1.205722e+02 2.157733e+02
#> [841] 1.643172e+02 2.532609e+03 7.761262e+02 1.495781e+02 1.285079e+02
#> [846] 8.965357e+01 1.615550e+02 7.168068e+01 1.224054e+02 6.582757e+01
#> [851] 2.337056e+02 6.496156e+01 4.747713e+02 1.119015e+02 7.382246e+01
#> [856] 7.424878e+01 3.409592e+02 1.562314e+06 1.842300e+02 6.403375e+01
#> [861] 6.818002e+01 7.426231e+01 6.499770e+01 6.971435e+01 1.024475e+02
#> [866] 1.070345e+02 8.363625e+01 2.725342e+02 9.977361e+01 7.119264e+01
#> [871] 9.219679e+02 2.439390e+02 6.928054e+01 6.082265e+01 1.184518e+02
#> [876] 2.185940e+02 1.747809e+02 6.688551e+01 7.081434e+01 6.579341e+01
#> [881] 8.708891e+01 8.421004e+01 7.415997e+01 8.450406e+01 8.216443e+01
#> [886] 6.477609e+01 6.377165e+01 3.073020e+02 7.486406e+01 8.284328e+01
#> [891] 7.927487e+01 7.222010e+01 1.191394e+02 1.099068e+02 2.663087e+02
#> [896] 1.077231e+02 8.511155e+01 1.287698e+02 7.203135e+01 1.252677e+02
#> [901] 1.103159e+02 9.150984e+01 7.418781e+01 1.010868e+02 1.030331e+02
#> [906] 7.240335e+02 6.429677e+01 7.668686e+01 1.135892e+02 7.461102e+01
#> [911] 6.880593e+01 8.660400e+01 3.943866e+02 1.126769e+02 2.400539e+03
#> [916] 1.051064e+02 6.890446e+01 6.402680e+01 1.009559e+02 2.106453e+02
#> [921] 8.223901e+02 1.430304e+03 1.603218e+02 8.457824e+01 8.526568e+01
#> [926] 1.695698e+02 6.477049e+01 7.834650e+01 1.546572e+02 8.188970e+01
#> [931] 6.695854e+01 7.191631e+01 6.336947e+01 7.661768e+01 7.436798e+01
#> [936] 1.389268e+02 5.013567e+02 1.527089e+02 8.714944e+01 7.577238e+01
#> [941] 2.631227e+02 1.083238e+02 9.473158e+01 7.838267e+01 1.436712e+02
#> [946] 1.944142e+02 6.756685e+01 8.301404e+01 8.171365e+01 2.199446e+02
#> [951] 6.994482e+01 6.484412e+01 7.840281e+01 1.223328e+02 1.461355e+02
#> [956] 1.519540e+05 1.720863e+03 6.978114e+01 6.902791e+01 7.450780e+01
#> [961] 7.275311e+01 5.451895e+02 1.056436e+02 8.652676e+01 7.175402e+01
#> [966] 7.489085e+01 2.242256e+02 1.591803e+02 1.098188e+02 7.139546e+01
#> [971] 1.071567e+02 7.240830e+01 8.124014e+01 7.033155e+01 1.615390e+02
#> [976] 1.406582e+02 1.491478e+02 6.978217e+01 6.914256e+01 6.959815e+01
#> [981] 6.524605e+01 1.029824e+02 9.533460e+01 1.162713e+02 6.570532e+02
#> [986] 6.418485e+01 6.753631e+01 6.638922e+01 7.258883e+01 8.717808e+01
#> [991] 9.107560e+01 6.650324e+01 2.128055e+02 6.472440e+01 8.099975e+01
#> [996] 6.718579e+01 7.002967e+01 7.063852e+01 9.828024e+01 8.441626e+01
doseFunction(
model,
alpha2 = samples@data$alpha2,
beta = samples@data$beta,
grade = 2L
)(x = 0.25)
#> [1] 61.08100 62.29098 2356.22065 72.32862 58.28887
#> [6] 74.79160 59.35864 62.94879 86.62611 58.87231
#> [11] 91.83389 85.39065 63.10887 57.85366 58.92140
#> [16] 68.04239 62.16223 57.72847 62.27789 54.69438
#> [21] 59.75792 56.62605 60.39671 60.52880 61.94425
#> [26] 53.52051 50.88591 59.51032 67.71141 67.85491
#> [31] 66.56328 75.89371 56.91408 62.11852 54.59148
#> [36] 60.39182 49.62354 61.00253 57.88572 87.79192
#> [41] 76.73953 60.50293 70.93297 118.10356 62.78856
#> [46] 53.70962 70.86817 65.47367 67.39473 72.19507
#> [51] 68.71276 129.26540 89.53180 67.43814 65.06243
#> [56] 61.53395 91.62348 73.12801 57.88077 79.64608
#> [61] 465.88584 62.16860 67.46002 59.16017 56.78452
#> [66] 64.99996 80.28496 74.63145 102.92039 52.25250
#> [71] 67.95220 135.48885 76.97631 130.45689 66.29687
#> [76] 56.67327 62.11758 86.39531 55.69880 63.85406
#> [81] 135.80263 71.99302 59.84180 61.65266 114.81193
#> [86] 71.98142 68.65283 74.70747 76.06877 97.36187
#> [91] 69.69893 57.34365 66.20089 59.13266 64.42450
#> [96] 65.50753 59.73033 79.08746 52.83364 1790.59040
#> [101] 53.73441 69.34586 73.47613 173.90690 61.61088
#> [106] 62.54576 69.81254 91.73745 162.78013 178.77125
#> [111] 77.50014 49.58902 61.74283 71.36458 98.78450
#> [116] 60.33244 97.48345 60.27020 219.52497 63.54698
#> [121] 64.67122 62.28462 63.13807 69.75704 73.11506
#> [126] 77.01405 55.69097 89.74344 84.32646 74.43006
#> [131] 55.41671 64.76834 56.05973 143.02396 117.37478
#> [136] 76.89265 93.13613 53.14483 56.45978 56.39357
#> [141] 69.65891 199.47578 55.98482 74.16786 70.55159
#> [146] 64.56670 60.17895 80.55830 67.47581 68.08584
#> [151] 79.92459 87.37771 120.38134 84.79044 76.60114
#> [156] 66.57059 50.31297 62.13660 178.72330 73.74377
#> [161] 84.21755 56.54998 55.43577 63.48429 88.30328
#> [166] 91.13677 125.01044 58.04696 69.98179 49.30608
#> [171] 66.79088 59.53986 78.70978 65.09822 56.41847
#> [176] 73.97903 97.44015 380.97064 61.64581 58.07620
#> [181] 60.45610 63.09521 168.62175 97.98601 55.97403
#> [186] 99.04591 49.56957 115.02284 571.34068 51.69408
#> [191] 127.05248 80.99301 71.11237 57.68820 225.56653
#> [196] 59.93169 56.91888 68.21239 60.78337 68.04325
#> [201] 59.58625 61.11208 61.09727 90.37592 60.18243
#> [206] 61.28164 61.63712 63.41514 54.37704 58.85548
#> [211] 455.74899 65.73591 54.56856 56.61811 78.48533
#> [216] 45.53020 52.19277 76.17375 136.25017 84.77642
#> [221] 57.89452 56.88641 63.43799 58.81249 60.19894
#> [226] 72.30441 86.82859 78.57239 63.85924 97.70370
#> [231] 55.79472 117.00638 79.21934 59.94028 95.91019
#> [236] 80.45092 77.79783 74.11009 53.58855 62.64850
#> [241] 60.23789 61.31816 112.83884 195.00158 73.91248
#> [246] 79.32687 59.92168 58.71383 71.66608 79.33233
#> [251] 130.36906 97.98425 78.60348 62.81726 59.21488
#> [256] 68.11838 103.86902 73.50870 57.26869 54.58172
#> [261] 52.60286 79.45358 52.58811 122.18044 53.80471
#> [266] 55.27902 82.97086 96.03035 55.45457 76.85216
#> [271] 76.90842 67.54799 121.37921 67.27203 60.42328
#> [276] 63.15571 55.12349 57.93519 65.45048 76.17412
#> [281] 62.81690 64.22104 75.14493 53.87690 77.09381
#> [286] 61.58607 59.13814 64.48749 57.53051 172.96458
#> [291] 58.84086 80.53426 66.53095 88.59013 71.33340
#> [296] 60.52385 97.27153 367.80081 134.84578 58.48717
#> [301] 75.81576 65.75247 62.90614 55.35557 63.84231
#> [306] 98.33595 84.88529 66.90610 76.91857 68.46775
#> [311] 77.39054 256.67065 64.95079 58.21746 74.26075
#> [316] 122.97316 63.27375 60.28218 91.52841 52.52688
#> [321] 45.42063 57.57504 75.59099 88.43642 65.53793
#> [326] 64.69494 76.29768 60.52613 70.19953 63.84978
#> [331] 52.46607 59.58620 56.48989 63.92479 60.43171
#> [336] 55.23018 162.04085 59.86762 82.69550 58.96208
#> [341] 69.48918 58.73150 87.39772 63.27988 89.46739
#> [346] 73.32195 59.16608 421.63120 58.50830 68.72937
#> [351] 95.45801 67.85200 91.67151 56.63442 90.97178
#> [356] 58.30022 80.69412 63.28525 346.65163 127.46974
#> [361] 65.06252 63.34470 64.27614 78.96590 91.63868
#> [366] 70.19033 52.12706 68.13839 99.09690 51.11851
#> [371] 55.35275 97.33867 57.95028 128.60631 98.89128
#> [376] 59.34033 69.81908 52.21771 126.95853 1141.37482
#> [381] 68.43605 76.05988 69.83342 66.17083 86.89422
#> [386] 205.82803 133.83031 187.17896 64.81152 44.65569
#> [391] 166.57145 58.97239 67.26975 69.80944 83.79976
#> [396] 63.89147 111.74531 76.16972 60.35172 65.06964
#> [401] 59.87228 60.39405 61.52071 57.58124 65.58403
#> [406] 100.35284 71.30791 61.22407 59.49754 63.16256
#> [411] 59.67186 67.58905 62.21140 67.89865 69.59155
#> [416] 59.91636 64.63877 69.97075 132.35967 66.26370
#> [421] 67.08476 59.60357 60.38398 75.62618 54.02015
#> [426] 70.14879 63.69428 53.81651 90.69363 98.69213
#> [431] 62.61705 62.73145 60.60734 58.11160 67.36500
#> [436] 74.85330 55.74881 61.04283 58.92187 61.50550
#> [441] 80.19943 103.28309 115.90279 796.25305 75.51174
#> [446] 162.57353 83.46711 62.44644 67.29815 64.65025
#> [451] 56.30577 64.25510 59.10955 86.72119 77.04104
#> [456] 52.40649 58.41241 51.77015 57.22661 55.75381
#> [461] 54.59883 62.70907 56.63676 63.10571 58.38679
#> [466] 152.19884 78.70994 72.97645 59.64718 62.32703
#> [471] 54.35683 57.22731 67.41235 65.45023 60.31427
#> [476] 86.22689 85.11519 109.90902 52.11413 65.26739
#> [481] 58.22779 59.33329 68.79781 65.60303 62.26562
#> [486] 56.92928 57.80962 61.49966 66.84691 63.62869
#> [491] 58.69225 59.20122 62.65286 125.84176 53.95272
#> [496] 61.92287 89.16984 64.85920 61.37272 75.64143
#> [501] 76.03663 64.22305 56.20439 56.55714 60.98649
#> [506] 56.32103 63.29939 57.86155 58.89549 63.02455
#> [511] 61.02760 66.85858 82.52539 72.68608 78.22264
#> [516] 61.70015 58.45264 72.60717 64.58831 56.16593
#> [521] 66.45167 57.27331 57.91219 59.40754 93.89975
#> [526] 63.12671 68.80883 102.62730 67.52673 58.48963
#> [531] 72.67495 55.66646 62.72475 60.03801 92.25448
#> [536] 76.99146 64.24743 64.46529 72.17570 55.79802
#> [541] 64.75536 69.86541 58.27326 69.68346 63.29205
#> [546] 2520.41960 130.74725 64.98434 74.84949 61.44960
#> [551] 68.89111 60.63990 57.70754 87.84794 55.08125
#> [556] 58.58615 67.56210 60.13065 58.43811 60.71788
#> [561] 94.53880 58.24363 64.77537 155.24244 57.98855
#> [566] 57.88173 668.10024 65.76409 57.59547 79.54412
#> [571] 124.02820 63.99236 61.42455 58.61295 56.04277
#> [576] 68.07988 77.98887 81.88178 73.44029 64.50256
#> [581] 65.19236 59.34142 60.54963 56.00571 75.21556
#> [586] 62.58214 57.28653 64.84208 92.04146 69.07745
#> [591] 52.76798 58.14627 90.55416 58.07792 117.09201
#> [596] 59.19536 64.64901 63.39586 88.19784 138.13525
#> [601] 62.08266 58.15082 77.68040 66.97460 63.59482
#> [606] 181.64608 64.87266 76.61458 60.00749 57.93023
#> [611] 2330.26861 71.39056 56.65739 57.26085 67.18603
#> [616] 56.69625 64.27178 68.02544 112.39076 64.37627
#> [621] 85.25221 60.76792 62.92096 54.87106 74.50763
#> [626] 67.44851 64.61275 63.63610 67.23917 56.74643
#> [631] 96.43539 52.39040 66.64578 74.36111 124.91971
#> [636] 70.44565 104.50188 109.93690 402.19291 60.59828
#> [641] 125.27617 64.28665 126.05668 54.02439 63.07201
#> [646] 61.08647 64.43233 61.74395 78.13275 62.88196
#> [651] 72.69387 98.14423 61.34235 60.48626 58.31563
#> [656] 60.80778 60.30719 57.80779 62.27146 59.69152
#> [661] 66.08535 61.67519 87.71646 69.61711 64.96076
#> [666] 61.57603 68.73875 63.84176 64.68608 61.39703
#> [671] 57.76530 55.69630 55.91002 75.72533 53.49283
#> [676] 89.07314 61.62535 53.03241 52.19541 89.25369
#> [681] 61.02140 100.60019 54.30727 161.22499 208.13821
#> [686] 174.20099 66.49770 80.80717 65.99059 57.49574
#> [691] 64.45716 60.01207 58.26244 60.55026 58.19310
#> [696] 78.37888 124.70767 60.12604 139.70998 55.89978
#> [701] 60.38103 90.78991 81.31184 61.94311 54.14750
#> [706] 65.42680 69.39416 71.08502 57.58545 59.10386
#> [711] 63.17327 100.98143 58.63117 60.03809 59.36876
#> [716] 75.42900 61.29786 201.53553 56.81184 83.58003
#> [721] 88.89785 69.21768 59.57720 61.58865 52.52218
#> [726] 127.34669 79.00857 64.89842 87.67311 56.53849
#> [731] 58.17319 56.70580 61.73418 100.69790 108.45159
#> [736] 60.02688 59.50855 60.18946 58.39352 73.63229
#> [741] 57.29523 58.26811 83.15466 52.74852 57.85371
#> [746] 85.06813 301.99801 134.35156 59.29692 57.14842
#> [751] 57.88373 104.77193 55.99638 88.44482 63.41197
#> [756] 1590.45475 146.24427 83.30913 61.43872 207.53677
#> [761] 123.14590 61.99321 65.56703 77.97139 53.41425
#> [766] 64.57049 100.37285 50.32805 74.12900 80.83929
#> [771] 62.58623 65.79094 135.33403 64.42144 58.12886
#> [776] 62.97625 60.97151 55.49273 131.56716 541.78394
#> [781] 68.68970 57.66216 65.50377 61.61873 88.99547
#> [786] 74.02569 70.45126 91.46434 66.71149 64.37102
#> [791] 87.41715 106.57425 59.16050 79.91373 62.30650
#> [796] 61.68212 70.19701 81.78097 71.79327 66.00771
#> [801] 73.07180 62.62332 60.16871 94.19858 163.85535
#> [806] 66.75176 82.15622 282.42482 81.29047 58.29093
#> [811] 63.79732 71.57178 119.13622 61.23546 65.48433
#> [816] 65.31045 141.93445 70.88340 69.87558 81.70710
#> [821] 99.64518 64.51334 56.09308 64.02667 59.57846
#> [826] 112.86541 67.24724 68.81670 66.44365 60.24856
#> [831] 73.09153 90.56698 79.31639 56.82367 61.35175
#> [836] 160019.80623 94.40235 59.68513 80.17305 106.88812
#> [841] 75.12708 251.75442 158.91092 74.95595 60.81649
#> [846] 71.57901 75.10836 63.63385 81.04028 59.34567
#> [851] 146.47330 59.62819 145.22647 87.52895 64.00863
#> [856] 62.53791 121.90885 14036.36797 101.89404 58.93454
#> [861] 60.28326 63.81888 57.65897 60.54015 78.16703
#> [866] 75.25180 59.85945 107.19185 67.62128 57.23618
#> [871] 269.09259 137.05537 59.97401 56.23009 83.70276
#> [876] 121.73893 71.71347 59.66458 57.28111 54.84079
#> [881] 72.75248 56.38991 58.11126 62.87075 59.09017
#> [886] 56.11707 58.18926 162.37963 63.44906 68.74929
#> [891] 68.00951 63.59841 68.30167 65.82200 69.27882
#> [896] 72.80528 53.95653 69.25089 60.62041 68.40143
#> [901] 78.91234 69.94904 64.20078 70.59836 75.50634
#> [906] 183.87399 55.47510 61.13929 78.50650 61.68485
#> [911] 61.23385 67.76580 156.81099 54.19878 541.98892
#> [916] 71.60963 56.15502 56.98920 69.30086 117.70733
#> [921] 147.06604 343.52341 90.77724 54.05169 67.49005
#> [926] 82.99543 58.22144 63.45649 77.43407 66.59963
#> [931] 54.21869 62.08075 55.59251 66.83760 61.97439
#> [936] 97.85230 166.56819 110.02621 61.76693 53.54844
#> [941] 88.49396 56.12563 69.67839 61.69356 83.91956
#> [946] 123.25313 60.73295 68.60133 65.43058 94.14106
#> [951] 58.33678 58.53402 63.17047 85.93014 78.29230
#> [956] 6532.90948 105.32298 56.14132 57.85018 64.40839
#> [961] 62.93660 149.50483 77.49391 67.29442 60.13104
#> [966] 61.23162 123.28982 93.78594 70.16690 58.58835
#> [971] 72.03812 51.16720 65.25137 59.09269 87.63349
#> [976] 71.72238 81.42506 60.02916 61.35884 59.37729
#> [981] 57.19125 54.37720 76.67099 75.55840 205.04046
#> [986] 57.06185 59.82791 57.39208 59.91060 63.00882
#> [991] 67.26693 58.40903 120.09357 58.11795 60.81948
#> [996] 56.83358 58.98511 61.96270 74.58539 63.67324
