![[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] 6.863119e+01 8.116661e+01 8.056249e+01 9.789948e+01 8.534848e+01
#> [6] 2.678572e+02 6.738993e+01 6.783162e+01 6.597161e+01 7.979540e+01
#> [11] 6.101148e+01 6.490478e+01 7.356286e+01 1.036092e+02 2.953936e+02
#> [16] 7.641425e+01 7.518301e+01 1.729182e+02 7.509025e+02 8.008204e+01
#> [21] 7.538115e+01 8.706357e+01 1.068381e+02 1.113716e+02 6.601910e+01
#> [26] 9.195059e+01 1.100421e+02 1.321310e+02 1.087687e+02 1.075418e+02
#> [31] 1.187392e+02 6.508624e+01 3.769746e+02 6.982654e+01 4.799053e+02
#> [36] 4.917576e+02 1.108341e+02 8.221024e+01 7.041046e+01 7.173671e+01
#> [41] 8.715487e+01 1.189479e+02 1.933385e+02 1.056049e+02 1.055013e+02
#> [46] 8.244467e+03 7.609084e+01 1.897277e+08 2.993327e+02 2.247446e+02
#> [51] 1.910890e+02 2.191096e+02 7.630036e+01 7.310496e+01 7.552398e+01
#> [56] 1.464232e+02 8.310554e+01 7.030816e+01 6.493470e+01 6.915696e+01
#> [61] 8.280732e+01 1.123547e+02 6.812175e+01 7.824933e+01 1.777793e+04
#> [66] 2.342798e+02 7.697767e+01 7.119191e+01 6.361567e+01 2.674276e+05
#> [71] 7.871941e+01 7.283810e+01 1.654118e+02 7.207963e+01 2.101412e+02
#> [76] 1.024619e+02 1.097905e+02 7.579264e+01 9.406576e+01 2.128580e+02
#> [81] 1.888990e+02 7.700890e+01 1.004256e+02 9.401758e+01 5.526996e+02
#> [86] 7.878857e+01 7.897325e+01 8.490546e+01 1.005343e+02 7.099758e+01
#> [91] 2.248149e+02 1.264135e+02 2.440144e+02 9.391836e+01 6.739373e+01
#> [96] 6.677058e+01 1.031211e+02 9.578381e+01 1.167648e+02 1.028726e+02
#> [101] 7.370738e+01 1.137985e+02 1.097165e+02 1.015942e+02 7.650982e+01
#> [106] 9.746415e+01 2.206803e+02 5.847198e+02 6.309031e+01 6.041101e+01
#> [111] 6.254185e+01 1.071290e+02 1.028460e+02 8.429482e+01 1.128053e+02
#> [116] 6.543004e+01 9.755975e+01 7.090529e+01 9.495814e+01 7.856989e+01
#> [121] 1.191789e+02 8.256751e+01 7.940584e+01 8.366351e+01 7.705593e+01
#> [126] 8.165142e+01 7.761611e+01 8.834370e+01 8.711546e+01 1.073845e+02
#> [131] 8.112719e+01 9.747541e+01 9.792188e+01 8.912032e+01 6.729034e+02
#> [136] 6.509773e+01 6.497258e+01 6.629010e+01 6.430836e+01 8.527523e+01
#> [141] 1.280482e+03 1.123418e+02 2.963844e+02 1.465246e+03 3.295531e+02
#> [146] 1.018387e+02 1.361760e+02 3.356283e+03 1.158408e+02 7.254534e+01
#> [151] 7.952873e+01 7.533749e+01 6.541019e+01 7.866723e+01 2.949978e+02
#> [156] 6.643902e+01 7.133941e+01 1.074616e+03 2.681038e+02 1.154006e+02
#> [161] 1.891370e+02 8.433050e+01 8.599579e+01 9.406086e+01 6.972453e+01
#> [166] 2.201286e+02 1.055451e+02 1.086471e+03 1.088281e+02 2.268617e+02
#> [171] 1.048295e+02 5.333503e+02 9.746024e+01 9.023170e+01 7.971353e+01
#> [176] 1.081663e+02 1.098626e+02 6.461381e+01 1.255767e+02 2.852689e+02
#> [181] 1.095153e+02 1.727817e+02 1.660298e+03 5.067975e+02 6.372396e+01
#> [186] 6.603619e+01 8.011739e+01 8.469302e+01 8.206754e+01 8.861976e+01
#> [191] 8.594571e+01 1.138767e+02 7.737286e+01 1.594029e+02 1.019056e+02
#> [196] 7.610599e+01 9.401592e+01 1.031511e+02 7.121362e+01 8.557943e+01
#> [201] 8.651063e+01 1.621558e+02 6.630499e+01 1.504705e+02 7.687791e+01
#> [206] 7.534067e+01 9.846970e+01 8.779108e+01 1.081890e+02 6.475687e+01
#> [211] 6.585070e+01 6.230957e+02 3.866587e+02 4.033499e+02 1.665928e+02
#> [216] 7.500890e+01 6.895175e+01 6.697686e+01 6.958822e+01 9.649781e+01
#> [221] 7.379089e+03 6.251283e+03 6.302059e+01 6.726595e+01 3.371056e+02
#> [226] 1.050481e+02 6.988950e+01 7.109168e+01 6.396744e+01 2.019911e+02
#> [231] 7.459846e+01 1.104971e+02 5.345787e+02 1.349021e+02 9.380338e+01
#> [236] 8.681688e+01 9.543806e+01 8.204008e+01 1.374033e+02 7.204146e+01
#> [241] 7.386019e+01 4.076000e+03 1.710310e+02 2.596148e+02 2.044435e+02
#> [246] 1.412947e+02 6.499148e+01 8.465452e+01 8.313368e+01 1.420499e+02
#> [251] 9.432491e+01 8.825978e+01 7.404074e+01 1.154572e+02 3.877733e+02
#> [256] 1.398333e+02 1.167390e+02 7.515473e+01 1.047774e+02 6.914315e+01
#> [261] 9.832944e+01 7.099664e+02 7.407362e+01 1.151386e+02 3.277829e+02
#> [266] 7.510682e+01 8.844355e+01 9.112918e+01 4.012495e+02 6.759980e+01
#> [271] 7.236422e+01 1.727356e+02 7.627426e+01 7.669468e+01 6.394698e+01
#> [276] 1.037491e+02 7.835586e+01 1.354382e+03 8.029726e+01 9.647742e+01
#> [281] 8.188433e+01 7.608671e+01 7.600847e+01 2.505113e+02 4.283842e+02
#> [286] 2.301868e+02 7.867362e+01 9.793039e+01 1.540598e+02 8.115365e+01
#> [291] 6.640335e+01 7.241811e+01 7.117124e+01 6.756699e+01 6.353679e+01
#> [296] 1.379612e+02 9.912454e+01 1.897997e+06 3.899179e+02 1.110288e+02
#> [301] 9.696956e+01 6.633689e+01 7.287626e+01 6.919301e+01 6.802847e+01
#> [306] 1.212590e+02 1.078110e+03 6.460236e+01 6.760269e+01 1.124944e+02
#> [311] 8.386192e+02 1.109992e+02 6.176203e+01 1.593487e+02 1.043087e+02
#> [316] 7.650192e+01 8.955977e+01 2.779115e+02 7.128664e+01 7.190467e+01
#> [321] 7.538114e+01 2.056240e+02 4.758770e+02 8.401293e+01 7.462376e+01
#> [326] 6.314747e+01 1.704012e+02 7.315725e+01 1.333767e+02 1.169819e+02
#> [331] 1.348590e+02 3.208815e+02 1.108367e+02 2.174954e+02 1.458923e+03
#> [336] 2.848618e+03 6.240895e+01 8.912641e+01 6.779911e+01 6.373593e+01
#> [341] 6.678535e+01 6.892995e+01 6.701902e+01 7.703065e+01 8.471032e+01
#> [346] 4.221691e+02 6.937891e+01 7.011754e+01 9.069088e+01 7.527446e+01
#> [351] 9.376321e+01 9.682008e+01 6.792352e+01 7.371881e+01 7.455791e+01
#> [356] 6.926318e+01 6.876952e+01 6.979467e+01 5.081327e+03 3.376048e+02
#> [361] 1.795979e+02 1.518086e+03 8.065944e+01 1.774769e+02 1.574771e+02
#> [366] 8.665073e+01 6.510329e+01 9.806045e+01 8.122127e+01 8.634598e+01
#> [371] 9.177892e+01 9.033193e+01 6.880431e+01 8.723118e+01 8.407930e+01
#> [376] 6.515591e+01 6.516670e+01 1.112161e+02 9.161017e+01 6.195087e+01
#> [381] 8.414042e+01 9.671242e+01 7.544125e+01 6.973529e+01 8.577355e+01
#> [386] 8.242515e+01 1.261303e+02 6.515539e+01 9.394061e+01 8.237292e+01
#> [391] 2.105002e+02 6.963006e+01 1.099280e+02 7.903145e+01 1.636415e+02
#> [396] 7.253078e+01 8.051194e+01 6.653826e+01 7.303589e+01 7.812364e+01
#> [401] 8.753461e+01 1.148172e+02 3.481980e+02 6.785035e+01 1.707207e+02
#> [406] 1.545542e+02 8.894823e+01 6.946014e+01 7.049759e+01 9.742566e+01
#> [411] 7.713015e+01 4.369845e+02 8.783002e+01 6.519502e+01 6.646855e+01
#> [416] 1.706083e+02 6.249586e+01 7.374863e+01 7.457338e+01 8.491453e+01
#> [421] 7.033181e+01 4.649399e+02 2.631749e+03 1.571649e+02 1.303743e+02
#> [426] 7.384252e+01 7.625808e+01 3.150230e+03 7.703714e+01 8.817127e+01
#> [431] 8.329830e+01 1.488709e+02 8.092882e+01 2.195390e+02 1.950814e+02
#> [436] 7.190709e+01 4.379528e+02 1.936949e+02 7.935058e+01 8.038224e+03
#> [441] 9.495040e+01 7.272532e+01 6.719669e+01 1.457668e+02 2.049129e+02
#> [446] 8.340166e+01 4.417870e+02 3.144428e+02 7.362543e+01 7.806813e+01
#> [451] 7.536117e+01 3.527829e+02 1.348185e+02 7.311342e+01 9.503790e+01
#> [456] 7.612406e+01 9.266145e+01 6.613383e+01 8.883253e+01 8.418497e+01
#> [461] 4.231658e+02 1.343237e+02 3.813198e+02 2.687097e+02 9.858336e+01
#> [466] 1.264012e+02 1.213024e+02 9.397154e+01 1.347792e+02 9.554354e+01
#> [471] 2.006616e+03 1.004983e+02 5.543423e+02 6.867048e+01 6.665569e+01
#> [476] 7.223370e+01 6.788865e+01 9.369559e+01 1.444863e+02 9.260302e+01
#> [481] 7.028251e+01 2.667701e+02 6.440834e+01 7.298212e+01 1.156219e+02
#> [486] 2.191051e+02 8.409238e+01 2.034342e+02 6.885151e+01 6.401232e+01
#> [491] 7.145446e+01 1.197754e+02 1.262810e+02 7.238272e+01 6.582092e+01
#> [496] 1.396776e+02 6.865850e+01 1.244021e+02 1.606505e+02 8.133920e+02
#> [501] 9.054005e+01 7.910614e+01 7.648381e+01 8.146617e+01 7.616001e+01
#> [506] 8.966769e+01 8.480908e+01 8.647065e+01 1.008361e+02 1.339230e+02
#> [511] 1.074907e+03 1.072800e+03 1.278154e+02 7.217422e+01 1.529426e+02
#> [516] 9.643087e+01 1.204713e+02 8.981628e+01 7.444387e+02 8.119268e+01
#> [521] 9.438414e+01 1.159718e+02 8.296080e+01 1.694591e+02 6.619703e+01
#> [526] 7.626400e+01 9.221442e+01 1.353197e+02 7.028337e+01 8.148824e+01
#> [531] 1.072162e+02 7.173070e+01 9.225522e+01 1.258733e+02 8.478351e+01
#> [536] 1.157092e+02 5.663063e+02 7.278682e+01 7.464873e+01 6.559363e+01
#> [541] 8.998527e+01 2.308871e+02 2.166078e+02 8.990380e+03 7.272763e+01
#> [546] 9.634719e+01 1.118779e+02 6.745597e+01 7.167676e+01 1.285884e+02
#> [551] 7.919587e+01 6.855630e+01 8.808747e+01 1.593792e+02 3.001772e+02
#> [556] 1.275087e+02 9.057436e+01 8.260512e+01 4.680043e+02 3.452882e+02
#> [561] 9.996796e+01 1.636494e+02 2.749755e+03 4.167638e+02 7.915324e+01
#> [566] 3.245422e+02 1.173980e+02 1.342094e+02 2.412295e+03 1.096875e+02
#> [571] 9.248560e+01 1.067177e+02 7.709290e+01 1.544300e+02 4.527380e+02
#> [576] 7.808773e+01 9.168781e+01 1.331734e+02 7.875204e+01 7.314221e+01
#> [581] 6.440411e+01 6.637119e+01 7.598782e+01 7.738490e+01 6.683651e+01
#> [586] 6.720131e+01 1.132797e+02 3.979561e+02 7.924473e+01 6.892228e+01
#> [591] 7.111751e+01 7.160101e+01 1.056294e+02 1.210371e+02 1.050409e+02
#> [596] 6.884478e+01 9.181473e+01 1.722773e+03 3.881408e+05 2.206406e+02
#> [601] 1.445176e+02 6.282944e+01 7.181495e+01 6.158483e+01 1.026998e+02
#> [606] 6.959020e+01 1.879726e+02 1.223820e+02 1.864778e+02 6.543710e+01
#> [611] 8.090240e+01 1.012825e+02 1.931371e+02 2.709645e+02 8.850734e+01
#> [616] 6.628211e+01 9.537689e+01 7.769271e+01 7.114765e+01 8.046449e+01
#> [621] 2.943527e+02 1.480363e+02 1.267315e+02 6.286296e+01 7.127156e+01
#> [626] 1.445043e+02 8.306674e+01 6.875288e+01 1.042674e+02 6.740990e+01
#> [631] 2.519673e+03 2.006037e+02 2.539836e+03 7.182437e+01 1.147980e+02
#> [636] 7.734792e+01 9.750280e+01 1.751412e+02 9.891710e+01 6.175355e+01
#> [641] 7.429380e+01 8.662062e+01 2.560888e+03 3.274342e+02 1.635815e+02
#> [646] 8.306169e+01 1.429634e+02 7.865297e+01 7.560380e+01 6.787971e+01
#> [651] 8.617378e+01 6.950403e+01 7.369486e+01 3.084263e+02 7.287791e+01
#> [656] 9.520887e+01 1.197034e+02 3.903502e+02 3.219608e+02 6.489241e+01
#> [661] 9.157519e+01 7.875287e+01 6.424203e+01 2.762173e+02 7.636982e+01
#> [666] 2.716783e+02 6.777983e+01 3.436319e+02 9.218003e+01 1.027561e+02
#> [671] 9.254154e+01 1.050791e+02 7.606029e+01 7.051364e+01 8.506869e+01
#> [676] 1.448016e+02 3.892085e+03 1.041145e+02 7.245251e+01 1.114854e+02
#> [681] 7.751408e+01 1.061718e+02 1.542165e+02 9.084298e+01 2.489366e+02
#> [686] 6.950599e+01 6.784495e+01 2.081396e+03 6.276503e+01 7.504489e+01
#> [691] 1.108720e+02 7.353598e+01 9.623938e+01 9.192983e+01 7.525566e+01
#> [696] 1.963056e+02 6.656999e+01 6.285665e+01 7.989342e+01 9.960900e+01
#> [701] 7.523505e+01 5.988821e+01 7.938734e+01 6.467968e+01 1.055786e+02
#> [706] 1.081696e+02 7.658512e+01 7.989574e+01 5.500081e+02 7.990116e+02
#> [711] 6.326832e+01 7.068126e+01 7.239226e+01 8.785534e+01 7.939494e+01
#> [716] 9.009276e+02 6.192584e+01 6.984718e+01 1.249329e+02 1.590367e+02
#> [721] 8.859319e+01 8.141721e+01 7.484271e+01 7.051894e+01 7.903447e+01
#> [726] 8.273569e+01 8.218595e+01 1.296236e+03 1.112448e+02 1.998473e+02
#> [731] 8.166887e+01 8.368459e+01 7.631556e+01 9.098099e+01 6.857739e+01
#> [736] 1.464967e+02 7.408598e+01 9.215864e+04 1.186250e+02 6.802885e+01
#> [741] 7.090224e+01 8.774427e+01 7.141234e+01 1.160248e+02 9.522165e+01
#> [746] 1.215722e+02 7.892442e+01 1.251766e+02 7.415319e+01 1.021224e+02
#> [751] 1.000021e+02 1.214444e+02 7.914003e+01 7.067446e+01 6.718126e+01
#> [756] 1.017487e+02 8.868887e+01 8.490291e+01 6.517257e+01 7.256369e+01
#> [761] 7.526944e+01 7.878392e+01 7.330822e+01 6.625327e+01 8.550908e+01
#> [766] 6.872783e+01 9.520965e+01 8.087439e+01 1.042940e+02 3.350582e+02
#> [771] 1.167985e+02 8.301989e+01 1.018829e+02 3.790633e+03 1.283581e+02
#> [776] 1.041076e+02 1.151004e+02 7.470536e+01 2.802698e+02 1.914148e+02
#> [781] 1.179634e+02 1.351732e+02 6.260611e+01 8.106740e+01 6.731071e+01
#> [786] 8.367251e+01 8.954540e+01 7.609938e+01 7.540383e+01 9.535766e+01
#> [791] 1.137041e+02 3.319684e+03 1.906876e+02 7.473540e+01 7.096240e+01
#> [796] 6.921167e+01 4.749498e+04 9.418912e+01 1.656455e+02 8.018405e+01
#> [801] 1.592935e+02 1.367945e+02 1.722319e+02 2.327503e+02 9.515355e+02
#> [806] 5.520550e+08 8.555988e+01 7.125170e+01 4.386670e+02 7.203641e+01
#> [811] 7.045487e+02 1.477558e+02 9.061863e+01 6.848631e+02 1.104674e+02
#> [816] 1.010634e+02 8.094513e+01 1.968067e+02 9.808213e+01 7.523181e+01
#> [821] 6.842317e+01 6.527728e+01 6.371655e+01 8.507693e+01 2.861158e+02
#> [826] 5.899384e+02 6.936261e+01 9.961503e+01 1.043655e+02 1.431205e+02
#> [831] 1.311962e+02 1.996892e+03 6.639011e+01 7.655823e+01 7.324987e+01
#> [836] 8.712555e+01 7.262816e+01 1.900751e+02 1.087406e+02 3.419524e+02
#> [841] 9.357703e+01 1.033752e+03 1.628053e+02 6.850040e+02 6.196050e+01
#> [846] 2.865391e+02 3.789756e+02 7.252477e+01 6.467840e+01 7.131820e+01
#> [851] 8.662796e+02 2.842670e+02 1.920678e+02 1.374056e+02 6.198306e+01
#> [856] 1.226868e+03 2.752790e+02 8.051503e+01 1.056401e+02 7.301205e+01
#> [861] 9.854078e+01 8.846965e+01 9.421429e+01 2.805613e+02 1.154994e+02
#> [866] 7.544221e+01 9.282312e+01 7.651930e+01 1.082688e+02 8.575569e+01
#> [871] 1.071053e+02 7.881980e+01 1.611371e+02 1.045410e+02 1.720019e+02
#> [876] 2.259735e+02 4.490344e+04 1.317575e+02 8.187952e+01 8.742081e+01
#> [881] 2.232051e+05 6.286987e+01 2.354845e+02 1.301306e+02 7.971891e+01
#> [886] 6.459716e+01 6.366038e+01 8.038876e+01 6.965270e+01 9.353525e+01
#> [891] 7.359911e+01 9.236067e+01 2.559779e+02 2.566091e+02 6.952524e+01
#> [896] 1.801627e+02 1.105394e+02 9.015558e+01 9.448445e+01 1.001194e+02
#> [901] 8.998466e+01 7.038827e+01 6.613056e+01 6.186742e+01 6.482438e+01
#> [906] 8.167913e+01 8.643001e+01 1.085597e+02 7.335001e+01 1.396356e+02
#> [911] 4.335258e+02 8.366468e+01 7.296484e+01 8.174949e+01 1.340715e+02
#> [916] 1.262500e+02 6.422706e+01 7.242131e+01 6.782525e+01 7.205114e+01
#> [921] 2.370359e+03 7.678517e+01 7.223055e+01 8.650834e+01 1.521336e+02
#> [926] 7.084374e+01 7.315062e+01 9.197118e+01 1.095355e+02 9.123985e+01
#> [931] 6.719559e+01 5.972967e+01 1.140784e+02 1.273098e+02 2.727800e+02
#> [936] 8.202840e+01 1.087311e+02 1.389584e+02 7.032388e+01 1.399602e+02
#> [941] 9.046938e+01 8.985585e+01 8.860348e+01 5.765629e+01 1.356644e+02
#> [946] 1.534744e+02 1.772329e+02 1.773449e+04 6.606185e+01 1.812197e+02
#> [951] 6.594749e+01 2.168474e+02 9.795431e+01 1.348424e+02 1.198979e+02
#> [956] 9.249838e+01 1.100122e+02 1.231466e+02 9.552964e+01 7.668488e+01
#> [961] 9.405453e+01 8.131538e+01 6.965433e+01 8.753658e+01 7.266801e+01
#> [966] 8.079918e+01 8.835043e+01 6.559166e+02 9.776671e+01 7.627110e+01
#> [971] 7.140600e+01 6.712965e+01 8.554244e+01 2.814290e+02 1.827258e+02
#> [976] 6.431341e+01 6.570195e+01 8.890461e+01 7.598419e+01 8.566936e+01
#> [981] 6.995112e+01 1.044533e+02 2.220732e+02 1.499402e+02 1.003811e+02
#> [986] 7.267337e+01 7.248421e+01 6.321886e+01 2.089981e+02 2.468789e+02
#> [991] 7.415663e+01 7.329942e+01 6.616572e+01 6.956440e+01 7.043345e+01
#> [996] 1.409199e+02 1.226911e+02 1.344526e+02 7.105119e+01 6.926599e+01
doseFunction(
model,
alpha2 = samples@data$alpha2,
beta = samples@data$beta,
grade = 2L
)(x = 0.25)
#> [1] 57.05540 62.01914 59.70097 69.48327 63.21585
#> [6] 124.30456 57.90506 56.29301 60.33189 68.37809
#> [11] 55.93087 59.93655 65.46681 82.10900 108.92375
#> [16] 60.92785 55.86868 103.39013 226.04663 62.12191
#> [21] 59.67476 64.69439 77.51612 56.20524 56.10040
#> [26] 68.57742 67.68771 80.37839 74.40353 68.32990
#> [31] 68.36671 58.48372 105.60943 61.33852 205.42455
#> [36] 181.00208 56.38052 69.90108 61.38820 61.46890
#> [41] 66.52352 76.23785 91.68306 62.38282 70.61320
#> [46] 228.49874 61.54678 49930.11497 134.44334 120.22209
#> [51] 102.48024 105.48569 59.08237 58.58665 61.86050
#> [56] 87.57549 58.36716 63.14302 59.35080 59.78849
#> [61] 62.41735 84.22787 59.42155 66.47790 1685.73725
#> [66] 115.81316 64.94502 59.80848 56.54789 15052.39950
#> [71] 59.23666 64.90377 101.86099 57.29748 123.44319
#> [76] 70.57910 69.67229 59.43731 66.95357 144.29660
#> [81] 126.44713 67.92586 84.55807 73.20002 228.65180
#> [86] 68.12717 63.76320 62.62133 76.54783 55.37740
#> [91] 92.18220 74.37393 123.61356 65.54744 55.85148
#> [96] 58.45410 69.36964 59.52916 84.15387 73.72290
#> [101] 62.29945 76.58600 68.78436 55.89575 59.62358
#> [106] 73.62010 127.97643 80.25661 54.75463 52.86771
#> [111] 54.41815 76.59857 69.53080 61.49784 83.63240
#> [116] 58.14978 77.40589 62.94396 78.52351 65.31212
#> [121] 91.16152 69.46741 66.17395 56.03394 65.59785
#> [126] 65.10750 60.79897 58.61871 57.66692 78.65819
#> [131] 64.04243 72.42604 71.93599 65.10118 206.55588
#> [136] 57.21414 60.20127 60.80090 59.28051 73.60182
#> [141] 172.36945 74.31109 118.92084 227.36487 67.02825
#> [146] 74.30572 64.63282 311.67181 81.29588 57.25256
#> [151] 65.12566 65.09702 55.77872 63.94325 62.17483
#> [156] 55.48716 56.07081 364.72051 111.54646 65.50679
#> [161] 99.88824 56.00236 65.01591 80.51543 55.56536
#> [166] 63.02371 65.03408 265.87508 73.07812 96.02056
#> [171] 56.97429 154.93798 67.60605 69.08209 61.67322
#> [176] 69.83343 72.20218 55.68038 67.26843 82.45581
#> [181] 61.65557 109.78664 162.02583 81.59874 57.15872
#> [186] 61.44439 70.70851 69.70577 67.27975 63.76475
#> [191] 63.93274 65.35923 63.91418 109.73913 75.06251
#> [196] 63.46118 63.26780 75.25339 58.85819 69.10862
#> [201] 65.04828 56.84477 57.88785 97.79888 66.47498
#> [206] 55.97845 67.06298 54.31063 57.49836 55.90128
#> [211] 58.89097 183.37066 119.47043 116.68157 88.71228
#> [216] 56.68529 59.13606 58.59920 61.61857 72.21728
#> [221] 1032.24352 1300.74133 54.92310 57.13538 129.50892
#> [226] 76.72837 59.20759 59.39582 57.42185 90.69131
#> [231] 59.22265 84.19687 138.31452 88.45124 74.28728
#> [236] 68.57788 63.37500 64.00311 87.28765 62.52088
#> [241] 63.65854 211.96951 89.58291 123.89026 85.48311
#> [246] 77.54253 56.69969 66.43542 63.51960 92.98195
#> [251] 72.39249 50.26196 58.53601 71.72809 141.06637
#> [256] 96.05612 77.13161 61.83395 76.89025 59.74257
#> [261] 70.44175 102.49233 62.10122 79.19451 140.40898
#> [266] 64.07345 60.47188 69.43995 124.44581 55.86371
#> [271] 60.19050 87.17378 64.65805 65.30007 58.77487
#> [276] 80.29043 67.62029 257.12539 59.76458 67.72510
#> [281] 68.56366 63.81895 55.43621 117.23265 149.26516
#> [286] 115.70735 59.02744 66.27659 70.63304 62.34289
#> [291] 59.66546 62.57331 63.72309 62.22129 56.59271
#> [296] 84.36204 71.00463 18797.99696 107.72592 71.31496
#> [301] 64.51895 56.25817 62.41457 58.42289 56.56929
#> [306] 83.86527 375.31027 59.48230 60.68132 87.11163
#> [311] 188.18511 82.24685 51.65632 80.18371 53.71577
#> [316] 59.08435 63.92476 124.74231 57.01466 65.34251
#> [321] 62.67872 69.30635 173.00171 68.69598 65.90628
#> [326] 54.48212 75.71298 56.16500 92.76980 77.57085
#> [331] 84.31983 86.06823 64.45034 92.59416 299.90593
#> [336] 446.25027 55.87005 66.25877 60.36949 58.13583
#> [341] 57.67141 59.53308 57.84529 59.90597 64.58893
#> [346] 120.95864 59.24986 60.50328 65.03632 58.94739
#> [351] 63.81828 69.38293 60.02522 63.26298 62.08887
#> [356] 58.55711 57.90619 62.19389 886.31636 156.99337
#> [361] 107.19965 476.95438 57.54408 92.57743 102.19599
#> [366] 61.61816 54.34948 68.57535 63.10764 61.16623
#> [371] 71.07254 62.95987 60.45909 62.45195 69.41931
#> [376] 59.57447 60.37949 78.89303 64.18576 54.77950
#> [381] 67.97035 74.50699 64.19579 61.96112 73.43455
#> [386] 68.35060 67.63862 55.46799 68.07997 61.39545
#> [391] 96.47968 56.58803 82.34718 54.02671 74.36163
#> [396] 59.78169 63.93731 59.98779 65.09060 62.93826
#> [401] 66.59720 62.04598 112.08979 55.64262 80.00029
#> [406] 97.27637 73.12065 59.24677 60.27564 70.39197
#> [411] 63.07378 135.18968 66.50840 59.22135 59.33121
#> [416] 114.46771 56.91085 61.01272 59.35657 70.55107
#> [421] 60.59853 222.14740 332.64348 91.91009 73.21196
#> [426] 51.92418 57.43875 388.99093 66.38980 70.82763
#> [431] 63.36192 85.05416 68.28209 121.16845 95.60242
#> [436] 59.56244 193.09898 82.72432 52.92408 566.62488
#> [441] 63.55968 61.62426 56.76532 94.68196 110.16726
#> [446] 57.39405 145.58751 158.90531 61.72134 63.48090
#> [451] 62.63791 148.13429 85.98091 59.73462 73.01087
#> [456] 65.67132 54.13842 43.05945 63.45852 64.30189
#> [461] 134.34880 76.29324 107.15428 76.18609 65.38867
#> [466] 85.29856 72.98508 64.73216 69.14889 67.32538
#> [471] 337.70177 73.43555 127.24428 58.11641 57.98371
#> [476] 61.27160 61.11026 67.05802 90.06926 70.39858
#> [481] 56.02561 96.47302 58.94272 60.68960 82.67751
#> [486] 97.31159 68.18883 95.56677 57.00115 56.25924
#> [491] 58.53824 77.39836 80.33238 59.59423 54.72441
#> [496] 100.25174 56.31186 79.55493 67.82302 347.52583
#> [501] 64.91463 64.34596 59.26836 66.44448 65.16403
#> [506] 75.10588 63.87368 69.21318 76.50517 90.21060
#> [511] 318.61750 204.98681 86.83210 56.15796 86.92759
#> [516] 58.86700 69.56301 67.67538 280.12514 64.69482
#> [521] 62.68063 55.06202 59.76454 114.69453 55.02604
#> [526] 63.72239 66.36419 80.54503 60.74866 48.97052
#> [531] 61.58577 60.89328 72.14719 76.45473 52.39204
#> [536] 60.89642 202.56988 60.52652 61.40768 59.30692
#> [541] 72.62506 107.43223 124.53732 1347.74253 64.85896
#> [546] 77.64086 84.95406 60.01200 60.54615 72.82145
#> [551] 65.04283 55.95486 61.57575 77.09058 169.62913
#> [556] 97.05323 75.94590 48.42896 211.71815 116.30773
#> [561] 63.84633 86.32139 452.12101 105.79955 59.67326
#> [566] 101.17075 87.28673 93.32327 720.69880 63.02476
#> [571] 69.46396 75.16752 61.69392 79.76267 109.02976
#> [576] 58.00651 62.36794 80.67534 57.65228 51.90461
#> [581] 55.06354 54.87066 63.12900 65.59359 57.93484
#> [586] 57.62381 83.12623 142.96501 51.10512 55.48194
#> [591] 58.07188 56.62231 61.11645 70.83363 63.14252
#> [596] 56.66641 69.86622 118.90426 9956.19732 97.12058
#> [601] 72.74871 49.34650 61.99404 55.48552 77.53756
#> [606] 60.57960 79.34048 69.92741 121.04864 59.04208
#> [611] 67.35803 79.02913 112.11719 136.98205 64.32422
#> [616] 56.59503 72.18682 54.27189 54.91009 64.80001
#> [621] 143.57360 79.52799 95.08438 57.48344 63.21456
#> [626] 71.51098 59.43965 59.47938 73.50151 57.47076
#> [631] 257.10850 92.55838 446.29958 60.13219 77.26604
#> [636] 56.89106 64.47540 84.79790 67.24719 53.60305
#> [641] 64.01063 61.61453 368.23970 140.83389 97.69133
#> [646] 71.16052 83.77124 53.45700 63.52552 57.75161
#> [651] 67.51207 62.14367 67.35147 162.57446 61.24639
#> [656] 71.69431 71.72009 128.64850 117.08637 55.68677
#> [661] 68.13678 60.34354 55.48293 155.86851 57.62396
#> [666] 95.67746 58.23330 125.73329 68.40499 74.73141
#> [671] 70.56849 76.35123 64.58468 60.61643 70.33993
#> [676] 83.26184 719.98816 77.70003 62.59825 75.57611
#> [681] 60.44589 76.25519 72.99517 73.79323 95.10264
#> [686] 56.44750 58.73566 517.81438 57.36784 64.23394
#> [691] 84.33627 63.72016 74.60311 52.17545 57.66128
#> [696] 79.19303 52.69297 51.77211 58.45101 75.91708
#> [701] 64.55262 55.29961 64.32134 59.43493 77.59163
#> [706] 74.12274 51.84946 64.97622 272.14714 357.85396
#> [711] 57.84689 61.78659 58.03689 58.85212 59.20059
#> [716] 229.33904 56.90922 62.66325 73.73680 71.86774
#> [721] 68.26473 54.68402 57.86251 55.77585 63.61930
#> [726] 61.71548 59.95984 474.89264 60.55916 84.18226
#> [731] 63.64455 66.99419 58.67549 60.03090 57.63025
#> [736] 83.43715 60.72681 3491.27551 73.22717 58.75138
#> [741] 58.65532 68.85312 59.97490 74.36807 74.64076
#> [746] 86.00697 56.86414 75.10851 57.27906 74.65779
#> [751] 67.63290 68.08433 56.17762 62.16834 56.90031
#> [756] 70.86209 61.92477 59.49171 54.95021 63.60216
#> [761] 59.72436 64.43064 60.30882 56.79391 70.15504
#> [766] 57.20236 71.07285 59.13057 80.28925 89.22530
#> [771] 73.37140 58.71070 55.88590 286.23219 86.87967
#> [776] 77.91593 76.88365 64.56518 152.26907 115.75591
#> [781] 86.50415 61.12366 56.84232 67.30814 60.81253
#> [786] 62.53656 65.46130 63.24944 63.96236 74.97683
#> [791] 86.20082 194.49548 72.50446 57.47765 62.10216
#> [796] 59.39263 3337.23392 65.80480 92.54526 65.84234
#> [801] 79.11054 74.00157 81.60137 84.68956 168.27486
#> [806] 136714.06395 66.20701 60.86919 222.18759 59.71297
#> [811] 155.24995 96.15730 71.94867 202.77096 62.80953
#> [816] 65.22954 63.24212 81.46928 60.23191 54.05734
#> [821] 56.10839 56.80603 54.63955 64.60460 124.85492
#> [826] 187.02229 55.20947 56.43808 70.09341 78.38085
#> [831] 71.90385 208.59664 55.99588 63.92411 62.16742
#> [836] 66.62243 62.54384 107.62102 77.76479 186.39139
#> [841] 75.49887 263.49410 78.19576 209.81073 54.97538
#> [846] 107.75601 149.66802 54.38408 58.83521 58.22070
#> [851] 150.55619 135.30744 83.02524 72.41780 54.89246
#> [856] 415.70074 134.09557 55.18755 77.72801 57.46586
#> [861] 61.46504 61.17708 60.20759 127.63897 62.14311
#> [866] 67.90532 67.92281 58.43708 69.54094 59.87256
#> [871] 63.64648 61.05081 83.98707 68.02033 100.52750
#> [876] 100.85399 3275.13439 62.18356 44.84374 69.31195
#> [881] 4625.12686 55.53790 104.63944 77.52814 59.40358
#> [886] 57.29479 56.35034 64.58131 56.37068 69.66525
#> [891] 57.58874 63.47706 120.21428 52.44860 57.74437
#> [896] 83.04214 59.53324 70.47106 72.32835 69.93763
#> [901] 65.57942 58.78877 56.02073 56.15745 58.98973
#> [906] 66.43599 60.61302 71.77457 52.04192 93.22382
#> [911] 216.39119 72.74974 66.54963 72.06902 76.73707
#> [916] 83.07843 56.06265 60.92065 60.70743 58.43848
#> [921] 231.75943 59.38280 55.46105 63.18345 82.04947
#> [926] 61.64480 58.90770 74.88781 85.05816 70.56220
#> [931] 56.81473 54.20467 81.27298 63.13511 99.08068
#> [936] 59.36346 72.80798 81.45996 58.83279 106.36775
#> [941] 64.82175 67.63484 68.32047 48.71505 71.72655
#> [946] 82.84170 64.52899 509.96987 55.41048 84.15217
#> [951] 58.85054 73.47948 60.62408 58.84093 79.93317
#> [956] 58.43600 72.30568 88.52210 74.04258 60.14418
#> [961] 70.15891 66.13467 61.62581 71.41087 60.28276
#> [966] 61.87608 64.97917 200.80123 66.81941 64.53946
#> [971] 59.83038 56.54937 67.50758 139.28131 125.75029
#> [976] 57.85754 60.37467 75.93556 60.06177 68.43588
#> [981] 61.90718 74.01790 116.57315 78.63102 73.55524
#> [986] 62.93535 62.52797 56.35421 124.86577 92.91045
#> [991] 61.07063 62.86547 59.89587 61.70511 62.39320
#> [996] 74.23043 78.49486 68.30901 60.01006 62.11635