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] 8.088542e+01 6.526374e+01 6.423585e+01 7.016567e+01 6.614313e+01
#> [6] 8.761114e+01 1.145361e+02 7.429762e+01 6.396176e+01 6.825156e+01
#> [11] 9.726337e+01 7.079898e+01 7.054038e+01 8.281558e+01 8.861650e+01
#> [16] 7.637561e+01 6.911731e+01 1.586582e+02 7.379398e+01 8.811966e+01
#> [21] 1.512752e+02 2.056158e+04 3.214743e+10 9.136274e+02 2.731080e+02
#> [26] 5.195823e+05 2.262167e+03 8.839906e+01 1.061657e+02 1.625132e+02
#> [31] 1.222100e+02 2.559043e+02 3.006484e+02 7.410065e+03 3.559648e+06
#> [36] 6.920985e+01 1.339907e+03 9.745870e+02 6.508967e+01 6.660966e+01
#> [41] 1.166486e+02 1.378422e+02 1.491551e+03 2.718303e+02 9.424560e+01
#> [46] 7.189019e+01 1.687848e+02 3.419021e+02 9.845500e+01 1.041935e+02
#> [51] 1.182198e+02 7.973810e+01 1.165271e+02 1.312016e+02 2.256728e+03
#> [56] 7.159261e+01 8.728874e+01 6.017251e+01 7.526322e+01 6.947818e+01
#> [61] 6.806674e+01 7.089915e+01 6.201432e+01 7.216540e+01 6.525281e+01
#> [66] 7.797341e+01 7.123089e+01 1.034336e+02 1.060941e+02 8.595497e+01
#> [71] 8.357203e+01 8.157153e+01 7.069744e+01 7.246675e+01 1.452543e+02
#> [76] 8.653516e+01 8.391325e+01 2.082984e+02 4.984184e+02 1.569491e+02
#> [81] 1.060313e+02 7.577691e+01 7.687584e+01 6.984732e+01 7.636875e+01
#> [86] 1.064553e+02 9.836703e+01 8.759670e+01 7.192436e+01 1.146047e+02
#> [91] 8.237717e+01 7.863289e+03 8.061916e+01 6.271298e+01 9.093194e+01
#> [96] 1.486078e+02 7.057372e+01 7.437486e+01 7.681231e+01 1.106089e+02
#> [101] 9.414640e+01 1.006084e+10 2.168571e+02 1.017828e+02 2.856925e+02
#> [106] 6.426427e+01 1.109010e+02 1.327511e+02 1.088988e+02 6.911849e+01
#> [111] 6.406160e+01 7.039526e+01 2.658643e+02 8.467138e+01 1.126238e+02
#> [116] 7.223356e+04 4.032884e+02 2.035563e+02 7.353260e+01 2.258434e+02
#> [121] 3.206213e+02 7.878215e+01 8.033806e+01 1.064639e+02 1.474816e+02
#> [126] 1.158840e+02 8.835304e+01 7.622925e+01 1.236275e+02 1.245449e+02
#> [131] 9.320678e+01 6.385945e+01 1.136671e+02 1.629016e+02 6.393654e+01
#> [136] 7.844311e+01 8.997814e+01 8.005474e+01 1.117432e+02 8.589448e+01
#> [141] 9.016836e+01 6.456690e+01 7.329743e+01 6.641966e+01 9.072655e+01
#> [146] 1.104810e+02 7.902965e+01 9.178003e+01 6.618946e+01 6.847861e+01
#> [151] 1.233603e+02 6.649200e+01 1.677785e+02 1.001228e+02 8.327287e+01
#> [156] 1.004387e+02 1.492086e+02 2.428524e+07 1.040857e+02 1.207005e+02
#> [161] 8.142220e+01 1.061712e+02 7.279288e+01 1.155229e+02 1.963836e+02
#> [166] 7.085766e+01 2.756753e+04 5.766689e+02 1.329835e+02 8.217253e+01
#> [171] 1.016916e+02 1.081678e+02 1.936538e+02 3.230871e+02 6.482128e+01
#> [176] 2.395914e+02 7.496587e+01 1.977601e+02 8.317286e+01 1.207636e+02
#> [181] 1.124979e+02 8.978052e+01 7.533701e+01 9.530781e+01 2.049538e+02
#> [186] 1.640943e+03 1.071631e+02 9.827097e+01 1.574924e+03 6.352970e+01
#> [191] 8.369980e+01 1.189850e+02 1.033914e+02 1.013421e+02 8.280190e+01
#> [196] 1.077433e+02 6.773889e+01 1.210067e+02 7.258562e+01 8.317524e+01
#> [201] 7.092547e+01 1.568899e+02 1.076447e+02 6.016841e+01 7.056310e+01
#> [206] 8.940077e+01 7.897180e+01 9.582460e+01 1.013782e+02 1.006924e+02
#> [211] 7.133396e+01 8.008434e+01 8.715422e+01 7.081370e+01 7.628385e+01
#> [216] 3.808454e+02 7.161218e+01 6.818916e+01 6.961031e+01 6.411388e+01
#> [221] 1.030737e+02 1.138472e+02 6.471096e+01 6.545447e+01 6.667873e+01
#> [226] 6.878787e+01 6.808600e+01 1.451078e+02 7.984118e+01 7.192185e+01
#> [231] 1.622707e+03 6.714005e+01 6.196531e+01 7.566720e+01 7.337564e+01
#> [236] 6.866165e+01 8.302885e+01 3.194738e+02 1.339659e+02 2.279043e+02
#> [241] 1.000110e+02 8.505729e+01 1.170796e+02 8.575316e+01 8.480683e+01
#> [246] 9.153033e+01 1.073862e+02 1.012009e+02 1.038167e+02 1.467297e+02
#> [251] 7.523727e+01 7.989088e+01 7.446240e+01 6.279254e+01 1.311414e+02
#> [256] 7.319487e+01 1.665424e+02 7.163812e+01 1.162350e+02 9.811734e+01
#> [261] 9.965999e+01 1.863009e+02 2.184171e+02 2.038123e+02 8.317629e+01
#> [266] 6.579152e+01 8.486933e+01 8.013723e+01 9.345607e+01 7.148557e+01
#> [271] 9.088733e+02 1.290383e+02 7.521956e+01 1.141700e+02 6.713624e+01
#> [276] 1.074941e+02 1.166605e+02 6.196515e+01 9.001161e+01 8.464399e+01
#> [281] 6.568346e+01 8.601216e+01 9.600900e+01 8.276927e+01 7.515835e+01
#> [286] 1.021220e+02 1.000419e+02 9.205427e+01 7.448245e+01 1.069771e+02
#> [291] 7.200893e+01 6.898342e+01 6.285374e+01 5.250267e+02 6.847436e+02
#> [296] 1.482851e+02 8.200383e+01 7.667577e+01 7.005604e+01 7.276443e+03
#> [301] 1.048703e+02 6.686914e+01 1.333323e+02 8.158693e+01 1.175960e+02
#> [306] 6.319126e+01 7.129765e+01 8.521615e+01 2.021536e+02 8.279574e+01
#> [311] 8.072133e+01 8.918237e+01 6.840172e+01 8.903920e+01 1.244429e+02
#> [316] 3.081308e+02 8.471292e+01 9.076682e+01 6.965790e+01 6.517536e+01
#> [321] 1.358930e+02 7.243737e+01 6.256306e+01 6.320912e+01 9.120422e+01
#> [326] 7.481246e+01 7.030171e+01 7.741814e+01 2.252524e+02 1.202619e+02
#> [331] 7.774213e+01 7.113016e+01 1.030735e+02 8.943945e+01 6.848151e+01
#> [336] 2.194785e+02 1.395307e+02 8.258515e+01 1.043869e+02 1.720057e+02
#> [341] 8.685380e+01 9.655929e+01 6.861166e+01 1.411704e+02 8.339419e+01
#> [346] 9.538082e+01 7.519590e+01 8.234728e+01 7.624971e+01 1.278358e+02
#> [351] 7.097985e+01 4.133046e+03 7.517932e+01 1.169836e+02 1.003165e+02
#> [356] 1.257126e+03 2.985218e+02 9.244488e+01 6.842171e+01 7.123411e+01
#> [361] 9.020031e+01 1.026662e+02 7.375358e+01 1.424299e+02 1.566847e+02
#> [366] 1.964082e+02 1.290521e+02 9.752801e+02 1.703121e+02 1.097898e+02
#> [371] 6.947439e+01 7.754376e+01 1.545993e+02 7.332026e+01 8.226448e+01
#> [376] 6.460463e+01 9.180804e+01 3.012396e+02 7.202216e+01 2.849393e+02
#> [381] 4.297046e+03 9.677714e+01 7.149062e+01 1.001670e+02 1.290418e+02
#> [386] 7.038702e+01 9.783089e+01 8.659499e+01 9.777447e+01 7.099699e+01
#> [391] 9.598649e+01 7.860775e+01 9.419715e+01 1.414536e+02 6.470866e+01
#> [396] 1.744057e+02 6.539909e+01 6.590788e+01 9.953592e+01 1.049379e+02
#> [401] 6.725448e+01 8.690188e+01 7.437727e+01 9.167790e+01 9.137210e+01
#> [406] 9.392796e+01 7.694345e+01 1.211297e+03 5.000005e+02 7.100205e+02
#> [411] 7.688942e+01 6.575839e+01 6.928852e+01 6.779403e+01 1.015726e+02
#> [416] 9.538999e+01 8.522947e+01 8.776541e+01 7.588500e+01 6.919900e+01
#> [421] 8.997735e+01 7.778939e+01 1.842165e+02 7.466355e+01 8.812537e+01
#> [426] 6.321990e+01 1.342933e+02 7.538897e+01 6.523901e+01 1.098847e+02
#> [431] 8.375490e+01 1.239173e+02 2.481778e+02 2.057136e+02 8.122328e+01
#> [436] 6.532841e+01 1.043425e+02 9.073271e+01 7.252578e+01 1.759468e+02
#> [441] 3.248242e+02 1.355709e+02 8.018290e+01 8.360545e+01 1.028031e+02
#> [446] 1.687665e+02 8.249054e+01 2.026285e+02 2.942998e+02 7.039236e+01
#> [451] 6.803133e+01 9.069613e+01 3.576144e+02 1.183284e+02 3.647710e+02
#> [456] 6.625066e+01 7.045611e+01 1.556851e+02 7.099060e+01 9.650658e+01
#> [461] 6.959717e+01 1.472496e+02 7.403648e+01 8.421320e+01 2.274865e+02
#> [466] 7.756567e+01 7.468978e+01 1.492484e+02 1.201862e+02 6.801384e+01
#> [471] 2.670283e+02 2.610197e+02 1.124759e+02 6.395525e+01 6.641779e+01
#> [476] 6.690665e+01 9.446393e+01 9.567069e+01 7.284431e+01 1.247192e+02
#> [481] 7.132906e+01 9.715046e+01 8.166699e+01 1.118384e+02 1.335539e+02
#> [486] 2.246664e+02 7.632151e+01 8.742365e+01 7.411577e+01 2.443616e+02
#> [491] 1.783839e+02 5.133813e+02 1.254514e+02 8.642498e+01 1.483744e+02
#> [496] 1.053331e+02 7.823920e+01 1.120050e+02 7.772191e+01 2.814318e+02
#> [501] 1.368995e+02 7.959693e+02 6.006473e+01 6.974552e+01 7.664509e+01
#> [506] 6.725283e+01 6.175519e+01 7.819322e+01 8.100981e+01 8.669894e+01
#> [511] 8.529298e+01 4.862115e+02 1.381202e+02 1.025416e+02 1.411296e+02
#> [516] 6.635686e+01 6.069072e+01 7.022877e+01 7.603134e+01 8.956779e+01
#> [521] 2.812019e+02 7.403702e+07 1.583679e+02 5.998672e+01 6.093937e+01
#> [526] 6.285422e+01 2.227718e+02 8.204639e+01 9.794419e+01 1.629792e+03
#> [531] 2.823503e+07 3.668208e+02 1.254327e+02 1.546661e+02 6.475651e+02
#> [536] 1.664657e+02 6.986434e+01 1.059454e+02 4.440766e+02 1.664972e+02
#> [541] 1.048651e+02 3.884826e+02 7.832691e+01 1.400968e+02 1.043294e+02
#> [546] 1.031505e+02 7.423222e+01 6.648473e+01 6.907511e+01 6.934951e+01
#> [551] 1.399645e+02 6.552120e+01 7.859328e+01 7.305187e+01 8.243819e+01
#> [556] 1.101014e+02 7.608384e+01 6.459829e+01 7.825657e+01 8.092199e+01
#> [561] 1.094467e+02 2.169310e+02 6.543055e+01 8.137445e+01 1.140804e+02
#> [566] 1.102613e+02 2.185377e+03 2.378135e+02 7.407370e+01 1.601313e+02
#> [571] 6.679986e+01 7.952797e+01 7.308536e+01 6.612690e+01 7.237889e+01
#> [576] 5.755777e+02 8.805095e+01 1.497962e+02 5.055716e+02 1.563249e+02
#> [581] 8.568513e+01 6.885901e+01 9.847047e+01 7.070179e+01 6.954561e+01
#> [586] 8.484543e+01 1.444394e+02 6.552037e+01 6.838071e+01 7.981895e+01
#> [591] 7.471682e+01 6.547934e+01 6.820763e+01 6.895235e+01 8.610207e+01
#> [596] 1.367009e+02 1.244027e+02 1.080531e+02 1.066460e+02 1.577227e+02
#> [601] 3.366446e+02 1.417715e+02 9.481499e+01 1.181444e+02 7.486610e+01
#> [606] 1.578282e+02 7.005654e+01 2.585732e+02 1.084200e+03 1.273020e+03
#> [611] 7.544211e+01 1.156579e+02 9.719510e+01 6.997276e+01 6.579416e+01
#> [616] 9.059253e+01 6.757229e+01 1.287509e+02 1.046236e+02 7.671404e+01
#> [621] 7.229989e+01 8.747087e+01 7.109991e+01 1.898008e+02 1.572576e+04
#> [626] 6.768648e+01 6.974518e+01 7.908378e+01 7.536001e+01 6.520141e+01
#> [631] 5.376811e+03 7.895137e+01 6.225099e+01 8.443776e+01 9.466849e+01
#> [636] 7.849867e+01 2.843609e+02 9.080355e+01 6.504999e+01 6.551225e+01
#> [641] 7.372758e+01 6.612945e+01 7.011906e+01 1.139473e+02 9.385382e+01
#> [646] 1.607191e+02 7.392768e+01 1.177255e+02 6.571213e+01 6.386204e+01
#> [651] 1.121163e+02 6.703203e+01 8.731690e+01 2.272942e+03 2.089011e+02
#> [656] 7.436321e+01 7.736530e+01 1.712112e+02 8.774588e+01 7.000507e+01
#> [661] 6.975126e+01 7.842163e+01 6.999336e+01 9.550701e+01 1.572171e+02
#> [666] 7.240199e+01 1.394834e+03 9.441970e+02 3.420200e+04 1.261084e+02
#> [671] 7.856632e+01 2.049698e+02 8.109916e+01 1.304591e+02 1.177794e+02
#> [676] 7.183116e+01 6.641804e+01 2.132352e+02 1.120780e+02 3.608176e+04
#> [681] 1.900220e+02 3.314045e+02 2.181125e+02 1.211884e+02 5.167765e+02
#> [686] 3.836358e+02 7.032074e+01 8.111829e+01 7.606011e+01 8.343025e+01
#> [691] 1.524740e+02 1.264647e+02 1.083021e+02 6.303124e+01 1.125398e+02
#> [696] 1.039291e+02 7.521294e+01 6.770861e+01 1.325283e+02 8.878317e+01
#> [701] 1.844457e+02 8.033475e+01 1.382150e+02 1.055921e+02 8.930472e+01
#> [706] 2.162167e+02 8.902811e+01 1.057327e+02 1.049892e+02 8.416657e+01
#> [711] 7.959387e+01 1.368745e+02 6.290304e+01 6.497264e+01 6.371146e+01
#> [716] 8.090950e+01 7.844284e+01 7.029241e+01 8.596877e+02 7.706287e+02
#> [721] 1.796842e+02 1.572853e+02 8.619983e+01 9.654256e+01 6.782840e+01
#> [726] 7.395249e+01 8.418346e+01 7.421779e+01 6.700635e+01 8.821537e+02
#> [731] 6.474911e+01 6.955940e+01 1.359929e+02 1.155041e+02 2.534116e+02
#> [736] 1.111942e+02 6.614335e+01 1.157930e+02 9.488274e+01 1.187667e+02
#> [741] 1.736418e+04 1.245895e+02 1.978702e+02 1.043824e+02 1.170108e+02
#> [746] 6.329889e+01 9.531069e+01 1.659528e+02 8.402844e+01 5.397445e+02
#> [751] 1.719917e+02 5.886950e+01 8.815996e+01 7.203714e+01 1.205991e+02
#> [756] 8.033831e+01 1.001797e+02 6.921391e+01 8.925519e+01 8.357989e+01
#> [761] 1.306076e+02 6.783411e+01 1.550608e+02 1.970712e+02 2.515793e+02
#> [766] 1.160161e+02 1.021616e+02 1.028732e+02 1.058645e+02 7.446612e+01
#> [771] 1.273889e+02 8.725158e+01 1.020112e+02 9.002625e+01 6.432880e+01
#> [776] 7.266572e+01 6.267997e+01 5.255554e+02 1.434106e+02 8.110795e+01
#> [781] 7.074942e+01 9.143808e+01 9.755454e+01 6.918048e+01 1.307014e+02
#> [786] 7.472079e+01 2.272234e+02 6.946162e+01 8.477158e+01 1.142240e+02
#> [791] 7.973984e+01 7.543839e+01 9.582552e+01 6.837608e+01 7.361497e+01
#> [796] 1.042837e+02 6.111504e+01 7.579299e+01 7.041419e+01 7.952779e+01
#> [801] 7.011206e+01 7.037126e+01 6.785802e+01 8.520693e+01 8.822714e+01
#> [806] 8.411601e+01 1.063448e+02 2.509015e+02 1.103272e+02 9.455631e+01
#> [811] 2.282753e+02 8.554791e+01 1.523099e+02 6.456991e+01 7.749745e+01
#> [816] 6.878925e+01 6.211147e+01 9.007987e+01 1.702704e+02 5.252699e+02
#> [821] 1.809934e+02 1.353517e+02 9.951104e+01 8.610426e+01 8.344266e+01
#> [826] 8.748858e+01 1.280949e+02 6.386345e+01 6.954526e+01 8.230779e+01
#> [831] 7.251234e+01 1.679340e+02 6.690190e+01 8.717731e+01 7.654338e+01
#> [836] 6.871352e+01 3.361437e+02 6.497166e+01 7.926209e+01 7.101330e+01
#> [841] 7.372509e+01 2.294193e+02 2.630182e+02 1.215807e+02 2.271131e+02
#> [846] 2.807315e+05 2.877368e+02 2.120457e+02 6.286190e+01 7.094109e+01
#> [851] 7.483188e+01 9.553520e+01 6.838069e+01 7.491676e+01 7.247958e+01
#> [856] 7.531143e+01 8.548063e+01 6.939428e+01 1.214798e+02 9.173216e+01
#> [861] 6.496346e+01 1.283491e+02 7.839318e+01 7.901542e+01 6.311125e+01
#> [866] 7.715881e+01 8.900737e+02 1.616745e+02 7.141122e+01 6.669253e+01
#> [871] 7.504125e+01 7.096415e+01 6.835251e+01 6.120517e+01 6.049157e+01
#> [876] 1.746721e+02 6.939971e+01 5.337589e+02 2.898942e+02 8.974972e+01
#> [881] 1.859449e+04 1.248245e+02 9.061592e+01 6.857166e+01 8.203598e+01
#> [886] 1.098412e+02 1.256170e+02 1.024984e+02 8.446904e+01 7.827825e+01
#> [891] 9.181739e+01 8.060172e+01 2.004689e+02 7.119009e+03 9.928365e+01
#> [896] 7.173359e+01 7.053766e+01 7.805213e+01 7.281677e+01 7.996823e+01
#> [901] 6.566792e+01 9.720937e+01 1.320696e+02 1.229359e+02 1.078188e+02
#> [906] 7.464026e+01 1.001904e+02 6.425219e+02 6.930247e+01 9.630992e+01
#> [911] 8.960073e+01 7.476841e+01 2.158500e+02 8.011825e+01 7.620360e+01
#> [916] 6.580836e+01 7.804021e+01 6.749600e+01 1.016986e+02 7.290773e+01
#> [921] 6.847846e+01 2.332101e+02 1.326364e+02 6.895802e+01 1.079530e+02
#> [926] 8.009026e+01 6.937666e+01 8.349575e+01 6.612901e+01 1.031605e+02
#> [931] 9.213425e+01 1.004542e+02 2.837423e+03 8.651951e+01 6.092642e+01
#> [936] 7.022698e+01 8.028926e+01 6.254550e+01 1.394862e+02 9.120961e+01
#> [941] 6.751041e+01 8.939930e+01 1.422072e+02 2.116331e+02 7.225448e+01
#> [946] 7.747772e+01 8.737380e+01 7.037257e+01 8.877674e+01 1.308382e+02
#> [951] 9.127290e+02 7.893247e+02 6.569689e+01 1.200407e+02 6.804761e+01
#> [956] 1.299389e+02 1.628372e+02 9.299679e+01 9.227680e+01 9.607588e+01
#> [961] 6.745105e+01 7.319091e+01 1.024963e+02 1.667325e+02 1.083688e+02
#> [966] 6.482559e+02 8.946189e+01 7.028713e+01 1.312928e+02 7.412931e+01
#> [971] 1.148463e+02 5.591560e+02 2.757608e+02 1.339870e+02 1.979073e+02
#> [976] 7.850708e+01 1.541854e+02 6.257279e+01 7.845535e+01 6.702158e+01
#> [981] 6.825997e+01 7.589213e+01 1.401152e+02 9.343483e+01 8.930004e+01
#> [986] 6.994579e+01 6.678670e+01 7.290033e+01 8.264838e+01 7.766656e+01
#> [991] 1.516577e+02 6.034953e+01 8.809673e+01 1.070284e+02 1.583198e+02
#> [996] 7.225976e+01 4.415269e+02 9.877248e+01 1.531919e+02 9.277334e+01
doseFunction(
model,
alpha2 = samples@data$alpha2,
beta = samples@data$beta,
grade = 2L
)(x = 0.25)
#> [1] 5.925999e+01 5.531300e+01 5.757795e+01 6.241956e+01 6.006592e+01
#> [6] 6.952153e+01 6.620278e+01 5.413703e+01 5.530973e+01 5.876131e+01
#> [11] 6.665219e+01 5.197982e+01 5.428308e+01 6.122297e+01 6.139928e+01
#> [16] 6.156058e+01 5.438875e+01 6.331477e+01 6.219429e+01 7.093272e+01
#> [21] 1.030714e+02 2.380767e+03 5.846617e+07 1.977158e+02 1.251535e+02
#> [26] 1.420717e+02 1.759989e+02 6.358976e+01 6.630972e+01 8.977205e+01
#> [31] 7.798844e+01 9.918325e+01 1.094710e+02 2.424146e+02 1.229974e+04
#> [36] 6.247246e+01 2.222286e+02 1.927348e+02 6.033875e+01 6.093451e+01
#> [41] 9.044765e+01 1.052997e+02 5.920084e+02 1.848245e+02 7.133188e+01
#> [46] 6.338017e+01 9.047176e+01 1.620996e+02 7.112794e+01 7.193943e+01
#> [51] 6.623888e+01 5.522173e+01 6.700597e+01 8.875116e+01 5.230368e+02
#> [56] 5.841493e+01 5.278651e+01 5.502572e+01 6.490782e+01 6.059248e+01
#> [61] 5.908302e+01 5.995270e+01 5.637005e+01 6.180411e+01 5.927403e+01
#> [66] 6.193376e+01 6.293153e+01 6.783429e+01 7.257693e+01 6.219758e+01
#> [71] 6.667685e+01 7.050440e+01 6.347017e+01 6.248917e+01 7.678850e+01
#> [76] 6.675224e+01 6.128292e+01 8.324474e+01 1.964282e+02 9.070072e+01
#> [81] 5.755098e+01 6.064838e+01 6.221818e+01 5.872740e+01 6.135591e+01
#> [86] 7.083921e+01 6.813520e+01 6.422972e+01 6.118911e+01 8.790513e+01
#> [91] 5.997962e+01 1.721951e+03 6.496347e+01 5.313675e+01 6.408342e+01
#> [96] 1.075847e+02 5.644565e+01 6.537906e+01 6.615002e+01 7.975852e+01
#> [101] 6.031103e+01 1.492777e+07 7.206497e+01 6.872091e+01 1.377792e+02
#> [106] 5.718768e+01 7.181093e+01 8.307849e+01 6.048505e+01 6.063607e+01
#> [111] 5.826860e+01 6.280562e+01 1.401992e+02 5.312371e+01 8.221577e+01
#> [116] 5.605997e+02 1.349387e+02 1.131972e+02 6.344227e+01 8.188382e+01
#> [121] 1.323929e+02 6.043183e+01 6.150978e+01 6.477632e+01 8.897522e+01
#> [126] 6.380011e+01 6.620688e+01 5.953769e+01 8.454486e+01 7.601462e+01
#> [131] 6.688170e+01 5.735713e+01 8.296696e+01 8.428273e+01 5.885041e+01
#> [136] 6.795608e+01 6.889966e+01 6.043837e+01 8.203806e+01 6.407763e+01
#> [141] 7.005452e+01 5.653859e+01 5.673094e+01 5.766764e+01 7.232964e+01
#> [146] 6.441725e+01 6.459880e+01 6.473617e+01 5.510201e+01 5.648377e+01
#> [151] 8.081666e+01 5.824293e+01 1.005133e+02 7.272083e+01 5.804155e+01
#> [156] 7.011859e+01 8.329819e+01 9.808536e+04 6.608738e+01 8.541884e+01
#> [161] 6.771562e+01 6.335705e+01 5.344768e+01 6.954527e+01 7.213937e+01
#> [166] 5.736855e+01 2.180706e+03 6.637841e+01 6.160940e+01 6.995454e+01
#> [171] 7.787638e+01 6.916036e+01 9.886831e+01 1.459570e+02 5.437401e+01
#> [176] 9.448113e+01 6.264875e+01 8.257761e+01 6.406032e+01 7.842364e+01
#> [181] 7.693698e+01 6.491393e+01 5.591294e+01 7.167586e+01 1.133229e+02
#> [186] 5.028257e+02 7.443433e+01 7.402399e+01 2.770608e+02 5.213684e+01
#> [191] 6.546315e+01 7.953096e+01 6.967340e+01 7.001715e+01 6.681571e+01
#> [196] 7.212755e+01 5.715309e+01 8.202926e+01 6.017629e+01 6.127350e+01
#> [201] 5.351081e+01 8.432281e+01 7.411437e+01 5.358135e+01 6.138732e+01
#> [206] 7.494281e+01 6.046580e+01 7.737293e+01 7.477434e+01 7.947313e+01
#> [211] 6.102070e+01 6.298819e+01 7.020264e+01 6.335712e+01 6.978381e+01
#> [216] 1.402981e+02 5.400788e+01 5.798256e+01 6.028319e+01 5.760676e+01
#> [221] 8.049233e+01 7.737539e+01 5.888053e+01 5.832421e+01 5.848106e+01
#> [226] 6.100337e+01 5.759535e+01 7.377037e+01 6.291907e+01 6.186806e+01
#> [231] 2.566241e+02 5.323519e+01 5.451338e+01 6.545852e+01 6.073211e+01
#> [236] 5.623731e+01 6.482497e+01 1.560979e+02 6.665067e+01 1.021882e+02
#> [241] 6.765248e+01 6.781585e+01 8.918029e+01 6.050308e+01 6.800536e+01
#> [246] 7.042467e+01 7.103741e+01 7.274749e+01 6.103529e+01 1.023255e+02
#> [251] 5.970299e+01 6.416927e+01 6.304258e+01 5.806523e+01 9.113304e+01
#> [256] 5.711377e+01 9.817738e+01 6.290796e+01 8.749409e+01 6.130928e+01
#> [261] 5.004077e+01 1.162352e+02 1.257128e+02 1.107053e+02 5.893841e+01
#> [266] 5.384485e+01 6.909759e+01 6.613584e+01 5.406579e+01 5.624904e+01
#> [271] 1.330870e+02 7.157112e+01 5.637626e+01 7.887567e+01 5.864633e+01
#> [276] 7.528479e+01 8.255261e+01 5.785977e+01 7.582025e+01 7.071142e+01
#> [281] 5.863349e+01 7.282418e+01 7.215456e+01 5.858658e+01 5.769198e+01
#> [286] 7.253055e+01 5.194912e+01 6.541036e+01 6.193444e+01 7.389431e+01
#> [291] 5.846365e+01 5.908411e+01 5.666328e+01 1.473721e+02 8.507580e+01
#> [296] 8.934755e+01 6.823548e+01 6.115225e+01 5.914501e+01 1.428730e+03
#> [301] 7.708992e+01 6.007561e+01 9.032299e+01 6.065484e+01 6.476470e+01
#> [306] 5.628246e+01 6.346866e+01 6.108431e+01 1.120308e+02 6.743984e+01
#> [311] 6.303131e+01 6.074384e+01 5.345669e+01 6.870173e+01 9.292199e+01
#> [316] 1.798641e+02 7.066952e+01 7.228572e+01 6.243629e+01 5.982018e+01
#> [321] 1.067181e+02 6.573551e+01 5.867966e+01 5.831897e+01 7.631543e+01
#> [326] 6.089884e+01 5.930713e+01 6.617060e+01 1.298918e+02 8.186458e+01
#> [331] 5.776012e+01 6.110486e+01 7.848297e+01 7.277774e+01 6.080560e+01
#> [336] 9.705790e+01 7.012282e+01 6.216462e+01 7.080178e+01 9.813772e+01
#> [341] 6.260192e+01 7.413979e+01 5.996880e+01 8.690276e+01 6.845416e+01
#> [346] 6.601430e+01 5.557224e+01 5.971935e+01 6.083196e+01 8.657432e+01
#> [351] 5.654607e+01 6.465489e+02 5.944430e+01 7.872979e+01 7.500291e+01
#> [356] 2.076589e+02 1.416172e+02 6.855590e+01 5.671892e+01 6.042894e+01
#> [361] 6.910907e+01 7.236936e+01 6.257249e+01 1.008598e+02 7.853487e+01
#> [366] 1.131125e+02 7.925006e+01 1.608436e+02 7.829993e+01 7.554794e+01
#> [371] 5.722196e+01 6.617468e+01 1.030007e+02 6.398481e+01 6.786560e+01
#> [376] 5.995478e+01 7.699708e+01 1.552301e+02 5.183626e+01 1.000762e+02
#> [381] 2.141539e+02 6.920143e+01 5.835959e+01 5.701992e+01 6.374865e+01
#> [386] 6.221050e+01 8.007697e+01 6.842070e+01 5.940632e+01 5.645321e+01
#> [391] 6.570648e+01 6.109174e+01 7.163749e+01 8.381254e+01 5.788787e+01
#> [396] 1.112426e+02 5.809426e+01 6.058703e+01 7.656549e+01 6.955369e+01
#> [401] 5.672806e+01 6.301195e+01 6.256980e+01 7.173117e+01 6.178571e+01
#> [406] 6.626945e+01 5.755895e+01 3.156902e+02 2.007367e+02 1.288028e+02
#> [411] 6.466365e+01 5.562166e+01 6.045585e+01 5.983036e+01 7.727624e+01
#> [416] 7.653876e+01 6.858683e+01 6.626794e+01 5.314592e+01 5.878044e+01
#> [421] 6.397041e+01 6.142082e+01 8.537902e+01 6.005171e+01 7.229292e+01
#> [426] 5.711050e+01 9.978001e+01 5.960089e+01 5.539988e+01 6.841609e+01
#> [431] 7.054980e+01 7.196136e+01 6.573115e+01 1.045859e+02 5.862926e+01
#> [436] 5.761452e+01 7.499506e+01 6.445616e+01 5.999241e+01 9.589931e+01
#> [441] 1.747023e+02 9.505710e+01 6.597796e+01 6.119122e+01 7.624678e+01
#> [446] 9.376032e+01 5.289557e+01 9.275655e+01 6.555993e+01 5.913013e+01
#> [451] 5.702510e+01 6.371019e+01 1.288302e+02 5.911758e+01 1.157185e+02
#> [456] 5.793934e+01 6.129267e+01 7.982053e+01 5.598323e+01 6.790006e+01
#> [461] 5.998137e+01 8.844763e+01 6.031463e+01 6.631541e+01 1.090007e+02
#> [466] 5.927581e+01 5.248483e+01 9.627651e+01 7.894853e+01 5.846997e+01
#> [471] 1.390022e+02 1.293503e+02 8.064559e+01 5.689126e+01 5.814462e+01
#> [476] 5.924703e+01 7.348757e+01 6.824474e+01 4.941825e+01 5.065970e+01
#> [481] 5.408972e+01 6.704526e+01 5.226991e+01 7.169470e+01 8.646030e+01
#> [486] 8.534495e+01 6.361962e+01 6.744758e+01 6.123843e+01 9.587553e+01
#> [491] 8.716355e+01 1.154201e+02 7.115728e+01 5.483179e+01 7.594228e+01
#> [496] 7.166490e+01 6.000705e+01 7.980841e+01 6.473579e+01 1.362544e+02
#> [501] 7.948554e+01 2.565728e+02 5.228536e+01 6.354964e+01 6.477680e+01
#> [506] 6.102254e+01 5.618169e+01 6.642519e+01 6.690714e+01 7.079687e+01
#> [511] 7.050391e+01 1.534186e+02 8.527529e+01 7.537359e+01 8.898290e+01
#> [516] 5.371020e+01 5.646430e+01 6.332646e+01 6.312005e+01 7.411350e+01
#> [521] 9.786201e+01 4.406054e+04 9.169280e+01 5.328488e+01 5.522532e+01
#> [526] 5.715348e+01 9.673338e+01 6.241564e+01 6.629980e+01 3.402373e+02
#> [531] 2.029432e+05 1.099463e+02 5.881373e+01 8.282073e+01 2.382367e+02
#> [536] 6.381668e+01 5.932809e+01 7.213184e+01 1.242430e+02 6.281185e+01
#> [541] 7.667187e+01 1.214731e+02 6.075627e+01 7.628732e+01 7.196453e+01
#> [546] 6.141203e+01 5.862212e+01 5.552354e+01 5.603789e+01 5.714377e+01
#> [551] 6.873717e+01 5.590693e+01 6.301428e+01 6.135309e+01 6.502996e+01
#> [556] 7.292677e+01 6.518037e+01 5.908589e+01 6.744190e+01 6.653392e+01
#> [561] 7.967147e+01 1.127113e+02 5.382597e+01 6.497792e+01 7.231788e+01
#> [566] 7.207154e+01 2.814263e+02 7.835333e+01 5.660342e+01 8.490592e+01
#> [571] 6.042019e+01 5.796723e+01 5.977916e+01 5.552054e+01 6.064995e+01
#> [576] 1.373238e+02 7.104018e+01 7.553972e+01 1.261899e+02 7.196330e+01
#> [581] 7.311982e+01 5.991836e+01 6.117449e+01 5.904820e+01 6.012511e+01
#> [586] 6.393423e+01 7.777784e+01 5.948239e+01 6.089999e+01 6.624836e+01
#> [591] 6.304928e+01 5.325522e+01 5.907758e+01 6.090495e+01 6.501884e+01
#> [596] 8.123435e+01 7.884402e+01 7.229090e+01 7.804611e+01 1.085383e+02
#> [601] 1.618358e+02 9.286781e+01 6.021284e+01 8.020074e+01 5.241559e+01
#> [606] 8.420908e+01 5.856341e+01 1.100217e+02 1.542661e+02 3.214220e+02
#> [611] 5.261496e+01 6.868463e+01 7.186517e+01 5.969738e+01 5.847243e+01
#> [616] 6.152501e+01 5.983180e+01 8.740719e+01 7.540338e+01 6.223978e+01
#> [621] 5.940524e+01 6.624094e+01 6.218525e+01 1.161805e+02 7.683708e+02
#> [626] 5.560138e+01 5.950961e+01 6.105095e+01 6.179932e+01 5.801704e+01
#> [631] 2.761702e+02 5.416631e+01 5.379108e+01 6.576074e+01 6.261029e+01
#> [636] 6.391515e+01 1.205187e+02 6.373108e+01 5.534959e+01 5.987562e+01
#> [641] 6.214496e+01 5.886530e+01 6.163170e+01 8.177808e+01 6.686396e+01
#> [646] 9.825547e+01 6.084521e+01 6.417349e+01 5.747037e+01 5.472231e+01
#> [651] 7.682555e+01 6.063066e+01 6.998713e+01 2.899332e+02 9.034729e+01
#> [656] 6.269816e+01 6.005459e+01 9.880928e+01 7.024082e+01 6.229387e+01
#> [661] 5.917588e+01 6.448105e+01 5.994281e+01 6.547345e+01 9.207563e+01
#> [666] 5.850614e+01 1.541997e+02 2.249990e+02 3.079079e+03 8.717565e+01
#> [671] 6.329659e+01 9.805770e+01 6.556458e+01 7.698731e+01 7.490594e+01
#> [676] 5.753314e+01 5.439222e+01 7.661286e+01 6.521755e+01 2.814785e+03
#> [681] 5.753753e+01 1.152582e+02 8.433249e+01 7.850917e+01 2.080579e+02
#> [686] 1.727568e+02 5.266015e+01 6.710204e+01 6.405604e+01 6.427720e+01
#> [691] 8.569056e+01 8.114634e+01 7.511584e+01 5.420978e+01 6.497548e+01
#> [696] 6.435341e+01 6.466515e+01 5.953541e+01 7.175888e+01 6.650158e+01
#> [701] 9.710741e+01 6.414676e+01 7.518516e+01 6.805202e+01 6.670789e+01
#> [706] 8.775067e+01 5.562802e+01 7.026742e+01 6.222763e+01 6.671833e+01
#> [711] 6.423482e+01 8.866427e+01 5.618518e+01 5.724844e+01 5.700937e+01
#> [716] 6.528024e+01 6.197675e+01 5.731708e+01 1.948225e+02 2.499434e+02
#> [721] 1.014445e+02 9.925354e+01 5.428000e+01 6.988520e+01 5.293694e+01
#> [726] 6.134146e+01 6.203640e+01 6.109492e+01 5.617471e+01 1.870969e+02
#> [731] 5.185216e+01 5.905997e+01 7.855297e+01 6.729018e+01 1.075215e+02
#> [736] 7.075266e+01 5.662456e+01 6.646977e+01 6.579456e+01 7.782568e+01
#> [741] 1.546116e+03 9.161376e+01 1.202621e+02 7.557432e+01 7.585201e+01
#> [746] 5.609479e+01 7.992952e+01 6.153705e+01 5.479292e+01 1.150992e+02
#> [751] 8.876304e+01 5.344199e+01 7.081530e+01 4.954859e+01 6.315657e+01
#> [756] 6.821839e+01 7.075281e+01 5.926949e+01 6.495832e+01 6.960130e+01
#> [761] 8.282852e+01 5.281443e+01 6.578474e+01 7.358812e+01 6.964384e+01
#> [766] 5.192793e+01 7.525206e+01 7.519384e+01 8.074878e+01 5.766556e+01
#> [771] 8.602789e+01 6.955279e+01 6.757387e+01 6.240710e+01 5.648755e+01
#> [776] 5.711954e+01 5.535268e+01 1.685317e+02 9.510778e+01 6.334537e+01
#> [781] 5.818247e+01 6.343283e+01 6.753514e+01 5.688799e+01 7.991068e+01
#> [786] 5.787062e+01 1.186647e+02 5.821252e+01 6.529323e+01 7.225885e+01
#> [791] 6.093333e+01 5.564508e+01 6.768882e+01 5.547205e+01 5.542377e+01
#> [796] 6.723807e+01 5.481761e+01 6.565259e+01 6.200705e+01 6.254698e+01
#> [801] 5.967568e+01 5.980686e+01 6.045464e+01 6.642233e+01 7.325786e+01
#> [806] 6.533046e+01 7.059961e+01 1.093179e+02 8.254054e+01 6.111267e+01
#> [811] 1.057113e+02 6.218005e+01 8.045206e+01 5.511948e+01 6.840295e+01
#> [816] 5.852141e+01 5.575714e+01 6.496614e+01 8.590522e+01 1.120092e+02
#> [821] 1.127239e+02 6.831783e+01 5.946455e+01 6.051781e+01 6.164526e+01
#> [826] 6.234749e+01 7.209662e+01 5.591092e+01 5.968686e+01 6.675258e+01
#> [831] 5.738323e+01 8.053653e+01 5.954092e+01 7.207861e+01 6.587337e+01
#> [836] 5.862071e+01 7.459712e+01 5.678540e+01 6.276934e+01 6.083662e+01
#> [841] 6.400852e+01 1.369590e+02 1.567997e+02 8.079597e+01 1.091060e+02
#> [846] 1.206943e+03 1.375685e+02 1.013156e+02 5.626902e+01 6.035941e+01
#> [851] 6.122133e+01 6.576108e+01 6.076685e+01 6.534866e+01 6.086371e+01
#> [856] 6.223082e+01 5.539022e+01 5.082659e+01 8.853785e+01 5.016195e+01
#> [861] 5.949811e+01 9.209960e+01 5.961611e+01 6.090262e+01 5.862822e+01
#> [866] 6.446316e+01 1.683468e+02 6.942663e+01 6.024661e+01 5.800444e+01
#> [871] 5.792102e+01 5.903047e+01 5.905064e+01 5.758922e+01 5.536758e+01
#> [876] 1.038374e+02 5.657780e+01 1.797931e+02 1.205173e+02 5.026042e+01
#> [881] 1.352633e+03 5.651420e+01 5.939997e+01 6.026428e+01 6.148543e+01
#> [886] 7.363359e+01 6.466744e+01 7.523982e+01 6.833530e+01 6.593705e+01
#> [891] 7.452247e+01 6.874107e+01 6.569660e+01 4.839410e+02 6.376657e+01
#> [896] 5.813968e+01 5.918713e+01 6.865285e+01 5.877639e+01 6.473897e+01
#> [901] 5.783480e+01 6.974691e+01 9.309600e+01 9.036875e+01 6.402417e+01
#> [906] 6.163400e+01 6.576398e+01 2.830947e+02 5.684054e+01 6.567669e+01
#> [911] 6.233642e+01 5.970481e+01 8.257443e+01 6.470069e+01 6.343986e+01
#> [916] 6.039908e+01 6.910450e+01 5.551620e+01 6.055000e+01 5.760997e+01
#> [921] 5.509072e+01 7.023425e+01 8.012904e+01 5.440013e+01 6.306826e+01
#> [926] 5.999752e+01 5.803511e+01 6.498109e+01 5.718575e+01 6.617442e+01
#> [931] 6.042074e+01 6.736858e+01 3.856940e+02 6.813931e+01 5.128529e+01
#> [936] 6.188917e+01 6.471570e+01 5.217348e+01 7.190115e+01 6.325551e+01
#> [941] 5.517755e+01 7.560331e+01 1.025642e+02 1.316912e+02 6.299246e+01
#> [946] 5.942899e+01 6.616008e+01 6.086578e+01 6.898343e+01 5.952806e+01
#> [951] 3.032688e+02 2.666652e+02 5.675221e+01 7.286205e+01 5.785451e+01
#> [956] 8.164655e+01 8.706421e+01 7.186908e+01 6.213114e+01 7.138507e+01
#> [961] 6.231327e+01 6.503073e+01 8.055056e+01 9.128891e+01 7.940133e+01
#> [966] 1.224274e+02 6.848340e+01 6.013670e+01 6.865575e+01 6.015158e+01
#> [971] 7.059110e+01 8.624808e+01 8.522626e+01 9.027253e+01 7.775198e+01
#> [976] 5.222805e+01 9.874144e+01 5.168036e+01 6.110357e+01 6.080223e+01
#> [981] 6.018341e+01 6.585143e+01 9.651483e+01 6.617380e+01 7.003466e+01
#> [986] 5.747238e+01 5.603502e+01 6.206126e+01 6.271644e+01 6.393178e+01
#> [991] 7.828052e+01 5.242859e+01 6.047430e+01 7.302376e+01 7.533527e+01
#> [996] 5.843593e+01 1.363329e+02 7.232135e+01 8.744138e+01 7.187755e+01