![[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] 7.080373e+01 8.869017e+01 6.826271e+01 8.652246e+01 9.438530e+01
#> [6] 8.849656e+01 7.725880e+01 7.015755e+01 1.853950e+02 1.533189e+02
#> [11] 1.549599e+02 1.042256e+02 8.648315e+01 6.832477e+01 7.026453e+01
#> [16] 6.711153e+01 7.047141e+01 1.135047e+02 7.701359e+01 7.531818e+01
#> [21] 1.111193e+02 5.325105e+03 8.114825e+01 9.195174e+01 2.069354e+02
#> [26] 3.180295e+02 9.239122e+01 6.454455e+01 1.618167e+02 2.935268e+02
#> [31] 6.499695e+01 1.656648e+02 4.733388e+02 1.277961e+02 3.246174e+02
#> [36] 9.662241e+01 1.239712e+02 2.102986e+02 2.462972e+02 7.146283e+01
#> [41] 7.005813e+01 6.651181e+01 7.725281e+01 1.001734e+02 8.700998e+01
#> [46] 6.979705e+01 6.370524e+01 1.716023e+02 6.834694e+01 9.869346e+01
#> [51] 9.994834e+01 5.144271e+02 1.098049e+03 1.005869e+02 1.044762e+02
#> [56] 6.594323e+01 1.707231e+02 1.498309e+02 6.737420e+01 4.969413e+02
#> [61] 1.057542e+02 1.544415e+02 1.159851e+02 7.448722e+01 7.317701e+01
#> [66] 1.156236e+02 7.144963e+01 9.740249e+01 8.176045e+01 6.958094e+01
#> [71] 8.011556e+01 8.843357e+01 1.648046e+02 7.360082e+01 8.212367e+01
#> [76] 3.396053e+03 7.767187e+01 7.778140e+01 6.622053e+01 1.061040e+02
#> [81] 7.170560e+01 6.426231e+01 7.717380e+01 7.414281e+01 1.078321e+02
#> [86] 6.679103e+01 8.474478e+01 8.432614e+01 8.400564e+01 7.618026e+01
#> [91] 6.280649e+01 8.608626e+02 8.089791e+01 7.721010e+01 8.305988e+01
#> [96] 8.559886e+01 8.143953e+01 1.285841e+02 3.009148e+03 9.355490e+02
#> [101] 1.169358e+03 6.920428e+01 9.368075e+01 7.636284e+01 1.312499e+02
#> [106] 8.339120e+01 6.299982e+01 6.275811e+01 8.040172e+01 5.661254e+02
#> [111] 2.306698e+02 9.102658e+01 6.596301e+01 2.032138e+02 8.515927e+01
#> [116] 5.949524e+01 7.859407e+01 6.046866e+01 8.396098e+01 6.976440e+01
#> [121] 6.755457e+01 1.062759e+02 9.537299e+01 1.113485e+02 9.804086e+01
#> [126] 1.116427e+02 1.203020e+02 1.644178e+02 8.783137e+02 2.210210e+07
#> [131] 5.790370e+16 5.735747e+01 5.933205e+01 2.997946e+02 1.224545e+02
#> [136] 1.163486e+02 6.411928e+01 6.458452e+01 6.468556e+01 1.006059e+02
#> [141] 2.229803e+02 6.580993e+01 6.584075e+01 7.397061e+01 6.329348e+01
#> [146] 7.321346e+02 9.524413e+01 1.449947e+02 6.415246e+01 6.888483e+01
#> [151] 1.073408e+02 7.548056e+01 7.003456e+01 1.104684e+02 1.637964e+02
#> [156] 6.859453e+01 1.108599e+02 7.731009e+01 2.118965e+02 1.336403e+02
#> [161] 8.883917e+01 1.153287e+02 6.951755e+01 8.624468e+01 8.251995e+01
#> [166] 6.241349e+01 4.940146e+02 1.041619e+02 7.971660e+02 7.069244e+01
#> [171] 1.093092e+02 8.669449e+01 8.099084e+01 8.234517e+01 8.096888e+01
#> [176] 4.328644e+02 9.359848e+01 4.008721e+02 7.125616e+01 7.702014e+01
#> [181] 6.876755e+01 6.756449e+01 6.622925e+01 6.675843e+01 6.380368e+01
#> [186] 8.486697e+01 6.751697e+01 6.628286e+01 8.893774e+01 8.080068e+01
#> [191] 1.802529e+02 7.352449e+01 1.030096e+02 8.980106e+01 8.091721e+01
#> [196] 1.301382e+02 8.653436e+01 9.253504e+01 8.179193e+01 6.992625e+01
#> [201] 2.404800e+02 6.973121e+01 9.212565e+02 5.433574e+02 3.836474e+05
#> [206] 3.739871e+02 4.424357e+02 1.054813e+02 2.708514e+02 7.356122e+01
#> [211] 7.274315e+01 1.272016e+02 1.724205e+02 8.289583e+01 1.700425e+02
#> [216] 7.168741e+01 7.365104e+02 1.822831e+03 8.209159e+01 1.785696e+02
#> [221] 1.141271e+02 6.744544e+01 6.383731e+01 8.221351e+01 6.969209e+01
#> [226] 7.554678e+01 8.958681e+01 8.102322e+01 1.508509e+02 7.718123e+01
#> [231] 1.489550e+02 1.654308e+02 9.382194e+01 9.166300e+01 1.364271e+02
#> [236] 1.204410e+02 8.443299e+01 7.456492e+01 1.561681e+02 8.024182e+01
#> [241] 1.209536e+02 7.299942e+01 4.412105e+03 1.859984e+05 1.347029e+03
#> [246] 8.319888e+01 7.650008e+01 9.826034e+01 6.151675e+01 6.261038e+01
#> [251] 5.903775e+04 1.868391e+02 6.823563e+01 9.698837e+01 1.348526e+02
#> [256] 6.880504e+02 8.583503e+01 8.177284e+01 1.214271e+02 8.240909e+01
#> [261] 6.887280e+01 9.725592e+01 7.017772e+01 1.084451e+02 9.259631e+01
#> [266] 6.251178e+01 2.626915e+02 8.082243e+01 1.002321e+02 1.432679e+02
#> [271] 1.174266e+02 1.944184e+02 1.264443e+03 8.194017e+01 1.307588e+02
#> [276] 1.019272e+03 7.332166e+01 8.190009e+01 5.983636e+01 1.142391e+03
#> [281] 1.235655e+02 8.461419e+01 8.774764e+01 6.623606e+01 6.312449e+01
#> [286] 7.668650e+01 7.218623e+01 8.530803e+01 4.197552e+02 8.722194e+01
#> [291] 1.206761e+02 9.830668e+01 1.447785e+02 7.097495e+01 4.055397e+02
#> [296] 7.835919e+01 2.129655e+02 6.630051e+01 2.604919e+02 8.358436e+02
#> [301] 1.151435e+02 7.912024e+01 6.908726e+01 1.865760e+02 7.872637e+01
#> [306] 7.206476e+01 2.441711e+02 7.778255e+01 7.484652e+01 2.361489e+02
#> [311] 2.113328e+02 9.261242e+01 7.232261e+01 6.772571e+01 1.005285e+02
#> [316] 6.189542e+01 2.088908e+03 1.005847e+03 6.483798e+01 1.562611e+02
#> [321] 8.232033e+01 6.750728e+01 6.948226e+01 1.110913e+02 1.224193e+02
#> [326] 1.597188e+02 8.994374e+01 1.321877e+02 1.001859e+02 7.955371e+01
#> [331] 7.628024e+01 1.135599e+02 8.830703e+01 8.948912e+01 7.914120e+01
#> [336] 6.839270e+01 7.163580e+01 7.907070e+01 9.238591e+01 5.033159e+02
#> [341] 7.080261e+01 8.803860e+01 7.508917e+01 7.375954e+01 3.973685e+02
#> [346] 8.379329e+01 1.051583e+02 8.794600e+01 7.282710e+01 7.944481e+01
#> [351] 9.972225e+01 1.861184e+02 9.437849e+01 8.435755e+01 8.371607e+01
#> [356] 1.285274e+02 8.803436e+01 8.141054e+01 6.451136e+01 9.420516e+01
#> [361] 6.446578e+01 8.817432e+01 9.524834e+01 6.671626e+01 1.393826e+02
#> [366] 1.179170e+02 9.610201e+01 6.736998e+01 1.546215e+02 6.463013e+01
#> [371] 1.082617e+02 1.996595e+02 7.931295e+01 7.109182e+01 1.135544e+02
#> [376] 1.307140e+02 1.203257e+02 1.076731e+02 1.511483e+02 7.425646e+01
#> [381] 1.340762e+02 1.635907e+02 6.242227e+02 6.106150e+01 1.051107e+02
#> [386] 6.961496e+01 1.845896e+02 1.186831e+02 9.704387e+01 3.241685e+02
#> [391] 6.312128e+01 1.324489e+02 1.359651e+02 1.079971e+02 9.002245e+01
#> [396] 9.297222e+01 6.928741e+01 8.770350e+01 8.567541e+01 2.665536e+02
#> [401] 1.134907e+02 9.633479e+01 1.013547e+02 8.547180e+01 1.160386e+02
#> [406] 9.132395e+01 8.498602e+01 6.520493e+01 1.081714e+02 6.761959e+01
#> [411] 9.213911e+01 9.898408e+01 1.183442e+02 7.236266e+01 6.579321e+02
#> [416] 9.834864e+01 1.556218e+02 6.718956e+01 9.378610e+01 9.704570e+01
#> [421] 8.119936e+01 7.100323e+01 6.962202e+01 3.409030e+02 7.972245e+02
#> [426] 7.240034e+01 1.238534e+02 1.373807e+02 6.911355e+01 3.832404e+02
#> [431] 3.070907e+02 6.785645e+01 9.810862e+01 8.381970e+01 7.689528e+01
#> [436] 8.393696e+01 1.096022e+02 1.010972e+02 7.939766e+01 2.231283e+02
#> [441] 1.227283e+02 1.524557e+02 9.306487e+01 1.104831e+02 9.054346e+01
#> [446] 8.709765e+01 6.968040e+01 9.522185e+01 9.952597e+01 8.278318e+01
#> [451] 8.891432e+01 7.486152e+01 1.059077e+02 1.069412e+02 1.425180e+02
#> [456] 1.031885e+02 9.871287e+01 6.500301e+01 7.283140e+01 1.052339e+02
#> [461] 7.298702e+01 7.718431e+01 9.721298e+01 8.382140e+01 7.646985e+01
#> [466] 1.206252e+02 3.084256e+02 5.057670e+02 6.996368e+01 9.615772e+01
#> [471] 8.443681e+01 4.944628e+02 6.606759e+01 6.250158e+01 1.280098e+02
#> [476] 9.031527e+01 8.358314e+01 7.594941e+01 9.439930e+01 7.278323e+01
#> [481] 8.153576e+01 7.022139e+01 1.711542e+02 3.450865e+02 7.180716e+01
#> [486] 9.224463e+01 1.527854e+02 6.200662e+02 9.060886e+01 9.274011e+01
#> [491] 6.541440e+01 7.088151e+01 6.877015e+01 9.226011e+01 9.208050e+01
#> [496] 1.094597e+03 1.051389e+02 6.558925e+01 7.488503e+01 7.145018e+01
#> [501] 6.837389e+01 1.010167e+02 9.284466e+01 9.134591e+01 7.309717e+01
#> [506] 1.057246e+02 2.057342e+02 1.956570e+02 3.243494e+02 5.849482e+01
#> [511] 2.168689e+02 7.686117e+01 6.600928e+01 9.241012e+01 7.708233e+01
#> [516] 1.414479e+02 7.171286e+01 1.082333e+02 9.280099e+01 7.334825e+01
#> [521] 7.329875e+01 1.467893e+02 7.642881e+01 8.684087e+01 9.618643e+01
#> [526] 6.793744e+01 7.485310e+01 6.037891e+02 6.329520e+01 1.311588e+02
#> [531] 6.594908e+01 9.597002e+01 1.101979e+02 7.359247e+01 6.154913e+01
#> [536] 6.416661e+01 6.529807e+01 7.290708e+01 7.422009e+01 1.165511e+02
#> [541] 6.879729e+01 6.888894e+01 8.122377e+01 8.411232e+01 9.277840e+01
#> [546] 6.977409e+01 6.519274e+01 7.634647e+01 2.564997e+02 1.724141e+02
#> [551] 7.263520e+01 6.319601e+01 1.968091e+02 6.744684e+01 8.042339e+01
#> [556] 1.228694e+02 1.212477e+02 8.541299e+01 1.156952e+02 6.835787e+01
#> [561] 7.874250e+01 7.929232e+02 6.828885e+01 6.719223e+01 6.431038e+01
#> [566] 1.310418e+02 1.831795e+02 1.199027e+02 6.803581e+01 8.127535e+01
#> [571] 2.239397e+02 6.436874e+01 7.370015e+01 7.189574e+01 3.716446e+02
#> [576] 7.199662e+01 2.378999e+02 2.910788e+02 7.421295e+01 7.028130e+01
#> [581] 7.395942e+01 1.296644e+02 1.140061e+02 8.045585e+01 6.463671e+01
#> [586] 7.301813e+01 7.929446e+01 8.892335e+01 8.833284e+01 8.496247e+01
#> [591] 9.138437e+01 7.667838e+01 7.111607e+01 9.824532e+01 8.421099e+01
#> [596] 1.779100e+02 8.140908e+01 2.162850e+03 4.254387e+02 1.061267e+02
#> [601] 1.043147e+02 7.312830e+01 2.513683e+02 1.173098e+02 1.004112e+02
#> [606] 9.152816e+01 6.162333e+01 6.355552e+01 9.149468e+01 1.747794e+02
#> [611] 6.728764e+01 9.720841e+01 7.765021e+01 1.075010e+03 3.422578e+02
#> [616] 5.501185e+02 1.215428e+02 8.094406e+01 6.277525e+01 7.868128e+01
#> [621] 6.439811e+01 1.081592e+02 6.696037e+01 4.634805e+02 7.532939e+01
#> [626] 6.415673e+01 1.289944e+02 7.157203e+01 9.347852e+01 8.513858e+01
#> [631] 1.416240e+02 1.454837e+02 1.207042e+02 6.795068e+01 1.394039e+02
#> [636] 9.988029e+01 9.191588e+01 3.957813e+02 8.552327e+01 7.022859e+01
#> [641] 6.085379e+02 6.335488e+01 7.754012e+01 1.351838e+02 9.227322e+01
#> [646] 9.406209e+01 5.328534e+02 7.505968e+01 5.892984e+03 1.428882e+02
#> [651] 1.232354e+02 6.868248e+01 4.607836e+02 1.079427e+02 9.126482e+01
#> [656] 1.243022e+02 8.657949e+01 9.896312e+01 6.430357e+01 9.723377e+01
#> [661] 7.633011e+01 1.168820e+02 7.509169e+01 8.471215e+01 7.361017e+01
#> [666] 6.607420e+01 7.754331e+01 7.014349e+01 7.594156e+01 6.933339e+01
#> [671] 1.556048e+02 7.183256e+01 8.080586e+02 7.302664e+01 8.989750e+01
#> [676] 1.118797e+02 1.049823e+02 6.420610e+01 8.010083e+01 7.914769e+01
#> [681] 7.983078e+01 9.970571e+01 1.769640e+02 6.762008e+01 1.344429e+02
#> [686] 4.143117e+02 3.340528e+02 5.812376e+02 7.385271e+01 9.128983e+01
#> [691] 5.817779e+01 8.138680e+01 1.378784e+02 7.762039e+01 1.193120e+02
#> [696] 3.761379e+03 1.268079e+02 8.552298e+01 6.982507e+01 7.108308e+01
#> [701] 1.567912e+02 3.793025e+02 8.346938e+01 7.167144e+01 1.072412e+02
#> [706] 6.206258e+01 6.383938e+01 7.107492e+01 1.140922e+02 1.880253e+02
#> [711] 8.313203e+01 8.179370e+01 9.417642e+01 1.097300e+02 6.498974e+01
#> [716] 8.386881e+01 1.590998e+02 7.270361e+01 7.282797e+01 6.660656e+01
#> [721] 7.687201e+01 7.181251e+01 1.066514e+02 7.981393e+01 9.861499e+01
#> [726] 8.461828e+01 7.306072e+01 6.911591e+01 7.891616e+01 7.676902e+01
#> [731] 1.252535e+02 7.250764e+01 1.068000e+02 6.182558e+02 1.324589e+03
#> [736] 1.051026e+02 1.283177e+02 6.435290e+01 6.533600e+03 9.697180e+01
#> [741] 1.437599e+02 8.984445e+01 2.271745e+02 1.272172e+02 7.173062e+01
#> [746] 7.735469e+01 7.861470e+01 2.021744e+02 6.836152e+01 8.102373e+01
#> [751] 1.096932e+02 1.123534e+02 1.499710e+02 1.259325e+02 7.730725e+01
#> [756] 6.810615e+01 8.564360e+01 7.018214e+01 2.094798e+02 6.151967e+01
#> [761] 9.773949e+01 1.744641e+02 9.668649e+01 2.161777e+02 7.349222e+01
#> [766] 4.215189e+06 8.098725e+01 8.528646e+01 8.923471e+01 1.298428e+02
#> [771] 7.959045e+01 6.621836e+01 1.226081e+02 9.457958e+01 1.065011e+02
#> [776] 1.136964e+02 5.256708e+02 4.224456e+03 1.021866e+02 1.076393e+02
#> [781] 8.246603e+01 6.149743e+01 7.151445e+01 1.109373e+02 9.148234e+01
#> [786] 9.341318e+01 6.684756e+01 1.677410e+02 7.464737e+01 1.310716e+02
#> [791] 7.125264e+01 9.464228e+01 8.642893e+02 2.924889e+02 1.215353e+02
#> [796] 8.789521e+01 4.190930e+02 6.670377e+03 1.209057e+02 6.472241e+01
#> [801] 9.090742e+02 2.399694e+02 2.009381e+02 7.488196e+01 6.926315e+01
#> [806] 1.270365e+02 1.197500e+02 1.101151e+02 1.612831e+02 2.031770e+02
#> [811] 3.150430e+03 4.193188e+04 7.394384e+01 7.832724e+01 9.672687e+01
#> [816] 6.990044e+01 6.796629e+01 7.606493e+01 8.014270e+01 9.708789e+01
#> [821] 8.050228e+01 2.041358e+02 6.680361e+01 1.250073e+02 6.732998e+01
#> [826] 8.305400e+01 9.686738e+01 9.671434e+01 9.550897e+01 1.708425e+02
#> [831] 7.290106e+01 1.412661e+02 7.239418e+01 6.769802e+01 6.702094e+01
#> [836] 8.262084e+01 7.631506e+01 9.576352e+01 9.626390e+01 7.100879e+01
#> [841] 1.683827e+02 1.968528e+02 1.091948e+02 6.893132e+01 1.122485e+02
#> [846] 6.940957e+01 1.049040e+02 8.646472e+01 9.409937e+01 5.740371e+02
#> [851] 1.538176e+02 7.635436e+01 7.130124e+01 5.599310e+02 1.416745e+02
#> [856] 2.708828e+02 8.094121e+01 2.338548e+02 2.032082e+02 1.251206e+02
#> [861] 5.566722e+02 6.306990e+01 9.303924e+01 2.586983e+02 1.033192e+02
#> [866] 7.160617e+01 8.570982e+01 7.878432e+01 2.259220e+02 1.897470e+02
#> [871] 8.907341e+01 6.033301e+01 6.077997e+01 8.679013e+01 1.308061e+02
#> [876] 1.008972e+02 1.316271e+02 6.636462e+01 7.359404e+01 6.867608e+01
#> [881] 6.882456e+01 7.210904e+01 8.547990e+01 1.275685e+02 7.609679e+01
#> [886] 1.434754e+02 1.319543e+02 6.753776e+02 1.438336e+02 1.059685e+02
#> [891] 7.434555e+01 1.484855e+02 3.209014e+02 1.084178e+02 1.425097e+02
#> [896] 1.290789e+02 1.196934e+02 7.205979e+01 1.069111e+02 7.598156e+01
#> [901] 8.922933e+01 8.917758e+01 5.574900e+02 1.733795e+02 6.940298e+01
#> [906] 2.409340e+04 1.805156e+02 6.725660e+02 1.362450e+02 6.620286e+01
#> [911] 1.341174e+02 9.992282e+01 9.078104e+01 7.591859e+01 7.209540e+01
#> [916] 1.064012e+02 6.388338e+01 7.131527e+01 7.318407e+01 9.837827e+01
#> [921] 1.773939e+02 7.920454e+01 6.179327e+01 7.733973e+01 8.372419e+01
#> [926] 8.467738e+01 1.393527e+02 1.060012e+02 2.690027e+03 6.272503e+01
#> [931] 6.551602e+03 9.309938e+01 2.176880e+02 9.415349e+01 9.340863e+01
#> [936] 1.061922e+02 2.117716e+02 7.630357e+01 7.264863e+01 7.965024e+01
#> [941] 1.308737e+02 6.575503e+01 3.671621e+02 7.384812e+01 1.054404e+02
#> [946] 1.183468e+02 7.017922e+01 1.153814e+02 7.265650e+01 3.557745e+02
#> [951] 1.045064e+02 7.186329e+01 1.702759e+02 6.812012e+01 9.537514e+01
#> [956] 1.260389e+02 7.598242e+01 1.121719e+02 9.445213e+01 8.741228e+01
#> [961] 7.343500e+01 2.470439e+02 1.983948e+03 2.191546e+02 1.038706e+02
#> [966] 2.090950e+02 8.053951e+01 1.454525e+02 7.490694e+01 9.885491e+01
#> [971] 8.159932e+01 7.095687e+01 7.177422e+01 1.181418e+02 1.784391e+02
#> [976] 1.040599e+02 7.078397e+01 1.058582e+02 8.092210e+01 8.486704e+01
#> [981] 6.320881e+01 6.652482e+01 6.536560e+01 9.660820e+01 2.477924e+02
#> [986] 1.183795e+02 3.865511e+02 8.007697e+01 1.368023e+02 1.444333e+02
#> [991] 6.676408e+01 6.653719e+01 7.212327e+01 6.475771e+01 6.911713e+01
#> [996] 6.460444e+01 7.637223e+01 7.118992e+01 1.831014e+02 1.825801e+02
doseFunction(
model,
alpha2 = samples@data$alpha2,
beta = samples@data$beta,
grade = 2L
)(x = 0.25)
#> [1] 5.750097e+01 6.905208e+01 5.862501e+01 6.885496e+01 6.520952e+01
#> [6] 6.637755e+01 5.357084e+01 5.838230e+01 1.077903e+02 8.105612e+01
#> [11] 8.693553e+01 8.230005e+01 6.634015e+01 5.885940e+01 6.354538e+01
#> [16] 6.047404e+01 6.310401e+01 7.568066e+01 6.268394e+01 6.196988e+01
#> [21] 8.039279e+01 2.959367e+02 6.585953e+01 7.194233e+01 9.933099e+01
#> [26] 7.687623e+01 6.195436e+01 5.631347e+01 9.824385e+01 7.107889e+01
#> [31] 5.455070e+01 8.052008e+01 1.813078e+02 8.752797e+01 1.692821e+02
#> [36] 7.721929e+01 8.735083e+01 8.638192e+01 9.702456e+01 5.685073e+01
#> [41] 5.959636e+01 6.081750e+01 6.980339e+01 8.361843e+01 7.606354e+01
#> [46] 6.450013e+01 5.915231e+01 1.142470e+02 5.633502e+01 6.861042e+01
#> [51] 6.304729e+01 1.591535e+02 3.537118e+02 7.641842e+01 6.947640e+01
#> [56] 5.622497e+01 8.842784e+01 8.693770e+01 5.684784e+01 1.789337e+02
#> [61] 8.030936e+01 5.813058e+01 6.951979e+01 5.302321e+01 6.443318e+01
#> [66] 8.138100e+01 5.996056e+01 7.295915e+01 6.876640e+01 6.244566e+01
#> [71] 6.960852e+01 6.902911e+01 7.773111e+01 5.736576e+01 6.147777e+01
#> [76] 4.237784e+02 5.457965e+01 6.393175e+01 6.051663e+01 7.804569e+01
#> [81] 6.109401e+01 5.582600e+01 6.700049e+01 6.513953e+01 8.284439e+01
#> [86] 5.778340e+01 6.772665e+01 6.382714e+01 5.527143e+01 5.631327e+01
#> [91] 5.577624e+01 1.133662e+02 6.691689e+01 6.203440e+01 6.710856e+01
#> [96] 6.634164e+01 5.975095e+01 6.299064e+01 3.821601e+02 2.673367e+02
#> [101] 2.489155e+02 6.021796e+01 7.503771e+01 6.128699e+01 9.347284e+01
#> [106] 6.308048e+01 5.048895e+01 5.532036e+01 6.414678e+01 1.100397e+02
#> [111] 1.154406e+02 6.788817e+01 5.753956e+01 1.291788e+02 7.154790e+01
#> [116] 5.459128e+01 6.705737e+01 5.532418e+01 7.316256e+01 5.307521e+01
#> [121] 5.157658e+01 6.391491e+01 6.264460e+01 5.947099e+01 7.612110e+01
#> [126] 8.160572e+01 8.364741e+01 4.356164e+01 1.540200e+02 3.236052e+03
#> [131] 5.721972e+12 5.425904e+01 5.262036e+01 1.193213e+02 5.378763e+01
#> [136] 5.460135e+01 5.513124e+01 5.722996e+01 5.787290e+01 7.042859e+01
#> [141] 8.920809e+01 5.731733e+01 5.931575e+01 5.906067e+01 5.790732e+01
#> [146] 2.763133e+02 6.538199e+01 8.341276e+01 5.508806e+01 5.988585e+01
#> [151] 7.793528e+01 6.052786e+01 6.005056e+01 8.037592e+01 9.483147e+01
#> [156] 5.756935e+01 8.166239e+01 6.635603e+01 1.288594e+02 8.256890e+01
#> [161] 6.373137e+01 7.227723e+01 5.623132e+01 6.897019e+01 6.356584e+01
#> [166] 5.817687e+01 2.064250e+02 7.576383e+01 1.217704e+02 5.406967e+01
#> [171] 7.613063e+01 6.193654e+01 6.017034e+01 6.349290e+01 6.299444e+01
#> [176] 1.988023e+02 5.981381e+01 1.926221e+02 6.125417e+01 6.469143e+01
#> [181] 6.208979e+01 6.158480e+01 5.862977e+01 6.016554e+01 5.906104e+01
#> [186] 7.215837e+01 6.086439e+01 6.034682e+01 6.342286e+01 6.565834e+01
#> [191] 1.161460e+02 5.853882e+01 6.990522e+01 6.967307e+01 6.303966e+01
#> [196] 9.315612e+01 6.882870e+01 7.323890e+01 6.772883e+01 5.701863e+01
#> [201] 1.104822e+02 6.079885e+01 1.294221e+02 1.857351e+02 3.584668e+03
#> [206] 1.169704e+02 2.091859e+02 7.032765e+01 1.062474e+02 6.336259e+01
#> [211] 6.213970e+01 7.866411e+01 9.185105e+01 6.672604e+01 9.999390e+01
#> [216] 5.260637e+01 1.405897e+02 3.209846e+02 5.422841e+01 1.122330e+02
#> [221] 8.198962e+01 5.942900e+01 5.824050e+01 6.784618e+01 6.118888e+01
#> [226] 5.953962e+01 6.945671e+01 6.597037e+01 8.874377e+01 6.406565e+01
#> [231] 8.176084e+01 8.673664e+01 5.096989e+01 5.059370e+01 8.083176e+01
#> [236] 7.848705e+01 6.235952e+01 5.895508e+01 8.491261e+01 4.745196e+01
#> [241] 7.607717e+01 5.850461e+01 2.981216e+02 3.997736e+03 1.801361e+02
#> [246] 6.277926e+01 6.380878e+01 8.004415e+01 5.642236e+01 5.789521e+01
#> [251] 8.218123e+03 1.161653e+02 5.915919e+01 7.266077e+01 7.845765e+01
#> [256] 1.082994e+02 5.609470e+01 6.096276e+01 8.283546e+01 6.423925e+01
#> [261] 5.961018e+01 7.360403e+01 6.242545e+01 8.045823e+01 6.703592e+01
#> [266] 5.572262e+01 1.536867e+02 6.068540e+01 7.220526e+01 7.590921e+01
#> [271] 8.397684e+01 1.020441e+02 3.090788e+02 5.372883e+01 8.206074e+01
#> [276] 2.792821e+02 5.914369e+01 6.702638e+01 5.356754e+01 1.468503e+02
#> [281] 8.196111e+01 5.602156e+01 6.208409e+01 5.931343e+01 5.853496e+01
#> [286] 6.917117e+01 6.466075e+01 6.677391e+01 1.105263e+02 6.395769e+01
#> [291] 8.484975e+01 6.847120e+01 8.807686e+01 5.933227e+01 1.671751e+02
#> [296] 5.858392e+01 1.044228e+02 5.703247e+01 1.221762e+02 1.292848e+02
#> [301] 7.637533e+01 5.772842e+01 5.710264e+01 1.187423e+02 6.503128e+01
#> [306] 5.894587e+01 1.262185e+02 6.023692e+01 5.300626e+01 5.976167e+01
#> [311] 8.256590e+01 7.099000e+01 5.963683e+01 5.624443e+01 7.196488e+01
#> [316] 5.413996e+01 6.312567e+02 2.194663e+02 5.794579e+01 1.132280e+02
#> [321] 6.281211e+01 5.702606e+01 5.535786e+01 7.454047e+01 5.747279e+01
#> [326] 8.593648e+01 6.674347e+01 8.479430e+01 6.307731e+01 6.618486e+01
#> [331] 5.699988e+01 7.897156e+01 5.056340e+01 7.307241e+01 6.288871e+01
#> [336] 5.552070e+01 6.220103e+01 6.001643e+01 6.272634e+01 1.211552e+02
#> [341] 5.567733e+01 7.033356e+01 6.188317e+01 5.909151e+01 1.946712e+02
#> [346] 6.537661e+01 7.656115e+01 7.176574e+01 5.919716e+01 5.935046e+01
#> [351] 6.837558e+01 9.765118e+01 6.996080e+01 6.434804e+01 6.487180e+01
#> [356] 8.932988e+01 7.155190e+01 6.717840e+01 5.607058e+01 6.959849e+01
#> [361] 5.843289e+01 6.721603e+01 7.936749e+01 5.359054e+01 9.598356e+01
#> [366] 8.670401e+01 7.613126e+01 6.096026e+01 9.232166e+01 5.907531e+01
#> [371] 8.652763e+01 1.105997e+02 5.184737e+01 6.095270e+01 7.940152e+01
#> [376] 9.521642e+01 7.699566e+01 6.467568e+01 6.732404e+01 6.125307e+01
#> [381] 7.631363e+01 1.061808e+02 2.797844e+02 5.392420e+01 6.814049e+01
#> [386] 6.170547e+01 9.028274e+01 7.810378e+01 7.517360e+01 1.530453e+02
#> [391] 5.258466e+01 9.150526e+01 8.319141e+01 6.697851e+01 5.572553e+01
#> [396] 6.397352e+01 5.884703e+01 6.627617e+01 7.060459e+01 1.034284e+02
#> [401] 6.566559e+01 6.907494e+01 5.858697e+01 5.456375e+01 8.682779e+01
#> [406] 6.904841e+01 6.246490e+01 5.840820e+01 8.446572e+01 5.844473e+01
#> [411] 6.689178e+01 6.857630e+01 8.011337e+01 6.080802e+01 1.923931e+02
#> [416] 7.410663e+01 6.882474e+01 5.484044e+01 7.493267e+01 6.395681e+01
#> [421] 6.545522e+01 5.846296e+01 6.110477e+01 8.380800e+01 1.824443e+02
#> [426] 6.233328e+01 8.528734e+01 7.431514e+01 5.842225e+01 1.904776e+02
#> [431] 1.236007e+02 5.753035e+01 7.264715e+01 6.057745e+01 6.161960e+01
#> [436] 5.961766e+01 7.458851e+01 6.927417e+01 6.268693e+01 1.296176e+02
#> [441] 7.581639e+01 5.519139e+01 6.272567e+01 7.348062e+01 7.133513e+01
#> [446] 6.152738e+01 5.241004e+01 7.379881e+01 6.057305e+01 6.367085e+01
#> [451] 7.028220e+01 6.106569e+01 6.786944e+01 6.370742e+01 8.101769e+01
#> [456] 7.075098e+01 6.838220e+01 5.546765e+01 6.284137e+01 8.304857e+01
#> [461] 6.149052e+01 6.103249e+01 6.959264e+01 6.422249e+01 5.166704e+01
#> [466] 7.691257e+01 1.350813e+02 1.697115e+02 6.034808e+01 7.242581e+01
#> [471] 6.760800e+01 1.595088e+02 5.812589e+01 5.527136e+01 8.503597e+01
#> [476] 5.751863e+01 6.264421e+01 5.847832e+01 7.137924e+01 6.022808e+01
#> [481] 6.426898e+01 5.548703e+01 8.887730e+01 1.499267e+02 5.848785e+01
#> [486] 7.250561e+01 5.493098e+01 1.633643e+02 6.299769e+01 5.465516e+01
#> [491] 5.732905e+01 5.911121e+01 5.997854e+01 6.726479e+01 6.633765e+01
#> [496] 1.398905e+02 7.140999e+01 5.499966e+01 6.205207e+01 5.674596e+01
#> [501] 6.006350e+01 7.757661e+01 6.803341e+01 6.911387e+01 6.455849e+01
#> [506] 7.760388e+01 6.149605e+01 5.014916e+01 6.604677e+01 5.463821e+01
#> [511] 1.284675e+02 6.703143e+01 5.946769e+01 6.600371e+01 6.142595e+01
#> [516] 7.041582e+01 5.642078e+01 7.111442e+01 6.811972e+01 5.987003e+01
#> [521] 6.281118e+01 7.797972e+01 6.123902e+01 6.402142e+01 5.793047e+01
#> [526] 6.032378e+01 6.079351e+01 1.610428e+02 5.382315e+01 1.016792e+02
#> [531] 5.692385e+01 6.997736e+01 8.513830e+01 6.096042e+01 5.759411e+01
#> [536] 6.006668e+01 6.012965e+01 6.258192e+01 6.301089e+01 6.888164e+01
#> [541] 6.161139e+01 6.119056e+01 6.073656e+01 6.287306e+01 7.313245e+01
#> [546] 5.901962e+01 5.790582e+01 6.089762e+01 7.353354e+01 8.274251e+01
#> [551] 6.197797e+01 5.703830e+01 1.184290e+02 5.587002e+01 6.411801e+01
#> [556] 7.781720e+01 8.498770e+01 6.229426e+01 7.624219e+01 5.699558e+01
#> [561] 5.991152e+01 2.491967e+02 5.787961e+01 5.955861e+01 5.888264e+01
#> [566] 8.793017e+01 1.088813e+02 7.304647e+01 6.031935e+01 6.686681e+01
#> [571] 1.074772e+02 5.512886e+01 6.490915e+01 6.490141e+01 1.714213e+02
#> [576] 6.217875e+01 1.407718e+02 1.288619e+02 6.313451e+01 4.830473e+01
#> [581] 5.990966e+01 8.195203e+01 6.269339e+01 6.777617e+01 5.737133e+01
#> [586] 6.255997e+01 6.358700e+01 7.251910e+01 6.423286e+01 6.563285e+01
#> [591] 5.513490e+01 5.668064e+01 5.249410e+01 7.315384e+01 6.432068e+01
#> [596] 1.163150e+02 5.923221e+01 4.034793e+02 8.147045e+01 7.285701e+01
#> [601] 7.537619e+01 5.549707e+01 1.262162e+02 7.289351e+01 7.245572e+01
#> [606] 7.145297e+01 5.633946e+01 5.863188e+01 7.771951e+01 1.278983e+02
#> [611] 5.993464e+01 7.859682e+01 6.474780e+01 1.504349e+02 1.041091e+02
#> [616] 2.289313e+02 7.631924e+01 6.411005e+01 5.479528e+01 5.944361e+01
#> [621] 5.747267e+01 8.191409e+01 5.562160e+01 1.557454e+02 6.554888e+01
#> [626] 5.629440e+01 8.545506e+01 6.384987e+01 7.538669e+01 6.972898e+01
#> [631] 9.021436e+01 8.223099e+01 7.261653e+01 5.992410e+01 1.016117e+02
#> [636] 7.105754e+01 7.153801e+01 1.808860e+02 6.815650e+01 5.351887e+01
#> [641] 1.645597e+02 5.348598e+01 5.860693e+01 9.045381e+01 7.088789e+01
#> [646] 7.426481e+01 1.088567e+02 5.640502e+01 6.060057e+02 7.800105e+01
#> [651] 9.188928e+01 6.045414e+01 1.354000e+02 6.899801e+01 6.737710e+01
#> [656] 6.779093e+01 6.124962e+01 7.133217e+01 5.552298e+01 8.063020e+01
#> [661] 6.544122e+01 7.679980e+01 6.061970e+01 7.047319e+01 6.126070e+01
#> [666] 5.912267e+01 6.759162e+01 6.097802e+01 6.155754e+01 5.852715e+01
#> [671] 8.563178e+01 5.962259e+01 2.746880e+02 5.509911e+01 6.381037e+01
#> [676] 7.585073e+01 7.972812e+01 5.827009e+01 6.442808e+01 5.992282e+01
#> [681] 6.363550e+01 7.389441e+01 9.440320e+01 5.764496e+01 9.281924e+01
#> [686] 1.575996e+02 1.238437e+02 1.942443e+02 5.697248e+01 6.214615e+01
#> [691] 5.274364e+01 6.512554e+01 9.443650e+01 6.662099e+01 7.096578e+01
#> [696] 2.385999e+02 6.223836e+01 5.747956e+01 5.364952e+01 6.067507e+01
#> [701] 9.474818e+01 1.480974e+02 6.494326e+01 5.840294e+01 7.415821e+01
#> [706] 5.473271e+01 5.531656e+01 5.978899e+01 7.254197e+01 9.254633e+01
#> [711] 6.815390e+01 6.561659e+01 6.973229e+01 7.146608e+01 5.678258e+01
#> [716] 6.756751e+01 1.000899e+02 6.205087e+01 5.910482e+01 5.622632e+01
#> [721] 6.481393e+01 6.247297e+01 7.875964e+01 6.353698e+01 5.882103e+01
#> [726] 6.740945e+01 6.533255e+01 6.105863e+01 5.949728e+01 5.699122e+01
#> [731] 6.836259e+01 5.735757e+01 7.280021e+01 2.192948e+02 2.437654e+02
#> [736] 7.648987e+01 6.954471e+01 5.587050e+01 1.004727e+03 5.423707e+01
#> [741] 8.660078e+01 6.683764e+01 8.001460e+01 8.084246e+01 6.195741e+01
#> [746] 6.288628e+01 5.805661e+01 1.069402e+02 6.068072e+01 6.581179e+01
#> [751] 6.922926e+01 8.512874e+01 1.045481e+02 7.930431e+01 5.921181e+01
#> [756] 5.732343e+01 6.402986e+01 5.727978e+01 1.187945e+02 5.673201e+01
#> [761] 7.992325e+01 9.016495e+01 7.252729e+01 9.024616e+01 5.385752e+01
#> [766] 2.886912e+04 6.668946e+01 6.759817e+01 6.772143e+01 8.281392e+01
#> [771] 5.844038e+01 5.526250e+01 7.441725e+01 6.453319e+01 6.994236e+01
#> [776] 8.631616e+01 1.984850e+02 4.643732e+02 7.135136e+01 7.363858e+01
#> [781] 6.477644e+01 5.106014e+01 6.050062e+01 7.864924e+01 6.966797e+01
#> [786] 7.560702e+01 5.957048e+01 1.022695e+02 6.287538e+01 9.759531e+01
#> [791] 6.178473e+01 6.340216e+01 2.604234e+02 8.408478e+01 8.113818e+01
#> [796] 6.254835e+01 1.070116e+02 1.731928e+02 7.909591e+01 5.824447e+01
#> [801] 2.158975e+02 1.134890e+02 1.311763e+02 5.881559e+01 5.988864e+01
#> [806] 8.287383e+01 6.571618e+01 7.939008e+01 1.081381e+02 1.203288e+02
#> [811] 2.344186e+02 2.510803e+03 5.220885e+01 5.733688e+01 6.064438e+01
#> [816] 5.788124e+01 5.656110e+01 5.908213e+01 6.000842e+01 6.664463e+01
#> [821] 6.819686e+01 1.039624e+02 5.844696e+01 7.692556e+01 5.950796e+01
#> [826] 6.680795e+01 7.268625e+01 6.749457e+01 6.240402e+01 8.829538e+01
#> [831] 5.625133e+01 8.535767e+01 5.916132e+01 5.964405e+01 5.843651e+01
#> [836] 6.110212e+01 5.586095e+01 6.819820e+01 6.935335e+01 5.853570e+01
#> [841] 7.827505e+01 8.873868e+01 6.788946e+01 5.961429e+01 8.037891e+01
#> [846] 6.142994e+01 6.870439e+01 6.664920e+01 7.633698e+01 2.401817e+02
#> [851] 9.742335e+01 5.018914e+01 5.299219e+01 9.246830e+01 8.413306e+01
#> [856] 1.142412e+02 6.550325e+01 6.574441e+01 1.021532e+02 7.219405e+01
#> [861] 1.785410e+02 5.602427e+01 7.450641e+01 1.339497e+02 8.047676e+01
#> [866] 5.871965e+01 6.851750e+01 6.399181e+01 9.317775e+01 7.281695e+01
#> [871] 6.916503e+01 5.672222e+01 5.655160e+01 7.118559e+01 8.682213e+01
#> [876] 6.937283e+01 9.173368e+01 5.705931e+01 6.372039e+01 5.666968e+01
#> [881] 5.378696e+01 6.160077e+01 6.513749e+01 7.831385e+01 5.961948e+01
#> [886] 8.546231e+01 8.553542e+01 2.559814e+02 8.405223e+01 7.045117e+01
#> [891] 5.570059e+01 8.475761e+01 1.631746e+02 5.976643e+01 7.879860e+01
#> [896] 7.897003e+01 7.768190e+01 6.048705e+01 8.570936e+01 6.581326e+01
#> [901] 7.183760e+01 7.045504e+01 1.504294e+02 1.004762e+02 5.691336e+01
#> [906] 1.250857e+03 1.125932e+02 2.205095e+02 7.991354e+01 5.916582e+01
#> [911] 8.905635e+01 6.746187e+01 6.619988e+01 6.155888e+01 5.894087e+01
#> [916] 7.346204e+01 5.083576e+01 6.014504e+01 5.876669e+01 6.990824e+01
#> [921] 7.344997e+01 6.150953e+01 5.685950e+01 6.489094e+01 6.307866e+01
#> [926] 6.171204e+01 7.213487e+01 7.341534e+01 3.854946e+02 5.436348e+01
#> [931] 5.579868e+02 6.794449e+01 1.112147e+02 7.207768e+01 6.338127e+01
#> [936] 5.461117e+01 7.847798e+01 5.906013e+01 5.961738e+01 6.251685e+01
#> [941] 8.092244e+01 5.859269e+01 8.965371e+01 6.292588e+01 6.994320e+01
#> [946] 6.642988e+01 5.316046e+01 6.918807e+01 5.861160e+01 1.428548e+02
#> [951] 7.416649e+01 5.923850e+01 9.436687e+01 5.909886e+01 6.388417e+01
#> [956] 8.183171e+01 6.078213e+01 7.583277e+01 7.131146e+01 7.348533e+01
#> [961] 6.229990e+01 1.241047e+02 4.972022e+02 9.815868e+01 7.621049e+01
#> [966] 8.610499e+01 6.321171e+01 8.536104e+01 6.043922e+01 6.925706e+01
#> [971] 6.289344e+01 6.202164e+01 6.320443e+01 7.597630e+01 7.111529e+01
#> [976] 6.526563e+01 5.955804e+01 7.440796e+01 5.557118e+01 6.563411e+01
#> [981] 5.436456e+01 5.947912e+01 5.699367e+01 5.144810e+01 8.804387e+01
#> [986] 7.960783e+01 9.662003e+01 5.348039e+01 6.049097e+01 7.552954e+01
#> [991] 5.904407e+01 5.962193e+01 6.544423e+01 5.862421e+01 6.083524e+01
#> [996] 5.396469e+01 5.900474e+01 6.258204e+01 1.029437e+02 9.733766e+01