![[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.370369e+01 6.055127e+01 6.969381e+01 1.164514e+02 6.723728e+01
#> [6] 2.043206e+02 7.610252e+01 8.158925e+01 7.028845e+01 7.747797e+01
#> [11] 1.189751e+02 7.159700e+01 6.820170e+01 8.637069e+01 8.102388e+01
#> [16] 6.263771e+01 4.086723e+02 2.376041e+02 1.316137e+02 7.185203e+02
#> [21] 6.578344e+01 6.766479e+01 1.471767e+02 1.239311e+02 8.335937e+01
#> [26] 7.124242e+01 8.677207e+01 6.530434e+01 6.388405e+01 8.511301e+01
#> [31] 6.817775e+01 7.199619e+01 1.073230e+02 6.656268e+01 7.437964e+01
#> [36] 7.528113e+01 9.095867e+01 6.829357e+01 6.478423e+01 9.110424e+01
#> [41] 1.629316e+04 1.787105e+02 9.941770e+01 1.373125e+02 8.095628e+01
#> [46] 8.035596e+01 2.130805e+02 1.486018e+02 6.563164e+01 8.232522e+01
#> [51] 6.353354e+01 6.747252e+01 8.629198e+01 4.280563e+02 7.824737e+01
#> [56] 1.044323e+02 7.316779e+01 7.319286e+01 1.658020e+02 9.572126e+01
#> [61] 1.794161e+02 7.712741e+01 1.321730e+02 1.802867e+02 1.033349e+02
#> [66] 1.776331e+02 1.017137e+02 1.326281e+02 4.740773e+02 6.882946e+01
#> [71] 6.550557e+01 6.704436e+01 6.728959e+01 7.427924e+01 8.823195e+01
#> [76] 7.613214e+01 7.026174e+01 1.222989e+02 6.187105e+01 1.111455e+02
#> [81] 1.235300e+04 1.870637e+02 9.464957e+01 7.004060e+01 7.321243e+01
#> [86] 1.081438e+02 1.444494e+02 1.013422e+02 6.634750e+01 9.527541e+01
#> [91] 7.853948e+01 7.886952e+01 6.654932e+01 2.087260e+02 8.676820e+01
#> [96] 2.185624e+02 6.448017e+01 7.194732e+01 1.071211e+02 6.698090e+01
#> [101] 7.425319e+01 6.793091e+01 6.614520e+01 1.007409e+02 1.037640e+04
#> [106] 7.775529e+04 6.544985e+01 8.049329e+01 1.064547e+02 8.873321e+01
#> [111] 7.821632e+01 6.226540e+01 3.873416e+03 6.530803e+02 8.760295e+01
#> [116] 8.023075e+01 8.132754e+01 7.502389e+01 1.193711e+02 6.233840e+02
#> [121] 7.454590e+01 2.576713e+02 1.707654e+02 4.279917e+02 3.752841e+02
#> [126] 2.527776e+02 9.427537e+01 1.265800e+02 1.588293e+02 9.631382e+01
#> [131] 7.430140e+01 8.972831e+01 7.725351e+01 7.140485e+01 7.621364e+01
#> [136] 7.413127e+01 1.458368e+02 2.194858e+02 7.376971e+01 5.142961e+02
#> [141] 6.165784e+01 1.578658e+02 3.360199e+02 8.272435e+01 1.542712e+02
#> [146] 4.619730e+02 6.765847e+01 8.294101e+01 8.060547e+01 8.064307e+01
#> [151] 8.488935e+01 6.602246e+01 1.177756e+02 8.046040e+01 1.313677e+02
#> [156] 1.135793e+02 7.674264e+01 1.241546e+02 1.250081e+02 2.591956e+04
#> [161] 7.217868e+01 1.483390e+02 1.034319e+02 1.398888e+02 6.398817e+01
#> [166] 8.403750e+01 8.515737e+01 2.534995e+03 6.884791e+01 6.102389e+01
#> [171] 6.800072e+01 7.359521e+01 6.769908e+01 7.386165e+01 6.627767e+01
#> [176] 1.093867e+02 9.799857e+01 8.532653e+01 1.048574e+02 6.920301e+01
#> [181] 1.640238e+02 1.411724e+02 2.424515e+02 1.941445e+02 6.089809e+01
#> [186] 9.481491e+01 7.529179e+01 1.013032e+02 7.802725e+01 7.672878e+01
#> [191] 7.862916e+01 8.185321e+01 7.222355e+01 5.529903e+02 8.605042e+01
#> [196] 9.080548e+01 7.372670e+01 9.356604e+01 8.044755e+01 1.896880e+02
#> [201] 5.106996e+02 8.108672e+01 1.079151e+02 1.047444e+02 2.008188e+02
#> [206] 9.889409e+01 9.936313e+01 7.701277e+01 6.215923e+01 6.971067e+01
#> [211] 1.541043e+02 7.692404e+01 1.320328e+02 1.267810e+02 2.179191e+02
#> [216] 1.432966e+02 9.052692e+01 1.480675e+02 1.981981e+02 7.445157e+01
#> [221] 6.273952e+01 6.332071e+01 1.077985e+02 1.026197e+03 8.320448e+01
#> [226] 8.061830e+01 8.865798e+01 9.117559e+01 6.797638e+01 3.153032e+02
#> [231] 6.868911e+01 6.878915e+01 7.964858e+01 1.191873e+02 6.294805e+01
#> [236] 8.409352e+01 1.776022e+02 1.009446e+02 9.905639e+01 1.954121e+02
#> [241] 1.933515e+02 7.177281e+01 8.549089e+01 9.036416e+01 7.290382e+01
#> [246] 1.069246e+02 1.635410e+02 1.772790e+02 7.334260e+01 1.796537e+02
#> [251] 6.653722e+01 1.120516e+02 9.776632e+01 6.830616e+01 7.193411e+01
#> [256] 6.906589e+01 6.951321e+01 7.677961e+01 1.593296e+02 6.735126e+01
#> [261] 1.523611e+02 8.680742e+01 1.009808e+02 7.227355e+01 1.005399e+02
#> [266] 1.371011e+02 7.163852e+01 1.324202e+02 2.060219e+02 7.805162e+01
#> [271] 8.482876e+01 8.925995e+01 8.484759e+01 7.134084e+01 3.209668e+02
#> [276] 1.033305e+02 9.379515e+01 7.971095e+01 7.025159e+02 1.122628e+02
#> [281] 9.006055e+01 1.841829e+02 8.427940e+01 1.430422e+02 6.024331e+01
#> [286] 1.086265e+02 2.799550e+02 1.491301e+02 7.641934e+01 6.803598e+01
#> [291] 8.744305e+01 8.672342e+01 1.127313e+02 2.928462e+02 1.003663e+02
#> [296] 1.241426e+03 6.403505e+01 7.563844e+01 7.475683e+01 1.271142e+02
#> [301] 3.193466e+02 7.364144e+01 9.457220e+01 6.672635e+01 7.188740e+01
#> [306] 1.121113e+02 1.306955e+02 7.348784e+01 7.582563e+01 6.673611e+01
#> [311] 9.511216e+01 1.258363e+02 8.278150e+01 7.881604e+01 8.666232e+01
#> [316] 7.088976e+01 6.365958e+01 6.368255e+01 7.876325e+01 9.546760e+01
#> [321] 2.793219e+02 6.114423e+02 1.609042e+02 1.642850e+02 8.211674e+01
#> [326] 6.838460e+01 6.299042e+01 1.386567e+02 7.440844e+01 6.400438e+01
#> [331] 6.915647e+01 6.882259e+01 1.342515e+02 9.157162e+01 1.054421e+02
#> [336] 7.478739e+01 1.120031e+02 1.469926e+02 9.047918e+01 6.478185e+01
#> [341] 8.274253e+01 1.518143e+02 4.641449e+02 6.008938e+02 1.351111e+02
#> [346] 3.001177e+02 6.369775e+01 6.665518e+01 9.060321e+01 9.959693e+01
#> [351] 7.513984e+03 6.358438e+01 6.997910e+01 9.248979e+01 1.090664e+02
#> [356] 1.618646e+02 9.594822e+01 1.438621e+02 8.436467e+01 7.237965e+01
#> [361] 9.545086e+01 7.710928e+01 6.933441e+01 7.351461e+01 7.923765e+01
#> [366] 7.281107e+01 9.238844e+01 7.494748e+01 6.559723e+01 6.606445e+01
#> [371] 7.053210e+01 1.074581e+02 6.553066e+01 1.209750e+02 6.550213e+01
#> [376] 2.103793e+02 8.569859e+01 9.279102e+01 8.692035e+01 7.118311e+01
#> [381] 6.825801e+01 1.358730e+02 9.764596e+01 9.095680e+01 2.499532e+02
#> [386] 8.867027e+01 7.933315e+01 1.062002e+02 2.553688e+02 1.127493e+02
#> [391] 1.109887e+02 8.875948e+01 7.695360e+01 1.002900e+02 7.968297e+01
#> [396] 9.119783e+03 3.869282e+03 6.340961e+01 9.978527e+01 7.606216e+01
#> [401] 7.428696e+01 7.858328e+01 8.301785e+01 7.909646e+01 5.550156e+02
#> [406] 1.020347e+02 1.621758e+02 1.162906e+02 7.348391e+01 7.762209e+01
#> [411] 8.941607e+01 7.373954e+01 1.733698e+02 1.654466e+02 7.445033e+01
#> [416] 6.427330e+01 9.266698e+01 4.095797e+02 2.967585e+02 1.247514e+02
#> [421] 1.345823e+02 8.865088e+01 2.080006e+02 8.737227e+01 8.704367e+01
#> [426] 8.313590e+02 8.421118e+01 1.126776e+02 9.406691e+01 7.336947e+01
#> [431] 1.124543e+02 1.407048e+02 7.438773e+01 7.792795e+01 6.434922e+01
#> [436] 6.413951e+01 6.326128e+01 1.311885e+02 2.092539e+02 1.550159e+02
#> [441] 1.094454e+02 9.695159e+01 8.501740e+01 7.520028e+01 1.249539e+02
#> [446] 8.206712e+01 1.057463e+02 3.400713e+03 7.471089e+01 8.221086e+01
#> [451] 8.632240e+01 1.588589e+02 6.498728e+01 7.351614e+01 1.272440e+02
#> [456] 9.279556e+01 1.239735e+04 1.061877e+02 8.091102e+01 7.622992e+01
#> [461] 8.038180e+01 1.082422e+02 7.062056e+01 1.055952e+03 3.903488e+02
#> [466] 6.974790e+01 7.746239e+01 1.530159e+02 7.328205e+01 6.995241e+01
#> [471] 2.194382e+04 1.738286e+02 9.369805e+01 6.228288e+01 3.939948e+02
#> [476] 7.535517e+01 6.629765e+01 1.286271e+02 9.701798e+01 6.448999e+01
#> [481] 3.153097e+02 1.485937e+02 1.601579e+02 8.387793e+01 3.726474e+03
#> [486] 6.959176e+01 7.870005e+01 8.460763e+01 9.964884e+01 9.298735e+01
#> [491] 9.377271e+01 7.575257e+01 1.116184e+02 6.995624e+01 7.297400e+01
#> [496] 6.972201e+01 1.169644e+02 8.900269e+01 1.879884e+02 5.871811e+02
#> [501] 1.107664e+02 5.669012e+02 1.991215e+03 6.628741e+01 1.123121e+02
#> [506] 1.096129e+02 8.388252e+01 4.755757e+03 1.355007e+02 6.683014e+01
#> [511] 6.285552e+01 1.039306e+02 7.627062e+01 1.495127e+02 1.394653e+02
#> [516] 6.503802e+01 9.453027e+01 6.605618e+01 1.454253e+02 6.021952e+01
#> [521] 6.879566e+01 7.268108e+01 7.408939e+01 6.712340e+01 6.535568e+01
#> [526] 6.376601e+01 6.847368e+01 1.933821e+02 1.040330e+02 7.093250e+01
#> [531] 7.178419e+01 9.684658e+01 1.032055e+02 3.544813e+02 7.592692e+01
#> [536] 7.780171e+01 7.023316e+01 7.194570e+01 7.290083e+01 7.958638e+01
#> [541] 8.753001e+01 6.953334e+01 6.422297e+01 8.263547e+01 2.125645e+02
#> [546] 1.481267e+03 1.167990e+02 2.743550e+02 1.016145e+02 7.312530e+01
#> [551] 1.541578e+02 9.188637e+01 7.779927e+01 1.521777e+02 9.100651e+01
#> [556] 6.391481e+01 6.630820e+01 8.006381e+01 9.504647e+01 2.783739e+02
#> [561] 7.986867e+01 9.968483e+01 7.123604e+01 8.891627e+01 9.457350e+01
#> [566] 2.927524e+02 6.389904e+01 8.443159e+01 6.411346e+01 2.401022e+02
#> [571] 6.984822e+01 1.399080e+02 1.957645e+02 7.918461e+01 7.532761e+01
#> [576] 2.302723e+03 3.241993e+02 2.502446e+02 6.464030e+01 1.166152e+02
#> [581] 9.283216e+01 1.400791e+02 1.067689e+02 6.997859e+02 7.148021e+01
#> [586] 1.778092e+02 4.501757e+02 1.010936e+02 1.903938e+02 1.196858e+02
#> [591] 6.743121e+01 8.880291e+01 7.256334e+01 8.538961e+01 3.178041e+02
#> [596] 7.415893e+01 6.839867e+01 6.366232e+01 6.577944e+01 9.317026e+01
#> [601] 2.377667e+02 9.459172e+01 7.771471e+01 6.452013e+01 7.662401e+01
#> [606] 1.596407e+02 6.529948e+01 7.233729e+01 1.039572e+02 6.266085e+01
#> [611] 6.891553e+01 1.205933e+02 8.545039e+01 2.362041e+02 1.797014e+02
#> [616] 1.066654e+02 8.113106e+01 6.602509e+01 7.701105e+01 7.662839e+01
#> [621] 1.140711e+02 6.727597e+01 1.224153e+02 1.190970e+02 8.157721e+01
#> [626] 6.509050e+01 1.011844e+02 4.273782e+03 3.427698e+02 1.077175e+02
#> [631] 3.347207e+02 1.116990e+02 2.825026e+02 9.066836e+01 9.247750e+01
#> [636] 6.888844e+01 6.562190e+01 1.054090e+02 8.390976e+01 7.330743e+01
#> [641] 7.532052e+01 7.199190e+01 3.256642e+02 6.768642e+01 8.631005e+01
#> [646] 7.758034e+01 1.245903e+02 7.172553e+01 7.480859e+01 1.339823e+02
#> [651] 7.425270e+01 1.326228e+02 7.826929e+01 7.874822e+01 1.760267e+02
#> [656] 1.124349e+03 2.342495e+02 6.178783e+02 1.048672e+02 1.895919e+02
#> [661] 4.193931e+02 6.888879e+01 2.608789e+02 6.858526e+01 8.507112e+01
#> [666] 9.185240e+01 9.711687e+01 6.905835e+01 8.346337e+01 6.585471e+01
#> [671] 2.400036e+02 1.864164e+02 1.429981e+02 9.874841e+01 7.296311e+01
#> [676] 9.582359e+01 7.900506e+01 1.278841e+02 1.943396e+02 8.582632e+01
#> [681] 8.262586e+01 8.221664e+01 6.357545e+01 6.614219e+01 7.927577e+01
#> [686] 2.029932e+02 7.513505e+01 6.812150e+01 6.797501e+01 7.082870e+01
#> [691] 7.599214e+01 1.265920e+03 6.513714e+01 6.735323e+01 7.262804e+01
#> [696] 8.351224e+01 6.366545e+01 6.295053e+01 7.428151e+01 9.414432e+01
#> [701] 7.855677e+01 6.503375e+01 1.063186e+02 7.549865e+01 6.062621e+03
#> [706] 9.109878e+01 4.978575e+02 1.464028e+02 8.053595e+01 6.999327e+01
#> [711] 1.175129e+02 1.049332e+02 7.471292e+01 7.583735e+01 6.997854e+01
#> [716] 7.169143e+01 7.596138e+01 6.289410e+01 6.620689e+01 5.424040e+02
#> [721] 8.652249e+01 2.478132e+02 1.591496e+02 7.624780e+01 3.444174e+02
#> [726] 1.233293e+03 7.995746e+01 2.488893e+02 6.787874e+01 4.324901e+02
#> [731] 1.833053e+02 6.880706e+01 7.421547e+01 8.143613e+01 7.518872e+01
#> [736] 7.585591e+01 6.729665e+01 2.393058e+02 1.130735e+02 8.240144e+01
#> [741] 8.216180e+01 1.050421e+02 1.072964e+02 6.353301e+01 8.030909e+01
#> [746] 6.943195e+01 7.620572e+01 8.081369e+01 3.113134e+02 2.331739e+02
#> [751] 7.186105e+01 3.155068e+02 8.189499e+01 8.553883e+01 7.821451e+02
#> [756] 7.363789e+01 1.183915e+02 9.216304e+01 8.011883e+01 7.071089e+01
#> [761] 8.363134e+01 6.893598e+01 2.042377e+02 1.995752e+02 6.032777e+01
#> [766] 2.296490e+10 1.248471e+02 6.453849e+01 1.019601e+02 8.498932e+01
#> [771] 9.338817e+02 8.338579e+01 2.389444e+02 6.435680e+01 6.353852e+01
#> [776] 8.520268e+01 7.646317e+01 3.012118e+02 7.087076e+01 6.090041e+02
#> [781] 6.525100e+02 1.828299e+02 7.149438e+01 7.427321e+01 2.275811e+02
#> [786] 9.091988e+01 7.854080e+01 6.972555e+01 8.828870e+01 7.858510e+01
#> [791] 9.985648e+01 8.150360e+01 8.577738e+07 2.161383e+02 1.003014e+02
#> [796] 1.212311e+02 1.151080e+02 8.911291e+01 6.450184e+02 5.634160e+02
#> [801] 8.455496e+01 7.975737e+01 2.288208e+02 7.195917e+04 1.376153e+02
#> [806] 1.146723e+02 7.827034e+01 6.629297e+01 7.610775e+01 6.584395e+01
#> [811] 6.430640e+01 8.952747e+01 1.911759e+02 9.302617e+01 7.054486e+01
#> [816] 7.220903e+01 6.751065e+01 7.591865e+01 3.934905e+02 1.332126e+02
#> [821] 6.722384e+01 7.592692e+01 7.544336e+01 1.490952e+02 1.639376e+02
#> [826] 7.264964e+01 6.669608e+01 1.601587e+02 8.306304e+01 1.181713e+02
#> [831] 8.506602e+01 7.959227e+01 6.194703e+01 6.669634e+01 1.011104e+02
#> [836] 9.732026e+01 4.301960e+02 3.356125e+02 8.158489e+01 1.345717e+02
#> [841] 1.371344e+02 5.030539e+02 7.501657e+01 1.200381e+02 6.414525e+01
#> [846] 7.365507e+01 1.249718e+02 8.683337e+01 1.394565e+02 2.747287e+02
#> [851] 1.908482e+02 8.242638e+01 1.190327e+02 9.020514e+02 1.034860e+02
#> [856] 9.696147e+01 9.051082e+01 6.919265e+01 8.312210e+01 1.214804e+02
#> [861] 7.927919e+01 1.987779e+02 7.996334e+01 1.297768e+02 6.439551e+01
#> [866] 7.719699e+01 6.390930e+01 7.266057e+01 6.501351e+01 7.084823e+01
#> [871] 6.791386e+01 7.502128e+01 1.095797e+02 7.059091e+01 6.634407e+01
#> [876] 1.107026e+02 2.879732e+02 1.771248e+02 7.244605e+01 9.899116e+01
#> [881] 7.200360e+01 7.273564e+01 7.906424e+01 6.800770e+01 2.213677e+02
#> [886] 1.353323e+02 1.412710e+02 6.623808e+01 6.317137e+01 6.331576e+01
#> [891] 7.100745e+01 3.530822e+02 6.959795e+01 8.611192e+01 7.708093e+01
#> [896] 6.797390e+01 7.516295e+01 9.271096e+01 9.079689e+01 1.076867e+02
#> [901] 7.138311e+01 7.371607e+01 8.165104e+01 8.843797e+01 2.332272e+02
#> [906] 1.413925e+02 1.160438e+02 9.228293e+01 7.135404e+01 8.626775e+01
#> [911] 8.035915e+01 7.200545e+01 6.497513e+01 7.121246e+01 3.258752e+02
#> [916] 1.626258e+02 1.202032e+02 9.249865e+01 1.116682e+02 3.418345e+02
#> [921] 8.305663e+01 1.768163e+03 1.211038e+02 7.022933e+01 1.217317e+02
#> [926] 2.016856e+02 7.899179e+01 7.100854e+01 7.913336e+01 1.077038e+02
#> [931] 9.325987e+01 7.482898e+01 7.135191e+01 7.172976e+01 8.768670e+01
#> [936] 9.764819e+01 7.675717e+01 7.116334e+01 7.487822e+01 6.992135e+01
#> [941] 2.370313e+02 9.798444e+01 8.732435e+01 8.647838e+01 9.550572e+01
#> [946] 4.396487e+02 2.284038e+02 1.823814e+02 6.488183e+01 8.748792e+01
#> [951] 1.098143e+02 9.364918e+01 3.318803e+02 8.518061e+01 1.321051e+02
#> [956] 9.404640e+02 2.438181e+02 2.225839e+04 2.283462e+02 3.305370e+02
#> [961] 7.571102e+01 7.111102e+01 7.039212e+01 7.253577e+01 6.016033e+01
#> [966] 8.790544e+01 6.647356e+01 2.912427e+02 6.874865e+01 7.992426e+01
#> [971] 9.719235e+01 1.336616e+02 1.669701e+02 1.938523e+02 6.934959e+01
#> [976] 6.818261e+01 6.213073e+01 7.908396e+08 2.042000e+03 1.648304e+03
#> [981] 5.887806e+01 6.989235e+01 7.699012e+01 8.249938e+01 1.049416e+02
#> [986] 8.576837e+03 1.282666e+03 7.308584e+01 6.282152e+01 6.753564e+01
#> [991] 6.754764e+01 1.080569e+02 9.699905e+01 8.928469e+01 8.537602e+01
#> [996] 7.363167e+01 1.003584e+02 7.182766e+01 6.384294e+01 1.702603e+02
doseFunction(
model,
alpha2 = samples@data$alpha2,
beta = samples@data$beta,
grade = 2L
)(x = 0.25)
#> [1] 6.284675e+01 5.602673e+01 6.225219e+01 8.500584e+01 5.931581e+01
#> [6] 8.771303e+01 6.402897e+01 6.541277e+01 5.551774e+01 6.733564e+01
#> [11] 8.433028e+01 5.923367e+01 5.952844e+01 6.920077e+01 6.543085e+01
#> [16] 5.543550e+01 1.175073e+02 1.160951e+02 8.557743e+01 1.828307e+02
#> [21] 5.871492e+01 5.919842e+01 1.072320e+02 8.281327e+01 7.083957e+01
#> [26] 5.927108e+01 6.157741e+01 5.717078e+01 5.673383e+01 7.011075e+01
#> [31] 5.452507e+01 5.925126e+01 8.412721e+01 5.436223e+01 6.379363e+01
#> [36] 6.303622e+01 4.694319e+01 5.934239e+01 5.468581e+01 6.411294e+01
#> [41] 1.122425e+03 7.342736e+01 6.829482e+01 7.300599e+01 6.626724e+01
#> [46] 6.661925e+01 1.003349e+02 9.874263e+01 5.779629e+01 6.659259e+01
#> [51] 5.876422e+01 6.169638e+01 7.227907e+01 1.872140e+02 6.125099e+01
#> [56] 8.412804e+01 6.067699e+01 6.234317e+01 8.447593e+01 6.133857e+01
#> [61] 9.431611e+01 6.254286e+01 8.654062e+01 9.571270e+01 6.065137e+01
#> [66] 8.163978e+01 6.897847e+01 7.747609e+01 1.090459e+02 6.225340e+01
#> [71] 5.931698e+01 6.063117e+01 6.023236e+01 6.614048e+01 6.602565e+01
#> [76] 6.363130e+01 6.053563e+01 7.139202e+01 4.757660e+01 8.881908e+01
#> [81] 2.728423e+03 1.007926e+02 6.375875e+01 5.712506e+01 6.150889e+01
#> [86] 8.327997e+01 9.646623e+01 7.005363e+01 5.864321e+01 7.029119e+01
#> [91] 5.717956e+01 5.767788e+01 5.501058e+01 1.156867e+02 5.538455e+01
#> [96] 9.881800e+01 5.497109e+01 5.529358e+01 7.695412e+01 6.027191e+01
#> [101] 6.080590e+01 5.557280e+01 5.914542e+01 7.861746e+01 8.680400e+02
#> [106] 6.762439e+03 5.552207e+01 6.498791e+01 7.500372e+01 6.929589e+01
#> [111] 6.225613e+01 5.151446e+01 5.481214e+01 1.963399e+02 4.947330e+01
#> [116] 5.177034e+01 6.292152e+01 6.378204e+01 8.490629e+01 2.536477e+02
#> [121] 6.322296e+01 1.181959e+02 7.556753e+01 1.485380e+02 8.347573e+01
#> [126] 8.516258e+01 6.422612e+01 6.267620e+01 7.829269e+01 6.893641e+01
#> [131] 6.333139e+01 7.171342e+01 6.634605e+01 5.873208e+01 6.218159e+01
#> [136] 6.069653e+01 8.957725e+01 8.311634e+01 5.777280e+01 1.411801e+02
#> [141] 5.338938e+01 1.070214e+02 1.516243e+02 6.848532e+01 8.935282e+01
#> [146] 1.332479e+02 6.157410e+01 6.923818e+01 6.616776e+01 6.608256e+01
#> [151] 6.805224e+01 5.473546e+01 7.334573e+01 6.204140e+01 8.242405e+01
#> [156] 7.410013e+01 6.691778e+01 8.465659e+01 7.706116e+01 3.531512e+02
#> [161] 6.066390e+01 1.016303e+02 5.012539e+01 7.245771e+01 5.486008e+01
#> [166] 6.805878e+01 6.356545e+01 1.471350e+02 5.588875e+01 5.514769e+01
#> [171] 6.120596e+01 6.439495e+01 5.847542e+01 6.301019e+01 5.858321e+01
#> [176] 8.288675e+01 7.925131e+01 7.080014e+01 7.266549e+01 5.888194e+01
#> [181] 8.827947e+01 6.025443e+01 1.488792e+02 1.170538e+02 5.344810e+01
#> [186] 7.350895e+01 6.208181e+01 7.483916e+01 6.109486e+01 6.454994e+01
#> [191] 6.272190e+01 6.267163e+01 5.917615e+01 2.497784e+02 6.707523e+01
#> [196] 6.384057e+01 5.761308e+01 7.348987e+01 5.465529e+01 8.394264e+01
#> [201] 1.194822e+02 6.074518e+01 5.700813e+01 6.922099e+01 8.364610e+01
#> [206] 5.604932e+01 6.457984e+01 6.077628e+01 5.422106e+01 5.851593e+01
#> [211] 9.541338e+01 5.879076e+01 7.228610e+01 8.178370e+01 1.178612e+02
#> [216] 8.875689e+01 6.561374e+01 8.528647e+01 1.234303e+02 5.852997e+01
#> [221] 5.575419e+01 5.396386e+01 5.397085e+01 1.673458e+02 4.594632e+01
#> [226] 6.450038e+01 5.844703e+01 6.802845e+01 5.931691e+01 1.509781e+02
#> [231] 5.880905e+01 5.963193e+01 6.048161e+01 7.214508e+01 5.347984e+01
#> [236] 6.805121e+01 9.414178e+01 7.928663e+01 7.844682e+01 1.163442e+02
#> [241] 6.416398e+01 6.189287e+01 6.622367e+01 6.974967e+01 6.235309e+01
#> [246] 7.799014e+01 7.912645e+01 7.348428e+01 6.022989e+01 1.056944e+02
#> [251] 5.566567e+01 7.563240e+01 6.868661e+01 5.549923e+01 5.820682e+01
#> [256] 5.781747e+01 5.973794e+01 6.490156e+01 9.788327e+01 5.614631e+01
#> [261] 8.001588e+01 6.453332e+01 7.384796e+01 6.002053e+01 7.041291e+01
#> [266] 5.757462e+01 6.245739e+01 9.394140e+01 1.010875e+02 4.905733e+01
#> [271] 6.700926e+01 7.152225e+01 6.240851e+01 5.921732e+01 1.253746e+02
#> [276] 7.525625e+01 6.985151e+01 6.278375e+01 2.537366e+02 5.977254e+01
#> [281] 5.009159e+01 7.024591e+01 4.877699e+01 8.919252e+01 5.425765e+01
#> [286] 7.041059e+01 1.459850e+02 9.322189e+01 5.972711e+01 5.459998e+01
#> [291] 6.803654e+01 6.697178e+01 7.325588e+01 9.189301e+01 6.642386e+01
#> [296] 1.901664e+02 5.594668e+01 6.404988e+01 6.443305e+01 8.403524e+01
#> [301] 1.551530e+02 5.399876e+01 6.860247e+01 5.886451e+01 5.938332e+01
#> [306] 6.837437e+01 7.201812e+01 5.956711e+01 5.938877e+01 5.766849e+01
#> [311] 7.254181e+01 8.404268e+01 6.401881e+01 6.646289e+01 6.998284e+01
#> [316] 6.219955e+01 5.546734e+01 5.439170e+01 5.948802e+01 5.502866e+01
#> [321] 1.198461e+02 1.083844e+02 6.324685e+01 9.753932e+01 5.651036e+01
#> [326] 5.584931e+01 5.535426e+01 7.520782e+01 5.839108e+01 5.550637e+01
#> [331] 5.824763e+01 5.850586e+01 9.983698e+01 6.994729e+01 5.948825e+01
#> [336] 5.535270e+01 6.162829e+01 9.271133e+01 6.502170e+01 5.673952e+01
#> [341] 6.938094e+01 9.062808e+01 1.182692e+02 1.769848e+02 8.148665e+01
#> [346] 8.566144e+01 5.774613e+01 5.986382e+01 7.105243e+01 7.357183e+01
#> [351] 9.494463e+02 5.571026e+01 6.006963e+01 7.213447e+01 7.194063e+01
#> [356] 7.673079e+01 7.082660e+01 7.511808e+01 6.541503e+01 5.831468e+01
#> [361] 7.282705e+01 6.412270e+01 5.758000e+01 5.580713e+01 6.737467e+01
#> [366] 5.929194e+01 6.934220e+01 6.492447e+01 5.780019e+01 5.829506e+01
#> [371] 6.018723e+01 5.747189e+01 5.558866e+01 7.903874e+01 5.894561e+01
#> [376] 1.041994e+02 6.556653e+01 6.712370e+01 6.697241e+01 5.852957e+01
#> [381] 6.188798e+01 8.511944e+01 6.858296e+01 6.309971e+01 1.232391e+02
#> [386] 6.120646e+01 6.642082e+01 7.875050e+01 1.137564e+02 7.145525e+01
#> [391] 7.764227e+01 6.087232e+01 6.440245e+01 7.539500e+01 6.034443e+01
#> [396] 4.830750e+02 4.462336e+02 5.557975e+01 5.330667e+01 5.533474e+01
#> [401] 6.243549e+01 6.120684e+01 6.893110e+01 6.646499e+01 2.240279e+02
#> [406] 7.381332e+01 9.515978e+01 7.194802e+01 5.633185e+01 6.588868e+01
#> [411] 6.671955e+01 6.390408e+01 1.083881e+02 1.061042e+02 6.421331e+01
#> [416] 5.768904e+01 6.654075e+01 1.211414e+02 9.357061e+01 7.098552e+01
#> [421] 6.129528e+01 5.927659e+01 1.249770e+02 6.165077e+01 6.300279e+01
#> [426] 1.607730e+02 6.196873e+01 8.041183e+01 6.459466e+01 5.881942e+01
#> [431] 7.659900e+01 8.534808e+01 6.064593e+01 6.453567e+01 5.657132e+01
#> [436] 5.799564e+01 5.786025e+01 9.751714e+01 9.195683e+01 8.215470e+01
#> [441] 7.062012e+01 6.340320e+01 6.678898e+01 6.007194e+01 6.122903e+01
#> [446] 6.839843e+01 7.066364e+01 4.728936e+02 6.223642e+01 6.314392e+01
#> [451] 6.027378e+01 8.694285e+01 5.262251e+01 5.629528e+01 7.772927e+01
#> [456] 7.089797e+01 8.932549e+02 8.554367e+01 6.514807e+01 6.363172e+01
#> [461] 6.323682e+01 7.851454e+01 5.931854e+01 2.273237e+02 9.145770e+01
#> [466] 5.440632e+01 6.069127e+01 7.778940e+01 6.276512e+01 6.112539e+01
#> [471] 1.025549e+03 8.247239e+01 6.496676e+01 5.572071e+01 1.720948e+02
#> [476] 6.109847e+01 5.802331e+01 7.649142e+01 7.764314e+01 5.572808e+01
#> [481] 1.650670e+02 9.728838e+01 8.312669e+01 5.701327e+01 4.343104e+02
#> [486] 5.641003e+01 5.951516e+01 6.300281e+01 6.122630e+01 6.668446e+01
#> [491] 6.619052e+01 6.151670e+01 7.685368e+01 4.908355e+01 5.342770e+01
#> [496] 6.003879e+01 7.336990e+01 6.489406e+01 1.165544e+02 2.074599e+02
#> [501] 7.720927e+01 1.837456e+02 4.021323e+02 5.921987e+01 6.347952e+01
#> [506] 7.820690e+01 6.907715e+01 1.062932e+03 7.011417e+01 5.754407e+01
#> [511] 5.704486e+01 7.927784e+01 6.452533e+01 7.036939e+01 8.610421e+01
#> [516] 5.837519e+01 6.766702e+01 5.783322e+01 8.052238e+01 4.920005e+01
#> [521] 6.015390e+01 6.003598e+01 5.874781e+01 5.628679e+01 5.791545e+01
#> [526] 5.800887e+01 5.962725e+01 1.180124e+02 7.956346e+01 6.115811e+01
#> [531] 6.143575e+01 6.736933e+01 6.383161e+01 1.474462e+02 4.546585e+01
#> [536] 4.794481e+01 5.984461e+01 6.109000e+01 6.166675e+01 6.334402e+01
#> [541] 6.121008e+01 5.730999e+01 5.784134e+01 6.688143e+01 1.150771e+02
#> [546] 4.533638e+02 8.637669e+01 1.131348e+02 7.217106e+01 6.311060e+01
#> [551] 9.357437e+01 6.199222e+01 5.816302e+01 7.967828e+01 6.758922e+01
#> [556] 5.579120e+01 6.062635e+01 6.803244e+01 6.413683e+01 1.193452e+02
#> [561] 6.113263e+01 6.641246e+01 6.064231e+01 6.296981e+01 6.798409e+01
#> [566] 1.121214e+02 5.340080e+01 6.693647e+01 5.807102e+01 1.466261e+02
#> [571] 6.232730e+01 1.052575e+02 1.101852e+02 6.338576e+01 6.471779e+01
#> [576] 5.692721e+02 1.458718e+02 1.260175e+02 5.251814e+01 7.287250e+01
#> [581] 6.616682e+01 8.998214e+01 6.331391e+01 2.355378e+02 6.060652e+01
#> [586] 1.217254e+02 2.318354e+02 8.045950e+01 9.552680e+01 8.638191e+01
#> [591] 5.442274e+01 6.481930e+01 5.877825e+01 6.392316e+01 1.243568e+02
#> [596] 6.240908e+01 5.879682e+01 5.732729e+01 5.938769e+01 6.957140e+01
#> [601] 7.927490e+01 7.311161e+01 6.564228e+01 5.591986e+01 6.307551e+01
#> [606] 7.939464e+01 5.796776e+01 5.761099e+01 6.820201e+01 5.603531e+01
#> [611] 6.049304e+01 8.170771e+01 6.865367e+01 8.291155e+01 1.049610e+02
#> [616] 6.657090e+01 5.533335e+01 5.385931e+01 6.234565e+01 6.186309e+01
#> [621] 5.524424e+01 5.470118e+01 5.330383e+01 7.489639e+01 6.236123e+01
#> [626] 5.694227e+01 7.043807e+01 5.380833e+02 1.371127e+02 8.252588e+01
#> [631] 1.392005e+02 6.486122e+01 1.398249e+02 6.557427e+01 7.148959e+01
#> [636] 6.288692e+01 5.954625e+01 8.086029e+01 6.505459e+01 5.388893e+01
#> [641] 5.769852e+01 6.000615e+01 1.269655e+02 5.649893e+01 6.234138e+01
#> [646] 5.692364e+01 7.455334e+01 6.369095e+01 6.607883e+01 9.693657e+01
#> [651] 5.983413e+01 8.290385e+01 6.142842e+01 6.445079e+01 1.101869e+02
#> [656] 3.083142e+02 5.620322e+01 1.879679e+02 6.822832e+01 1.239892e+02
#> [661] 1.604419e+02 5.739056e+01 1.188770e+02 6.153835e+01 6.201309e+01
#> [666] 6.737926e+01 5.395585e+01 5.936371e+01 6.653947e+01 5.666396e+01
#> [671] 1.248487e+02 6.262012e+01 8.162392e+01 5.471372e+01 5.213021e+01
#> [676] 7.281528e+01 6.494437e+01 8.795981e+01 1.001890e+02 6.998794e+01
#> [681] 6.507577e+01 6.551400e+01 5.559477e+01 5.759232e+01 6.832967e+01
#> [686] 1.056555e+02 6.308216e+01 6.082309e+01 6.043474e+01 6.251917e+01
#> [691] 6.668644e+01 3.143869e+02 5.654465e+01 5.950389e+01 5.990273e+01
#> [696] 6.997072e+01 5.868099e+01 5.822704e+01 6.406744e+01 7.309243e+01
#> [701] 6.562430e+01 5.861306e+01 7.705654e+01 5.896220e+01 6.702619e+02
#> [706] 6.349251e+01 1.338463e+02 6.130319e+01 5.731753e+01 5.978644e+01
#> [711] 8.395828e+01 7.872027e+01 6.175238e+01 5.847790e+01 6.153202e+01
#> [716] 6.096698e+01 5.920852e+01 5.654186e+01 5.885202e+01 1.583027e+02
#> [721] 6.211749e+01 1.336027e+02 7.762633e+01 5.803356e+01 1.221151e+02
#> [726] 1.257819e+02 5.984931e+01 1.266058e+02 6.121576e+01 1.291686e+02
#> [731] 7.653346e+01 5.161140e+01 6.277175e+01 6.744701e+01 6.130666e+01
#> [736] 5.835029e+01 5.603716e+01 8.079116e+01 6.644541e+01 5.466266e+01
#> [741] 5.942141e+01 5.955137e+01 7.120081e+01 5.535751e+01 6.676734e+01
#> [746] 5.692609e+01 5.633360e+01 6.152526e+01 1.178849e+02 1.107901e+02
#> [751] 6.022024e+01 9.205374e+01 6.460775e+01 6.822438e+01 2.716118e+02
#> [756] 5.723442e+01 8.022996e+01 7.076004e+01 6.276955e+01 5.753016e+01
#> [761] 5.933384e+01 5.833579e+01 9.856446e+01 9.975981e+01 5.220351e+01
#> [766] 4.505885e+06 8.673611e+01 5.724859e+01 7.145501e+01 6.981619e+01
#> [771] 3.988090e+02 6.627850e+01 1.485209e+02 5.962670e+01 5.817805e+01
#> [776] 6.935981e+01 6.402295e+01 1.567090e+02 6.233374e+01 1.484544e+02
#> [781] 1.065750e+02 7.216894e+01 5.499361e+01 6.299327e+01 1.183532e+02
#> [786] 6.724475e+01 6.103542e+01 5.635456e+01 5.574856e+01 5.942213e+01
#> [791] 8.161507e+01 6.995181e+01 6.789855e+03 9.289776e+01 7.079852e+01
#> [796] 8.769512e+01 6.040852e+01 6.539246e+01 1.009331e+02 1.143487e+02
#> [801] 6.912337e+01 6.467536e+01 9.997120e+01 6.972334e+03 9.900607e+01
#> [806] 8.553364e+01 6.596243e+01 5.921941e+01 6.562024e+01 5.941776e+01
#> [811] 5.631043e+01 7.125960e+01 8.854773e+01 6.674454e+01 5.702922e+01
#> [816] 5.993268e+01 5.983552e+01 6.253892e+01 1.474099e+02 8.983351e+01
#> [821] 6.209798e+01 6.352181e+01 5.714221e+01 9.796971e+01 9.197350e+01
#> [826] 5.557744e+01 5.843686e+01 1.022713e+02 6.860103e+01 8.115768e+01
#> [831] 6.716346e+01 5.980791e+01 5.242140e+01 5.568146e+01 8.184612e+01
#> [836] 7.647459e+01 9.098775e+01 1.340049e+02 6.553271e+01 6.783620e+01
#> [841] 8.807580e+01 1.415547e+02 6.068924e+01 7.946282e+01 5.863542e+01
#> [846] 6.571536e+01 8.842104e+01 6.426029e+01 9.413960e+01 1.024497e+02
#> [851] 9.629099e+01 5.424358e+01 8.534521e+01 2.847754e+02 8.215544e+01
#> [856] 7.403897e+01 5.964315e+01 5.981258e+01 7.259382e+01 8.195547e+01
#> [861] 6.642307e+01 1.179652e+02 5.610905e+01 8.908214e+01 5.747899e+01
#> [866] 6.320855e+01 5.843977e+01 6.341494e+01 5.778683e+01 6.012958e+01
#> [871] 6.018635e+01 6.125889e+01 7.497140e+01 6.199093e+01 5.695342e+01
#> [876] 7.391720e+01 1.073800e+02 9.195539e+01 6.114031e+01 7.771552e+01
#> [881] 6.131382e+01 5.667855e+01 5.085838e+01 5.882337e+01 1.050369e+02
#> [886] 8.289525e+01 8.234100e+01 5.834962e+01 5.786792e+01 5.841183e+01
#> [891] 6.160002e+01 1.445254e+02 6.019507e+01 5.383452e+01 5.157729e+01
#> [896] 5.576267e+01 6.028785e+01 7.443362e+01 6.922282e+01 7.951403e+01
#> [901] 6.244787e+01 6.276495e+01 6.712066e+01 6.033312e+01 9.931468e+01
#> [906] 9.328545e+01 8.944370e+01 6.489107e+01 6.142655e+01 7.333934e+01
#> [911] 6.982047e+01 6.392118e+01 5.754022e+01 6.243743e+01 1.489919e+02
#> [916] 6.001109e+01 6.199570e+01 6.211264e+01 7.672561e+01 1.230857e+02
#> [921] 6.019210e+01 2.147024e+02 6.850198e+01 5.259524e+01 7.198791e+01
#> [926] 1.126093e+02 6.617526e+01 5.944048e+01 6.362497e+01 7.874233e+01
#> [931] 6.265174e+01 5.416141e+01 5.970097e+01 6.064403e+01 6.574640e+01
#> [936] 7.060996e+01 6.142746e+01 5.902118e+01 6.366716e+01 5.932739e+01
#> [941] 9.302623e+01 7.131007e+01 6.566337e+01 6.699624e+01 7.917179e+01
#> [946] 1.676305e+02 8.643376e+01 8.125469e+01 5.641026e+01 6.641329e+01
#> [951] 6.651158e+01 5.778651e+01 1.279509e+02 6.492168e+01 7.731963e+01
#> [956] 2.992049e+02 8.337278e+01 2.533900e+03 1.343526e+02 1.577914e+02
#> [961] 6.641383e+01 6.057139e+01 5.790116e+01 5.817339e+01 5.481596e+01
#> [966] 6.292377e+01 5.828072e+01 1.559352e+02 6.014549e+01 6.482905e+01
#> [971] 5.304182e+01 5.517062e+01 8.620572e+01 1.106819e+02 5.428896e+01
#> [976] 5.626599e+01 5.557180e+01 1.186092e+05 3.584533e+02 3.372703e+02
#> [981] 5.463643e+01 6.331552e+01 6.757253e+01 6.880236e+01 7.555305e+01
#> [986] 6.435170e+02 2.489779e+02 5.791797e+01 5.643494e+01 6.138789e+01
#> [991] 6.092209e+01 7.597493e+01 6.931168e+01 6.616360e+01 5.578966e+01
#> [996] 6.190371e+01 5.792169e+01 6.082329e+01 5.911423e+01 7.323485e+01