![[Experimental]](figures/lifecycle-experimental.svg)
A function that returns a dose() method that computes the dose reaching a
specific target value of a given independent variable, based on the model
specific parameters.
Usage
doseFunction(model, ...)
# S4 method for class 'GeneralModel'
doseFunction(model, ...)
# S4 method for class 'ModelPseudo'
doseFunction(model, ...)
# S4 method for class 'LogisticLogNormalOrdinal'
doseFunction(model, grade, ...)Value
A dose() method that computes doses.
Functions
doseFunction(GeneralModel): Return a dose function for aGeneralModel.doseFunction(ModelPseudo): Return a dose function for aModelPseudo.doseFunction(LogisticLogNormalOrdinal): Return a grade-specific dose function for aLogisticLogNormalOrdinalmodel.
Examples
my_model <- LogisticLogNormal(
mean = c(-0.85, 1),
cov = matrix(c(1, -0.5, -0.5, 1), nrow = 2),
ref_dose = 50
)
dose_fun <- doseFunction(my_model, alpha0 = 2, alpha1 = 3)
dose_fun(0.6)
#> [1] 29.38579
data_ordinal <- .DefaultDataOrdinal()
model <- .DefaultLogisticLogNormalOrdinal()
options <- .DefaultMcmcOptions()
suppressWarnings({
samples <- mcmc(data_ordinal, model, options)
})
doseFunction(
model,
alpha1 = samples@data$alpha2,
beta = samples@data$beta,
grade = 1L
)(x = 0.75)
#> [1] 7.755929e+01 6.996370e+02 5.725308e+07 1.656955e+02 6.516023e+01
#> [6] 7.413838e+01 1.399318e+02 1.232457e+02 6.803173e+02 7.883601e+01
#> [11] 1.370648e+03 1.773238e+02 2.304939e+02 9.064345e+01 1.021758e+02
#> [16] 3.866093e+02 6.947058e+01 8.400707e+01 1.051445e+02 3.317974e+02
#> [21] 7.137159e+01 1.033091e+02 2.236438e+02 8.016906e+01 1.620204e+02
#> [26] 1.210168e+03 6.549486e+01 7.601734e+01 8.547325e+01 6.067081e+01
#> [31] 9.050137e+01 8.715702e+01 8.293274e+01 6.700995e+01 8.139737e+01
#> [36] 1.636348e+02 1.202443e+02 8.155723e+02 1.252209e+02 6.933839e+01
#> [41] 9.748678e+01 9.092237e+01 8.252676e+02 1.004513e+02 8.335494e+01
#> [46] 8.753762e+01 7.688506e+01 3.339215e+02 7.100382e+01 1.012440e+02
#> [51] 6.384330e+01 1.812466e+02 8.257373e+01 1.124074e+02 1.207791e+02
#> [56] 1.113132e+02 9.501156e+01 9.365412e+01 8.474314e+01 6.689317e+01
#> [61] 1.774184e+02 9.072007e+01 7.273457e+01 6.737107e+01 7.691422e+01
#> [66] 1.029761e+02 5.906697e+01 1.135336e+02 2.209807e+02 6.663264e+01
#> [71] 4.490817e+02 6.513204e+02 1.730141e+02 8.722845e+01 1.146514e+02
#> [76] 1.016102e+02 8.871598e+01 6.585374e+01 1.593122e+02 1.114516e+02
#> [81] 1.017666e+02 6.835733e+01 7.113511e+01 8.188902e+01 8.354113e+01
#> [86] 5.524529e+02 1.129563e+02 6.982860e+01 7.636226e+01 8.951482e+01
#> [91] 8.150722e+01 8.320432e+01 6.836526e+01 6.828259e+01 6.523702e+01
#> [96] 6.947797e+01 7.546774e+01 8.138368e+01 1.514030e+02 7.328154e+01
#> [101] 2.849101e+02 1.100629e+02 8.412475e+01 1.442212e+02 8.249565e+01
#> [106] 7.469390e+01 6.471710e+01 1.250827e+02 9.566486e+01 7.632170e+01
#> [111] 1.146996e+02 7.755109e+01 7.726490e+01 6.381607e+01 8.346619e+01
#> [116] 9.283981e+01 1.620270e+02 9.774633e+01 1.106650e+02 6.671115e+01
#> [121] 1.013247e+02 6.627484e+01 7.764271e+01 7.142917e+01 1.082966e+02
#> [126] 8.355200e+01 9.805848e+01 7.061525e+01 9.181940e+01 9.595875e+01
#> [131] 8.436190e+01 1.116569e+02 7.556749e+01 7.608308e+01 1.130120e+02
#> [136] 9.776162e+01 3.499458e+02 1.745952e+02 7.281616e+01 6.654850e+01
#> [141] 1.183672e+02 1.803716e+03 7.935052e+02 1.778837e+03 3.277340e+02
#> [146] 7.234828e+01 7.190692e+01 7.553593e+01 6.924013e+01 8.487983e+01
#> [151] 8.669880e+01 1.322947e+02 8.328140e+01 1.074778e+02 7.624564e+01
#> [156] 9.376449e+01 8.378032e+01 2.051144e+02 1.192229e+04 7.686878e+01
#> [161] 1.241328e+02 6.405424e+01 1.736390e+02 6.801449e+01 6.293235e+01
#> [166] 1.475303e+02 1.141615e+02 6.896697e+01 8.041630e+01 7.658224e+01
#> [171] 7.907541e+01 6.798558e+01 1.383558e+02 3.692265e+02 9.196066e+01
#> [176] 7.129977e+01 8.660812e+01 6.581340e+01 6.979650e+01 7.236611e+01
#> [181] 2.649987e+02 8.396851e+01 1.222795e+02 7.327662e+01 1.070103e+02
#> [186] 1.229588e+02 1.194704e+03 8.416287e+01 1.917775e+02 8.609099e+01
#> [191] 8.367619e+01 6.684816e+01 1.645296e+02 7.213886e+01 2.905361e+02
#> [196] 9.187445e+01 6.804082e+01 3.692423e+02 9.353390e+01 2.365328e+04
#> [201] 2.030722e+04 4.978420e+03 1.378615e+03 7.326583e+01 7.053973e+01
#> [206] 8.917509e+01 9.795325e+01 9.894596e+01 8.869887e+01 6.665540e+01
#> [211] 6.834271e+01 9.575927e+01 1.768444e+02 9.815822e+01 2.470363e+02
#> [216] 8.834152e+01 6.675450e+01 6.206838e+01 2.190126e+02 1.440922e+02
#> [221] 1.035323e+02 8.574565e+01 6.842046e+01 8.168974e+01 1.502256e+02
#> [226] 1.338577e+03 1.955385e+02 7.518432e+01 7.445417e+01 1.213906e+02
#> [231] 8.090722e+01 7.214108e+01 6.724983e+01 7.220886e+01 7.538896e+01
#> [236] 7.319441e+01 7.602084e+01 6.372340e+01 7.515266e+01 8.419992e+01
#> [241] 7.997951e+01 1.259370e+02 7.426054e+01 3.223419e+02 7.053902e+01
#> [246] 7.284319e+01 3.522945e+02 1.795419e+02 9.330929e+01 7.454422e+01
#> [251] 7.654890e+01 8.087088e+01 7.138944e+01 6.874303e+01 1.479368e+02
#> [256] 8.073198e+01 1.238766e+02 8.594888e+01 7.824154e+01 1.135797e+02
#> [261] 1.256255e+02 7.131189e+01 6.753816e+01 9.138320e+01 1.042689e+02
#> [266] 1.978399e+02 1.075138e+02 1.963463e+02 2.348412e+02 8.298067e+01
#> [271] 1.460940e+02 9.406053e+01 8.118789e+01 6.988968e+01 6.834992e+01
#> [276] 6.837035e+01 7.016318e+01 8.554923e+01 1.155906e+02 1.057621e+02
#> [281] 1.184886e+02 7.326422e+01 7.837215e+01 1.898922e+02 1.109165e+02
#> [286] 1.218145e+02 2.172679e+03 1.025801e+02 2.097489e+02 1.202113e+02
#> [291] 1.462806e+02 7.853246e+01 7.336539e+01 6.875637e+01 1.049983e+02
#> [296] 7.717876e+01 8.349701e+01 8.634228e+01 9.172397e+01 1.672268e+02
#> [301] 7.893117e+01 1.003499e+02 7.311746e+01 9.628949e+01 7.697100e+07
#> [306] 6.555496e+02 4.145140e+04 2.392423e+03 9.465724e+01 2.459552e+02
#> [311] 4.554903e+02 1.187715e+03 6.344320e+01 7.850786e+01 1.598158e+02
#> [316] 2.351589e+03 4.139764e+03 7.228559e+01 6.854639e+01 6.884527e+01
#> [321] 7.031511e+01 1.351779e+02 1.097230e+02 6.750111e+01 1.566312e+02
#> [326] 1.236591e+02 6.590697e+01 6.243238e+01 7.540668e+01 7.775888e+01
#> [331] 6.362364e+01 6.597962e+01 2.383163e+02 1.394399e+02 6.710068e+01
#> [336] 7.903786e+01 2.149727e+03 3.630421e+02 8.780885e+01 1.526464e+02
#> [341] 1.376035e+02 1.266077e+02 1.445760e+02 2.294345e+02 1.897117e+02
#> [346] 1.276311e+02 8.955668e+01 1.071632e+02 9.058861e+01 7.001957e+01
#> [351] 1.794627e+02 8.187953e+01 3.409076e+02 3.194887e+02 1.308857e+02
#> [356] 6.652870e+01 1.115219e+02 1.073888e+02 6.333571e+01 1.476781e+02
#> [361] 1.361811e+02 1.161108e+02 8.781503e+01 4.040696e+02 2.675634e+02
#> [366] 7.941096e+01 8.303467e+01 1.340302e+02 8.362516e+01 6.721087e+01
#> [371] 6.083330e+01 8.536053e+01 1.191110e+04 8.115799e+01 6.750281e+01
#> [376] 1.365311e+02 6.549052e+01 8.078836e+01 9.565361e+01 7.385114e+01
#> [381] 6.227369e+01 1.169840e+03 6.803339e+01 8.281766e+01 1.799354e+03
#> [386] 1.246543e+02 1.014408e+02 8.689642e+01 7.799935e+01 1.491838e+02
#> [391] 1.314447e+02 3.540004e+02 1.031419e+02 2.140600e+02 1.388510e+02
#> [396] 1.205604e+03 2.230839e+02 7.233689e+01 1.949476e+03 4.139861e+02
#> [401] 9.553839e+01 1.508024e+02 9.391840e+01 6.983859e+02 8.206557e+01
#> [406] 1.846597e+02 1.285562e+02 8.888195e+01 1.105748e+02 1.331365e+02
#> [411] 7.371614e+01 6.955332e+01 7.246365e+01 7.859235e+01 8.780115e+01
#> [416] 3.266885e+02 7.029390e+01 6.576433e+01 7.906939e+01 7.517689e+01
#> [421] 3.667748e+02 2.356238e+02 3.084150e+03 7.268920e+01 1.453440e+02
#> [426] 1.432684e+02 8.756453e+01 8.099217e+02 1.089898e+02 1.941580e+02
#> [431] 1.747727e+02 1.053358e+03 1.046343e+02 6.455757e+01 1.005118e+02
#> [436] 7.041195e+01 1.164379e+02 9.740773e+01 8.470093e+01 1.731405e+02
#> [441] 6.340388e+01 6.830475e+01 1.414558e+02 7.016861e+01 2.394718e+03
#> [446] 9.427140e+02 1.170855e+02 9.293176e+01 8.636542e+01 1.552382e+02
#> [451] 3.356919e+05 6.751948e+01 7.062864e+01 3.450254e+02 7.033036e+01
#> [456] 4.477983e+02 7.577009e+02 6.695810e+01 7.471243e+01 9.878970e+01
#> [461] 8.435523e+01 7.155932e+01 2.224048e+02 1.610698e+02 8.909253e+01
#> [466] 1.999241e+02 7.903558e+01 1.023492e+02 2.047560e+02 1.082527e+02
#> [471] 9.175779e+01 1.218219e+02 7.936741e+01 8.034614e+01 6.952940e+01
#> [476] 1.115167e+03 4.199057e+03 6.722594e+01 1.541507e+02 2.191173e+02
#> [481] 1.126665e+02 1.804034e+02 1.763969e+02 7.669445e+01 6.052352e+01
#> [486] 4.258255e+02 6.744054e+01 7.037438e+01 2.423970e+02 6.023450e+01
#> [491] 7.653117e+01 7.660903e+01 9.106736e+01 6.459652e+01 9.599337e+01
#> [496] 2.741206e+02 1.469775e+02 1.034920e+02 1.107036e+02 6.267684e+01
#> [501] 7.586496e+01 7.053574e+01 6.902830e+01 6.590563e+01 8.876293e+01
#> [506] 7.230661e+01 1.709936e+02 1.263481e+02 1.660910e+02 7.084792e+01
#> [511] 6.433036e+01 2.599510e+02 1.024089e+02 9.996307e+01 3.370894e+02
#> [516] 3.178988e+02 7.906013e+01 1.215280e+02 2.963286e+02 6.612693e+01
#> [521] 1.500469e+02 1.409612e+02 7.510564e+01 6.988018e+01 7.046510e+01
#> [526] 7.492049e+01 7.048705e+01 8.590468e+01 7.085488e+01 8.544675e+01
#> [531] 6.598726e+01 8.238883e+01 7.475416e+01 1.069779e+02 8.320145e+01
#> [536] 8.633510e+01 7.745097e+01 1.090403e+02 9.845573e+01 1.612207e+02
#> [541] 8.597568e+01 6.971924e+01 1.236380e+02 1.683861e+02 8.150622e+01
#> [546] 7.386158e+01 9.182422e+02 8.840544e+01 9.489147e+01 1.781787e+02
#> [551] 9.243473e+01 6.803807e+01 7.093955e+01 9.164124e+01 6.787640e+02
#> [556] 2.141217e+02 8.675743e+01 9.084848e+01 1.090290e+02 7.523221e+01
#> [561] 2.511509e+02 6.574240e+01 1.023467e+02 5.770973e+02 8.183855e+01
#> [566] 1.107769e+02 7.350361e+01 6.505074e+01 8.447632e+01 7.211885e+01
#> [571] 6.423208e+01 8.804785e+01 1.514036e+02 8.393495e+01 9.051838e+01
#> [576] 8.976464e+01 7.400180e+01 1.035880e+02 2.771819e+02 8.173340e+01
#> [581] 4.444460e+02 7.754323e+01 7.004366e+01 7.872994e+01 1.674390e+02
#> [586] 7.656699e+01 7.201414e+01 7.442343e+01 3.285474e+02 1.283454e+02
#> [591] 1.478539e+02 2.043813e+02 6.439155e+01 7.574430e+01 4.202622e+02
#> [596] 7.658141e+01 7.322757e+01 9.719749e+01 1.049263e+02 7.569105e+01
#> [601] 7.093981e+01 1.128788e+02 8.578600e+01 3.706361e+02 9.965616e+01
#> [606] 8.425333e+01 6.842939e+01 6.600601e+01 1.456629e+02 2.459820e+02
#> [611] 9.315049e+01 8.460834e+01 1.414815e+02 1.364857e+02 9.780004e+01
#> [616] 1.347868e+02 8.819658e+01 7.188884e+01 8.560393e+01 2.098256e+02
#> [621] 1.205564e+03 1.121074e+02 1.419368e+03 2.378835e+02 1.214426e+02
#> [626] 7.401344e+01 5.441387e+02 7.027863e+01 8.938394e+01 1.788554e+02
#> [631] 6.639942e+01 6.201647e+01 6.814495e+01 2.857483e+02 5.295877e+02
#> [636] 1.622763e+02 1.139645e+02 9.949806e+01 7.507269e+01 7.585440e+01
#> [641] 9.851310e+01 1.300587e+02 8.702431e+01 9.828091e+01 4.130512e+03
#> [646] 1.171964e+02 1.080272e+02 8.611161e+01 8.315485e+01 7.066955e+01
#> [651] 8.532149e+01 8.855927e+01 1.463237e+02 9.144930e+01 1.315272e+02
#> [656] 8.489015e+01 1.164365e+02 6.956076e+01 6.795010e+01 8.428981e+01
#> [661] 1.618600e+02 7.701195e+01 9.687204e+01 3.486812e+02 8.427282e+01
#> [666] 3.446460e+02 3.464456e+03 4.226655e+03 1.147776e+02 8.521076e+01
#> [671] 1.442221e+02 1.180480e+02 2.713282e+02 1.861624e+02 8.379403e+01
#> [676] 1.916216e+02 3.949343e+06 6.531020e+01 1.653194e+02 1.491308e+02
#> [681] 7.763166e+01 1.115359e+02 2.286153e+02 1.061211e+02 1.652328e+02
#> [686] 4.779375e+04 2.928474e+02 1.406120e+02 1.004763e+02 8.737194e+01
#> [691] 7.398795e+01 2.118268e+04 6.692493e+01 1.301198e+02 1.029233e+02
#> [696] 7.970165e+01 6.513771e+01 9.065489e+01 1.209128e+02 1.917026e+02
#> [701] 7.658017e+01 7.318371e+01 1.089437e+02 7.872798e+01 9.556935e+01
#> [706] 1.165408e+02 7.436587e+01 7.128896e+01 8.517396e+01 1.168621e+02
#> [711] 9.935459e+01 1.754764e+02 1.166812e+02 1.473287e+02 3.184695e+02
#> [716] 8.286542e+01 1.062720e+02 6.697096e+01 9.760822e+01 1.971763e+04
#> [721] 7.463506e+01 6.574144e+01 8.264064e+02 9.833273e+01 2.298474e+02
#> [726] 1.200299e+02 1.543941e+02 6.649189e+01 1.221380e+02 7.064580e+01
#> [731] 7.220738e+01 9.407627e+01 7.545549e+01 6.466065e+01 6.620126e+01
#> [736] 8.963985e+01 1.000585e+02 6.364855e+01 1.021786e+02 1.213981e+02
#> [741] 3.280906e+02 1.494220e+02 7.411276e+01 7.655835e+01 8.426229e+01
#> [746] 3.732170e+02 1.099601e+02 1.370087e+02 6.375830e+01 6.965612e+01
#> [751] 7.856772e+01 1.095848e+02 8.113121e+01 1.414320e+02 3.382107e+02
#> [756] 6.794917e+01 6.922956e+01 1.077964e+02 1.194705e+02 1.192266e+02
#> [761] 6.511631e+01 8.226911e+01 8.686823e+01 8.538254e+01 1.192338e+02
#> [766] 1.285820e+02 1.334451e+02 6.892645e+01 6.158930e+01 7.089579e+01
#> [771] 1.669910e+02 6.643389e+01 1.514538e+02 1.558725e+02 2.281694e+02
#> [776] 6.343399e+01 7.824056e+01 6.759664e+01 1.089919e+02 7.376496e+01
#> [781] 1.062669e+02 1.218265e+02 1.027291e+02 1.098158e+02 9.254623e+01
#> [786] 7.470570e+01 7.569800e+01 1.085491e+02 6.982428e+01 7.970768e+01
#> [791] 3.675670e+02 1.317406e+02 8.274503e+01 6.927386e+01 6.657868e+01
#> [796] 7.039529e+01 7.390326e+01 7.513820e+01 1.682906e+02 9.562371e+01
#> [801] 1.109689e+02 6.143151e+01 6.685928e+01 7.516752e+01 3.198047e+02
#> [806] 4.222747e+02 7.821107e+01 8.734598e+01 1.458767e+02 9.745860e+01
#> [811] 7.167319e+01 8.004871e+01 8.051489e+01 8.046090e+01 7.938592e+01
#> [816] 8.226903e+01 9.533458e+01 1.401321e+02 1.343866e+02 2.406958e+03
#> [821] 1.083950e+02 9.636402e+01 1.773844e+02 9.167947e+01 6.508319e+01
#> [826] 7.391932e+01 1.219206e+02 7.457024e+01 9.990312e+01 1.466385e+02
#> [831] 7.829133e+01 9.449687e+01 2.541342e+02 8.019315e+01 7.549543e+01
#> [836] 1.273666e+02 1.321362e+02 9.281506e+01 2.498413e+02 5.828845e+01
#> [841] 6.337453e+01 8.271950e+01 8.888069e+01 1.148932e+02 1.529819e+03
#> [846] 9.698089e+01 9.467458e+01 1.661179e+02 1.424656e+02 7.475519e+01
#> [851] 6.895912e+01 1.251548e+03 3.444020e+02 6.728956e+01 9.150103e+01
#> [856] 6.763872e+01 2.092820e+02 9.744233e+01 9.819795e+01 7.188185e+01
#> [861] 7.889500e+01 2.946701e+02 1.390942e+02 8.855860e+01 6.836339e+03
#> [866] 6.696240e+01 8.172564e+01 1.004473e+02 1.143542e+02 6.635023e+01
#> [871] 6.431310e+01 7.156088e+01 8.188743e+01 7.898918e+01 8.659964e+01
#> [876] 8.734315e+01 6.907560e+01 6.325830e+01 6.549028e+01 1.006242e+02
#> [881] 6.985731e+01 1.126351e+02 7.187074e+01 7.877607e+01 7.300089e+01
#> [886] 1.292140e+02 1.158691e+02 6.983663e+01 1.280581e+02 8.061487e+01
#> [891] 5.125817e+02 7.857530e+01 9.917526e+01 8.428676e+01 6.901463e+01
#> [896] 6.448278e+01 6.773594e+01 6.800693e+03 2.315101e+02 9.340853e+01
#> [901] 7.616004e+01 1.325923e+02 1.490825e+02 7.158037e+01 8.561623e+01
#> [906] 6.615932e+01 3.805403e+02 3.354200e+02 8.640045e+01 1.026874e+02
#> [911] 7.362458e+01 3.300987e+02 6.544537e+01 1.465379e+02 4.264985e+02
#> [916] 4.866378e+03 2.552898e+02 1.127403e+02 7.626350e+01 7.631138e+01
#> [921] 7.164222e+01 6.987219e+01 6.420044e+01 2.223530e+02 6.614370e+01
#> [926] 6.262125e+01 1.119211e+02 1.385875e+02 1.406495e+02 1.261383e+02
#> [931] 6.992942e+01 7.881365e+01 7.259678e+01 7.294228e+01 9.119735e+01
#> [936] 1.903487e+02 1.025080e+02 8.396589e+01 1.663403e+02 2.604709e+02
#> [941] 8.093334e+01 5.867309e+02 9.110362e+01 2.438154e+02 8.135089e+01
#> [946] 7.045704e+01 7.688336e+01 3.825930e+02 1.234037e+02 1.278122e+02
#> [951] 9.066063e+01 1.208082e+02 1.613748e+02 1.150282e+02 7.404398e+01
#> [956] 7.660897e+01 1.317790e+02 4.168683e+02 1.477982e+02 1.237979e+02
#> [961] 2.066462e+02 6.936476e+01 8.753774e+01 7.625271e+01 8.824375e+01
#> [966] 1.065450e+02 9.781960e+01 9.748276e+01 7.493885e+01 9.177409e+01
#> [971] 2.055489e+02 6.867878e+01 1.001294e+02 1.132144e+02 7.991484e+01
#> [976] 2.300093e+02 7.675597e+01 4.511941e+02 1.405548e+02 2.419194e+02
#> [981] 2.198016e+04 9.768872e+03 7.875334e+01 9.452083e+01 9.413337e+01
#> [986] 1.186647e+02 1.352578e+02 1.606597e+02 7.742159e+01 8.876620e+01
#> [991] 8.218281e+01 2.018662e+03 1.019089e+02 2.278748e+02 1.105304e+02
#> [996] 9.690659e+01 6.800423e+01 2.581921e+02 6.628763e+01 1.411874e+02
doseFunction(
model,
alpha2 = samples@data$alpha2,
beta = samples@data$beta,
grade = 2L
)(x = 0.25)
#> [1] 63.89896 174.24906 163863.88563 71.94701 52.89997
#> [6] 60.98296 86.67584 88.62765 283.68881 67.96382
#> [11] 278.29012 75.47513 90.98739 66.58999 68.67148
#> [16] 153.75195 50.28664 64.99469 74.84494 98.97790
#> [21] 62.89304 81.44537 137.89727 40.83129 111.20471
#> [26] 155.09248 51.82313 61.03807 65.53926 54.37584
#> [31] 66.61880 65.68705 55.61158 59.14873 68.44656
#> [36] 100.01343 75.71642 186.85557 65.28151 59.01128
#> [41] 63.81576 57.91227 228.99817 75.28678 66.46374
#> [46] 64.18515 63.52678 158.63272 60.69847 70.15537
#> [51] 49.49539 81.42872 69.29350 63.24064 63.92814
#> [56] 71.54206 61.54979 74.65318 65.78086 58.83697
#> [61] 121.07742 71.83108 60.34507 60.00907 64.00783
#> [66] 76.87324 53.40675 78.13177 93.38742 56.18710
#> [71] 174.04150 213.77610 79.66901 61.25446 68.38921
#> [76] 71.07499 67.76092 58.17411 107.02144 66.90066
#> [81] 70.23369 56.33026 62.56101 67.66177 56.24303
#> [86] 143.29788 74.07303 61.50647 62.26484 73.20831
#> [91] 66.94172 69.17521 60.89066 62.05349 57.98787
#> [96] 58.99231 64.13596 67.31740 102.36054 61.86516
#> [101] 144.27684 74.82532 71.02847 103.42117 70.49653
#> [106] 60.09707 48.76112 82.37477 49.85716 62.75928
#> [111] 74.61122 62.76754 64.48059 54.75217 65.75686
#> [116] 62.17438 52.98926 75.39370 82.08731 49.78876
#> [121] 63.69372 56.16799 60.24018 58.02412 63.60076
#> [126] 62.14304 77.52603 62.60333 65.36303 63.28989
#> [131] 65.53476 84.31689 54.05933 63.68656 81.45838
#> [136] 68.42366 169.99473 95.06382 61.30338 58.68494
#> [141] 66.47036 269.57255 257.55600 313.17689 88.64078
#> [146] 57.86534 59.15402 63.16823 59.80182 74.06548
#> [151] 63.51534 90.62574 64.02477 78.30860 60.39034
#> [156] 67.31697 62.34022 111.28390 1463.23476 60.59816
#> [161] 70.92542 56.30552 81.77555 60.45563 57.52361
#> [166] 97.99436 63.93438 52.62450 67.68081 63.28463
#> [171] 66.91305 59.00098 63.01044 135.50155 61.10415
#> [176] 58.44565 71.65262 58.78301 60.55589 64.04509
#> [181] 136.31238 72.15468 65.56500 59.58932 70.97782
#> [186] 84.22404 264.75818 59.99997 106.72849 68.60115
#> [191] 57.10709 57.64286 108.60006 59.03343 113.07375
#> [196] 65.88548 59.37824 154.49987 66.14416 132.73551
#> [201] 116.16743 568.25723 101.34297 62.05878 61.61352
#> [206] 71.66059 77.38197 76.69463 67.50043 58.06134
#> [211] 57.97799 69.16157 63.15257 77.00760 99.72982
#> [216] 58.43187 54.57698 53.68729 116.84357 82.68464
#> [221] 69.03250 62.99311 59.37925 67.68130 91.28472
#> [226] 160.54862 93.50423 59.63715 49.82075 70.32338
#> [231] 61.80765 61.45041 57.08304 55.89916 61.45083
#> [236] 60.69840 64.11506 52.24889 56.23567 60.84464
#> [241] 58.53784 66.52921 60.45072 100.95211 62.80027
#> [246] 65.16561 169.76916 113.39164 70.49630 62.15962
#> [251] 63.13994 63.30219 63.33730 60.82014 65.61041
#> [256] 65.45563 83.86501 60.66861 60.58684 81.94454
#> [261] 83.17024 58.61423 56.96466 66.02360 68.17364
#> [266] 117.06510 51.61196 123.22997 124.13911 61.65597
#> [271] 70.30878 62.17647 68.94175 60.46837 59.18788
#> [276] 59.84297 61.43971 67.89256 82.07451 78.57788
#> [281] 83.48054 65.09595 66.92240 102.58287 69.91987
#> [286] 72.19679 451.50080 73.90554 95.78469 81.47532
#> [291] 88.85175 60.52168 60.82990 56.05720 74.54016
#> [296] 59.99189 63.85506 68.29583 68.70423 96.11396
#> [301] 61.62831 77.92275 58.26422 74.50562 1498.85982
#> [306] 123.24129 3059.04588 581.56132 73.33655 136.78386
#> [311] 173.03430 286.89032 57.21921 66.25882 115.34347
#> [316] 223.36277 471.28700 62.23965 60.16762 60.84261
#> [321] 60.78310 93.21713 65.50023 54.25428 96.55847
#> [326] 62.14559 57.98266 56.37270 65.14944 66.18884
#> [331] 57.96291 57.27339 128.82242 84.20645 52.33314
#> [336] 66.75133 266.71940 178.95709 53.65843 85.19112
#> [341] 93.39907 90.27214 72.66513 116.07814 86.94807
#> [346] 75.82760 63.91725 74.21965 63.17717 61.43533
#> [351] 111.20960 66.04445 96.93042 78.25386 97.06861
#> [356] 60.15189 88.28524 84.56621 56.58844 93.21373
#> [361] 91.68885 83.66538 61.31097 102.49335 91.30674
#> [366] 61.39331 59.30370 82.29952 67.87035 55.26662
#> [371] 53.89112 64.92656 1002.78546 60.79984 56.46447
#> [376] 97.45284 58.23881 64.31660 74.40264 63.39633
#> [381] 56.09532 444.88998 60.22230 65.58096 264.95932
#> [386] 78.00172 72.19799 69.83943 62.55052 82.57935
#> [391] 86.29926 84.76517 71.23671 77.18600 72.09672
#> [396] 327.40166 100.00939 54.49477 240.41372 76.22544
#> [401] 57.63121 105.86710 70.62579 173.81779 55.29615
#> [406] 88.41380 91.56128 69.30763 65.83073 80.50501
#> [411] 56.15446 58.31050 61.37133 65.71561 65.58340
#> [416] 137.04492 60.67273 58.83865 65.45069 65.30268
#> [421] 137.71139 114.52800 400.84285 60.57938 62.40469
#> [426] 82.53761 69.66102 224.73094 59.59409 75.69931
#> [431] 71.34875 167.94636 57.79603 51.26996 61.18752
#> [436] 53.96390 76.31049 61.20310 63.94584 86.58556
#> [441] 55.48423 55.70467 72.76910 59.18742 322.93124
#> [446] 286.95079 71.67245 71.79013 65.46412 85.62052
#> [451] 9623.35483 61.04438 63.96386 191.92491 63.58966
#> [456] 241.63800 352.35169 60.59744 67.20424 81.61257
#> [461] 69.26462 59.88396 95.14976 87.62157 57.70947
#> [466] 96.15126 55.15401 74.83395 125.30095 71.68677
#> [471] 74.15615 93.32956 66.68377 66.02606 60.26316
#> [476] 353.26373 601.72474 57.60981 102.54466 81.08459
#> [481] 56.26593 77.02823 76.91031 53.67242 55.07873
#> [486] 102.12167 59.36870 61.02949 103.42584 53.63168
#> [491] 66.48559 63.16620 54.40289 56.82159 75.08147
#> [496] 152.78519 75.28265 52.45364 75.12394 52.73065
#> [501] 64.62354 62.46780 60.49370 58.62069 68.82094
#> [506] 56.03184 110.67354 70.36833 91.64405 59.49300
#> [511] 58.03383 132.96306 69.45493 74.20300 92.74482
#> [516] 114.78625 63.02390 66.77956 110.90367 57.86273
#> [521] 96.73705 75.94209 64.68162 62.39778 61.01213
#> [526] 65.18858 61.85807 72.85714 60.57989 63.57865
#> [531] 58.47400 64.07915 57.70406 83.62522 70.69056
#> [536] 71.67871 61.51639 60.88511 63.21938 77.25586
#> [541] 63.45310 59.91249 85.88801 80.23021 61.37185
#> [546] 62.72467 234.92041 61.78721 71.40672 115.57743
#> [551] 65.06691 59.27248 58.77988 59.08147 123.12775
#> [556] 122.28891 68.99417 66.65226 81.37282 64.57031
#> [561] 85.84308 57.82756 72.27346 291.53635 62.27966
#> [566] 78.09866 58.14515 58.49715 70.14155 63.50281
#> [571] 58.42259 72.51213 101.07231 69.87652 73.47404
#> [576] 66.97622 59.41921 82.14862 107.72076 58.34827
#> [581] 117.10338 64.32094 60.24232 68.19862 80.18760
#> [586] 50.89370 52.50307 62.56300 107.59292 71.42807
#> [591] 66.14367 117.01866 57.14387 59.49689 142.52051
#> [596] 58.33293 62.97938 68.52528 67.51647 58.79972
#> [601] 62.78962 74.19316 67.19049 155.20017 72.36722
#> [606] 62.33834 58.75465 58.78325 105.10629 116.03798
#> [611] 61.59025 66.07250 91.54495 67.79743 68.82473
#> [616] 65.55534 56.12295 57.84804 62.60128 93.47733
#> [621] 301.18257 70.92856 280.98228 65.79982 79.75872
#> [626] 62.49948 136.78211 61.57569 64.38981 60.26505
#> [631] 58.63066 55.95237 60.40809 105.90589 163.68127
#> [636] 74.25075 80.57978 78.64138 64.88424 64.26029
#> [641] 71.22865 87.58898 71.49576 76.30292 650.76847
#> [646] 69.20087 78.96558 62.30876 61.70629 59.65179
#> [651] 62.77747 59.53944 61.12183 65.93730 76.08276
#> [656] 63.35989 75.91663 59.48934 60.06556 65.20601
#> [661] 79.62813 57.76707 74.90424 87.32378 61.50912
#> [666] 62.98859 578.97663 593.50277 79.01493 62.48350
#> [671] 59.94420 67.42674 139.70244 94.36937 60.86810
#> [676] 110.89542 33054.70927 54.81024 75.82426 84.67893
#> [681] 60.14343 78.65753 104.04445 67.23016 55.43261
#> [686] 8013.65967 179.22664 65.60297 76.70036 65.24542
#> [691] 58.73546 1250.18847 58.01515 92.51090 81.89897
#> [696] 67.96567 57.15180 67.80573 77.69297 95.90516
#> [701] 65.70262 62.67472 87.69204 67.56572 73.69460
#> [706] 87.58696 64.96437 60.27850 63.20433 82.01376
#> [711] 63.89736 66.87773 61.38964 77.37128 110.81595
#> [716] 62.35832 77.53756 52.70144 73.27743 2614.30554
#> [721] 63.56306 54.97263 183.19993 63.64813 92.71768
#> [726] 72.61617 76.28839 56.72029 77.06961 60.26194
#> [731] 58.21425 73.84613 64.87995 58.83692 60.38742
#> [736] 68.92175 80.74810 54.65759 75.67834 60.42278
#> [741] 92.02048 87.08998 61.69543 61.20593 66.00293
#> [746] 154.37218 72.65560 87.48127 55.28296 57.18751
#> [751] 62.86344 80.47206 63.47390 51.98019 105.62756
#> [756] 60.14770 58.71380 70.11974 62.74165 83.57940
#> [761] 58.02013 59.06578 62.65682 63.83979 66.71336
#> [766] 71.10770 87.11563 59.05407 56.00874 62.32302
#> [771] 98.61869 59.03919 84.16608 84.95658 81.46010
#> [776] 54.33138 60.15959 54.86831 63.77952 64.76113
#> [781] 68.76033 79.50022 74.64279 72.37811 68.47286
#> [786] 63.13906 50.80345 49.49731 58.73216 56.13312
#> [791] 124.17011 61.85762 55.28602 55.64297 57.72050
#> [796] 61.36595 62.46729 63.50426 112.81943 72.62151
#> [801] 86.68788 57.74968 59.48506 60.09364 128.02218
#> [806] 119.53037 66.21640 70.68509 87.10518 64.71308
#> [811] 60.24805 63.02259 62.80395 68.04117 63.34444
#> [816] 64.96013 65.54930 74.55824 76.02192 361.78353
#> [821] 85.13272 66.21503 82.95636 68.99869 59.30751
#> [826] 64.83664 77.42253 61.28406 75.33853 95.69471
#> [831] 62.54441 68.69195 131.19648 67.19247 62.29676
#> [836] 90.79246 81.53409 57.77147 146.77568 52.03787
#> [841] 57.00745 65.87370 64.52254 81.51251 287.80047
#> [846] 69.71518 56.90190 59.37275 81.83413 62.02557
#> [851] 60.30309 365.03515 86.78976 59.75594 70.98408
#> [856] 57.39446 65.32238 65.73857 75.36144 61.95552
#> [861] 55.83511 121.51930 67.19066 59.95380 1027.22214
#> [866] 60.46760 70.31253 62.06019 73.09653 56.93016
#> [871] 55.41151 63.40315 64.17158 62.44109 68.89003
#> [876] 71.40604 61.59248 56.72087 56.54516 75.80514
#> [881] 55.52159 70.85621 57.81020 64.66845 59.82490
#> [886] 77.26124 77.25207 61.92477 91.29541 61.43036
#> [891] 183.72055 50.61539 63.44642 63.79186 59.46249
#> [896] 57.43761 58.98847 974.31737 127.38660 62.35782
#> [901] 58.68751 92.05407 51.69294 60.80612 68.00178
#> [906] 58.24697 155.13365 125.04461 67.29909 78.27661
#> [911] 64.06023 149.11744 54.85703 87.93597 116.20076
#> [916] 734.34837 114.24182 72.70319 65.73175 66.62521
#> [921] 63.40109 60.40566 58.79324 112.81455 59.74196
#> [926] 54.42001 77.96819 79.80809 69.84513 82.53805
#> [931] 59.98809 62.50240 55.87436 55.85885 66.96288
#> [936] 116.90327 70.24252 64.33235 86.95875 101.44520
#> [941] 64.45808 258.57817 74.42779 78.33801 68.72363
#> [946] 62.26344 59.51922 84.31809 69.90434 78.40336
#> [951] 66.75043 81.61685 82.89873 73.86272 59.73659
#> [956] 57.21084 72.94773 109.24869 76.73530 90.14851
#> [961] 99.81726 58.86385 71.74132 62.25126 63.39200
#> [966] 75.35140 70.83161 68.72737 64.12834 62.67237
#> [971] 104.25276 60.27733 79.23571 75.52551 63.60267
#> [976] 123.06404 58.52361 166.05731 79.78348 114.62497
#> [981] 1258.47291 724.48268 59.41181 67.34813 65.25468
#> [986] 71.15085 90.18784 96.56078 61.80775 62.18615
#> [991] 55.17643 325.15859 62.15843 103.21725 76.53048
#> [996] 73.80327 58.69008 75.37047 57.67684 81.54834