![[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] 1.537783e+02 3.191120e+02 5.131526e+02 1.002257e+02 7.908856e+01
#> [6] 7.329636e+01 7.425253e+01 7.974919e+01 1.349595e+02 9.222898e+01
#> [11] 9.865042e+01 9.272385e+01 6.968617e+01 1.662612e+02 1.038210e+02
#> [16] 7.602771e+01 1.610981e+02 2.038920e+03 3.877217e+02 1.034791e+02
#> [21] 6.954767e+01 6.368914e+01 5.995594e+01 6.187090e+01 6.674952e+01
#> [26] 8.165682e+01 6.844388e+01 9.398843e+01 1.022576e+02 4.882925e+02
#> [31] 6.948932e+01 6.919945e+01 6.564865e+01 9.374500e+01 5.339089e+02
#> [36] 6.570380e+03 8.567471e+01 8.271169e+01 7.560143e+01 6.756809e+01
#> [41] 6.923754e+01 1.007834e+02 6.811873e+01 8.500441e+01 9.588043e+01
#> [46] 9.432386e+01 3.648102e+02 7.155909e+01 1.160990e+02 8.029892e+01
#> [51] 7.642666e+01 3.264990e+02 1.115659e+02 8.540878e+01 9.660861e+01
#> [56] 8.461954e+01 1.396728e+02 9.924042e+01 6.583189e+01 8.753898e+01
#> [61] 7.972212e+01 1.610037e+02 8.963029e+01 1.304387e+02 7.462416e+01
#> [66] 6.645334e+01 6.275072e+01 8.503303e+01 2.433415e+02 5.926935e+01
#> [71] 6.359563e+01 7.060168e+01 7.841869e+01 1.241044e+02 9.839752e+01
#> [76] 8.224060e+01 8.068142e+01 4.822254e+02 3.764073e+02 7.792342e+01
#> [81] 1.350084e+02 7.219689e+01 1.993716e+02 1.088347e+02 1.130098e+02
#> [86] 7.561923e+01 6.646325e+01 8.493935e+01 8.678742e+01 6.546099e+01
#> [91] 7.019029e+01 6.687610e+01 9.579908e+01 7.682972e+01 1.051975e+02
#> [96] 1.067303e+02 7.860867e+01 1.460199e+02 7.053835e+01 1.655395e+02
#> [101] 2.932280e+02 1.559929e+02 9.122439e+01 7.587648e+01 7.857198e+01
#> [106] 8.712052e+01 7.479648e+01 6.295694e+01 1.164826e+02 8.121985e+01
#> [111] 1.380817e+02 1.428807e+08 1.378287e+02 1.050884e+02 1.878645e+02
#> [116] 8.527982e+01 6.257645e+01 1.104880e+02 7.595543e+01 1.561877e+02
#> [121] 6.216290e+02 8.426231e+01 7.415965e+01 7.839033e+01 7.746712e+01
#> [126] 1.810640e+04 9.659307e+01 2.261543e+03 2.002967e+02 1.842622e+02
#> [131] 9.236004e+01 6.535142e+01 2.618173e+02 6.306028e+03 1.443506e+02
#> [136] 8.113654e+01 8.000682e+01 7.328763e+01 7.792699e+01 9.910200e+01
#> [141] 7.720870e+01 6.601775e+01 1.280207e+02 6.856908e+01 7.576134e+01
#> [146] 6.724666e+01 9.101179e+01 2.113412e+02 1.880404e+02 2.650786e+02
#> [151] 7.966924e+01 8.454333e+01 8.080787e+01 6.421219e+01 1.267155e+02
#> [156] 9.448800e+01 1.542686e+02 2.145052e+02 7.135624e+01 1.699678e+02
#> [161] 1.472185e+02 6.444633e+01 1.954292e+02 6.908416e+01 7.616584e+01
#> [166] 1.041252e+02 6.882806e+01 6.990041e+01 7.660330e+01 9.102849e+01
#> [171] 8.640297e+01 1.444974e+02 1.440392e+02 8.201388e+01 7.335079e+01
#> [176] 1.126103e+02 9.098581e+01 9.938722e+01 1.902272e+03 2.418675e+02
#> [181] 2.231587e+03 5.552480e+02 7.114813e+01 2.847264e+02 1.109871e+02
#> [186] 5.866769e+01 8.039731e+01 6.475808e+01 4.205587e+02 7.571154e+01
#> [191] 7.916493e+01 7.225215e+01 6.729163e+01 1.010067e+02 8.267532e+01
#> [196] 2.382368e+02 7.461983e+01 9.431829e+01 7.437535e+01 4.185706e+02
#> [201] 8.221811e+01 1.335520e+02 7.848183e+01 7.760195e+01 7.491576e+01
#> [206] 9.611888e+01 7.827308e+01 1.612691e+02 6.856147e+01 8.236323e+01
#> [211] 3.004521e+02 1.050804e+02 6.608112e+01 9.710053e+01 1.128713e+04
#> [216] 3.517210e+03 1.489661e+02 7.487586e+01 6.837848e+01 6.495737e+01
#> [221] 8.431643e+01 9.861909e+01 1.735476e+02 2.637275e+02 9.628609e+01
#> [226] 1.078814e+02 7.822860e+01 1.420956e+02 9.465565e+01 6.861270e+01
#> [231] 7.674735e+01 9.736349e+01 7.036005e+01 8.265923e+01 7.793648e+01
#> [236] 4.397721e+02 9.835314e+01 6.906536e+01 8.473322e+01 7.766430e+01
#> [241] 4.269974e+03 8.365780e+01 9.486597e+01 7.101071e+01 8.094832e+01
#> [246] 1.187255e+02 7.545452e+01 1.581186e+02 2.198975e+02 8.470045e+01
#> [251] 1.037586e+02 7.835686e+01 6.555898e+01 4.790431e+03 3.217711e+02
#> [256] 2.396764e+02 4.190503e+02 7.308758e+01 1.514158e+02 8.407325e+01
#> [261] 7.151838e+01 7.145052e+01 2.695898e+02 7.567935e+01 7.105985e+01
#> [266] 6.287197e+01 6.610178e+01 6.312840e+01 6.744417e+01 7.185650e+01
#> [271] 1.050821e+02 8.661789e+01 8.347781e+01 8.601059e+01 7.213954e+01
#> [276] 7.135440e+01 7.474252e+01 1.049617e+02 6.846446e+01 6.328681e+01
#> [281] 6.891013e+01 8.044303e+01 9.677843e+01 1.178305e+02 1.133302e+02
#> [286] 1.077351e+02 6.624500e+01 8.479506e+01 1.623983e+02 7.161879e+01
#> [291] 8.633709e+01 7.657114e+01 2.697017e+02 7.200826e+01 7.151154e+01
#> [296] 7.622299e+01 6.961336e+01 7.733121e+01 5.959752e+01 1.367039e+02
#> [301] 8.819365e+01 6.850689e+01 1.057262e+03 6.948883e+02 1.549133e+02
#> [306] 1.095803e+02 1.111056e+02 1.550419e+02 6.327079e+01 1.161713e+02
#> [311] 7.224486e+01 9.038058e+01 7.345721e+01 2.227076e+02 1.090236e+02
#> [316] 9.078690e+01 1.288641e+02 8.655262e+02 1.040129e+02 1.191589e+02
#> [321] 1.324840e+02 8.019275e+01 1.219320e+02 6.172603e+02 5.383210e+02
#> [326] 8.233905e+01 7.133030e+01 1.414773e+02 7.082119e+01 1.353690e+02
#> [331] 7.144229e+01 7.183440e+01 6.510504e+01 1.835498e+02 6.641271e+01
#> [336] 4.332955e+06 3.491809e+03 1.721031e+02 3.156933e+02 2.170179e+02
#> [341] 4.169854e+02 6.765147e+01 8.146432e+01 8.178888e+01 6.663059e+01
#> [346] 8.001934e+01 4.468825e+02 6.810493e+01 6.468137e+01 7.978344e+01
#> [351] 1.147241e+02 9.126161e+02 5.221081e+02 8.020946e+01 1.411673e+02
#> [356] 6.351666e+01 6.423911e+01 6.779735e+01 2.639991e+02 7.970651e+01
#> [361] 7.634813e+01 1.757456e+02 3.638827e+02 1.354959e+02 6.630216e+01
#> [366] 6.552501e+01 8.317048e+01 1.133785e+02 2.844910e+02 6.970984e+01
#> [371] 1.192107e+02 2.750613e+02 2.053883e+04 8.522463e+03 2.361193e+02
#> [376] 7.930694e+01 1.456654e+02 7.977349e+01 8.648151e+01 8.462235e+01
#> [381] 1.025030e+02 2.184815e+02 9.114859e+01 1.538511e+02 4.718565e+02
#> [386] 6.119551e+01 1.903374e+02 8.242634e+01 1.813859e+02 7.666667e+01
#> [391] 8.739498e+01 7.095067e+01 1.217282e+02 2.826944e+02 1.125466e+03
#> [396] 4.196329e+02 1.143952e+02 2.584505e+02 3.083107e+02 7.747078e+01
#> [401] 6.485893e+01 7.476262e+02 3.479889e+02 1.267449e+02 8.898902e+01
#> [406] 9.044240e+01 9.962366e+01 6.142672e+01 7.127757e+01 9.549119e+01
#> [411] 8.817238e+01 6.792712e+01 8.897497e+02 6.912164e+01 1.019380e+02
#> [416] 1.036436e+02 7.827023e+01 1.483853e+02 1.031133e+02 9.917119e+01
#> [421] 6.410440e+01 2.746079e+02 6.937386e+01 1.184974e+02 3.017127e+02
#> [426] 6.218883e+01 6.917686e+01 1.006599e+02 7.655078e+01 6.173360e+01
#> [431] 6.280067e+01 7.265254e+01 7.755781e+01 7.414042e+01 8.104226e+01
#> [436] 5.480176e+02 1.478089e+02 3.395247e+02 6.723947e+01 8.348355e+01
#> [441] 8.720543e+01 9.043413e+01 1.852051e+02 6.402468e+01 1.667386e+02
#> [446] 1.122609e+02 1.086781e+02 1.008716e+02 8.564630e+01 7.500296e+01
#> [451] 9.262753e+03 2.137657e+02 1.511819e+02 6.129021e+02 1.432631e+03
#> [456] 8.704829e+01 8.596507e+01 2.422775e+02 8.788937e+01 7.911922e+01
#> [461] 1.526048e+02 4.733251e+02 1.429829e+02 7.692566e+01 1.326465e+02
#> [466] 1.091438e+02 7.256597e+01 7.954769e+01 1.338075e+02 5.654091e+02
#> [471] 1.128080e+02 1.495249e+02 6.721594e+01 1.926350e+02 9.951786e+01
#> [476] 9.100249e+01 6.697212e+01 1.404042e+02 1.909659e+02 1.226667e+02
#> [481] 8.273953e+01 4.415973e+02 1.087997e+02 1.023640e+02 1.036550e+02
#> [486] 1.028953e+02 8.110428e+01 9.907209e+01 1.110297e+02 1.471976e+02
#> [491] 6.636260e+01 1.140099e+02 1.627142e+02 7.975052e+01 7.137977e+01
#> [496] 9.223536e+01 7.157735e+01 2.262896e+02 2.403520e+02 2.312532e+02
#> [501] 2.995545e+02 7.062126e+01 1.999303e+02 6.170016e+01 8.577505e+01
#> [506] 8.250780e+01 9.848783e+01 7.843186e+01 8.708347e+01 6.608402e+01
#> [511] 1.146933e+02 1.079493e+02 8.224897e+01 9.586710e+01 6.044946e+02
#> [516] 9.594114e+01 1.129976e+02 1.061457e+02 9.793825e+01 3.518713e+02
#> [521] 3.519689e+02 1.064227e+02 1.061906e+02 2.615024e+02 1.121256e+02
#> [526] 1.546002e+02 7.448878e+03 8.120659e+01 9.075855e+01 1.109298e+02
#> [531] 9.293525e+01 8.592456e+01 1.678854e+02 7.467861e+01 7.452073e+01
#> [536] 6.747808e+01 8.023552e+01 8.293063e+01 1.403664e+02 6.426990e+01
#> [541] 7.783573e+01 7.393723e+01 7.678748e+01 1.498351e+02 3.334040e+04
#> [546] 9.793840e+03 6.775214e+06 9.135795e+01 7.094709e+01 7.076397e+01
#> [551] 9.040903e+01 1.454297e+02 7.895440e+01 9.080486e+01 6.881568e+01
#> [556] 6.262893e+01 6.148264e+01 6.523398e+01 8.155703e+01 1.283421e+02
#> [561] 6.790516e+01 6.292289e+01 3.891483e+02 3.908112e+03 1.443885e+02
#> [566] 1.126638e+02 7.831187e+01 2.078414e+02 7.665345e+01 9.613071e+01
#> [571] 8.342614e+01 3.300804e+02 9.078440e+01 8.301670e+01 8.984839e+01
#> [576] 7.050744e+01 7.699295e+01 8.663294e+01 7.181331e+01 8.419010e+01
#> [581] 3.743204e+03 9.341508e+01 9.126851e+01 7.097095e+01 6.821629e+01
#> [586] 8.540778e+01 1.100102e+02 6.254995e+01 6.811210e+01 8.035183e+01
#> [591] 8.459516e+03 3.292221e+02 8.591303e+01 9.622928e+01 7.612678e+01
#> [596] 8.133756e+01 6.378495e+01 7.437404e+01 8.527616e+01 8.269324e+01
#> [601] 1.115918e+02 1.445903e+02 2.487809e+02 1.762039e+02 2.106090e+02
#> [606] 7.170749e+01 7.219604e+01 6.953462e+01 9.493513e+01 4.056326e+02
#> [611] 1.219197e+02 1.613211e+02 6.560251e+01 8.029952e+01 6.376593e+01
#> [616] 1.005724e+02 6.860565e+01 4.138859e+06 8.493843e+01 9.226757e+01
#> [621] 8.402084e+01 8.985535e+01 7.937501e+01 6.528136e+01 6.166532e+01
#> [626] 8.292170e+01 1.084628e+02 8.441953e+01 8.690720e+01 6.995477e+01
#> [631] 8.920109e+01 3.503281e+03 1.362494e+04 7.695189e+01 7.461872e+01
#> [636] 6.303793e+01 7.740799e+01 1.716436e+02 1.181221e+02 1.449381e+03
#> [641] 6.466393e+01 7.378083e+01 9.004621e+01 6.859850e+01 7.471622e+01
#> [646] 8.171030e+01 8.254913e+01 8.280631e+01 1.367560e+02 6.891488e+01
#> [651] 9.080573e+01 9.280431e+01 1.174619e+02 6.785525e+01 6.693317e+01
#> [656] 7.343006e+01 6.952791e+01 8.869595e+01 9.087756e+01 7.747896e+01
#> [661] 9.683892e+01 1.013253e+02 6.388044e+01 8.412116e+02 7.763523e+01
#> [666] 7.295698e+01 8.253559e+01 7.394480e+01 7.916496e+01 1.206350e+02
#> [671] 1.064549e+02 8.579058e+01 7.519810e+01 6.883828e+01 8.068489e+01
#> [676] 6.516656e+01 6.433321e+01 6.904394e+01 7.123737e+01 2.724179e+02
#> [681] 7.562237e+01 1.221664e+02 8.564403e+01 7.603154e+01 5.164394e+02
#> [686] 1.612958e+02 6.800943e+01 8.527040e+01 2.917551e+02 7.740563e+01
#> [691] 7.365079e+01 7.531264e+01 7.775541e+02 7.570914e+04 1.055111e+02
#> [696] 9.383480e+01 1.051867e+02 8.034568e+01 9.015042e+01 6.326505e+02
#> [701] 1.974759e+05 1.752150e+02 6.845538e+01 8.476984e+01 1.228388e+02
#> [706] 1.578718e+02 8.356185e+01 1.246693e+02 1.105476e+02 7.751631e+01
#> [711] 6.721544e+01 7.969411e+01 7.080109e+01 9.867443e+01 6.710005e+01
#> [716] 6.776111e+01 7.273388e+01 7.419897e+01 7.402002e+01 1.116302e+02
#> [721] 1.458377e+02 7.218606e+01 2.742362e+02 2.101873e+02 1.116535e+02
#> [726] 1.523225e+02 1.100888e+02 7.715737e+01 7.192159e+01 6.831379e+01
#> [731] 8.637678e+01 6.411071e+01 7.856623e+01 8.253623e+01 7.154356e+01
#> [736] 1.703401e+02 7.383909e+01 6.952655e+01 1.078290e+02 9.168547e+01
#> [741] 8.837902e+01 8.656848e+01 2.622640e+02 2.270256e+02 8.618440e+01
#> [746] 7.222992e+01 1.164911e+02 1.896900e+02 7.456965e+01 6.745826e+01
#> [751] 1.093502e+02 1.088939e+02 1.117060e+02 8.049226e+01 6.625817e+01
#> [756] 7.737818e+01 1.125527e+02 9.958521e+01 7.162767e+01 9.881495e+01
#> [761] 7.649538e+01 7.140882e+01 1.520987e+02 1.039380e+02 1.444404e+02
#> [766] 8.331480e+01 7.662536e+01 1.188379e+02 6.757126e+02 7.781620e+02
#> [771] 1.470658e+02 7.741835e+01 6.674377e+01 8.391996e+01 8.322496e+01
#> [776] 1.220243e+02 8.510298e+01 8.098844e+01 7.224080e+02 9.989203e+01
#> [781] 9.133762e+01 2.526352e+02 9.637137e+01 8.554321e+01 6.468308e+01
#> [786] 6.449240e+01 6.496131e+02 5.547120e+02 1.429153e+02 6.575223e+01
#> [791] 6.319379e+01 8.170011e+01 7.175070e+01 6.787790e+01 1.745875e+02
#> [796] 1.075473e+02 7.236871e+01 6.252572e+01 7.064477e+01 6.958618e+01
#> [801] 8.736665e+01 1.456276e+02 6.949537e+01 8.348101e+01 6.825959e+01
#> [806] 9.298699e+01 7.301924e+01 6.702569e+01 7.862756e+01 7.917350e+01
#> [811] 6.630900e+01 9.123171e+01 7.708051e+01 7.491504e+01 1.675518e+03
#> [816] 7.480943e+01 1.004395e+02 7.018690e+01 3.710667e+02 9.095449e+01
#> [821] 7.777309e+01 7.441704e+01 6.355173e+01 1.118767e+02 1.064054e+02
#> [826] 7.576689e+01 1.110425e+02 2.760443e+02 1.130859e+02 1.016257e+02
#> [831] 2.793161e+02 1.109493e+02 7.987699e+02 7.892659e+01 9.917600e+01
#> [836] 7.322760e+01 7.825249e+01 8.669691e+01 6.386397e+01 6.805182e+01
#> [841] 7.257094e+01 9.626095e+01 6.846074e+01 6.437053e+01 1.473702e+02
#> [846] 1.312446e+02 9.433376e+01 9.247332e+01 6.974883e+01 7.924124e+01
#> [851] 1.744268e+02 8.453280e+01 1.489933e+02 7.523732e+01 2.464860e+02
#> [856] 8.283566e+01 2.278400e+02 2.147752e+02 1.107868e+02 8.743971e+01
#> [861] 9.670011e+01 5.624932e+01 6.676239e+01 6.570401e+01 1.050478e+02
#> [866] 7.249758e+01 1.117220e+02 8.448167e+01 7.601958e+01 1.032203e+02
#> [871] 1.362108e+02 9.059278e+01 8.779972e+01 1.638730e+02 8.146934e+01
#> [876] 9.779535e+01 8.237123e+01 7.157483e+01 7.036787e+01 6.897462e+01
#> [881] 1.433247e+02 1.249213e+03 8.452791e+01 1.568271e+02 7.467987e+01
#> [886] 1.082523e+02 6.522581e+01 7.096105e+01 9.954226e+01 1.022756e+02
#> [891] 1.067466e+02 8.805563e+01 6.499261e+01 7.821764e+01 1.074411e+02
#> [896] 6.796964e+01 6.578678e+01 3.845417e+02 7.208906e+01 9.653230e+01
#> [901] 7.060190e+01 9.006884e+01 8.087568e+01 8.558269e+01 8.678287e+01
#> [906] 9.854343e+01 6.248932e+01 1.097733e+02 6.514923e+01 8.428065e+01
#> [911] 8.195997e+01 8.011441e+01 9.758374e+01 1.602574e+02 1.546258e+02
#> [916] 1.060849e+03 8.049835e+01 6.542364e+01 7.006885e+01 2.737713e+02
#> [921] 1.050304e+02 6.368374e+01 7.555667e+01 6.677563e+01 7.905942e+01
#> [926] 6.448084e+01 1.386302e+02 6.482150e+01 2.465661e+02 7.580557e+01
#> [931] 9.612771e+01 7.528251e+01 3.157618e+03 2.199286e+02 1.178008e+02
#> [936] 8.729750e+01 2.238354e+02 9.141089e+01 1.938818e+02 7.917603e+01
#> [941] 7.249046e+01 9.637967e+01 7.217722e+01 8.811037e+01 7.182535e+01
#> [946] 1.973840e+02 4.610331e+02 7.623799e+01 1.099291e+05 4.091602e+02
#> [951] 1.501179e+02 9.213141e+01 8.748350e+01 7.823812e+01 7.925851e+01
#> [956] 4.281334e+02 7.239190e+01 4.491903e+02 1.012692e+02 6.669525e+01
#> [961] 2.425264e+02 1.067889e+03 7.459795e+01 8.273849e+01 1.016183e+03
#> [966] 6.727679e+01 7.519345e+01 8.416539e+01 1.653009e+02 1.342866e+02
#> [971] 8.789616e+01 1.189751e+04 7.312429e+02 7.734952e+01 9.074041e+01
#> [976] 2.201670e+02 1.114534e+02 8.935951e+01 8.025281e+01 7.345371e+01
#> [981] 7.995736e+01 2.363368e+04 1.563067e+02 1.534054e+05 1.435804e+03
#> [986] 1.088992e+02 8.360003e+01 1.627972e+02 1.873471e+02 4.999570e+02
#> [991] 7.817809e+01 1.421013e+02 8.152020e+01 6.898854e+01 2.534048e+02
#> [996] 8.242656e+01 1.055111e+02 6.989453e+01 7.608029e+01 1.095577e+02
doseFunction(
model,
alpha2 = samples@data$alpha2,
beta = samples@data$beta,
grade = 2L
)(x = 0.25)
#> [1] 102.64369 117.15507 173.19490 71.22183 61.67253
#> [6] 51.76125 53.96554 63.03466 80.22609 63.76243
#> [11] 73.58703 73.13015 58.92376 92.05023 80.47926
#> [16] 61.15206 87.14735 204.45844 98.25725 79.94967
#> [21] 61.82569 57.57232 55.82457 57.79289 60.75425
#> [26] 67.45191 60.27317 60.39204 69.78866 214.43111
#> [31] 59.91844 58.08053 52.99112 67.13045 198.64808
#> [36] 1066.02196 70.52019 67.82385 58.37147 59.67157
#> [41] 59.71479 72.43338 58.22979 68.82723 77.98779
#> [46] 75.88972 163.45826 63.82877 70.37588 52.44708
#> [51] 66.60672 106.67466 79.79577 68.57640 67.81432
#> [56] 58.45054 86.32107 62.00457 53.33612 68.35645
#> [61] 62.85234 75.32664 60.44759 78.11796 62.65782
#> [66] 55.12764 55.75952 59.73150 113.78676 55.59465
#> [71] 58.63939 61.73216 61.90516 84.52092 59.88831
#> [76] 60.39791 65.98417 89.75259 108.29628 57.94949
#> [81] 81.42866 53.42056 98.71184 73.47647 68.78734
#> [86] 58.65497 56.80091 66.47641 70.76410 56.44036
#> [91] 60.88633 55.29109 72.47116 57.33011 74.16017
#> [96] 79.45878 61.39491 86.05020 57.33565 95.31114
#> [101] 98.54148 87.20335 65.55387 59.35118 53.08111
#> [106] 62.60026 60.46424 55.15029 80.98620 64.00429
#> [111] 86.64543 583330.39523 92.23741 66.33146 77.92762
#> [116] 66.60909 56.35123 83.48601 55.78520 63.96886
#> [121] 171.06556 60.17308 53.98550 64.36821 59.88044
#> [126] 1278.49906 66.86609 369.92133 102.66589 110.83670
#> [131] 66.60119 55.27425 113.43403 515.02278 78.23182
#> [136] 61.34911 63.60315 62.34070 62.95867 68.25069
#> [141] 63.31721 59.17604 82.80674 56.94383 63.05997
#> [146] 60.66742 66.11568 111.21927 100.54634 101.82062
#> [151] 61.20333 58.35351 66.30506 58.29126 86.79551
#> [156] 69.80188 86.68309 97.56564 56.27666 76.24091
#> [161] 90.23292 58.70386 84.79313 58.02103 62.67021
#> [166] 69.44579 57.70274 62.52495 66.42088 69.34593
#> [171] 69.44018 86.20205 76.14161 63.45192 60.12592
#> [176] 85.33929 66.75812 67.46416 363.97886 52.05881
#> [181] 260.90534 106.19171 61.10459 121.52344 73.94386
#> [186] 53.12384 63.93035 57.53067 172.31499 61.92645
#> [191] 64.54918 61.77687 58.45639 67.98989 61.33361
#> [196] 128.32194 51.53858 74.77654 60.87828 186.07438
#> [201] 63.67870 89.51992 63.83710 65.63240 62.98478
#> [206] 76.80463 67.02720 105.39871 57.98192 68.41600
#> [211] 78.74429 66.87972 59.03371 72.07895 512.04610
#> [216] 527.34267 61.72034 57.66817 51.65687 54.59062
#> [221] 65.55367 68.57742 64.94444 101.27983 62.09949
#> [226] 75.00316 65.64658 75.45972 64.26092 58.71435
#> [231] 60.64752 68.61954 50.44323 65.15543 65.59982
#> [236] 168.32801 34.63065 62.01075 64.24251 58.20151
#> [241] 1199.93360 70.05334 76.90459 57.32991 58.71591
#> [246] 77.05295 50.03685 87.38992 110.22228 68.07181
#> [251] 85.03653 58.64464 57.29685 574.32262 125.50043
#> [256] 115.32960 125.65734 63.70923 114.91535 65.59379
#> [261] 62.11496 60.81248 117.25892 60.71397 58.92219
#> [266] 57.57641 59.50211 57.67368 57.30444 59.88072
#> [271] 70.10827 68.02256 69.24475 69.34729 62.47230
#> [276] 62.61721 59.96655 58.19032 60.51459 57.09245
#> [281] 59.18287 62.68461 70.22948 58.67018 73.06364
#> [286] 76.75359 56.72083 62.15662 106.19745 50.01209
#> [291] 65.94034 59.84028 115.55993 59.76624 59.80392
#> [296] 61.51105 59.70895 66.47765 53.22405 85.84855
#> [301] 61.89391 59.30138 327.23517 201.33780 81.64246
#> [306] 70.27999 71.37221 100.30354 51.85703 80.61376
#> [311] 62.33767 64.52661 61.96634 108.28248 78.41229
#> [316] 71.32386 82.83253 93.97822 55.89803 67.99137
#> [321] 75.68695 60.16611 84.50705 208.43749 200.66392
#> [326] 71.76437 64.25572 93.36202 55.65083 73.21321
#> [331] 58.78337 59.94618 58.92496 102.48526 60.41929
#> [336] 32842.35545 774.63335 103.98650 147.34554 84.16288
#> [341] 110.28235 51.93071 65.64018 62.94937 59.94887
#> [346] 70.27771 228.41589 57.75038 56.58256 63.19398
#> [351] 73.28156 211.84036 109.50067 59.33453 88.02060
#> [356] 55.10647 54.46606 59.52296 95.17231 66.67123
#> [361] 50.38366 75.21190 162.52121 93.12424 57.39925
#> [366] 57.15909 64.35865 76.95983 107.10155 52.86701
#> [371] 69.14476 105.24343 898.42097 1226.96525 94.14393
#> [376] 64.11909 77.73841 65.82313 60.20080 64.35403
#> [381] 63.38451 126.72795 67.16914 67.96743 169.72040
#> [386] 55.41071 90.40657 60.51417 105.96171 62.88169
#> [391] 69.33600 60.45933 75.73265 145.40859 385.78884
#> [396] 182.07637 68.88086 128.91416 125.09362 49.24027
#> [401] 54.04964 266.82495 90.50089 65.78705 60.49315
#> [406] 69.44370 62.49375 56.44730 60.96843 72.03055
#> [411] 62.14351 59.01584 220.42987 59.31791 67.91697
#> [416] 71.44140 56.20381 95.88954 80.69451 75.19822
#> [421] 56.03598 107.02590 57.88150 71.06076 58.28784
#> [426] 55.51634 60.88377 65.32816 61.36522 57.51216
#> [431] 54.71899 60.23930 64.64825 62.48245 68.10224
#> [436] 223.79674 60.49043 124.03360 59.06314 70.81779
#> [441] 70.20336 72.66042 89.23025 50.55186 63.56553
#> [446] 67.48590 64.86308 63.28959 66.75815 56.92006
#> [451] 784.77043 95.92652 77.18273 151.57407 151.59065
#> [456] 69.61358 65.01740 124.77822 67.95576 64.02048
#> [461] 65.14213 148.08868 80.70572 61.24474 98.28179
#> [466] 68.51771 56.47508 55.52430 70.79601 169.86955
#> [471] 86.82142 100.17326 59.89060 99.95289 63.74128
#> [476] 74.06842 60.06584 98.64653 92.45776 73.23983
#> [481] 66.72151 155.24532 70.39523 61.89763 68.84419
#> [486] 63.83782 66.47980 73.05467 76.66094 68.33567
#> [491] 54.01796 51.97677 82.48467 61.42278 57.88739
#> [496] 63.95781 63.54764 107.99509 114.48363 98.01250
#> [501] 127.36295 59.16539 76.12958 54.56878 72.93252
#> [506] 68.39105 67.23322 60.05553 63.28001 56.44712
#> [511] 84.31393 58.02190 63.88303 72.86267 177.26979
#> [516] 59.27508 56.93599 63.39299 72.29407 152.65947
#> [521] 110.82922 75.22607 56.64446 135.50401 73.54896
#> [526] 67.73738 779.70658 50.90358 45.03372 63.82678
#> [531] 66.27305 66.60543 72.25932 59.85920 61.59684
#> [536] 57.92511 63.35185 61.97390 92.00351 56.09547
#> [541] 66.09566 60.03014 61.80157 77.63037 560.15841
#> [546] 204.63266 24410.95904 54.97601 60.67499 61.20376
#> [551] 70.96419 79.56822 54.16161 72.82189 60.81448
#> [556] 58.69142 57.07407 59.87637 62.25586 76.88968
#> [561] 59.40693 54.25418 112.25635 327.57861 71.86868
#> [566] 78.36211 61.69115 108.01892 58.71344 65.96318
#> [571] 53.76047 118.88774 63.15354 71.10452 65.45825
#> [576] 59.82619 65.84156 71.28477 61.37708 62.24624
#> [581] 465.12005 65.60075 57.17799 52.51153 55.71016
#> [586] 60.47077 75.02304 56.07263 60.29437 64.62481
#> [591] 1162.86949 145.56944 67.97830 69.44512 63.95339
#> [596] 67.05639 55.00526 61.32535 61.85883 62.29261
#> [601] 70.82342 77.82396 112.12719 82.92763 97.50664
#> [606] 60.99972 60.75157 59.14228 73.27538 153.19076
#> [611] 63.65019 70.65397 56.50763 65.06594 58.44312
#> [616] 68.30097 61.87411 11256.04233 60.79631 66.74048
#> [621] 67.06521 73.55340 66.27761 59.56722 56.81754
#> [626] 69.22605 79.96900 65.15832 65.77083 59.21888
#> [631] 58.90574 431.87923 1160.43382 58.59016 61.37734
#> [636] 54.34717 61.36791 107.59760 71.76883 346.52328
#> [641] 54.42422 61.44808 65.17083 60.21392 63.47998
#> [646] 64.98063 66.31664 62.04032 92.54480 56.37799
#> [651] 61.01801 70.12622 73.81050 59.72891 49.12552
#> [656] 54.41695 59.00862 61.79229 62.30962 59.15488
#> [661] 65.53218 68.57144 54.95796 266.93619 61.79263
#> [666] 62.20732 63.46008 63.79966 59.84707 79.53530
#> [671] 72.20933 57.70513 54.50636 58.98336 59.75528
#> [676] 56.74125 59.06686 61.40467 62.27067 106.64455
#> [681] 61.61639 81.71543 64.93234 51.02083 98.84146
#> [686] 87.92836 56.59757 65.36116 128.58211 63.21919
#> [691] 63.21984 64.96192 149.00075 1317.90998 77.75248
#> [696] 64.32297 61.18718 55.15525 71.64499 152.94763
#> [701] 943.79725 106.26555 58.51759 62.26148 68.72177
#> [706] 72.91587 66.85099 69.73934 64.06051 62.15312
#> [711] 58.33734 61.52697 63.21788 69.95530 58.33874
#> [716] 59.10050 58.11341 56.79614 59.78432 71.16729
#> [721] 100.14676 58.85124 124.50639 99.01077 67.01056
#> [726] 69.61985 63.96385 59.74917 55.28978 59.86175
#> [731] 69.10928 54.36436 58.04379 55.71503 49.98528
#> [736] 115.28863 61.38446 57.47932 67.56428 67.11823
#> [741] 64.06680 59.53140 116.77149 110.39058 68.87431
#> [746] 57.68716 86.26582 92.59445 64.52862 57.62615
#> [751] 77.88952 59.96913 75.50667 66.62536 53.60814
#> [756] 65.05977 80.36840 69.42765 59.39089 63.38251
#> [761] 45.69862 62.58999 104.43992 78.62556 87.19690
#> [766] 52.28402 56.21062 72.18257 187.71403 136.66773
#> [771] 85.20922 64.36050 58.81871 64.47100 64.66527
#> [776] 86.02787 70.80516 67.32202 152.57380 61.20179
#> [781] 60.77348 155.76064 69.33782 68.29032 58.12220
#> [786] 59.59533 241.95211 214.78713 73.51469 57.03656
#> [791] 54.46130 62.83439 62.73068 59.60977 113.70757
#> [796] 72.99601 61.63247 56.90100 59.54954 58.58472
#> [801] 73.63854 91.42246 56.26699 61.38014 57.41407
#> [806] 68.93574 63.56133 58.62119 66.49777 65.95159
#> [811] 59.59643 73.76122 64.33829 63.21773 168.91741
#> [816] 56.66568 50.34500 51.74341 152.25132 68.17938
#> [821] 65.00788 63.58285 59.27756 80.81162 78.06604
#> [826] 63.96718 86.50273 101.64278 69.66315 70.63747
#> [831] 85.36645 79.29999 224.57722 64.27694 67.49218
#> [836] 60.26747 61.88871 68.51350 56.58678 61.73896
#> [841] 65.16448 77.99449 55.81570 56.53078 102.48045
#> [846] 58.01394 64.38943 74.12219 60.38965 64.81798
#> [851] 102.48884 63.46198 60.24902 60.29630 109.82198
#> [856] 62.49742 115.96716 74.78765 63.13707 62.67471
#> [861] 63.94304 51.20583 60.81513 60.24805 81.43065
#> [866] 61.62811 74.45304 68.46488 57.39064 78.57339
#> [871] 76.55371 70.30705 65.83684 78.36806 68.49869
#> [876] 66.89287 62.49376 58.21181 60.71281 56.04679
#> [881] 100.75908 157.56224 58.78068 97.20108 64.18964
#> [886] 73.90016 57.12024 60.04668 68.41033 60.50597
#> [891] 65.94004 62.00971 57.90053 63.46403 56.98599
#> [896] 61.24836 60.96472 153.33503 62.77852 74.06279
#> [901] 57.56885 64.25024 61.27613 63.13779 68.97493
#> [906] 71.67432 55.14410 71.50192 57.27057 69.72635
#> [911] 60.84157 67.55207 72.30074 95.58968 74.06077
#> [916] 149.36219 66.81242 56.81936 55.72112 145.73561
#> [921] 69.83806 51.76446 61.96462 52.78672 65.93049
#> [926] 57.45269 95.98608 56.02760 132.76921 59.64547
#> [931] 64.47970 55.80072 407.83696 74.39698 58.40326
#> [936] 65.37188 119.11564 56.01760 96.17182 63.39680
#> [941] 58.36726 56.19453 61.17808 63.35078 57.72107
#> [946] 109.67599 147.42842 55.97115 561.68325 143.37667
#> [951] 99.72788 71.78199 70.22973 62.96938 57.87614
#> [956] 144.77825 60.78083 193.59334 77.57144 58.36027
#> [961] 132.91564 244.87158 58.72184 54.98495 316.20737
#> [966] 55.75083 63.07043 64.11393 77.01863 69.70815
#> [971] 72.44364 457.54887 312.07648 58.06091 52.15224
#> [976] 103.31203 71.02110 71.77110 58.70640 60.02753
#> [981] 65.37721 1242.95115 92.69890 15142.68483 524.08084
#> [986] 79.01682 62.13882 81.70263 100.92669 164.98817
#> [991] 61.57859 96.25022 60.58889 61.29700 130.37677
#> [996] 67.51707 51.30460 54.76777 59.15275 71.76000