Skip to contents

[Experimental]

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, ...)

Arguments

model

(GeneralModel or ModelPseudo)
the model.

...

model specific parameters.

grade

(integer)
the toxicity grade for which the dose function is required

Value

A dose() method that computes doses.

Functions

  • doseFunction(GeneralModel):

  • doseFunction(ModelPseudo):

  • doseFunction(LogisticLogNormalOrdinal):

See also

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] 6.503084e+01 8.498407e+01 1.722265e+05 1.084608e+02 9.078398e+01
#>    [6] 1.181502e+02 7.027493e+01 7.434530e+01 1.324969e+02 6.394657e+01
#>   [11] 1.124658e+02 1.091526e+02 7.120897e+01 6.327834e+01 6.402245e+01
#>   [16] 7.769486e+01 6.891442e+01 6.453608e+01 7.005443e+01 6.153396e+01
#>   [21] 8.923055e+01 6.450479e+01 7.760224e+01 6.985643e+01 8.563263e+01
#>   [26] 6.384704e+01 6.466711e+01 6.653645e+01 8.308099e+01 7.840831e+01
#>   [31] 8.379316e+01 1.238230e+02 6.395280e+01 7.695206e+01 5.856745e+01
#>   [36] 6.771935e+01 8.929541e+01 7.751109e+01 7.365970e+01 1.578591e+02
#>   [41] 1.031466e+02 7.131603e+01 8.747489e+01 2.472047e+02 1.084568e+02
#>   [46] 6.090196e+01 8.630686e+01 7.371421e+01 7.697519e+01 1.165463e+02
#>   [51] 8.345645e+01 5.639928e+02 2.101983e+02 9.786847e+01 8.185470e+01
#>   [56] 7.669657e+01 1.791204e+02 2.796758e+02 6.848033e+01 6.687694e+02
#>   [61] 1.917765e+03 7.379869e+01 9.362901e+01 7.127939e+01 1.003004e+02
#>   [66] 9.207894e+01 1.053104e+02 1.137647e+02 2.473878e+02 6.325475e+01
#>   [71] 8.353038e+01 5.700847e+02 1.056529e+02 2.220344e+02 7.537199e+01
#>   [76] 6.536868e+01 8.322058e+01 1.316995e+02 6.258951e+01 7.727648e+01
#>   [81] 5.580148e+02 1.222888e+02 7.684023e+01 6.940840e+01 1.829506e+02
#>   [86] 9.330323e+01 1.017972e+02 1.009770e+02 1.184292e+02 1.535989e+02
#>   [91] 9.198828e+01 7.276377e+01 1.093656e+02 6.975872e+01 8.332459e+01
#>   [96] 9.252075e+01 7.698194e+01 1.291220e+02 6.209793e+01 1.150519e+04
#>  [101] 6.729705e+01 1.327842e+02 9.013334e+01 5.972891e+02 7.441872e+01
#>  [106] 8.070157e+01 8.907183e+01 2.618728e+02 4.245492e+02 1.149803e+03
#>  [111] 1.848492e+02 6.718782e+01 6.955089e+01 8.917288e+01 1.532461e+02
#>  [116] 8.127902e+01 1.556842e+02 6.747357e+01 8.603749e+02 7.582119e+01
#>  [121] 7.830263e+01 7.855383e+01 8.011230e+01 1.113070e+02 1.057506e+02
#>  [126] 1.427723e+02 7.027705e+01 1.285299e+02 1.557013e+02 1.104287e+02
#>  [131] 6.629406e+01 7.818412e+01 7.646418e+01 2.477438e+02 2.919791e+02
#>  [136] 3.789521e+02 3.758356e+02 7.422668e+01 1.030223e+02 1.231875e+02
#>  [141] 1.136044e+02 5.720861e+02 6.529358e+01 1.152568e+02 1.436768e+02
#>  [146] 1.109442e+02 7.159707e+01 1.027552e+02 7.826992e+01 8.248965e+01
#>  [151] 1.136497e+02 1.415924e+02 3.731579e+02 2.335800e+02 2.285431e+02
#>  [156] 1.272584e+02 7.363271e+01 6.942581e+01 3.482417e+02 1.532114e+02
#>  [161] 1.271255e+02 6.524294e+01 8.122912e+01 8.158074e+01 1.222241e+02
#>  [166] 1.599526e+02 2.856903e+02 7.599746e+01 1.332819e+02 8.176940e+01
#>  [171] 9.259351e+01 7.028315e+01 9.820258e+01 8.488094e+01 6.883887e+01
#>  [176] 9.495413e+01 2.122723e+02 1.207067e+03 7.432917e+01 1.060834e+02
#>  [181] 7.410315e+01 9.485360e+01 2.862606e+02 1.277007e+02 6.184486e+01
#>  [186] 2.220644e+02 6.278293e+01 3.106339e+02 3.008494e+03 6.720637e+01
#>  [191] 1.999328e+02 1.885044e+02 9.685263e+01 7.010346e+01 6.258795e+03
#>  [196] 7.459464e+01 6.245622e+01 8.148639e+01 6.918182e+01 8.059190e+01
#>  [201] 6.728309e+01 7.093447e+01 6.739035e+01 1.347430e+02 9.137718e+01
#>  [206] 8.057651e+01 7.666685e+01 8.992839e+01 8.200693e+01 9.541110e+01
#>  [211] 5.816691e+03 1.228327e+02 6.181771e+01 6.586417e+01 1.161272e+02
#>  [216] 7.909117e+01 8.846124e+01 9.592866e+01 2.188370e+02 2.346454e+02
#>  [221] 7.509577e+01 6.835090e+01 7.621722e+01 9.874818e+01 7.094198e+01
#>  [226] 9.709109e+01 1.258394e+02 1.019503e+02 7.501154e+01 1.556508e+02
#>  [231] 8.438213e+01 2.819699e+02 1.418858e+02 7.023149e+01 1.952700e+02
#>  [236] 1.191799e+02 2.219140e+02 1.232989e+02 7.852397e+01 7.187521e+01
#>  [241] 7.363604e+01 7.369231e+01 2.086001e+02 6.744804e+02 1.400254e+02
#>  [246] 1.063645e+02 8.228697e+01 1.139842e+02 9.281450e+01 1.366660e+02
#>  [251] 2.320574e+02 2.307426e+02 1.089490e+02 9.650974e+01 7.888544e+01
#>  [256] 1.011646e+02 2.109782e+02 1.181977e+02 8.340524e+01 6.861892e+01
#>  [261] 6.539948e+01 1.338659e+02 6.326899e+01 2.574105e+02 8.496280e+01
#>  [266] 6.952776e+01 1.156693e+02 1.435487e+02 6.809930e+01 1.039647e+02
#>  [271] 1.134985e+02 1.028105e+02 2.173632e+02 8.080296e+01 7.245273e+01
#>  [276] 8.069801e+01 8.118066e+01 7.216292e+01 9.182271e+01 1.381751e+02
#>  [281] 1.149461e+02 8.718183e+01 1.804098e+02 6.727142e+01 1.048407e+02
#>  [286] 7.171278e+01 6.947606e+01 1.045978e+02 1.239224e+02 5.015247e+02
#>  [291] 6.855526e+01 1.352003e+02 1.021910e+02 1.375786e+02 8.951363e+01
#>  [296] 6.640493e+01 3.112651e+02 4.127177e+03 3.041731e+02 7.035787e+01
#>  [301] 1.015686e+02 7.364512e+01 6.796977e+01 5.943621e+01 8.923011e+01
#>  [306] 1.577525e+02 1.629621e+02 8.218035e+01 9.135446e+01 7.986531e+01
#>  [311] 9.670040e+01 4.853841e+02 8.192245e+01 6.720938e+01 1.068993e+02
#>  [316] 2.597964e+02 7.510108e+01 7.207820e+01 1.980493e+02 1.286661e+02
#>  [321] 2.136060e+02 6.841335e+01 1.177406e+02 1.791463e+02 9.074243e+01
#>  [326] 2.702403e+02 9.715915e+01 8.290010e+01 1.259079e+02 7.857410e+01
#>  [331] 6.045756e+01 7.301783e+01 6.358667e+01 8.649604e+01 6.980643e+01
#>  [336] 6.388502e+01 3.132787e+02 6.698287e+01 1.966805e+02 6.854826e+01
#>  [341] 9.058445e+01 6.831299e+01 1.327763e+02 7.936325e+01 2.039416e+02
#>  [346] 1.247671e+02 6.931938e+01 4.062341e+03 8.208494e+01 8.075341e+01
#>  [351] 1.516907e+02 9.224136e+01 1.251799e+02 6.336984e+01 1.611782e+02
#>  [356] 7.360378e+01 1.576472e+02 1.112493e+02 4.703930e+03 3.626039e+02
#>  [361] 9.220015e+01 8.836313e+01 8.002635e+01 1.367394e+02 3.389248e+02
#>  [366] 1.308288e+02 7.963813e+01 9.076880e+01 1.744151e+02 6.520522e+01
#>  [371] 6.929842e+01 1.457986e+02 6.658277e+01 2.697688e+02 1.663014e+02
#>  [376] 7.842516e+01 9.521462e+01 7.574806e+01 2.253561e+02 3.693896e+03
#>  [381] 1.020538e+02 1.230892e+02 1.490889e+02 1.063988e+02 1.374652e+02
#>  [386] 5.152618e+02 5.754561e+02 5.066471e+02 8.476937e+01 7.325551e+01
#>  [391] 2.257513e+03 6.965076e+01 8.062518e+01 8.458528e+01 1.338855e+02
#>  [396] 7.285312e+01 2.082680e+02 9.952664e+01 6.566586e+01 7.791174e+01
#>  [401] 7.779504e+01 6.864016e+01 7.248550e+01 6.423324e+01 7.424535e+01
#>  [406] 1.603833e+02 1.062326e+02 1.072339e+02 6.616223e+01 7.243314e+01
#>  [411] 6.690499e+01 8.163824e+01 7.277554e+01 8.352850e+01 1.179314e+02
#>  [416] 6.805935e+01 8.401525e+01 1.101491e+02 2.335976e+02 7.657881e+01
#>  [421] 7.606132e+01 6.527077e+01 6.875819e+01 1.021488e+02 6.264788e+01
#>  [426] 1.038294e+02 1.433259e+02 6.639518e+01 2.413614e+02 1.261039e+02
#>  [431] 6.956612e+01 7.645859e+01 6.974067e+01 6.546398e+01 8.177875e+01
#>  [436] 1.047113e+02 6.124354e+01 6.790119e+01 6.616305e+01 7.224999e+01
#>  [441] 1.157478e+02 1.415222e+02 1.811120e+02 6.160772e+03 1.752306e+02
#>  [446] 3.917627e+02 1.260762e+02 7.492817e+01 9.686728e+01 1.179086e+02
#>  [451] 7.978477e+01 7.517379e+01 7.211045e+01 1.784269e+02 1.611152e+02
#>  [456] 2.556675e+02 1.062776e+02 7.380160e+01 6.896856e+01 7.064532e+01
#>  [461] 6.270236e+01 7.097065e+01 6.343107e+01 9.220503e+01 1.650248e+02
#>  [466] 4.627322e+02 1.111316e+02 8.922099e+01 6.821423e+01 8.164998e+01
#>  [471] 6.391668e+01 6.365634e+01 8.325735e+01 7.684420e+01 9.749441e+01
#>  [476] 1.390125e+02 1.239318e+02 5.634085e+02 6.152209e+01 9.083780e+01
#>  [481] 6.802851e+01 7.634628e+01 8.627444e+01 7.977474e+01 8.794724e+01
#>  [486] 6.526277e+01 7.312176e+01 6.939773e+01 7.674708e+01 7.291890e+01
#>  [491] 6.434157e+01 6.454216e+01 7.392111e+01 2.493895e+02 6.819520e+01
#>  [496] 8.919941e+01 1.903033e+02 8.139815e+01 7.570364e+01 1.042115e+02
#>  [501] 1.372302e+02 9.100375e+01 7.689756e+01 8.163560e+01 7.762703e+01
#>  [506] 6.085078e+01 7.070630e+01 6.177928e+01 6.235916e+01 7.049648e+01
#>  [511] 7.107035e+01 8.209730e+01 1.268741e+02 9.644858e+01 1.073805e+02
#>  [516] 8.612022e+01 9.308050e+01 1.518829e+02 9.127490e+01 7.056902e+01
#>  [521] 8.097765e+01 6.473212e+01 6.442351e+01 7.042769e+01 1.689398e+02
#>  [526] 8.189825e+01 1.031744e+02 1.815059e+02 9.062112e+01 6.669275e+01
#>  [531] 8.708299e+01 6.476419e+01 7.790505e+01 8.322525e+01 1.976810e+02
#>  [536] 9.557449e+01 7.557802e+01 8.054657e+01 1.021680e+02 1.032896e+02
#>  [541] 8.950749e+01 9.795868e+01 6.694175e+01 1.458403e+02 1.139801e+02
#>  [546] 5.693004e+04 5.022792e+02 8.388899e+01 1.781849e+02 8.330739e+01
#>  [551] 1.085215e+02 9.429007e+01 7.018223e+01 1.291206e+02 6.226269e+01
#>  [556] 7.755749e+01 1.437325e+02 7.281637e+01 6.699461e+01 6.988115e+01
#>  [561] 1.528451e+02 7.107760e+01 9.871291e+01 4.452729e+02 6.532630e+01
#>  [566] 6.419435e+01 4.342461e+03 8.373562e+01 6.807697e+01 2.172637e+02
#>  [571] 3.746384e+02 1.161669e+02 8.790881e+01 6.777657e+01 6.517945e+01
#>  [576] 7.682197e+01 9.404811e+01 1.181210e+02 1.126046e+02 1.067075e+02
#>  [581] 1.148417e+02 8.749214e+01 9.033284e+01 6.845147e+01 9.704365e+01
#>  [586] 7.858482e+01 9.497891e+01 7.916923e+01 1.308223e+02 9.974410e+01
#>  [591] 7.656605e+01 8.610852e+01 1.247100e+02 6.745917e+01 1.639858e+02
#>  [596] 7.006606e+01 8.470136e+01 8.361961e+01 1.193119e+02 4.923638e+02
#>  [601] 7.217552e+01 6.662284e+01 9.513087e+01 9.631167e+01 7.387117e+01
#>  [606] 4.173352e+02 9.479093e+01 1.100325e+02 7.291247e+01 6.589898e+01
#>  [611] 2.615167e+04 9.960697e+01 6.984198e+01 6.914031e+01 9.320362e+01
#>  [616] 7.256220e+01 7.809848e+01 9.959519e+01 3.622340e+02 8.230438e+01
#>  [621] 1.642883e+02 8.762355e+01 1.107786e+02 7.106525e+01 1.104771e+02
#>  [626] 7.630980e+01 7.194540e+01 7.399865e+01 8.127346e+01 6.320093e+01
#>  [631] 1.497209e+02 6.258932e+01 7.940668e+01 1.022388e+02 2.286406e+02
#>  [636] 1.033732e+02 2.481919e+02 3.264940e+02 2.347604e+03 7.284574e+01
#>  [641] 3.242767e+02 8.211219e+01 2.429147e+02 7.488144e+01 1.025128e+02
#>  [646] 6.871923e+01 9.680284e+01 8.008450e+01 1.146113e+02 8.699601e+01
#>  [651] 2.137526e+02 2.440099e+02 7.529815e+01 7.097522e+01 7.353278e+01
#>  [656] 7.533854e+01 7.203241e+01 6.779574e+01 7.730931e+01 7.127327e+01
#>  [661] 9.319975e+01 8.695814e+01 2.291523e+02 8.296521e+01 7.513481e+01
#>  [666] 7.614569e+01 8.207810e+01 7.122129e+01 7.820029e+01 7.032229e+01
#>  [671] 6.639852e+01 6.953890e+01 6.190981e+01 9.036169e+01 6.142469e+01
#>  [676] 1.449840e+02 9.800269e+01 7.122956e+01 1.098650e+02 2.545892e+02
#>  [681] 7.618136e+01 3.679906e+02 6.299330e+01 5.085663e+02 1.069025e+03
#>  [686] 3.847088e+02 1.362140e+02 1.884074e+02 8.899823e+01 7.013242e+01
#>  [691] 7.645829e+01 6.967685e+01 6.502035e+01 6.837303e+01 7.270277e+01
#>  [696] 1.603585e+02 4.558941e+02 1.231059e+02 2.622426e+02 6.818516e+01
#>  [701] 7.356042e+01 2.960673e+02 1.266247e+02 8.165447e+01 7.789208e+01
#>  [706] 7.978521e+01 8.645060e+01 9.217071e+01 6.830740e+01 7.019757e+01
#>  [711] 8.311072e+01 1.355279e+02 6.405546e+01 6.552920e+01 6.497825e+01
#>  [716] 1.000158e+02 1.043084e+02 1.461556e+03 6.311240e+01 1.383562e+02
#>  [721] 1.415537e+02 1.042629e+02 8.413007e+01 9.601770e+01 6.108690e+01
#>  [726] 1.896029e+02 1.080678e+02 7.539670e+01 1.148878e+02 6.244713e+01
#>  [731] 6.541991e+01 6.357502e+01 7.176918e+01 1.857017e+02 2.408103e+02
#>  [736] 7.888159e+01 6.730535e+01 6.871393e+01 6.755808e+01 9.540366e+01
#>  [741] 6.644645e+01 7.545907e+01 1.729810e+02 6.873647e+01 7.037858e+01
#>  [746] 7.324532e+02 5.382044e+03 5.398474e+02 7.796392e+01 6.479670e+01
#>  [751] 6.670543e+01 1.660615e+02 1.002886e+02 1.712647e+02 7.962938e+01
#>  [756] 1.397672e+04 2.498252e+02 1.158853e+02 6.994027e+01 5.278646e+02
#>  [761] 2.647985e+02 7.701714e+01 8.533231e+01 1.093354e+02 8.072386e+01
#>  [766] 8.976489e+01 2.477205e+02 8.787242e+01 9.298816e+01 1.201584e+02
#>  [771] 8.330275e+01 9.374226e+01 2.487042e+02 7.353362e+01 6.603288e+01
#>  [776] 7.319286e+01 7.118204e+01 6.401546e+01 5.982924e+02 3.180768e+03
#>  [781] 8.138040e+01 6.466809e+01 7.979642e+01 8.437622e+01 1.179181e+02
#>  [786] 9.138436e+01 8.357337e+01 1.822999e+02 8.735049e+01 7.206720e+01
#>  [791] 1.110641e+02 4.418442e+02 7.807747e+01 1.312790e+02 8.393741e+01
#>  [796] 7.392201e+01 1.036455e+02 2.548565e+02 9.643718e+01 9.197963e+01
#>  [801] 9.447140e+01 8.025034e+01 6.710468e+01 1.407838e+02 3.274664e+02
#>  [806] 8.000429e+01 1.413720e+02 9.985223e+02 1.379333e+02 6.835176e+01
#>  [811] 7.719082e+01 8.473138e+01 2.676810e+02 8.757422e+01 9.038539e+01
#>  [816] 7.451720e+01 2.798526e+02 9.543949e+01 1.245684e+02 1.176589e+02
#>  [821] 1.856986e+02 8.984961e+01 6.986043e+01 8.460862e+01 8.768672e+01
#>  [826] 2.529359e+02 2.991581e+02 9.471875e+01 8.506682e+01 1.076897e+02
#>  [831] 1.279349e+02 1.291407e+02 1.229130e+02 6.983130e+01 7.905521e+01
#>  [836] 3.638210e+08 1.045398e+03 7.995731e+01 1.205722e+02 2.157733e+02
#>  [841] 1.643172e+02 2.532609e+03 7.761262e+02 1.495781e+02 1.285079e+02
#>  [846] 8.965357e+01 1.615550e+02 7.168068e+01 1.224054e+02 6.582757e+01
#>  [851] 2.337056e+02 6.496156e+01 4.747713e+02 1.119015e+02 7.382246e+01
#>  [856] 7.424878e+01 3.409592e+02 1.562314e+06 1.842300e+02 6.403375e+01
#>  [861] 6.818002e+01 7.426231e+01 6.499770e+01 6.971435e+01 1.024475e+02
#>  [866] 1.070345e+02 8.363625e+01 2.725342e+02 9.977361e+01 7.119264e+01
#>  [871] 9.219679e+02 2.439390e+02 6.928054e+01 6.082265e+01 1.184518e+02
#>  [876] 2.185940e+02 1.747809e+02 6.688551e+01 7.081434e+01 6.579341e+01
#>  [881] 8.708891e+01 8.421004e+01 7.415997e+01 8.450406e+01 8.216443e+01
#>  [886] 6.477609e+01 6.377165e+01 3.073020e+02 7.486406e+01 8.284328e+01
#>  [891] 7.927487e+01 7.222010e+01 1.191394e+02 1.099068e+02 2.663087e+02
#>  [896] 1.077231e+02 8.511155e+01 1.287698e+02 7.203135e+01 1.252677e+02
#>  [901] 1.103159e+02 9.150984e+01 7.418781e+01 1.010868e+02 1.030331e+02
#>  [906] 7.240335e+02 6.429677e+01 7.668686e+01 1.135892e+02 7.461102e+01
#>  [911] 6.880593e+01 8.660400e+01 3.943866e+02 1.126769e+02 2.400539e+03
#>  [916] 1.051064e+02 6.890446e+01 6.402680e+01 1.009559e+02 2.106453e+02
#>  [921] 8.223901e+02 1.430304e+03 1.603218e+02 8.457824e+01 8.526568e+01
#>  [926] 1.695698e+02 6.477049e+01 7.834650e+01 1.546572e+02 8.188970e+01
#>  [931] 6.695854e+01 7.191631e+01 6.336947e+01 7.661768e+01 7.436798e+01
#>  [936] 1.389268e+02 5.013567e+02 1.527089e+02 8.714944e+01 7.577238e+01
#>  [941] 2.631227e+02 1.083238e+02 9.473158e+01 7.838267e+01 1.436712e+02
#>  [946] 1.944142e+02 6.756685e+01 8.301404e+01 8.171365e+01 2.199446e+02
#>  [951] 6.994482e+01 6.484412e+01 7.840281e+01 1.223328e+02 1.461355e+02
#>  [956] 1.519540e+05 1.720863e+03 6.978114e+01 6.902791e+01 7.450780e+01
#>  [961] 7.275311e+01 5.451895e+02 1.056436e+02 8.652676e+01 7.175402e+01
#>  [966] 7.489085e+01 2.242256e+02 1.591803e+02 1.098188e+02 7.139546e+01
#>  [971] 1.071567e+02 7.240830e+01 8.124014e+01 7.033155e+01 1.615390e+02
#>  [976] 1.406582e+02 1.491478e+02 6.978217e+01 6.914256e+01 6.959815e+01
#>  [981] 6.524605e+01 1.029824e+02 9.533460e+01 1.162713e+02 6.570532e+02
#>  [986] 6.418485e+01 6.753631e+01 6.638922e+01 7.258883e+01 8.717808e+01
#>  [991] 9.107560e+01 6.650324e+01 2.128055e+02 6.472440e+01 8.099975e+01
#>  [996] 6.718579e+01 7.002967e+01 7.063852e+01 9.828024e+01 8.441626e+01
doseFunction(
  model,
  alpha2 = samples@data$alpha2,
  beta = samples@data$beta,
  grade = 2L
)(x = 0.25)
#>    [1]     61.08100     62.29098   2356.22065     72.32862     58.28887
#>    [6]     74.79160     59.35864     62.94879     86.62611     58.87231
#>   [11]     91.83389     85.39065     63.10887     57.85366     58.92140
#>   [16]     68.04239     62.16223     57.72847     62.27789     54.69438
#>   [21]     59.75792     56.62605     60.39671     60.52880     61.94425
#>   [26]     53.52051     50.88591     59.51032     67.71141     67.85491
#>   [31]     66.56328     75.89371     56.91408     62.11852     54.59148
#>   [36]     60.39182     49.62354     61.00253     57.88572     87.79192
#>   [41]     76.73953     60.50293     70.93297    118.10356     62.78856
#>   [46]     53.70962     70.86817     65.47367     67.39473     72.19507
#>   [51]     68.71276    129.26540     89.53180     67.43814     65.06243
#>   [56]     61.53395     91.62348     73.12801     57.88077     79.64608
#>   [61]    465.88584     62.16860     67.46002     59.16017     56.78452
#>   [66]     64.99996     80.28496     74.63145    102.92039     52.25250
#>   [71]     67.95220    135.48885     76.97631    130.45689     66.29687
#>   [76]     56.67327     62.11758     86.39531     55.69880     63.85406
#>   [81]    135.80263     71.99302     59.84180     61.65266    114.81193
#>   [86]     71.98142     68.65283     74.70747     76.06877     97.36187
#>   [91]     69.69893     57.34365     66.20089     59.13266     64.42450
#>   [96]     65.50753     59.73033     79.08746     52.83364   1790.59040
#>  [101]     53.73441     69.34586     73.47613    173.90690     61.61088
#>  [106]     62.54576     69.81254     91.73745    162.78013    178.77125
#>  [111]     77.50014     49.58902     61.74283     71.36458     98.78450
#>  [116]     60.33244     97.48345     60.27020    219.52497     63.54698
#>  [121]     64.67122     62.28462     63.13807     69.75704     73.11506
#>  [126]     77.01405     55.69097     89.74344     84.32646     74.43006
#>  [131]     55.41671     64.76834     56.05973    143.02396    117.37478
#>  [136]     76.89265     93.13613     53.14483     56.45978     56.39357
#>  [141]     69.65891    199.47578     55.98482     74.16786     70.55159
#>  [146]     64.56670     60.17895     80.55830     67.47581     68.08584
#>  [151]     79.92459     87.37771    120.38134     84.79044     76.60114
#>  [156]     66.57059     50.31297     62.13660    178.72330     73.74377
#>  [161]     84.21755     56.54998     55.43577     63.48429     88.30328
#>  [166]     91.13677    125.01044     58.04696     69.98179     49.30608
#>  [171]     66.79088     59.53986     78.70978     65.09822     56.41847
#>  [176]     73.97903     97.44015    380.97064     61.64581     58.07620
#>  [181]     60.45610     63.09521    168.62175     97.98601     55.97403
#>  [186]     99.04591     49.56957    115.02284    571.34068     51.69408
#>  [191]    127.05248     80.99301     71.11237     57.68820    225.56653
#>  [196]     59.93169     56.91888     68.21239     60.78337     68.04325
#>  [201]     59.58625     61.11208     61.09727     90.37592     60.18243
#>  [206]     61.28164     61.63712     63.41514     54.37704     58.85548
#>  [211]    455.74899     65.73591     54.56856     56.61811     78.48533
#>  [216]     45.53020     52.19277     76.17375    136.25017     84.77642
#>  [221]     57.89452     56.88641     63.43799     58.81249     60.19894
#>  [226]     72.30441     86.82859     78.57239     63.85924     97.70370
#>  [231]     55.79472    117.00638     79.21934     59.94028     95.91019
#>  [236]     80.45092     77.79783     74.11009     53.58855     62.64850
#>  [241]     60.23789     61.31816    112.83884    195.00158     73.91248
#>  [246]     79.32687     59.92168     58.71383     71.66608     79.33233
#>  [251]    130.36906     97.98425     78.60348     62.81726     59.21488
#>  [256]     68.11838    103.86902     73.50870     57.26869     54.58172
#>  [261]     52.60286     79.45358     52.58811    122.18044     53.80471
#>  [266]     55.27902     82.97086     96.03035     55.45457     76.85216
#>  [271]     76.90842     67.54799    121.37921     67.27203     60.42328
#>  [276]     63.15571     55.12349     57.93519     65.45048     76.17412
#>  [281]     62.81690     64.22104     75.14493     53.87690     77.09381
#>  [286]     61.58607     59.13814     64.48749     57.53051    172.96458
#>  [291]     58.84086     80.53426     66.53095     88.59013     71.33340
#>  [296]     60.52385     97.27153    367.80081    134.84578     58.48717
#>  [301]     75.81576     65.75247     62.90614     55.35557     63.84231
#>  [306]     98.33595     84.88529     66.90610     76.91857     68.46775
#>  [311]     77.39054    256.67065     64.95079     58.21746     74.26075
#>  [316]    122.97316     63.27375     60.28218     91.52841     52.52688
#>  [321]     45.42063     57.57504     75.59099     88.43642     65.53793
#>  [326]     64.69494     76.29768     60.52613     70.19953     63.84978
#>  [331]     52.46607     59.58620     56.48989     63.92479     60.43171
#>  [336]     55.23018    162.04085     59.86762     82.69550     58.96208
#>  [341]     69.48918     58.73150     87.39772     63.27988     89.46739
#>  [346]     73.32195     59.16608    421.63120     58.50830     68.72937
#>  [351]     95.45801     67.85200     91.67151     56.63442     90.97178
#>  [356]     58.30022     80.69412     63.28525    346.65163    127.46974
#>  [361]     65.06252     63.34470     64.27614     78.96590     91.63868
#>  [366]     70.19033     52.12706     68.13839     99.09690     51.11851
#>  [371]     55.35275     97.33867     57.95028    128.60631     98.89128
#>  [376]     59.34033     69.81908     52.21771    126.95853   1141.37482
#>  [381]     68.43605     76.05988     69.83342     66.17083     86.89422
#>  [386]    205.82803    133.83031    187.17896     64.81152     44.65569
#>  [391]    166.57145     58.97239     67.26975     69.80944     83.79976
#>  [396]     63.89147    111.74531     76.16972     60.35172     65.06964
#>  [401]     59.87228     60.39405     61.52071     57.58124     65.58403
#>  [406]    100.35284     71.30791     61.22407     59.49754     63.16256
#>  [411]     59.67186     67.58905     62.21140     67.89865     69.59155
#>  [416]     59.91636     64.63877     69.97075    132.35967     66.26370
#>  [421]     67.08476     59.60357     60.38398     75.62618     54.02015
#>  [426]     70.14879     63.69428     53.81651     90.69363     98.69213
#>  [431]     62.61705     62.73145     60.60734     58.11160     67.36500
#>  [436]     74.85330     55.74881     61.04283     58.92187     61.50550
#>  [441]     80.19943    103.28309    115.90279    796.25305     75.51174
#>  [446]    162.57353     83.46711     62.44644     67.29815     64.65025
#>  [451]     56.30577     64.25510     59.10955     86.72119     77.04104
#>  [456]     52.40649     58.41241     51.77015     57.22661     55.75381
#>  [461]     54.59883     62.70907     56.63676     63.10571     58.38679
#>  [466]    152.19884     78.70994     72.97645     59.64718     62.32703
#>  [471]     54.35683     57.22731     67.41235     65.45023     60.31427
#>  [476]     86.22689     85.11519    109.90902     52.11413     65.26739
#>  [481]     58.22779     59.33329     68.79781     65.60303     62.26562
#>  [486]     56.92928     57.80962     61.49966     66.84691     63.62869
#>  [491]     58.69225     59.20122     62.65286    125.84176     53.95272
#>  [496]     61.92287     89.16984     64.85920     61.37272     75.64143
#>  [501]     76.03663     64.22305     56.20439     56.55714     60.98649
#>  [506]     56.32103     63.29939     57.86155     58.89549     63.02455
#>  [511]     61.02760     66.85858     82.52539     72.68608     78.22264
#>  [516]     61.70015     58.45264     72.60717     64.58831     56.16593
#>  [521]     66.45167     57.27331     57.91219     59.40754     93.89975
#>  [526]     63.12671     68.80883    102.62730     67.52673     58.48963
#>  [531]     72.67495     55.66646     62.72475     60.03801     92.25448
#>  [536]     76.99146     64.24743     64.46529     72.17570     55.79802
#>  [541]     64.75536     69.86541     58.27326     69.68346     63.29205
#>  [546]   2520.41960    130.74725     64.98434     74.84949     61.44960
#>  [551]     68.89111     60.63990     57.70754     87.84794     55.08125
#>  [556]     58.58615     67.56210     60.13065     58.43811     60.71788
#>  [561]     94.53880     58.24363     64.77537    155.24244     57.98855
#>  [566]     57.88173    668.10024     65.76409     57.59547     79.54412
#>  [571]    124.02820     63.99236     61.42455     58.61295     56.04277
#>  [576]     68.07988     77.98887     81.88178     73.44029     64.50256
#>  [581]     65.19236     59.34142     60.54963     56.00571     75.21556
#>  [586]     62.58214     57.28653     64.84208     92.04146     69.07745
#>  [591]     52.76798     58.14627     90.55416     58.07792    117.09201
#>  [596]     59.19536     64.64901     63.39586     88.19784    138.13525
#>  [601]     62.08266     58.15082     77.68040     66.97460     63.59482
#>  [606]    181.64608     64.87266     76.61458     60.00749     57.93023
#>  [611]   2330.26861     71.39056     56.65739     57.26085     67.18603
#>  [616]     56.69625     64.27178     68.02544    112.39076     64.37627
#>  [621]     85.25221     60.76792     62.92096     54.87106     74.50763
#>  [626]     67.44851     64.61275     63.63610     67.23917     56.74643
#>  [631]     96.43539     52.39040     66.64578     74.36111    124.91971
#>  [636]     70.44565    104.50188    109.93690    402.19291     60.59828
#>  [641]    125.27617     64.28665    126.05668     54.02439     63.07201
#>  [646]     61.08647     64.43233     61.74395     78.13275     62.88196
#>  [651]     72.69387     98.14423     61.34235     60.48626     58.31563
#>  [656]     60.80778     60.30719     57.80779     62.27146     59.69152
#>  [661]     66.08535     61.67519     87.71646     69.61711     64.96076
#>  [666]     61.57603     68.73875     63.84176     64.68608     61.39703
#>  [671]     57.76530     55.69630     55.91002     75.72533     53.49283
#>  [676]     89.07314     61.62535     53.03241     52.19541     89.25369
#>  [681]     61.02140    100.60019     54.30727    161.22499    208.13821
#>  [686]    174.20099     66.49770     80.80717     65.99059     57.49574
#>  [691]     64.45716     60.01207     58.26244     60.55026     58.19310
#>  [696]     78.37888    124.70767     60.12604    139.70998     55.89978
#>  [701]     60.38103     90.78991     81.31184     61.94311     54.14750
#>  [706]     65.42680     69.39416     71.08502     57.58545     59.10386
#>  [711]     63.17327    100.98143     58.63117     60.03809     59.36876
#>  [716]     75.42900     61.29786    201.53553     56.81184     83.58003
#>  [721]     88.89785     69.21768     59.57720     61.58865     52.52218
#>  [726]    127.34669     79.00857     64.89842     87.67311     56.53849
#>  [731]     58.17319     56.70580     61.73418    100.69790    108.45159
#>  [736]     60.02688     59.50855     60.18946     58.39352     73.63229
#>  [741]     57.29523     58.26811     83.15466     52.74852     57.85371
#>  [746]     85.06813    301.99801    134.35156     59.29692     57.14842
#>  [751]     57.88373    104.77193     55.99638     88.44482     63.41197
#>  [756]   1590.45475    146.24427     83.30913     61.43872    207.53677
#>  [761]    123.14590     61.99321     65.56703     77.97139     53.41425
#>  [766]     64.57049    100.37285     50.32805     74.12900     80.83929
#>  [771]     62.58623     65.79094    135.33403     64.42144     58.12886
#>  [776]     62.97625     60.97151     55.49273    131.56716    541.78394
#>  [781]     68.68970     57.66216     65.50377     61.61873     88.99547
#>  [786]     74.02569     70.45126     91.46434     66.71149     64.37102
#>  [791]     87.41715    106.57425     59.16050     79.91373     62.30650
#>  [796]     61.68212     70.19701     81.78097     71.79327     66.00771
#>  [801]     73.07180     62.62332     60.16871     94.19858    163.85535
#>  [806]     66.75176     82.15622    282.42482     81.29047     58.29093
#>  [811]     63.79732     71.57178    119.13622     61.23546     65.48433
#>  [816]     65.31045    141.93445     70.88340     69.87558     81.70710
#>  [821]     99.64518     64.51334     56.09308     64.02667     59.57846
#>  [826]    112.86541     67.24724     68.81670     66.44365     60.24856
#>  [831]     73.09153     90.56698     79.31639     56.82367     61.35175
#>  [836] 160019.80623     94.40235     59.68513     80.17305    106.88812
#>  [841]     75.12708    251.75442    158.91092     74.95595     60.81649
#>  [846]     71.57901     75.10836     63.63385     81.04028     59.34567
#>  [851]    146.47330     59.62819    145.22647     87.52895     64.00863
#>  [856]     62.53791    121.90885  14036.36797    101.89404     58.93454
#>  [861]     60.28326     63.81888     57.65897     60.54015     78.16703
#>  [866]     75.25180     59.85945    107.19185     67.62128     57.23618
#>  [871]    269.09259    137.05537     59.97401     56.23009     83.70276
#>  [876]    121.73893     71.71347     59.66458     57.28111     54.84079
#>  [881]     72.75248     56.38991     58.11126     62.87075     59.09017
#>  [886]     56.11707     58.18926    162.37963     63.44906     68.74929
#>  [891]     68.00951     63.59841     68.30167     65.82200     69.27882
#>  [896]     72.80528     53.95653     69.25089     60.62041     68.40143
#>  [901]     78.91234     69.94904     64.20078     70.59836     75.50634
#>  [906]    183.87399     55.47510     61.13929     78.50650     61.68485
#>  [911]     61.23385     67.76580    156.81099     54.19878    541.98892
#>  [916]     71.60963     56.15502     56.98920     69.30086    117.70733
#>  [921]    147.06604    343.52341     90.77724     54.05169     67.49005
#>  [926]     82.99543     58.22144     63.45649     77.43407     66.59963
#>  [931]     54.21869     62.08075     55.59251     66.83760     61.97439
#>  [936]     97.85230    166.56819    110.02621     61.76693     53.54844
#>  [941]     88.49396     56.12563     69.67839     61.69356     83.91956
#>  [946]    123.25313     60.73295     68.60133     65.43058     94.14106
#>  [951]     58.33678     58.53402     63.17047     85.93014     78.29230
#>  [956]   6532.90948    105.32298     56.14132     57.85018     64.40839
#>  [961]     62.93660    149.50483     77.49391     67.29442     60.13104
#>  [966]     61.23162    123.28982     93.78594     70.16690     58.58835
#>  [971]     72.03812     51.16720     65.25137     59.09269     87.63349
#>  [976]     71.72238     81.42506     60.02916     61.35884     59.37729
#>  [981]     57.19125     54.37720     76.67099     75.55840    205.04046
#>  [986]     57.06185     59.82791     57.39208     59.91060     63.00882
#>  [991]     67.26693     58.40903    120.09357     58.11795     60.81948
#>  [996]     56.83358     58.98511     61.96270     74.58539     63.67324