Optimizer functions for gradient and likelihood boosting with bamlss. In each boosting iteration the function selects the model term with the largest contribution to the log-likelihood, AIC or BIC.

## Gradient boosting optimizer.
opt_boost(x, y, family, weights = NULL,
  offset = NULL, nu = 0.1, nu.adapt = TRUE, df = 4, maxit = 400,
  mstop = NULL, maxq = NULL, qsel.splitfactor = FALSE,
  verbose = TRUE, digits = 4, flush = TRUE,
  eps = .Machine$double.eps^0.25,
  nback = NULL, plot = TRUE, initialize = TRUE,
  stop.criterion = NULL, select.type = 1,
  force.stop = TRUE, hatmatrix = !is.null(stop.criterion),
  reverse.edf = FALSE, approx.edf = FALSE,
  always = FALSE, ...)

boost(x, y, family, weights = NULL,
  offset = NULL, nu = 0.1, nu.adapt = TRUE, df = 4, maxit = 400,
  mstop = NULL, maxq = NULL, qsel.splitfactor = FALSE,
  verbose = TRUE, digits = 4, flush = TRUE,
  eps = .Machine$double.eps^0.25,
  nback = NULL, plot = TRUE, initialize = TRUE,
  stop.criterion = NULL, select.type = 1,
  force.stop = TRUE, hatmatrix = !is.null(stop.criterion),
  reverse.edf = FALSE, approx.edf = FALSE,
  always = FALSE, ...)

## Modified likelihood based boosting.
opt_boostm(x, y, family, offset = NULL,
  nu = 0.1, df = 3, maxit = 400, mstop = NULL,
  verbose = TRUE, digits = 4, flush = TRUE,
  eps = .Machine$double.eps^0.25, plot = TRUE,
  initialize = TRUE, stop.criterion = "BIC",
  force.stop = !is.null(stop.criterion),
  do.optim = TRUE, always = FALSE, ...)

boostm(x, y, family, offset = NULL,
  nu = 0.1, df = 3, maxit = 400, mstop = NULL,
  verbose = TRUE, digits = 4, flush = TRUE,
  eps = .Machine$double.eps^0.25, plot = TRUE,
  initialize = TRUE, stop.criterion = "BIC",
  force.stop = !is.null(stop.criterion),
  do.optim = TRUE, always = FALSE, ...)

## Boosting summary extractor.
boost_summary(object, ...)

## Plot all boosting paths.
boost_plot(x, which = c("loglik", "loglik.contrib", "parameters",
  "aic", "bic", "user"), intercept = TRUE, spar = TRUE, mstop = NULL,
  name = NULL, drop = NULL, labels = NULL, color = NULL, ...)

## Boosting summary printing and plotting.
# S3 method for boost_summary
print(x, summary = TRUE, plot = TRUE,
  which = c("loglik", "loglik.contrib"), intercept = TRUE,
  spar = TRUE, ...)
# S3 method for boost_summary
plot(x, ...)

## Model frame for out-of-sample selection.
boost_frame(formula, train, test, family = "gaussian", ...)

Arguments

x

For function opt_boost() the x list, as returned from function bamlss.frame, holding all model matrices and other information that is used for fitting the model. For the plotting function the corresponding bamlss object fitted with the opt_boost() optimizer.

y

The model response, as returned from function bamlss.frame.

family

A bamlss family object, see family.bamlss.

weights

Prior weights on the data, as returned from function bamlss.frame.

offset

Can be used to supply model offsets for use in fitting, returned from function bamlss.frame.

nu

Numeric, between [0, 1], controls the step size, i.e., the amount that should be added to model term parameters.

nu.adapt

Logical. If set to TRUE (default) step size nu is divided by 2, if current boosting iteration did not improve the loglikelihood.

df

Integer, defines the initial degrees of freedom that should be assigned to each smooth model term. May also be a named vector, the names must match the model term labels, e.g., as provided in summary.bamlss.

maxit

Integer, the maximum number of boosting iterations.

mstop

For convenience, overwrites maxit.

maxq

Integer, defines the maximum number of selected base-learners. The algorithm stops if this numer is exceeded.

qsel.splitfactor

Logical, if set to TRUE dummy variables of categorical predictors are counted individually.

name

Character, the name of the coefficient (group) that should be plotted. Note that the string provided in name will be removed from the labels on the 4th axis.

drop

Character, the name of the coefficient (group) that should not be plotted.

labels

A character string of labels that should be used on the 4 axis.

color

Colors or color function that creates colors for the (group) paths.

verbose

Print information during runtime of the algorithm.

digits

Set the digits for printing when verbose = TRUE.

flush

use flush.console for displaying the current output in the console.

eps

The tolerance used as stopping mechanism, see argument nback.

nback

Integer. If nback is not NULL, then the algorithm stops if the the change in the log-likelihood of the last nback iterations is smaller or equal to eps. If maxit = NULL the maximum number of iterations is set to 10000.

plot

Should the boosting summary be printed and plotted?

initialize

Logical, should intercepts be initialized?

stop.criterion

Character, selects the information criterion that should be used to determine the optimum number of boosting iterations. Either "AIC" or "BIC" is possible. Note that this feature requires to compute hat-matrices for each distributional parameter, therefore, the routine may be slow and computer storage intensive.

select.type

Should model terms be selected by the log-likelihood contribution, select.type = 1, or by the corresponding stop.criterion, select.type = 2.

force.stop

Logical, should the algorithm stop if the information criterion increases?

do.optim

Logical. Should smoothing parameters be optimized in each boosting iteration?

hatmatrix

Logical, if set to TRUE the hat-matrices for each distributional parameter will be computed. The hat-matrices are used to determine the effective (equivalent) degrees of freedom in each boosting iteration, i.e., it is possible to compute information criteria like the AIC or BIC for selecting the optimum number of boosting iterations.

reverse.edf

Logical. Instead of computing degrees of freedom with hat-matrices, the actual smoothing parameters are reverse engineered to compute the corresponding actual smoother matrix. Note that this option is still experimental.

approx.edf

Logical. Another experimental and fast approximation of the degrees of freedom.

always

Logical or character. Should the intercepts forced to be updated in each boosting iteration? If always = TRUE each intercept of each distributional parameter is updated, if always = "best" only the intercept corresponding to the distributional of the best fitting model term is updated.

object

A bamlss object that was fitted using opt_boost().

summary

Should the summary be printed?

which

Which of the three provided plots should be created?

intercept

Should the coefficient paths of intercepts be dropped in the plot?

spar

Should graphical parmeters be set with par?

formula

See bamlss.frame.

train, test

Data frames used for training and testing the model..

...

For function opt_boost(), arguments passed to bamlss.engine.setup. for function boost_summary() arguments passed to function print.boost_summary().

Value

For function boost_summary() a list containing information on selection frequencies etc. For function opt_boost() and opt_boostm() a list containing the following objects:

fitted.values

A named list of the fitted values based on the last boosting iteration of the modeled parameters of the selected distribution.

parameters

A matrix, each row corresponds to the parameter values of one boosting iteration.

boost_summary

The boosting summary which can be printed and plotted.

WARNINGS

The function does not take care of variable scaling for the linear parts! This must be done by the user, e.g., one option is to use argument scale.d in function bamlss.frame, which uses scale.

Function opt_boost() does not select the optimum stopping iteration! The modified likelihood based algorithm implemented in function opt_boostm() is still experimental!

See also

Examples

if (FALSE) ## Simulate data.
set.seed(123)
d <- GAMart()

## Estimate model.
f <- num ~ x1 + x2 + x3 + lon + lat +
  s(x1) + s(x2) + s(x3) + s(lon) + s(lat) + te(lon,lat)

b <- bamlss(f, data = d, optimizer = opt_boost,
  sampler = FALSE, scale.d = TRUE, nu = 0.01,
  maxit = 1000, plot = FALSE)
#> 
#> logLik -181.739 eps 0.0123 iteration    2 qsel 1
#> logLik -180.556 eps 0.0122 iteration    3 qsel 1
#> logLik -179.396 eps 0.0122 iteration    4 qsel 1
#> logLik -178.260 eps 0.0122 iteration    5 qsel 1
#> logLik -177.145 eps 0.0123 iteration    6 qsel 1
#> logLik -176.053 eps 0.0123 iteration    7 qsel 1
#> logLik -174.983 eps 0.0125 iteration    8 qsel 1
#> logLik -173.933 eps 0.0126 iteration    9 qsel 1
#> logLik -172.905 eps 0.0129 iteration   10 qsel 1
#> logLik -171.896 eps 0.0132 iteration   11 qsel 1
#> logLik -170.908 eps 0.0136 iteration   12 qsel 1
#> logLik -169.939 eps 0.0141 iteration   13 qsel 1
#> logLik -168.989 eps 0.0147 iteration   14 qsel 1
#> logLik -168.018 eps 0.0138 iteration   15 qsel 2
#> logLik -167.087 eps 0.0157 iteration   16 qsel 2
#> logLik -166.134 eps 0.0146 iteration   17 qsel 2
#> logLik -165.220 eps 0.0176 iteration   18 qsel 2
#> logLik -164.286 eps 0.0203 iteration   19 qsel 2
#> logLik -163.390 eps 0.0223 iteration   20 qsel 2
#> logLik -162.474 eps 0.0231 iteration   21 qsel 2
#> logLik -161.595 eps 0.1738 iteration   22 qsel 2
#> logLik -160.697 eps 0.0237 iteration   23 qsel 2
#> logLik -159.835 eps 0.0786 iteration   24 qsel 2
#> logLik -158.953 eps 0.0207 iteration   25 qsel 2
#> logLik -158.108 eps 0.0292 iteration   26 qsel 2
#> logLik -157.244 eps 0.0249 iteration   27 qsel 2
#> logLik -156.415 eps 0.0497 iteration   28 qsel 2
#> logLik -155.567 eps 0.0285 iteration   29 qsel 2
#> logLik -154.754 eps 0.0354 iteration   30 qsel 2
#> logLik -153.922 eps 0.0373 iteration   31 qsel 2
#> logLik -153.125 eps 0.0446 iteration   32 qsel 2
#> logLik -152.309 eps 0.0276 iteration   33 qsel 2
#> logLik -151.527 eps 0.0355 iteration   34 qsel 2
#> logLik -150.727 eps 0.0346 iteration   35 qsel 2
#> logLik -149.960 eps 0.0582 iteration   36 qsel 2
#> logLik -149.176 eps 0.0391 iteration   37 qsel 2
#> logLik -148.423 eps 0.0557 iteration   38 qsel 2
#> logLik -147.654 eps 0.0988 iteration   39 qsel 2
#> logLik -146.916 eps 0.0856 iteration   40 qsel 2
#> logLik -146.162 eps 0.0470 iteration   41 qsel 2
#> logLik -145.438 eps 0.0415 iteration   42 qsel 2
#> logLik -144.698 eps 0.0734 iteration   43 qsel 2
#> logLik -143.988 eps 0.0699 iteration   44 qsel 2
#> logLik -143.262 eps 0.4393 iteration   45 qsel 2
#> logLik -142.566 eps 0.0535 iteration   46 qsel 2
#> logLik -141.854 eps 0.0478 iteration   47 qsel 2
#> logLik -141.171 eps 0.0869 iteration   48 qsel 2
#> logLik -140.473 eps 0.0449 iteration   49 qsel 2
#> logLik -139.803 eps 0.0610 iteration   50 qsel 2
#> logLik -139.118 eps 0.1395 iteration   51 qsel 2
#> logLik -138.461 eps 0.0702 iteration   52 qsel 2
#> logLik -137.790 eps 0.0446 iteration   53 qsel 2
#> logLik -137.145 eps 0.0364 iteration   54 qsel 2
#> logLik -136.487 eps 0.0337 iteration   55 qsel 2
#> logLik -135.855 eps 0.0329 iteration   56 qsel 2
#> logLik -135.209 eps 0.0245 iteration   57 qsel 2
#> logLik -134.589 eps 0.0234 iteration   58 qsel 2
#> logLik -133.955 eps 0.0323 iteration   59 qsel 2
#> logLik -133.347 eps 0.0339 iteration   60 qsel 2
#> logLik -132.726 eps 0.0256 iteration   61 qsel 2
#> logLik -132.129 eps 0.0234 iteration   62 qsel 2
#> logLik -131.520 eps 0.0300 iteration   63 qsel 2
#> logLik -130.853 eps 0.0553 iteration   64 qsel 3
#> logLik -130.268 eps 0.0280 iteration   65 qsel 3
#> logLik -129.671 eps 0.0327 iteration   66 qsel 3
#> logLik -129.017 eps 0.0293 iteration   67 qsel 3
#> logLik -128.443 eps 0.0255 iteration   68 qsel 3
#> logLik -127.859 eps 0.0414 iteration   69 qsel 3
#> logLik -127.218 eps 0.0252 iteration   70 qsel 3
#> logLik -126.656 eps 0.0258 iteration   71 qsel 3
#> logLik -126.083 eps 0.0365 iteration   72 qsel 3
#> logLik -125.455 eps 0.0448 iteration   73 qsel 3
#> logLik -124.904 eps 0.0654 iteration   74 qsel 3
#> logLik -124.343 eps 0.0241 iteration   75 qsel 3
#> logLik -123.728 eps 0.0333 iteration   76 qsel 3
#> logLik -123.188 eps 0.0256 iteration   77 qsel 3
#> logLik -122.638 eps 0.0335 iteration   78 qsel 3
#> logLik -122.035 eps 0.0458 iteration   79 qsel 3
#> logLik -121.506 eps 0.0365 iteration   80 qsel 3
#> logLik -120.967 eps 0.0203 iteration   81 qsel 3
#> logLik -120.448 eps 0.0174 iteration   82 qsel 3
#> logLik -119.857 eps 0.0268 iteration   83 qsel 3
#> logLik -119.330 eps 0.0205 iteration   84 qsel 3
#> logLik -118.821 eps 0.0200 iteration   85 qsel 3
#> logLik -118.242 eps 0.0210 iteration   86 qsel 3
#> logLik -117.725 eps 0.0227 iteration   87 qsel 3
#> logLik -117.226 eps 0.0150 iteration   88 qsel 3
#> logLik -116.659 eps 0.0312 iteration   89 qsel 3
#> logLik -116.153 eps 0.0155 iteration   90 qsel 3
#> logLik -115.664 eps 0.0157 iteration   91 qsel 3
#> logLik -115.108 eps 0.0170 iteration   92 qsel 3
#> logLik -114.612 eps 0.0155 iteration   93 qsel 3
#> logLik -114.132 eps 0.0129 iteration   94 qsel 3
#> logLik -113.587 eps 0.0271 iteration   95 qsel 3
#> logLik -113.102 eps 0.0142 iteration   96 qsel 3
#> logLik -112.632 eps 0.0123 iteration   97 qsel 3
#> logLik -112.155 eps 0.0139 iteration   98 qsel 3
#> logLik -111.622 eps 0.0148 iteration   99 qsel 3
#> logLik -111.161 eps 0.0118 iteration  100 qsel 3
#> logLik -110.694 eps 0.0148 iteration  101 qsel 3
#> logLik -110.171 eps 0.0653 iteration  102 qsel 3
#> logLik -109.719 eps 0.0129 iteration  103 qsel 3
#> logLik -109.262 eps 0.0202 iteration  104 qsel 3
#> logLik -108.750 eps 0.0232 iteration  105 qsel 3
#> logLik -108.307 eps 0.0133 iteration  106 qsel 3
#> logLik -107.859 eps 0.0328 iteration  107 qsel 3
#> logLik -107.357 eps 0.0154 iteration  108 qsel 3
#> logLik -106.922 eps 0.0129 iteration  109 qsel 3
#> logLik -106.484 eps 0.0150 iteration  110 qsel 3
#> logLik -105.991 eps 0.0142 iteration  111 qsel 3
#> logLik -105.565 eps 0.0105 iteration  112 qsel 3
#> logLik -105.136 eps 0.0149 iteration  113 qsel 3
#> logLik -104.718 eps 0.0104 iteration  114 qsel 3
#> logLik -104.236 eps 0.0138 iteration  115 qsel 3
#> logLik -103.815 eps 0.0203 iteration  116 qsel 3
#> logLik -103.405 eps 0.0108 iteration  117 qsel 3
#> logLik -102.932 eps 0.0153 iteration  118 qsel 3
#> logLik -102.520 eps 0.0313 iteration  119 qsel 3
#> logLik -102.118 eps 0.0125 iteration  120 qsel 3
#> logLik -101.655 eps 0.0174 iteration  121 qsel 3
#> logLik -101.251 eps 0.0179 iteration  122 qsel 3
#> logLik -100.857 eps 0.0124 iteration  123 qsel 3
#> logLik -100.403 eps 0.0178 iteration  124 qsel 3
#> logLik -100.007 eps 0.0192 iteration  125 qsel 3
#> logLik -99.6211 eps 0.0181 iteration  126 qsel 3
#> logLik -99.1757 eps 0.0225 iteration  127 qsel 3
#> logLik -98.7880 eps 0.0260 iteration  128 qsel 3
#> logLik -98.4094 eps 0.0175 iteration  129 qsel 3
#> logLik -97.9728 eps 0.0268 iteration  130 qsel 3
#> logLik -97.5930 eps 0.0186 iteration  131 qsel 3
#> logLik -97.2218 eps 0.0200 iteration  132 qsel 3
#> logLik -96.7937 eps 0.0226 iteration  133 qsel 3
#> logLik -96.4216 eps 0.0201 iteration  134 qsel 3
#> logLik -96.0555 eps 0.0010 iteration  135 qsel 3
#> logLik -95.6900 eps 0.0157 iteration  136 qsel 3
#> logLik -95.2686 eps 0.0221 iteration  137 qsel 3
#> logLik -94.9016 eps 0.0010 iteration  138 qsel 3
#> logLik -94.5413 eps 0.0009 iteration  139 qsel 3
#> logLik -94.1801 eps 0.1458 iteration  140 qsel 3
#> logLik -93.8106 eps 0.0606 iteration  141 qsel 3
#> logLik -93.3939 eps 0.0201 iteration  142 qsel 3
#> logLik -93.0292 eps 0.0009 iteration  143 qsel 3
#> logLik -92.6712 eps 0.0009 iteration  144 qsel 3
#> logLik -92.3140 eps 0.0171 iteration  145 qsel 3
#> logLik -91.9490 eps 0.0230 iteration  146 qsel 3
#> logLik -91.5369 eps 0.0143 iteration  147 qsel 3
#> logLik -91.1747 eps 0.0009 iteration  148 qsel 3
#> logLik -90.8189 eps 0.0009 iteration  149 qsel 3
#> logLik -90.4658 eps 0.0106 iteration  150 qsel 3
#> logLik -90.1051 eps 0.0193 iteration  151 qsel 3
#> logLik -89.6977 eps 0.0147 iteration  152 qsel 3
#> logLik -89.3380 eps 0.0009 iteration  153 qsel 3
#> logLik -88.9846 eps 0.0009 iteration  154 qsel 3
#> logLik -88.6355 eps 0.0100 iteration  155 qsel 3
#> logLik -88.2792 eps 0.0207 iteration  156 qsel 3
#> logLik -87.8763 eps 0.0395 iteration  157 qsel 3
#> logLik -87.5191 eps 0.0009 iteration  158 qsel 3
#> logLik -87.1683 eps 0.0009 iteration  159 qsel 3
#> logLik -86.8230 eps 0.0108 iteration  160 qsel 3
#> logLik -86.4710 eps 0.0390 iteration  161 qsel 3
#> logLik -86.1201 eps 0.0009 iteration  162 qsel 3
#> logLik -85.7755 eps 0.0009 iteration  163 qsel 3
#> logLik -85.4340 eps 0.0109 iteration  164 qsel 3
#> logLik -85.0326 eps 0.0579 iteration  165 qsel 3
#> logLik -84.6849 eps 0.0224 iteration  166 qsel 3
#> logLik -84.3365 eps 0.0009 iteration  167 qsel 3
#> logLik -83.9944 eps 0.0009 iteration  168 qsel 3
#> logLik -83.6568 eps 0.0291 iteration  169 qsel 3
#> logLik -83.2599 eps 0.0290 iteration  170 qsel 3
#> logLik -82.9164 eps 0.0178 iteration  171 qsel 3
#> logLik -82.5707 eps 0.0009 iteration  172 qsel 3
#> logLik -82.2313 eps 0.0009 iteration  173 qsel 3
#> logLik -81.8975 eps 0.0166 iteration  174 qsel 3
#> logLik -81.5578 eps 0.0302 iteration  175 qsel 3
#> logLik -81.1657 eps 0.0219 iteration  176 qsel 3
#> logLik -80.8227 eps 0.0009 iteration  177 qsel 3
#> logLik -80.4859 eps 0.0009 iteration  178 qsel 3
#> logLik -80.1552 eps 0.0009 iteration  179 qsel 3
#> logLik -79.8239 eps 0.0117 iteration  180 qsel 3
#> logLik -79.4871 eps 0.0190 iteration  181 qsel 3
#> logLik -79.0978 eps 0.0131 iteration  182 qsel 3
#> logLik -78.7637 eps 0.0009 iteration  183 qsel 3
#> logLik -78.4355 eps 0.0009 iteration  184 qsel 3
#> logLik -78.1081 eps 0.0123 iteration  185 qsel 3
#> logLik -77.7754 eps 0.0122 iteration  186 qsel 3
#> logLik -77.3906 eps 0.0181 iteration  187 qsel 3
#> logLik -77.0591 eps 0.0009 iteration  188 qsel 3
#> logLik -76.7337 eps 0.0009 iteration  189 qsel 3
#> logLik -76.4100 eps 0.0117 iteration  190 qsel 3
#> logLik -76.0814 eps 0.0170 iteration  191 qsel 3
#> logLik -75.7010 eps 0.0115 iteration  192 qsel 3
#> logLik -75.3723 eps 0.0009 iteration  193 qsel 3
#> logLik -75.0495 eps 0.0008 iteration  194 qsel 3
#> logLik -74.7296 eps 0.0086 iteration  195 qsel 3
#> logLik -74.4051 eps 0.0496 iteration  196 qsel 3
#> logLik -74.0290 eps 0.0111 iteration  197 qsel 3
#> logLik -73.7032 eps 0.0008 iteration  198 qsel 3
#> logLik -73.3832 eps 0.0008 iteration  199 qsel 3
#> logLik -73.0670 eps 0.0094 iteration  200 qsel 3
#> logLik -72.7465 eps 0.0111 iteration  201 qsel 3
#> logLik -72.3747 eps 0.0104 iteration  202 qsel 3
#> logLik -72.0517 eps 0.0008 iteration  203 qsel 3
#> logLik -71.7346 eps 0.0008 iteration  204 qsel 3
#> logLik -71.4221 eps 0.0085 iteration  205 qsel 3
#> logLik -71.1056 eps 0.0105 iteration  206 qsel 3
#> logLik -70.7380 eps 0.0126 iteration  207 qsel 3
#> logLik -70.4180 eps 0.0008 iteration  208 qsel 3
#> logLik -70.1038 eps 0.0008 iteration  209 qsel 3
#> logLik -69.7949 eps 0.0083 iteration  210 qsel 3
#> logLik -69.4823 eps 0.0111 iteration  211 qsel 3
#> logLik -69.1189 eps 0.0116 iteration  212 qsel 3
#> logLik -68.8019 eps 0.0008 iteration  213 qsel 3
#> logLik -68.4906 eps 0.0008 iteration  214 qsel 3
#> logLik -68.1850 eps 0.0008 iteration  215 qsel 3
#> logLik -67.8785 eps 0.0085 iteration  216 qsel 3
#> logLik -67.5687 eps 0.0174 iteration  217 qsel 3
#> logLik -67.2081 eps 0.0120 iteration  218 qsel 3
#> logLik -66.8997 eps 0.0008 iteration  219 qsel 3
#> logLik -66.5969 eps 0.0008 iteration  220 qsel 3
#> logLik -66.2941 eps 0.0095 iteration  221 qsel 3
#> logLik -65.9882 eps 0.0122 iteration  222 qsel 3
#> logLik -65.6858 eps 0.0008 iteration  223 qsel 3
#> logLik -65.3875 eps 0.0105 iteration  224 qsel 3
#> logLik -65.0298 eps 0.0213 iteration  225 qsel 3
#> logLik -64.7273 eps 0.0008 iteration  226 qsel 3
#> logLik -64.4253 eps 0.0134 iteration  227 qsel 3
#> logLik -64.1259 eps 0.0008 iteration  228 qsel 3
#> logLik -63.8312 eps 0.0119 iteration  229 qsel 3
#> logLik -63.4775 eps 0.0175 iteration  230 qsel 3
#> logLik -63.1780 eps 0.0008 iteration  231 qsel 3
#> logLik -62.8799 eps 0.0156 iteration  232 qsel 3
#> logLik -62.5834 eps 0.0008 iteration  233 qsel 3
#> logLik -62.2922 eps 0.0103 iteration  234 qsel 3
#> logLik -61.9426 eps 0.0730 iteration  235 qsel 3
#> logLik -61.6462 eps 0.0008 iteration  236 qsel 3
#> logLik -61.3518 eps 0.0231 iteration  237 qsel 3
#> logLik -61.0583 eps 0.0008 iteration  238 qsel 3
#> logLik -60.7702 eps 0.0008 iteration  239 qsel 3
#> logLik -60.4814 eps 0.0224 iteration  240 qsel 3
#> logLik -60.1345 eps 0.0128 iteration  241 qsel 3
#> logLik -59.8439 eps 0.0188 iteration  242 qsel 3
#> logLik -59.5535 eps 0.0008 iteration  243 qsel 3
#> logLik -59.2683 eps 0.0008 iteration  244 qsel 3
#> logLik -58.9829 eps 0.0102 iteration  245 qsel 3
#> logLik -58.6401 eps 0.0500 iteration  246 qsel 3
#> logLik -58.3532 eps 0.0169 iteration  247 qsel 3
#> logLik -58.0658 eps 0.0008 iteration  248 qsel 3
#> logLik -57.7836 eps 0.0008 iteration  249 qsel 3
#> logLik -57.5017 eps 0.0856 iteration  250 qsel 3
#> logLik -57.1628 eps 0.0155 iteration  251 qsel 3
#> logLik -56.8797 eps 0.0157 iteration  252 qsel 3
#> logLik -56.5953 eps 0.0008 iteration  253 qsel 3
#> logLik -56.3161 eps 0.0007 iteration  254 qsel 3
#> logLik -56.0376 eps 0.0120 iteration  255 qsel 3
#> logLik -55.7026 eps 0.0108 iteration  256 qsel 3
#> logLik -55.4231 eps 0.0112 iteration  257 qsel 3
#> logLik -55.1419 eps 0.0007 iteration  258 qsel 3
#> logLik -54.8657 eps 0.0007 iteration  259 qsel 3
#> logLik -54.5906 eps 0.0091 iteration  260 qsel 3
#> logLik -54.2595 eps 0.0111 iteration  261 qsel 3
#> logLik -53.9836 eps 0.0007 iteration  262 qsel 3
#> logLik -53.7067 eps 0.0552 iteration  263 qsel 3
#> logLik -53.4336 eps 0.0007 iteration  264 qsel 3
#> logLik -53.1618 eps 0.0101 iteration  265 qsel 3
#> logLik -52.8346 eps 0.0177 iteration  266 qsel 3
#> logLik -52.5617 eps 0.0007 iteration  267 qsel 3
#> logLik -52.2885 eps 0.0127 iteration  268 qsel 3
#> logLik -52.0184 eps 0.0007 iteration  269 qsel 3
#> logLik -51.7500 eps 0.0122 iteration  270 qsel 3
#> logLik -51.4265 eps 0.0129 iteration  271 qsel 3
#> logLik -51.1568 eps 0.0007 iteration  272 qsel 3
#> logLik -50.8920 eps 0.0007 iteration  273 qsel 3
#> logLik -50.6270 eps 0.0127 iteration  274 qsel 3
#> logLik -50.3561 eps 0.0161 iteration  275 qsel 3
#> logLik -50.0364 eps 0.0292 iteration  276 qsel 3
#> logLik -49.7698 eps 0.0007 iteration  277 qsel 3
#> logLik -49.5081 eps 0.0007 iteration  278 qsel 3
#> logLik -49.2463 eps 0.0142 iteration  279 qsel 3
#> logLik -48.9790 eps 0.0392 iteration  280 qsel 3
#> logLik -48.6630 eps 0.0128 iteration  281 qsel 3
#> logLik -48.3996 eps 0.0007 iteration  282 qsel 3
#> logLik -48.1409 eps 0.0007 iteration  283 qsel 3
#> logLik -47.8823 eps 0.0350 iteration  284 qsel 3
#> logLik -47.6186 eps 0.0144 iteration  285 qsel 3
#> logLik -47.3606 eps 0.0007 iteration  286 qsel 3
#> logLik -47.1061 eps 0.0218 iteration  287 qsel 3
#> logLik -46.7928 eps 0.0197 iteration  288 qsel 3
#> logLik -46.5334 eps 0.0152 iteration  289 qsel 3
#> logLik -46.2739 eps 0.0007 iteration  290 qsel 3
#> logLik -46.0191 eps 0.0007 iteration  291 qsel 3
#> logLik -45.7677 eps 0.0369 iteration  292 qsel 3
#> logLik -45.5114 eps 0.0158 iteration  293 qsel 3
#> logLik -45.2020 eps 0.0174 iteration  294 qsel 3
#> logLik -44.9456 eps 0.0007 iteration  295 qsel 3
#> logLik -44.6940 eps 0.0007 iteration  296 qsel 3
#> logLik -44.4457 eps 0.0123 iteration  297 qsel 3
#> logLik -44.1929 eps 0.0832 iteration  298 qsel 3
#> logLik -43.8870 eps 0.0248 iteration  299 qsel 3
#> logLik -43.6339 eps 0.0007 iteration  300 qsel 3
#> logLik -43.3854 eps 0.0007 iteration  301 qsel 3
#> logLik -43.1402 eps 0.0129 iteration  302 qsel 3
#> logLik -42.8908 eps 0.0230 iteration  303 qsel 3
#> logLik -42.5884 eps 0.0128 iteration  304 qsel 3
#> logLik -42.3386 eps 0.0007 iteration  305 qsel 3
#> logLik -42.0933 eps 0.0007 iteration  306 qsel 3
#> logLik -41.8511 eps 0.0142 iteration  307 qsel 3
#> logLik -41.6050 eps 0.0155 iteration  308 qsel 3
#> logLik -41.3062 eps 0.0118 iteration  309 qsel 3
#> logLik -41.0595 eps 0.0007 iteration  310 qsel 3
#> logLik -40.8174 eps 0.0007 iteration  311 qsel 3
#> logLik -40.5782 eps 0.0127 iteration  312 qsel 3
#> logLik -40.3354 eps 0.0205 iteration  313 qsel 3
#> logLik -40.0401 eps 0.0106 iteration  314 qsel 3
#> logLik -39.7967 eps 0.0007 iteration  315 qsel 3
#> logLik -39.5578 eps 0.0007 iteration  316 qsel 3
#> logLik -39.3216 eps 0.0085 iteration  317 qsel 3
#> logLik -39.0820 eps 0.0253 iteration  318 qsel 3
#> logLik -38.7901 eps 0.0148 iteration  319 qsel 3
#> logLik -38.5500 eps 0.0007 iteration  320 qsel 3
#> logLik -38.3142 eps 0.0006 iteration  321 qsel 3
#> logLik -38.0810 eps 0.0080 iteration  322 qsel 3
#> logLik -37.8447 eps 0.0137 iteration  323 qsel 3
#> logLik -37.5562 eps 0.0167 iteration  324 qsel 3
#> logLik -37.3193 eps 0.0006 iteration  325 qsel 3
#> logLik -37.0867 eps 0.0006 iteration  326 qsel 3
#> logLik -36.8565 eps 0.0079 iteration  327 qsel 3
#> logLik -36.6233 eps 0.0126 iteration  328 qsel 3
#> logLik -36.3382 eps 0.0156 iteration  329 qsel 3
#> logLik -36.1045 eps 0.0006 iteration  330 qsel 3
#> logLik -35.8751 eps 0.0006 iteration  331 qsel 3
#> logLik -35.6478 eps 0.0175 iteration  332 qsel 3
#> logLik -35.4178 eps 0.0153 iteration  333 qsel 3
#> logLik -35.1359 eps 0.0124 iteration  334 qsel 3
#> logLik -34.9055 eps 0.0006 iteration  335 qsel 3
#> logLik -34.6793 eps 0.0006 iteration  336 qsel 3
#> logLik -34.4549 eps 0.0073 iteration  337 qsel 3
#> logLik -34.2280 eps 0.0440 iteration  338 qsel 3
#> logLik -33.9495 eps 0.0136 iteration  339 qsel 3
#> logLik -33.7223 eps 0.0006 iteration  340 qsel 3
#> logLik -33.4993 eps 0.0006 iteration  341 qsel 3
#> logLik -33.2776 eps 0.0070 iteration  342 qsel 3
#> logLik -33.0539 eps 0.0127 iteration  343 qsel 3
#> logLik -32.8318 eps 0.0006 iteration  344 qsel 3
#> logLik -32.6137 eps 0.3435 iteration  345 qsel 3
#> logLik -32.3375 eps 0.0116 iteration  346 qsel 3
#> logLik -32.1162 eps 0.0006 iteration  347 qsel 3
#> logLik -31.8955 eps 0.0127 iteration  348 qsel 3
#> logLik -31.6766 eps 0.0006 iteration  349 qsel 3
#> logLik -31.4613 eps 0.0077 iteration  350 qsel 3
#> logLik -31.1883 eps 0.0157 iteration  351 qsel 3
#> logLik -30.9701 eps 0.0006 iteration  352 qsel 3
#> logLik -30.7525 eps 0.0142 iteration  353 qsel 3
#> logLik -30.5368 eps 0.0006 iteration  354 qsel 3
#> logLik -30.3242 eps 0.0085 iteration  355 qsel 3
#> logLik -30.0544 eps 0.0367 iteration  356 qsel 3
#> logLik -29.8395 eps 0.0006 iteration  357 qsel 3
#> logLik -29.6248 eps 0.0114 iteration  358 qsel 3
#> logLik -29.4123 eps 0.0006 iteration  359 qsel 3
#> logLik -29.2024 eps 0.0174 iteration  360 qsel 3
#> logLik -28.9358 eps 0.0100 iteration  361 qsel 3
#> logLik -28.7240 eps 0.0006 iteration  362 qsel 3
#> logLik -28.5122 eps 0.0415 iteration  363 qsel 3
#> logLik -28.3029 eps 0.0006 iteration  364 qsel 3
#> logLik -28.0957 eps 0.0084 iteration  365 qsel 3
#> logLik -27.8322 eps 0.0087 iteration  366 qsel 3
#> logLik -27.6236 eps 0.0006 iteration  367 qsel 3
#> logLik -27.4148 eps 0.6297 iteration  368 qsel 3
#> logLik -27.2086 eps 0.0006 iteration  369 qsel 3
#> logLik -27.0041 eps 0.0098 iteration  370 qsel 3
#> logLik -26.7437 eps 0.0097 iteration  371 qsel 3
#> logLik -26.5382 eps 0.0006 iteration  372 qsel 3
#> logLik -26.3323 eps 0.0099 iteration  373 qsel 3
#> logLik -26.1292 eps 0.0006 iteration  374 qsel 3
#> logLik -25.9274 eps 0.0083 iteration  375 qsel 3
#> logLik -25.6700 eps 0.0142 iteration  376 qsel 3
#> logLik -25.4676 eps 0.0231 iteration  377 qsel 3
#> logLik -25.2639 eps 0.0006 iteration  378 qsel 3
#> logLik -25.0639 eps 0.0006 iteration  379 qsel 3
#> logLik -24.8647 eps 0.0070 iteration  380 qsel 3
#> logLik -24.6103 eps 0.0116 iteration  381 qsel 3
#> logLik -24.4107 eps 0.0073 iteration  382 qsel 3
#> logLik -24.2102 eps 0.0006 iteration  383 qsel 3
#> logLik -24.0134 eps 0.0006 iteration  384 qsel 3
#> logLik -23.8167 eps 0.0071 iteration  385 qsel 3
#> logLik -23.5653 eps 0.0084 iteration  386 qsel 3
#> logLik -23.3685 eps 0.0070 iteration  387 qsel 3
#> logLik -23.1711 eps 0.0006 iteration  388 qsel 3
#> logLik -22.9774 eps 0.0006 iteration  389 qsel 3
#> logLik -22.7833 eps 0.0058 iteration  390 qsel 3
#> logLik -22.5349 eps 0.0130 iteration  391 qsel 3
#> logLik -22.3408 eps 0.0077 iteration  392 qsel 3
#> logLik -22.1465 eps 0.0006 iteration  393 qsel 3
#> logLik -21.9556 eps 0.0064 iteration  394 qsel 3
#> logLik -21.7107 eps 0.0076 iteration  395 qsel 3
#> logLik -21.5173 eps 0.0006 iteration  396 qsel 3
#> logLik -21.3259 eps 0.0060 iteration  397 qsel 3
#> logLik -21.1348 eps 0.0005 iteration  398 qsel 3
#> logLik -20.9464 eps 0.0066 iteration  399 qsel 3
#> logLik -20.7043 eps 0.0064 iteration  400 qsel 3
#> logLik -20.5140 eps 0.0005 iteration  401 qsel 3
#> logLik -20.3272 eps 0.0005 iteration  402 qsel 3
#> logLik -20.1414 eps 0.0058 iteration  403 qsel 3
#> logLik -19.9520 eps 0.0059 iteration  404 qsel 3
#> logLik -19.7128 eps 0.0062 iteration  405 qsel 3
#> logLik -19.5256 eps 0.0005 iteration  406 qsel 3
#> logLik -19.3419 eps 0.0005 iteration  407 qsel 3
#> logLik -19.1585 eps 0.0060 iteration  408 qsel 3
#> logLik -18.9718 eps 0.0058 iteration  409 qsel 3
#> logLik -18.7892 eps 0.0005 iteration  410 qsel 3
#> logLik -18.6087 eps 0.0056 iteration  411 qsel 3
#> logLik -18.3716 eps 0.0082 iteration  412 qsel 3
#> logLik -18.1881 eps 0.0066 iteration  413 qsel 3
#> logLik -18.0051 eps 0.0005 iteration  414 qsel 3
#> logLik -17.8255 eps 0.0005 iteration  415 qsel 3
#> logLik -17.6474 eps 0.0080 iteration  416 qsel 3
#> logLik -17.4663 eps 0.0061 iteration  417 qsel 3
#> logLik -17.2321 eps 0.0075 iteration  418 qsel 3
#> logLik -17.0522 eps 0.0005 iteration  419 qsel 3
#> logLik -16.8756 eps 0.0005 iteration  420 qsel 3
#> logLik -16.6999 eps 0.0047 iteration  421 qsel 3
#> logLik -16.5214 eps 0.0084 iteration  422 qsel 3
#> logLik -16.2898 eps 0.0054 iteration  423 qsel 3
#> logLik -16.1130 eps 0.0005 iteration  424 qsel 3
#> logLik -15.9394 eps 0.0005 iteration  425 qsel 3
#> logLik -15.7660 eps 0.0043 iteration  426 qsel 3
#> logLik -15.5900 eps 0.0053 iteration  427 qsel 3
#> logLik -15.3612 eps 0.0049 iteration  428 qsel 3
#> logLik -14.7152 eps 0.0089 iteration  429 qsel 4
#> logLik -14.5374 eps 0.0005 iteration  430 qsel 4
#> logLik -14.3629 eps 0.0005 iteration  431 qsel 4
#> logLik -14.1915 eps 0.0045 iteration  432 qsel 4
#> logLik -14.0179 eps 0.0052 iteration  433 qsel 4
#> logLik -13.7917 eps 0.0048 iteration  434 qsel 4
#> logLik -13.1438 eps 0.0108 iteration  435 qsel 4
#> logLik -12.9651 eps 0.0005 iteration  436 qsel 4
#> logLik -12.7897 eps 0.0005 iteration  437 qsel 4
#> logLik -12.6176 eps 0.0005 iteration  438 qsel 4
#> logLik -12.4476 eps 0.0044 iteration  439 qsel 4
#> logLik -12.2759 eps 0.0062 iteration  440 qsel 4
#> logLik -11.6245 eps 0.0088 iteration  441 qsel 4
#> logLik -11.4003 eps 0.0050 iteration  442 qsel 4
#> logLik -11.2240 eps 0.0005 iteration  443 qsel 4
#> logLik -11.0509 eps 0.0005 iteration  444 qsel 4
#> logLik -10.8811 eps 0.0005 iteration  445 qsel 4
#> logLik -10.7126 eps 0.0042 iteration  446 qsel 4
#> logLik -10.5428 eps 0.0049 iteration  447 qsel 4
#> logLik  -9.8880 eps 0.0096 iteration  448 qsel 4
#> logLik  -9.6657 eps 0.0053 iteration  449 qsel 4
#> logLik  -9.4918 eps 0.0005 iteration  450 qsel 4
#> logLik  -9.3210 eps 0.0005 iteration  451 qsel 4
#> logLik  -9.1535 eps 0.0005 iteration  452 qsel 4
#> logLik  -8.9864 eps 0.0045 iteration  453 qsel 4
#> logLik  -8.8185 eps 0.0049 iteration  454 qsel 4
#> logLik  -8.1606 eps 0.0172 iteration  455 qsel 4
#> logLik  -7.9403 eps 0.0096 iteration  456 qsel 4
#> logLik  -7.7686 eps 0.0005 iteration  457 qsel 4
#> logLik  -7.6001 eps 0.0005 iteration  458 qsel 4
#> logLik  -7.4348 eps 0.0005 iteration  459 qsel 4
#> logLik  -7.2691 eps 0.0075 iteration  460 qsel 4
#> logLik  -7.1030 eps 0.0057 iteration  461 qsel 4
#> logLik  -6.4424 eps 0.0196 iteration  462 qsel 4
#> logLik  -6.2239 eps 0.0078 iteration  463 qsel 4
#> logLik  -6.0545 eps 0.0005 iteration  464 qsel 4
#> logLik  -5.8882 eps 0.0005 iteration  465 qsel 4
#> logLik  -5.7245 eps 0.0178 iteration  466 qsel 4
#> logLik  -5.5603 eps 0.0005 iteration  467 qsel 4
#> logLik  -5.3960 eps 0.0087 iteration  468 qsel 4
#> logLik  -4.7329 eps 0.0110 iteration  469 qsel 4
#> logLik  -4.5669 eps 0.0005 iteration  470 qsel 4
#> logLik  -4.4040 eps 0.0005 iteration  471 qsel 4
#> logLik  -4.2417 eps 0.0050 iteration  472 qsel 4
#> logLik  -4.0240 eps 0.0063 iteration  473 qsel 4
#> logLik  -3.8619 eps 0.0147 iteration  474 qsel 4
#> logLik  -3.1984 eps 0.0097 iteration  475 qsel 4
#> logLik  -3.0315 eps 0.0005 iteration  476 qsel 4
#> logLik  -2.8677 eps 0.0005 iteration  477 qsel 4
#> logLik  -2.7070 eps 0.0005 iteration  478 qsel 4
#> logLik  -2.5460 eps 0.0061 iteration  479 qsel 4
#> logLik  -2.3302 eps 0.0069 iteration  480 qsel 4
#> logLik  -2.1702 eps 0.0005 iteration  481 qsel 4
#> logLik  -2.0120 eps 0.0061 iteration  482 qsel 4
#> logLik  -1.8511 eps 0.0061 iteration  483 qsel 4
#> logLik  -1.1836 eps 0.0133 iteration  484 qsel 4
#> logLik  -0.9709 eps 0.0070 iteration  485 qsel 4
#> logLik  -0.8070 eps 0.0005 iteration  486 qsel 4
#> logLik  -0.6461 eps 0.0005 iteration  487 qsel 4
#> logLik  -0.4883 eps 0.0005 iteration  488 qsel 4
#> logLik  -0.3313 eps 0.0110 iteration  489 qsel 4
#> logLik   0.3372 eps 0.0201 iteration  490 qsel 4
#> logLik   0.4969 eps 0.0005 iteration  491 qsel 4
#> logLik   0.6566 eps 0.0112 iteration  492 qsel 4
#> logLik   0.8680 eps 0.0115 iteration  493 qsel 4
#> logLik   1.0269 eps 0.0005 iteration  494 qsel 4
#> logLik   1.1827 eps 0.0005 iteration  495 qsel 4
#> logLik   1.3383 eps 0.0073 iteration  496 qsel 4
#> logLik   2.0080 eps 0.0087 iteration  497 qsel 4
#> logLik   2.1657 eps 0.0005 iteration  498 qsel 4
#> logLik   2.3236 eps 0.0050 iteration  499 qsel 4
#> logLik   2.5332 eps 0.0050 iteration  500 qsel 4
#> logLik   2.6900 eps 0.0005 iteration  501 qsel 4
#> logLik   2.8439 eps 0.0005 iteration  502 qsel 4
#> logLik   2.9981 eps 0.0043 iteration  503 qsel 4
#> logLik   3.6687 eps 0.0082 iteration  504 qsel 4
#> logLik   3.8245 eps 0.0046 iteration  505 qsel 4
#> logLik   4.0317 eps 0.0045 iteration  506 qsel 4
#> logLik   4.1895 eps 0.0005 iteration  507 qsel 4
#> logLik   4.3443 eps 0.0005 iteration  508 qsel 4
#> logLik   4.4968 eps 0.0038 iteration  509 qsel 4
#> logLik   5.1661 eps 0.0099 iteration  510 qsel 4
#> logLik   5.3227 eps 0.0005 iteration  511 qsel 4
#> logLik   5.4768 eps 0.0073 iteration  512 qsel 4
#> logLik   5.6821 eps 0.0050 iteration  513 qsel 4
#> logLik   5.8379 eps 0.0005 iteration  514 qsel 4
#> logLik   5.9907 eps 0.0005 iteration  515 qsel 4
#> logLik   6.1419 eps 0.0041 iteration  516 qsel 4
#> logLik   6.8113 eps 0.0119 iteration  517 qsel 4
#> logLik   6.9659 eps 0.0005 iteration  518 qsel 4
#> logLik   7.1183 eps 0.0045 iteration  519 qsel 4
#> logLik   7.2709 eps 0.0005 iteration  520 qsel 4
#> logLik   7.4206 eps 0.0004 iteration  521 qsel 4
#> logLik   7.5705 eps 0.0039 iteration  522 qsel 4
#> logLik   7.7752 eps 0.0071 iteration  523 qsel 4
#> logLik   8.4444 eps 0.0075 iteration  524 qsel 4
#> logLik   8.5970 eps 0.0005 iteration  525 qsel 4
#> logLik   8.7478 eps 0.0039 iteration  526 qsel 4
#> logLik   8.8984 eps 0.0004 iteration  527 qsel 4
#> logLik   9.0467 eps 0.0033 iteration  528 qsel 4
#> logLik   9.2490 eps 0.0041 iteration  529 qsel 4
#> logLik   9.9159 eps 0.0069 iteration  530 qsel 4
#> logLik  10.0694 eps 0.0005 iteration  531 qsel 4
#> logLik  10.2200 eps 0.0004 iteration  532 qsel 4
#> logLik  10.3691 eps 0.0037 iteration  533 qsel 4
#> logLik  10.5178 eps 0.0004 iteration  534 qsel 4
#> logLik  10.6648 eps 0.0032 iteration  535 qsel 4
#> logLik  10.8654 eps 0.0038 iteration  536 qsel 4
#> logLik  11.5314 eps 0.0067 iteration  537 qsel 4
#> logLik  11.6829 eps 0.0004 iteration  538 qsel 4
#> logLik  11.8316 eps 0.0004 iteration  539 qsel 4
#> logLik  11.9791 eps 0.0036 iteration  540 qsel 4
#> logLik  12.1258 eps 0.0004 iteration  541 qsel 4
#> logLik  12.2715 eps 0.0031 iteration  542 qsel 4
#> logLik  12.4703 eps 0.0038 iteration  543 qsel 4
#> logLik  13.1352 eps 0.0067 iteration  544 qsel 4
#> logLik  13.2847 eps 0.0004 iteration  545 qsel 4
#> logLik  13.4314 eps 0.0004 iteration  546 qsel 4
#> logLik  13.5755 eps 0.0030 iteration  547 qsel 4
#> logLik  13.7215 eps 0.0036 iteration  548 qsel 4
#> logLik  13.9180 eps 0.0037 iteration  549 qsel 4
#> logLik  14.0647 eps 0.0004 iteration  550 qsel 4
#> logLik  14.2086 eps 0.0004 iteration  551 qsel 4
#> logLik  14.3507 eps 0.0030 iteration  552 qsel 4
#> logLik  15.0160 eps 0.0068 iteration  553 qsel 4
#> logLik  15.1616 eps 0.0004 iteration  554 qsel 4
#> logLik  15.3060 eps 0.0036 iteration  555 qsel 4
#> logLik  15.5008 eps 0.0038 iteration  556 qsel 4
#> logLik  15.6456 eps 0.0004 iteration  557 qsel 4
#> logLik  15.7877 eps 0.0004 iteration  558 qsel 4
#> logLik  15.9285 eps 0.0030 iteration  559 qsel 4
#> logLik  16.5919 eps 0.0072 iteration  560 qsel 4
#> logLik  16.7356 eps 0.0004 iteration  561 qsel 4
#> logLik  16.8784 eps 0.0037 iteration  562 qsel 4
#> logLik  17.0716 eps 0.0039 iteration  563 qsel 4
#> logLik  17.2145 eps 0.0004 iteration  564 qsel 4
#> logLik  17.3546 eps 0.0004 iteration  565 qsel 4
#> logLik  17.4943 eps 0.0030 iteration  566 qsel 4
#> logLik  18.1555 eps 0.0079 iteration  567 qsel 4
#> logLik  18.2973 eps 0.0004 iteration  568 qsel 4
#> logLik  18.4386 eps 0.0039 iteration  569 qsel 4
#> logLik  18.5785 eps 0.0004 iteration  570 qsel 4
#> logLik  18.7167 eps 0.0031 iteration  571 qsel 4
#> logLik  18.9086 eps 0.0043 iteration  572 qsel 4
#> logLik  19.5656 eps 0.0099 iteration  573 qsel 4
#> logLik  19.7082 eps 0.0004 iteration  574 qsel 4
#> logLik  19.8482 eps 0.0004 iteration  575 qsel 4
#> logLik  19.9879 eps 0.0044 iteration  576 qsel 4
#> logLik  20.1260 eps 0.0004 iteration  577 qsel 4
#> logLik  20.2630 eps 0.0033 iteration  578 qsel 4
#> logLik  20.4532 eps 0.0051 iteration  579 qsel 4
#> logLik  21.1074 eps 0.0300 iteration  580 qsel 4
#> logLik  21.2482 eps 0.0004 iteration  581 qsel 4
#> logLik  21.3863 eps 0.0004 iteration  582 qsel 4
#> logLik  21.5245 eps 0.0068 iteration  583 qsel 4
#> logLik  21.6607 eps 0.0004 iteration  584 qsel 4
#> logLik  21.7965 eps 0.0038 iteration  585 qsel 4
#> logLik  21.9851 eps 0.0198 iteration  586 qsel 4
#> logLik  22.6363 eps 0.0139 iteration  587 qsel 4
#> logLik  22.7752 eps 0.0004 iteration  588 qsel 4
#> logLik  22.9115 eps 0.0101 iteration  589 qsel 4
#> logLik  23.0485 eps 0.0004 iteration  590 qsel 4
#> logLik  23.1829 eps 0.0004 iteration  591 qsel 4
#> logLik  23.3175 eps 0.0052 iteration  592 qsel 4
#> logLik  23.5044 eps 0.0089 iteration  593 qsel 4
#> logLik  24.1523 eps 0.0244 iteration  594 qsel 4
#> logLik  24.2893 eps 0.0004 iteration  595 qsel 4
#> logLik  24.4241 eps 0.0138 iteration  596 qsel 4
#> logLik  24.5593 eps 0.0004 iteration  597 qsel 4
#> logLik  24.6924 eps 0.0046 iteration  598 qsel 4
#> logLik  24.8773 eps 0.0139 iteration  599 qsel 4
#> logLik  25.5200 eps 0.0165 iteration  600 qsel 4
#> logLik  25.6577 eps 0.0004 iteration  601 qsel 4
#> logLik  25.7928 eps 0.0004 iteration  602 qsel 4
#> logLik  25.9261 eps 0.0052 iteration  603 qsel 4
#> logLik  26.0593 eps 0.0004 iteration  604 qsel 4
#> logLik  26.1913 eps 0.0039 iteration  605 qsel 4
#> logLik  26.3746 eps 0.0078 iteration  606 qsel 4
#> logLik  26.5070 eps 0.0004 iteration  607 qsel 4
#> logLik  26.6369 eps 0.0004 iteration  608 qsel 4
#> logLik  26.7669 eps 0.0039 iteration  609 qsel 4
#> logLik  26.8987 eps 0.0085 iteration  610 qsel 4
#> logLik  27.5415 eps 0.0103 iteration  611 qsel 4
#> logLik  27.7228 eps 0.0050 iteration  612 qsel 4
#> logLik  27.8558 eps 0.0004 iteration  613 qsel 4
#> logLik  27.9864 eps 0.0004 iteration  614 qsel 4
#> logLik  28.1150 eps 0.0035 iteration  615 qsel 4
#> logLik  28.2450 eps 0.0045 iteration  616 qsel 4
#> logLik  28.8823 eps 0.0090 iteration  617 qsel 4
#> logLik  29.0150 eps 0.0004 iteration  618 qsel 4
#> logLik  29.1452 eps 0.0004 iteration  619 qsel 4
#> logLik  29.2730 eps 0.0004 iteration  620 qsel 4
#> logLik  29.4005 eps 0.0036 iteration  621 qsel 4
#> logLik  29.5811 eps 0.0044 iteration  622 qsel 4
#> logLik  29.7097 eps 0.0043 iteration  623 qsel 4
#> logLik  30.3428 eps 0.0104 iteration  624 qsel 4
#> logLik  30.4736 eps 0.0004 iteration  625 qsel 4
#> logLik  30.6020 eps 0.0004 iteration  626 qsel 4
#> logLik  30.7281 eps 0.0044 iteration  627 qsel 4
#> logLik  30.8551 eps 0.0070 iteration  628 qsel 4
#> logLik  31.0337 eps 0.0048 iteration  629 qsel 4
#> logLik  31.6609 eps 0.0160 iteration  630 qsel 4
#> logLik  31.7923 eps 0.0004 iteration  631 qsel 4
#> logLik  31.9213 eps 0.0004 iteration  632 qsel 4
#> logLik  32.0478 eps 0.0004 iteration  633 qsel 4
#> logLik  32.1728 eps 0.0093 iteration  634 qsel 4
#> logLik  32.2984 eps 0.0078 iteration  635 qsel 4
#> logLik  32.4754 eps 0.0087 iteration  636 qsel 4
#> logLik  33.0982 eps 0.0162 iteration  637 qsel 4
#> logLik  33.2277 eps 0.0004 iteration  638 qsel 4
#> logLik  33.3547 eps 0.0004 iteration  639 qsel 4
#> logLik  33.4794 eps 0.0004 iteration  640 qsel 4
#> logLik  33.6033 eps 0.0048 iteration  641 qsel 4
#> logLik  33.7275 eps 0.4839 iteration  642 qsel 4
#> logLik  33.9030 eps 0.0076 iteration  643 qsel 4
#> logLik  34.5212 eps 0.0099 iteration  644 qsel 4
#> logLik  34.6487 eps 0.0004 iteration  645 qsel 4
#> logLik  34.7739 eps 0.0004 iteration  646 qsel 4
#> logLik  34.8967 eps 0.0004 iteration  647 qsel 4
#> logLik  35.0196 eps 0.0045 iteration  648 qsel 4
#> logLik  35.1424 eps 0.0065 iteration  649 qsel 4
#> logLik  35.7557 eps 0.0078 iteration  650 qsel 4
#> logLik  35.9297 eps 0.0083 iteration  651 qsel 4
#> logLik  36.0554 eps 0.0004 iteration  652 qsel 4
#> logLik  36.1787 eps 0.0004 iteration  653 qsel 4
#> logLik  36.3002 eps 0.0043 iteration  654 qsel 4
#> logLik  36.4218 eps 0.0004 iteration  655 qsel 4
#> logLik  36.5432 eps 0.0045 iteration  656 qsel 4
#> logLik  37.1517 eps 0.0074 iteration  657 qsel 4
#> logLik  37.2746 eps 0.0004 iteration  658 qsel 4
#> logLik  37.3952 eps 0.0004 iteration  659 qsel 4
#> logLik  37.5156 eps 0.0033 iteration  660 qsel 4
#> logLik  37.6889 eps 0.0091 iteration  661 qsel 4
#> logLik  37.8088 eps 0.0042 iteration  662 qsel 4
#> logLik  37.9291 eps 0.0004 iteration  663 qsel 4
#> logLik  38.0477 eps 0.0034 iteration  664 qsel 4
#> logLik  38.6513 eps 0.0089 iteration  665 qsel 4
#> logLik  38.8224 eps 0.0094 iteration  666 qsel 4
#> logLik  38.9449 eps 0.0004 iteration  667 qsel 4
#> logLik  39.0651 eps 0.0004 iteration  668 qsel 4
#> logLik  39.1836 eps 0.0043 iteration  669 qsel 4
#> logLik  39.3021 eps 0.0004 iteration  670 qsel 4
#> logLik  39.4196 eps 0.0105 iteration  671 qsel 4
#> logLik  40.0182 eps 0.0070 iteration  672 qsel 4
#> logLik  40.1878 eps 0.0048 iteration  673 qsel 4
#> logLik  40.3084 eps 0.0004 iteration  674 qsel 4
#> logLik  40.4267 eps 0.0004 iteration  675 qsel 4
#> logLik  40.5429 eps 0.0004 iteration  676 qsel 4
#> logLik  40.6592 eps 0.0051 iteration  677 qsel 4
#> logLik  40.7768 eps 0.0039 iteration  678 qsel 4
#> logLik  41.3702 eps 0.0066 iteration  679 qsel 4
#> logLik  41.5383 eps 0.0048 iteration  680 qsel 4
#> logLik  41.6571 eps 0.0004 iteration  681 qsel 4
#> logLik  41.7736 eps 0.0004 iteration  682 qsel 4
#> logLik  41.8886 eps 0.0047 iteration  683 qsel 4
#> logLik  42.0046 eps 0.0055 iteration  684 qsel 4
#> logLik  42.5915 eps 0.0089 iteration  685 qsel 4
#> logLik  42.7578 eps 0.0045 iteration  686 qsel 4
#> logLik  42.8769 eps 0.0004 iteration  687 qsel 4
#> logLik  42.9937 eps 0.0004 iteration  688 qsel 4
#> logLik  43.1084 eps 0.0004 iteration  689 qsel 4
#> logLik  43.2225 eps 0.0033 iteration  690 qsel 4
#> logLik  43.3371 eps 0.0041 iteration  691 qsel 4
#> logLik  43.9188 eps 0.0064 iteration  692 qsel 4
#> logLik  44.0837 eps 0.0040 iteration  693 qsel 4
#> logLik  44.2008 eps 0.0004 iteration  694 qsel 4
#> logLik  44.3158 eps 0.0004 iteration  695 qsel 4
#> logLik  44.4286 eps 0.0030 iteration  696 qsel 4
#> logLik  44.5419 eps 0.0004 iteration  697 qsel 4
#> logLik  44.6554 eps 0.0033 iteration  698 qsel 4
#> logLik  45.2317 eps 0.0065 iteration  699 qsel 4
#> logLik  45.3462 eps 0.0004 iteration  700 qsel 4
#> logLik  45.4585 eps 0.0003 iteration  701 qsel 4
#> logLik  45.5703 eps 0.0030 iteration  702 qsel 4
#> logLik  45.7345 eps 0.0039 iteration  703 qsel 4
#> logLik  45.8464 eps 0.0030 iteration  704 qsel 4
#> logLik  46.4161 eps 0.0080 iteration  705 qsel 4
#> logLik  46.5308 eps 0.0004 iteration  706 qsel 4
#> logLik  46.6434 eps 0.0003 iteration  707 qsel 4
#> logLik  46.7540 eps 0.0032 iteration  708 qsel 4
#> logLik  46.9165 eps 0.0048 iteration  709 qsel 4
#> logLik  47.0282 eps 0.0003 iteration  710 qsel 4
#> logLik  47.1389 eps 0.0031 iteration  711 qsel 4
#> logLik  47.2490 eps 0.0003 iteration  712 qsel 4
#> logLik  47.3581 eps 0.0031 iteration  713 qsel 4
#> logLik  47.9239 eps 0.0181 iteration  714 qsel 4
#> logLik  48.0846 eps 0.0058 iteration  715 qsel 4
#> logLik  48.1966 eps 0.0003 iteration  716 qsel 4
#> logLik  48.3065 eps 0.0003 iteration  717 qsel 4
#> logLik  48.4159 eps 0.0033 iteration  718 qsel 4
#> logLik  48.5242 eps 0.0003 iteration  719 qsel 4
#> logLik  48.6323 eps 0.0065 iteration  720 qsel 4
#> logLik  49.1927 eps 0.0276 iteration  721 qsel 4
#> logLik  49.3521 eps 0.0115 iteration  722 qsel 4
#> logLik  49.4622 eps 0.0003 iteration  723 qsel 4
#> logLik  49.5703 eps 0.0003 iteration  724 qsel 4
#> logLik  49.6784 eps 0.0078 iteration  725 qsel 4
#> logLik  49.7853 eps 0.0053 iteration  726 qsel 4
#> logLik  50.3391 eps 0.0088 iteration  727 qsel 4
#> logLik  50.4967 eps 0.0041 iteration  728 qsel 4
#> logLik  50.6070 eps 0.0003 iteration  729 qsel 4
#> logLik  50.7152 eps 0.0003 iteration  730 qsel 4
#> logLik  50.8220 eps 0.0029 iteration  731 qsel 4
#> logLik  50.9287 eps 0.0003 iteration  732 qsel 4
#> logLik  51.0346 eps 0.0033 iteration  733 qsel 4
#> logLik  51.5829 eps 0.0068 iteration  734 qsel 4
#> logLik  51.6906 eps 0.0003 iteration  735 qsel 4
#> logLik  51.7962 eps 0.0003 iteration  736 qsel 4
#> logLik  51.9018 eps 0.0028 iteration  737 qsel 4
#> logLik  52.0066 eps 0.0028 iteration  738 qsel 4
#> logLik  52.1636 eps 0.0039 iteration  739 qsel 4
#> logLik  52.7053 eps 0.0076 iteration  740 qsel 4
#> logLik  52.8131 eps 0.0003 iteration  741 qsel 4
#> logLik  52.9189 eps 0.0003 iteration  742 qsel 4
#> logLik  53.0227 eps 0.0003 iteration  743 qsel 4
#> logLik  53.1264 eps 0.0026 iteration  744 qsel 4
#> logLik  53.2309 eps 0.0029 iteration  745 qsel 4
#> logLik  53.3866 eps 0.0067 iteration  746 qsel 4
#> logLik  53.9228 eps 0.0068 iteration  747 qsel 4
#> logLik  54.0287 eps 0.0003 iteration  748 qsel 4
#> logLik  54.1326 eps 0.0003 iteration  749 qsel 4
#> logLik  54.2352 eps 0.0025 iteration  750 qsel 4
#> logLik  54.3383 eps 0.0031 iteration  751 qsel 4
#> logLik  54.4923 eps 0.0037 iteration  752 qsel 4
#> logLik  55.0220 eps 0.0066 iteration  753 qsel 4
#> logLik  55.1280 eps 0.0003 iteration  754 qsel 4
#> logLik  55.2320 eps 0.0003 iteration  755 qsel 4
#> logLik  55.3340 eps 0.0003 iteration  756 qsel 4
#> logLik  55.4357 eps 0.0024 iteration  757 qsel 4
#> logLik  55.5376 eps 0.0031 iteration  758 qsel 4
#> logLik  55.6904 eps 0.0036 iteration  759 qsel 4
#> logLik  55.7920 eps 0.0003 iteration  760 qsel 4
#> logLik  55.8921 eps 0.0023 iteration  761 qsel 4
#> logLik  56.4177 eps 0.0310 iteration  762 qsel 4
#> logLik  56.5203 eps 0.0003 iteration  763 qsel 4
#> logLik  56.6209 eps 0.0003 iteration  764 qsel 4
#> logLik  56.7217 eps 0.0032 iteration  765 qsel 4
#> logLik  56.8731 eps 0.0067 iteration  766 qsel 4
#> logLik  56.9730 eps 0.0003 iteration  767 qsel 4
#> logLik  57.0721 eps 0.0025 iteration  768 qsel 4
#> logLik  57.5924 eps 0.0059 iteration  769 qsel 4
#> logLik  57.6931 eps 0.0003 iteration  770 qsel 4
#> logLik  57.7925 eps 0.0026 iteration  771 qsel 4
#> logLik  57.9423 eps 0.0036 iteration  772 qsel 4
#> logLik  58.0422 eps 0.0003 iteration  773 qsel 4
#> logLik  58.1403 eps 0.0023 iteration  774 qsel 4
#> logLik  58.6540 eps 0.0052 iteration  775 qsel 4
#> logLik  58.7547 eps 0.0003 iteration  776 qsel 4
#> logLik  58.8536 eps 0.0003 iteration  777 qsel 4
#> logLik  58.9519 eps 0.0024 iteration  778 qsel 4
#> logLik  59.0492 eps 0.0003 iteration  779 qsel 4
#> logLik  59.1464 eps 0.0022 iteration  780 qsel 4
#> logLik  59.2953 eps 0.0032 iteration  781 qsel 4
#> logLik  59.8037 eps 0.0050 iteration  782 qsel 4
#> logLik  59.9025 eps 0.0003 iteration  783 qsel 4
#> logLik  59.9995 eps 0.0024 iteration  784 qsel 4
#> logLik  60.0969 eps 0.0003 iteration  785 qsel 4
#> logLik  60.1930 eps 0.0022 iteration  786 qsel 4
#> logLik  60.3403 eps 0.0031 iteration  787 qsel 4
#> logLik  60.8423 eps 0.0049 iteration  788 qsel 4
#> logLik  60.9411 eps 0.0003 iteration  789 qsel 4
#> logLik  61.0381 eps 0.0003 iteration  790 qsel 4
#> logLik  61.1340 eps 0.0024 iteration  791 qsel 4
#> logLik  61.2295 eps 0.0003 iteration  792 qsel 4
#> logLik  61.3247 eps 0.0022 iteration  793 qsel 4
#> logLik  61.4708 eps 0.0031 iteration  794 qsel 4
#> logLik  61.5656 eps 0.0003 iteration  795 qsel 4
#> logLik  61.6600 eps 0.0023 iteration  796 qsel 4
#> logLik  62.1579 eps 0.0048 iteration  797 qsel 4
#> logLik  62.2534 eps 0.0003 iteration  798 qsel 4
#> logLik  62.3475 eps 0.0022 iteration  799 qsel 4
#> logLik  62.4922 eps 0.0032 iteration  800 qsel 4
#> logLik  62.5869 eps 0.0003 iteration  801 qsel 4
#> logLik  62.6800 eps 0.0024 iteration  802 qsel 4
#> logLik  62.7733 eps 0.0003 iteration  803 qsel 4
#> logLik  62.8662 eps 0.0022 iteration  804 qsel 4
#> logLik  63.3590 eps 0.0049 iteration  805 qsel 4
#> logLik  63.5022 eps 0.0036 iteration  806 qsel 4
#> logLik  63.5968 eps 0.0003 iteration  807 qsel 4
#> logLik  63.6897 eps 0.0003 iteration  808 qsel 4
#> logLik  63.7814 eps 0.0023 iteration  809 qsel 4
#> logLik  63.8736 eps 0.0025 iteration  810 qsel 4
#> logLik  64.3602 eps 0.0052 iteration  811 qsel 4
#> logLik  64.5018 eps 0.0052 iteration  812 qsel 4
#> logLik  64.5964 eps 0.0003 iteration  813 qsel 4
#> logLik  64.6892 eps 0.0003 iteration  814 qsel 4
#> logLik  64.7802 eps 0.0003 iteration  815 qsel 4
#> logLik  64.8711 eps 0.0029 iteration  816 qsel 4
#> logLik  64.9622 eps 0.0028 iteration  817 qsel 4
#> logLik  65.4436 eps 0.0071 iteration  818 qsel 4
#> logLik  65.5356 eps 0.0003 iteration  819 qsel 4
#> logLik  65.6260 eps 0.0003 iteration  820 qsel 4
#> logLik  65.7159 eps 0.0041 iteration  821 qsel 4
#> logLik  65.8570 eps 0.0048 iteration  822 qsel 4
#> logLik  65.9468 eps 0.0043 iteration  823 qsel 4
#> logLik  66.4220 eps 0.0068 iteration  824 qsel 4
#> logLik  66.5140 eps 0.0003 iteration  825 qsel 4
#> logLik  66.6042 eps 0.0003 iteration  826 qsel 4
#> logLik  66.6931 eps 0.0029 iteration  827 qsel 4
#> logLik  66.7820 eps 0.0003 iteration  828 qsel 4
#> logLik  66.8709 eps 0.0035 iteration  829 qsel 4
#> logLik  67.0108 eps 0.0481 iteration  830 qsel 4
#> logLik  67.4809 eps 0.0090 iteration  831 qsel 4
#> logLik  67.5710 eps 0.0003 iteration  832 qsel 4
#> logLik  67.6594 eps 0.0003 iteration  833 qsel 4
#> logLik  67.7475 eps 0.0069 iteration  834 qsel 4
#> logLik  67.8352 eps 0.0031 iteration  835 qsel 4
#> logLik  67.9737 eps 0.0054 iteration  836 qsel 4
#> logLik  68.4377 eps 0.0167 iteration  837 qsel 4
#> logLik  68.5277 eps 0.0003 iteration  838 qsel 4
#> logLik  68.6159 eps 0.0003 iteration  839 qsel 4
#> logLik  68.7030 eps 0.0042 iteration  840 qsel 4
#> logLik  68.7899 eps 0.0003 iteration  841 qsel 4
#> logLik  68.8766 eps 0.0099 iteration  842 qsel 4
#> logLik  69.0140 eps 0.0082 iteration  843 qsel 4
#> logLik  69.4730 eps 0.0061 iteration  844 qsel 4
#> logLik  69.5611 eps 0.0003 iteration  845 qsel 4
#> logLik  69.6475 eps 0.0003 iteration  846 qsel 4
#> logLik  69.7338 eps 0.0030 iteration  847 qsel 4
#> logLik  69.8193 eps 0.0030 iteration  848 qsel 4
#> logLik  69.9553 eps 0.0037 iteration  849 qsel 4
#> logLik  70.0413 eps 0.0003 iteration  850 qsel 4
#> logLik  70.1262 eps 0.0028 iteration  851 qsel 4
#> logLik  70.5804 eps 0.0055 iteration  852 qsel 4
#> logLik  70.6670 eps 0.0003 iteration  853 qsel 4
#> logLik  70.7519 eps 0.0003 iteration  854 qsel 4
#> logLik  70.8365 eps 0.0028 iteration  855 qsel 4
#> logLik  70.9205 eps 0.0028 iteration  856 qsel 4
#> logLik  71.0553 eps 0.0033 iteration  857 qsel 4
#> logLik  71.5037 eps 0.0059 iteration  858 qsel 4
#> logLik  71.5900 eps 0.0003 iteration  859 qsel 4
#> logLik  71.6748 eps 0.0003 iteration  860 qsel 4
#> logLik  71.7582 eps 0.0032 iteration  861 qsel 4
#> logLik  71.8417 eps 0.0003 iteration  862 qsel 4
#> logLik  71.9249 eps 0.0029 iteration  863 qsel 4
#> logLik  72.0587 eps 0.0033 iteration  864 qsel 4
#> logLik  72.5021 eps 0.0099 iteration  865 qsel 4
#> logLik  72.5866 eps 0.0003 iteration  866 qsel 4
#> logLik  72.6695 eps 0.0003 iteration  867 qsel 4
#> logLik  72.7520 eps 0.0034 iteration  868 qsel 4
#> logLik  72.8342 eps 0.0047 iteration  869 qsel 4
#> logLik  72.9667 eps 0.0063 iteration  870 qsel 4
#> logLik  73.4044 eps 0.0052 iteration  871 qsel 4
#> logLik  73.4887 eps 0.0003 iteration  872 qsel 4
#> logLik  73.5714 eps 0.0003 iteration  873 qsel 4
#> logLik  73.6527 eps 0.0036 iteration  874 qsel 4
#> logLik  73.7342 eps 0.0003 iteration  875 qsel 4
#> logLik  73.8157 eps 0.0035 iteration  876 qsel 4
#> logLik  73.9470 eps 0.0030 iteration  877 qsel 4
#> logLik  74.3799 eps 0.0048 iteration  878 qsel 4
#> logLik  74.4624 eps 0.0003 iteration  879 qsel 4
#> logLik  74.5432 eps 0.0003 iteration  880 qsel 4
#> logLik  74.6238 eps 0.0027 iteration  881 qsel 4
#> logLik  74.7043 eps 0.0035 iteration  882 qsel 4
#> logLik  74.8343 eps 0.0029 iteration  883 qsel 4
#> logLik  75.2616 eps 0.0047 iteration  884 qsel 4
#> logLik  75.3438 eps 0.0003 iteration  885 qsel 4
#> logLik  75.4244 eps 0.0003 iteration  886 qsel 4
#> logLik  75.5041 eps 0.0026 iteration  887 qsel 4
#> logLik  75.5835 eps 0.0003 iteration  888 qsel 4
#> logLik  75.6630 eps 0.0060 iteration  889 qsel 4
#> logLik  75.7414 eps 0.0029 iteration  890 qsel 4
#> logLik  75.8704 eps 0.0028 iteration  891 qsel 4
#> logLik  76.2932 eps 0.0047 iteration  892 qsel 4
#> logLik  76.3738 eps 0.0003 iteration  893 qsel 4
#> logLik  76.4529 eps 0.0003 iteration  894 qsel 4
#> logLik  76.5314 eps 0.0068 iteration  895 qsel 4
#> logLik  76.6094 eps 0.0003 iteration  896 qsel 4
#> logLik  76.6870 eps 0.0039 iteration  897 qsel 4
#> logLik  76.8150 eps 0.0029 iteration  898 qsel 4
#> logLik  77.2331 eps 0.0051 iteration  899 qsel 4
#> logLik  77.3119 eps 0.0003 iteration  900 qsel 4
#> logLik  77.3893 eps 0.0036 iteration  901 qsel 4
#> logLik  77.4670 eps 0.0003 iteration  902 qsel 4
#> logLik  77.5438 eps 0.0047 iteration  903 qsel 4
#> logLik  77.6202 eps 0.0003 iteration  904 qsel 4
#> logLik  77.6965 eps 0.0030 iteration  905 qsel 4
#> logLik  78.1100 eps 0.0052 iteration  906 qsel 4
#> logLik  78.2369 eps 0.0028 iteration  907 qsel 4
#> logLik  78.3143 eps 0.0003 iteration  908 qsel 4
#> logLik  78.3902 eps 0.0030 iteration  909 qsel 4
#> logLik  78.4664 eps 0.0003 iteration  910 qsel 4
#> logLik  78.5416 eps 0.0034 iteration  911 qsel 4
#> logLik  78.9497 eps 0.0085 iteration  912 qsel 4
#> logLik  79.0754 eps 0.0029 iteration  913 qsel 4
#> logLik  79.1525 eps 0.0003 iteration  914 qsel 4
#> logLik  79.2281 eps 0.0003 iteration  915 qsel 4
#> logLik  79.3033 eps 0.0054 iteration  916 qsel 4
#> logLik  79.3778 eps 0.0003 iteration  917 qsel 4
#> logLik  79.4521 eps 0.0036 iteration  918 qsel 4
#> logLik  79.8557 eps 0.0050 iteration  919 qsel 4
#> logLik  79.9804 eps 0.0026 iteration  920 qsel 4
#> logLik  80.0557 eps 0.0003 iteration  921 qsel 4
#> logLik  80.1301 eps 0.0026 iteration  922 qsel 4
#> logLik  80.2042 eps 0.0003 iteration  923 qsel 4
#> logLik  80.2776 eps 0.0026 iteration  924 qsel 4
#> logLik  80.6759 eps 0.0048 iteration  925 qsel 4
#> logLik  80.7504 eps 0.0003 iteration  926 qsel 4
#> logLik  80.8240 eps 0.0023 iteration  927 qsel 4
#> logLik  80.9477 eps 0.0025 iteration  928 qsel 4
#> logLik  81.0215 eps 0.0003 iteration  929 qsel 4
#> logLik  81.0939 eps 0.0002 iteration  930 qsel 4
#> logLik  81.1664 eps 0.0024 iteration  931 qsel 4
#> logLik  81.2389 eps 0.0022 iteration  932 qsel 4
#> logLik  81.6330 eps 0.0049 iteration  933 qsel 4
#> logLik  81.7556 eps 0.0025 iteration  934 qsel 4
#> logLik  81.8290 eps 0.0002 iteration  935 qsel 4
#> logLik  81.9011 eps 0.0002 iteration  936 qsel 4
#> logLik  81.9726 eps 0.0023 iteration  937 qsel 4
#> logLik  82.0444 eps 0.0021 iteration  938 qsel 4
#> logLik  82.4334 eps 0.0052 iteration  939 qsel 4
#> logLik  82.5547 eps 0.0025 iteration  940 qsel 4
#> logLik  82.6278 eps 0.0002 iteration  941 qsel 4
#> logLik  82.6996 eps 0.0002 iteration  942 qsel 4
#> logLik  82.7704 eps 0.0021 iteration  943 qsel 4
#> logLik  82.8411 eps 0.0002 iteration  944 qsel 4
#> logLik  82.9118 eps 0.0024 iteration  945 qsel 4
#> logLik  83.2965 eps 0.0062 iteration  946 qsel 4
#> logLik  83.4169 eps 0.0027 iteration  947 qsel 4
#> logLik  83.4883 eps 0.0002 iteration  948 qsel 4
#> logLik  83.5583 eps 0.0021 iteration  949 qsel 4
#> logLik  83.6286 eps 0.0002 iteration  950 qsel 4
#> logLik  83.6984 eps 0.0025 iteration  951 qsel 4
#> logLik  84.0782 eps 0.0098 iteration  952 qsel 4
#> logLik  84.1974 eps 0.0034 iteration  953 qsel 4
#> logLik  84.2684 eps 0.0002 iteration  954 qsel 4
#> logLik  84.3381 eps 0.0002 iteration  955 qsel 4
#> logLik  84.4075 eps 0.0021 iteration  956 qsel 4
#> logLik  84.4763 eps 0.0036 iteration  957 qsel 4
#> logLik  84.8512 eps 0.0086 iteration  958 qsel 4
#> logLik  84.9214 eps 0.0002 iteration  959 qsel 4
#> logLik  84.9903 eps 0.0002 iteration  960 qsel 4
#> logLik  85.0589 eps 0.0024 iteration  961 qsel 4
#> logLik  85.1774 eps 0.0030 iteration  962 qsel 4
#> logLik  85.2456 eps 0.0002 iteration  963 qsel 4
#> logLik  85.3137 eps 0.0038 iteration  964 qsel 4
#> logLik  85.6844 eps 0.0075 iteration  965 qsel 4
#> logLik  85.7528 eps 0.0002 iteration  966 qsel 4
#> logLik  85.8207 eps 0.0027 iteration  967 qsel 4
#> logLik  85.9381 eps 0.0037 iteration  968 qsel 4
#> logLik  86.0060 eps 0.0002 iteration  969 qsel 4
#> logLik  86.0731 eps 0.0063 iteration  970 qsel 4
#> logLik  86.1399 eps 0.0055 iteration  971 qsel 4
#> logLik  86.5060 eps 0.0102 iteration  972 qsel 4
#> logLik  86.5743 eps 0.0002 iteration  973 qsel 4
#> logLik  86.6414 eps 0.0002 iteration  974 qsel 4
#> logLik  86.7076 eps 0.0029 iteration  975 qsel 4
#> logLik  86.8241 eps 0.0034 iteration  976 qsel 4
#> logLik  86.8905 eps 0.0002 iteration  977 qsel 4
#> logLik  86.9567 eps 0.0034 iteration  978 qsel 4
#> logLik  87.3187 eps 0.0066 iteration  979 qsel 4
#> logLik  87.3854 eps 0.0002 iteration  980 qsel 4
#> logLik  87.4507 eps 0.0002 iteration  981 qsel 4
#> logLik  87.5162 eps 0.0028 iteration  982 qsel 4
#> logLik  87.6318 eps 0.0030 iteration  983 qsel 4
#> logLik  87.6972 eps 0.0023 iteration  984 qsel 4
#> logLik  88.0546 eps 0.0063 iteration  985 qsel 4
#> logLik  88.1209 eps 0.0002 iteration  986 qsel 4
#> logLik  88.1859 eps 0.0002 iteration  987 qsel 4
#> logLik  88.2504 eps 0.0024 iteration  988 qsel 4
#> logLik  88.3649 eps 0.0029 iteration  989 qsel 4
#> logLik  88.4297 eps 0.0022 iteration  990 qsel 4
#> logLik  88.7825 eps 0.0171 iteration  991 qsel 4
#> logLik  88.8483 eps 0.0002 iteration  992 qsel 4
#> logLik  88.9129 eps 0.0002 iteration  993 qsel 4
#> logLik  88.9766 eps 0.0028 iteration  994 qsel 4
#> logLik  89.0406 eps 0.0025 iteration  995 qsel 4
#> logLik  89.1540 eps 0.0032 iteration  996 qsel 4
#> logLik  89.2182 eps 0.0002 iteration  997 qsel 4
#> logLik  89.2812 eps 0.0002 iteration  998 qsel 4
#> logLik  89.3439 eps 0.0023 iteration  999 qsel 4
#> logLik  89.6935 eps 0.0076 iteration 1000 qsel 4
#> 
#>  elapsed time:  2.88sec

## Plot estimated effects.
## plot(b)

## Print and plot the boosting summary.
boost_summary(b, plot = FALSE)
#> 
#> logLik. = 89.6935 -> at mstop = 1000 
#> ---
#>             mu % selected LogLik contrib.
#> s(x1)                21.6           69.70
#> s(x3)                19.3           55.76
#> te(lon,lat)           8.4           45.71
#> s(x2)                15.9           43.29
#> (Intercept)           0.0            0.00
#> x1                    0.0            0.00
#> x2                    0.0            0.00
#> x3                    0.0            0.00
#> lon                   0.0            0.00
#> lat                   0.0            0.00
#> s(lon)                0.0            0.00
#> s(lat)                0.0            0.00
#> ---
#>             sigma % selected LogLik contrib.
#> (Intercept)             34.7           58.18
#> 
## boost_plot(b, which = 1)
## boost_plot(b, which = 2)
## boost_plot(b, which = 3, name = "mu.s.te(lon,lat).")

## Extract estimated parameters for certain
## boosting iterations.
parameters(b, mstop = 1)
#>     mu.p.(Intercept)              mu.p.x1              mu.p.x2 
#>          -0.05841161           0.00000000           0.00000000 
#>              mu.p.x3             mu.p.lon             mu.p.lat 
#>           0.00000000           0.00000000           0.00000000 
#>        mu.s.s(x1).b1        mu.s.s(x1).b2        mu.s.s(x1).b3 
#>           0.00000000           0.00000000           0.00000000 
#>        mu.s.s(x1).b4        mu.s.s(x1).b5        mu.s.s(x1).b6 
#>           0.00000000           0.00000000           0.00000000 
#>        mu.s.s(x1).b7        mu.s.s(x1).b8        mu.s.s(x1).b9 
#>           0.00000000           0.00000000           0.00000000 
#>        mu.s.s(x2).b1        mu.s.s(x2).b2        mu.s.s(x2).b3 
#>           0.00000000           0.00000000           0.00000000 
#>        mu.s.s(x2).b4        mu.s.s(x2).b5        mu.s.s(x2).b6 
#>           0.00000000           0.00000000           0.00000000 
#>        mu.s.s(x2).b7        mu.s.s(x2).b8        mu.s.s(x2).b9 
#>           0.00000000           0.00000000           0.00000000 
#>        mu.s.s(x3).b1        mu.s.s(x3).b2        mu.s.s(x3).b3 
#>           0.00000000           0.00000000           0.00000000 
#>        mu.s.s(x3).b4        mu.s.s(x3).b5        mu.s.s(x3).b6 
#>           0.00000000           0.00000000           0.00000000 
#>        mu.s.s(x3).b7        mu.s.s(x3).b8        mu.s.s(x3).b9 
#>           0.00000000           0.00000000           0.00000000 
#>       mu.s.s(lon).b1       mu.s.s(lon).b2       mu.s.s(lon).b3 
#>           0.00000000           0.00000000           0.00000000 
#>       mu.s.s(lon).b4       mu.s.s(lon).b5       mu.s.s(lon).b6 
#>           0.00000000           0.00000000           0.00000000 
#>       mu.s.s(lon).b7       mu.s.s(lon).b8       mu.s.s(lon).b9 
#>           0.00000000           0.00000000           0.00000000 
#>       mu.s.s(lat).b1       mu.s.s(lat).b2       mu.s.s(lat).b3 
#>           0.00000000           0.00000000           0.00000000 
#>       mu.s.s(lat).b4       mu.s.s(lat).b5       mu.s.s(lat).b6 
#>           0.00000000           0.00000000           0.00000000 
#>       mu.s.s(lat).b7       mu.s.s(lat).b8       mu.s.s(lat).b9 
#>           0.00000000           0.00000000           0.00000000 
#>  mu.s.te(lon,lat).b1  mu.s.te(lon,lat).b2  mu.s.te(lon,lat).b3 
#>           0.00000000           0.00000000           0.00000000 
#>  mu.s.te(lon,lat).b4  mu.s.te(lon,lat).b5  mu.s.te(lon,lat).b6 
#>           0.00000000           0.00000000           0.00000000 
#>  mu.s.te(lon,lat).b7  mu.s.te(lon,lat).b8  mu.s.te(lon,lat).b9 
#>           0.00000000           0.00000000           0.00000000 
#> mu.s.te(lon,lat).b10 mu.s.te(lon,lat).b11 mu.s.te(lon,lat).b12 
#>           0.00000000           0.00000000           0.00000000 
#> mu.s.te(lon,lat).b13 mu.s.te(lon,lat).b14 mu.s.te(lon,lat).b15 
#>           0.00000000           0.00000000           0.00000000 
#> mu.s.te(lon,lat).b16 mu.s.te(lon,lat).b17 mu.s.te(lon,lat).b18 
#>           0.00000000           0.00000000           0.00000000 
#> mu.s.te(lon,lat).b19 mu.s.te(lon,lat).b20 mu.s.te(lon,lat).b21 
#>           0.00000000           0.00000000           0.00000000 
#> mu.s.te(lon,lat).b22 mu.s.te(lon,lat).b23 mu.s.te(lon,lat).b24 
#>           0.00000000           0.00000000           0.00000000 
#>  sigma.p.(Intercept) 
#>          -1.05304572 
parameters(b, mstop = 100)
#>     mu.p.(Intercept)              mu.p.x1              mu.p.x2 
#>         -0.058411614          0.000000000          0.000000000 
#>              mu.p.x3             mu.p.lon             mu.p.lat 
#>          0.000000000          0.000000000          0.000000000 
#>        mu.s.s(x1).b1        mu.s.s(x1).b2        mu.s.s(x1).b3 
#>          0.039268945         -0.001350402          0.015219138 
#>        mu.s.s(x1).b4        mu.s.s(x1).b5        mu.s.s(x1).b6 
#>          0.009161954          0.007104356          0.005601264 
#>        mu.s.s(x1).b7        mu.s.s(x1).b8        mu.s.s(x1).b9 
#>          0.001180029          0.001274986         -0.126444256 
#>        mu.s.s(x2).b1        mu.s.s(x2).b2        mu.s.s(x2).b3 
#>         -0.001271167          0.069043611          0.001003697 
#>        mu.s.s(x2).b4        mu.s.s(x2).b5        mu.s.s(x2).b6 
#>         -0.024173043          0.004058185         -0.021312378 
#>        mu.s.s(x2).b7        mu.s.s(x2).b8        mu.s.s(x2).b9 
#>         -0.007109845          0.055989599          0.007363287 
#>        mu.s.s(x3).b1        mu.s.s(x3).b2        mu.s.s(x3).b3 
#>          0.024911838         -0.029259307         -0.012201572 
#>        mu.s.s(x3).b4        mu.s.s(x3).b5        mu.s.s(x3).b6 
#>          0.022586312          0.016285732          0.035045844 
#>        mu.s.s(x3).b7        mu.s.s(x3).b8        mu.s.s(x3).b9 
#>          0.011050418         -0.160889392         -0.021359849 
#>       mu.s.s(lon).b1       mu.s.s(lon).b2       mu.s.s(lon).b3 
#>          0.000000000          0.000000000          0.000000000 
#>       mu.s.s(lon).b4       mu.s.s(lon).b5       mu.s.s(lon).b6 
#>          0.000000000          0.000000000          0.000000000 
#>       mu.s.s(lon).b7       mu.s.s(lon).b8       mu.s.s(lon).b9 
#>          0.000000000          0.000000000          0.000000000 
#>       mu.s.s(lat).b1       mu.s.s(lat).b2       mu.s.s(lat).b3 
#>          0.000000000          0.000000000          0.000000000 
#>       mu.s.s(lat).b4       mu.s.s(lat).b5       mu.s.s(lat).b6 
#>          0.000000000          0.000000000          0.000000000 
#>       mu.s.s(lat).b7       mu.s.s(lat).b8       mu.s.s(lat).b9 
#>          0.000000000          0.000000000          0.000000000 
#>  mu.s.te(lon,lat).b1  mu.s.te(lon,lat).b2  mu.s.te(lon,lat).b3 
#>          0.000000000          0.000000000          0.000000000 
#>  mu.s.te(lon,lat).b4  mu.s.te(lon,lat).b5  mu.s.te(lon,lat).b6 
#>          0.000000000          0.000000000          0.000000000 
#>  mu.s.te(lon,lat).b7  mu.s.te(lon,lat).b8  mu.s.te(lon,lat).b9 
#>          0.000000000          0.000000000          0.000000000 
#> mu.s.te(lon,lat).b10 mu.s.te(lon,lat).b11 mu.s.te(lon,lat).b12 
#>          0.000000000          0.000000000          0.000000000 
#> mu.s.te(lon,lat).b13 mu.s.te(lon,lat).b14 mu.s.te(lon,lat).b15 
#>          0.000000000          0.000000000          0.000000000 
#> mu.s.te(lon,lat).b16 mu.s.te(lon,lat).b17 mu.s.te(lon,lat).b18 
#>          0.000000000          0.000000000          0.000000000 
#> mu.s.te(lon,lat).b19 mu.s.te(lon,lat).b20 mu.s.te(lon,lat).b21 
#>          0.000000000          0.000000000          0.000000000 
#> mu.s.te(lon,lat).b22 mu.s.te(lon,lat).b23 mu.s.te(lon,lat).b24 
#>          0.000000000          0.000000000          0.000000000 
#>  sigma.p.(Intercept) 
#>         -1.053045720 

## Also works with predict().
head(do.call("cbind", predict(b, mstop = 1)))
#>               mu     sigma
#> [1,] -0.05841161 -1.053046
#> [2,] -0.05841161 -1.053046
#> [3,] -0.05841161 -1.053046
#> [4,] -0.05841161 -1.053046
#> [5,] -0.05841161 -1.053046
#> [6,] -0.05841161 -1.053046
head(do.call("cbind", predict(b, mstop = 100)))
#>               mu     sigma
#> [1,]  0.01910224 -1.053046
#> [2,]  0.02862474 -1.053046
#> [3,]  0.01147507 -1.053046
#> [4,] -0.02623836 -1.053046
#> [5,] -0.15807300 -1.053046
#> [6,] -0.01981386 -1.053046

## Another example using the modified likelihood
## boosting algorithm.
f <- list(
  num ~ x1 + x2 + x3 + lon + lat +
    s(x1) + s(x2) + s(x3) + s(lon) + s(lat) + te(lon,lat),
  sigma ~ x1 + x2 + x3 + lon + lat +
    s(x1) + s(x2) + s(x3) + s(lon) + s(lat) + te(lon,lat)
)

b <- bamlss(f, data = d, optimizer = opt_boostm,
  sampler = FALSE, scale.d = TRUE, nu = 0.05,
  maxit = 400, stop.criterion = "AIC", force.stop = FALSE)
#> 
#> AIC 342.6370 logLik -165.617 edf 5.70  eps 0.2729 iteration   2 qsel 1
#> AIC 317.4701 logLik -153.617 edf 5.11  eps 0.8732 iteration   3 qsel 1
#> AIC 299.7005 logLik -144.837 edf 5.01  eps 0.4551 iteration   4 qsel 1
#> AIC 286.6101 logLik -136.499 edf 6.80  eps 0.5352 iteration   5 qsel 2
#> AIC 273.5653 logLik -129.989 edf 6.79  eps 0.3087 iteration   6 qsel 2
#> AIC 264.4006 logLik -125.410 edf 6.78  eps 0.1677 iteration   7 qsel 2
#> AIC 255.5750 logLik -121.381 edf 6.40  eps 0.3306 iteration   8 qsel 2
#> AIC 247.9529 logLik -115.907 edf 8.06  eps 0.4191 iteration   9 qsel 3
#> AIC 240.3521 logLik -112.304 edf 7.87  eps 0.2329 iteration  10 qsel 3
#> AIC 233.1309 logLik -108.695 edf 7.86  eps 0.1806 iteration  11 qsel 3
#> AIC 226.5740 logLik -105.445 edf 7.84  eps 0.1467 iteration  12 qsel 3
#> AIC 220.0163 logLik -102.163 edf 7.84  eps 0.1087 iteration  13 qsel 3
#> AIC 214.1120 logLik -99.2019 edf 7.85  eps 0.0723 iteration  14 qsel 3
#> AIC 208.5310 logLik -96.4445 edf 7.82  eps 0.1589 iteration  15 qsel 3
#> AIC 203.2088 logLik -93.7769 edf 7.82  eps 0.2531 iteration  16 qsel 3
#> AIC 198.4152 logLik -91.3695 edf 7.83  eps 0.0531 iteration  17 qsel 3
#> AIC 194.0921 logLik -89.1967 edf 7.84  eps 0.0852 iteration  18 qsel 3
#> AIC 189.9517 logLik -87.1583 edf 7.81  eps 0.2007 iteration  19 qsel 3
#> AIC 186.0524 logLik -85.1994 edf 7.82  eps 0.0575 iteration  20 qsel 3
#> AIC 182.1352 logLik -83.4495 edf 7.61  eps 0.7721 iteration  21 qsel 3
#> AIC 178.5933 logLik -81.6707 edf 7.62  eps 0.0817 iteration  22 qsel 3
#> AIC 175.4069 logLik -80.0647 edf 7.63  eps 0.1155 iteration  23 qsel 3
#> AIC 172.3283 logLik -78.5459 edf 7.61  eps 0.1152 iteration  24 qsel 3
#> AIC 169.4546 logLik -77.0982 edf 7.62  eps 0.0430 iteration  25 qsel 3
#> AIC 166.8687 logLik -75.7909 edf 7.64  eps 0.1481 iteration  26 qsel 3
#> AIC 164.5381 logLik -74.6104 edf 7.65  eps 0.0552 iteration  27 qsel 3
#> AIC 162.2282 logLik -72.4554 edf 8.65  eps 0.1934 iteration  28 qsel 4
#> AIC 159.5092 logLik -71.0960 edf 8.65  eps 0.1042 iteration  29 qsel 4
#> AIC 157.2037 logLik -69.9550 edf 8.64  eps 0.0606 iteration  30 qsel 4
#> AIC 155.0331 logLik -68.8662 edf 8.65  eps 0.0353 iteration  31 qsel 4
#> AIC 153.0994 logLik -67.8827 edf 8.66  eps 0.0586 iteration  32 qsel 4
#> AIC 151.3325 logLik -66.9993 edf 8.66  eps 0.1045 iteration  33 qsel 4
#> AIC 149.5649 logLik -66.1016 edf 8.68  eps 0.0526 iteration  34 qsel 4
#> AIC 147.8498 logLik -65.2553 edf 8.66  eps 0.0779 iteration  35 qsel 4
#> AIC 146.2610 logLik -64.4449 edf 8.68  eps 0.0306 iteration  36 qsel 4
#> AIC 144.8354 logLik -63.7130 edf 8.70  eps 0.0260 iteration  37 qsel 4
#> AIC 143.5296 logLik -63.0639 edf 8.70  eps 0.0418 iteration  38 qsel 4
#> AIC 142.2431 logLik -62.4031 edf 8.71  eps 0.0288 iteration  39 qsel 4
#> AIC 141.0534 logLik -61.8083 edf 8.71  eps 0.0332 iteration  40 qsel 4
#> AIC 139.8818 logLik -61.2049 edf 8.73  eps 0.0274 iteration  41 qsel 4
#> AIC 138.8355 logLik -60.6600 edf 8.75  eps 0.0261 iteration  42 qsel 4
#> AIC 137.8056 logLik -60.1448 edf 8.75  eps 0.0308 iteration  43 qsel 4
#> AIC 136.8642 logLik -59.6532 edf 8.77  eps 0.0320 iteration  44 qsel 4
#> AIC 136.0137 logLik -59.2219 edf 8.78  eps 0.0332 iteration  45 qsel 4
#> AIC 135.1710 logLik -58.7790 edf 8.80  eps 0.0283 iteration  46 qsel 4
#> AIC 134.3783 logLik -58.3826 edf 8.80  eps 0.0297 iteration  47 qsel 4
#> AIC 133.6133 logLik -57.9792 edf 8.82  eps 0.0193 iteration  48 qsel 4
#> AIC 132.8894 logLik -57.6094 edf 8.83  eps 0.0265 iteration  49 qsel 4
#> AIC 132.2101 logLik -57.2458 edf 8.85  eps 0.0314 iteration  50 qsel 4
#> AIC 131.5717 logLik -56.9164 edf 8.86  eps 0.0276 iteration  51 qsel 4
#> AIC 130.9665 logLik -56.5890 edf 8.89  eps 1.5190 iteration  52 qsel 4
#> AIC 130.3910 logLik -56.2902 edf 8.90  eps 0.0340 iteration  53 qsel 4
#> AIC 129.8521 logLik -55.9961 edf 8.92  eps 0.0231 iteration  54 qsel 4
#> AIC 129.3243 logLik -55.7322 edf 8.92  eps 0.0614 iteration  55 qsel 4
#> AIC 128.8039 logLik -55.4594 edf 8.94  eps 0.0366 iteration  56 qsel 4
#> AIC 128.3140 logLik -55.1932 edf 8.96  eps 0.0395 iteration  57 qsel 4
#> AIC 127.8342 logLik -54.9403 edf 8.97  eps 0.0177 iteration  58 qsel 4
#> AIC 127.3884 logLik -54.7022 edf 8.99  eps 0.0140 iteration  59 qsel 4
#> AIC 126.9619 logLik -54.4641 edf 9.01  eps 0.0148 iteration  60 qsel 4
#> AIC 126.5465 logLik -54.2412 edf 9.03  eps 0.0135 iteration  61 qsel 4
#> AIC 126.1585 logLik -54.0303 edf 9.04  eps 0.0236 iteration  62 qsel 4
#> AIC 125.7861 logLik -53.8177 edf 9.07  eps 0.0138 iteration  63 qsel 4
#> AIC 125.4232 logLik -53.6192 edf 9.09  eps 0.0214 iteration  64 qsel 4
#> AIC 125.0770 logLik -53.4461 edf 9.09  eps 0.0188 iteration  65 qsel 4
#> AIC 124.7430 logLik -53.2592 edf 9.11  eps 0.0121 iteration  66 qsel 4
#> AIC 124.4093 logLik -53.0684 edf 9.13  eps 0.0213 iteration  67 qsel 4
#> AIC 124.0931 logLik -52.8906 edf 9.15  eps 0.0117 iteration  68 qsel 4
#> AIC 123.7950 logLik -52.7201 edf 9.17  eps 0.0106 iteration  69 qsel 4
#> AIC 123.5092 logLik -52.5501 edf 9.20  eps 0.0178 iteration  70 qsel 4
#> AIC 123.2277 logLik -52.3870 edf 9.22  eps 0.0157 iteration  71 qsel 4
#> AIC 122.9611 logLik -52.2284 edf 9.25  eps 0.0112 iteration  72 qsel 4
#> AIC 122.7058 logLik -52.0724 edf 9.28  eps 0.0332 iteration  73 qsel 4
#> AIC 122.4557 logLik -51.9212 edf 9.30  eps 0.0099 iteration  74 qsel 4
#> AIC 122.2085 logLik -51.7677 edf 9.33  eps 0.0361 iteration  75 qsel 4
#> AIC 121.9699 logLik -51.6123 edf 9.37  eps 0.0114 iteration  76 qsel 4
#> AIC 121.7355 logLik -51.4511 edf 9.41  eps 0.0094 iteration  77 qsel 4
#> AIC 121.4573 logLik -51.3129 edf 9.41  eps 0.0075 iteration  78 qsel 4
#> AIC 121.2151 logLik -51.1925 edf 9.41  eps 0.0063 iteration  79 qsel 4
#> AIC 120.9911 logLik -51.0805 edf 9.41  eps 0.0241 iteration  80 qsel 4
#> AIC 120.7589 logLik -50.9452 edf 9.43  eps 0.0148 iteration  81 qsel 4
#> AIC 120.5759 logLik -50.8250 edf 9.46  eps 0.0094 iteration  82 qsel 4
#> AIC 120.4204 logLik -50.7182 edf 9.49  eps 0.0078 iteration  83 qsel 4
#> AIC 120.2687 logLik -50.6423 edf 9.49  eps 0.0114 iteration  84 qsel 4
#> AIC 120.1288 logLik -50.5465 edf 9.51  eps 0.0088 iteration  85 qsel 4
#> AIC 120.0170 logLik -50.4615 edf 9.54  eps 0.0112 iteration  86 qsel 4
#> AIC 119.9137 logLik -50.4098 edf 9.54  eps 0.0098 iteration  87 qsel 4
#> AIC 119.8143 logLik -50.3337 edf 9.57  eps 0.0068 iteration  88 qsel 4
#> AIC 119.7236 logLik -49.2884 edf 10.5  eps 0.0092 iteration  89 qsel 5
#> AIC 118.3341 logLik -48.5937 edf 10.5  eps 0.0081 iteration  90 qsel 5
#> AIC 117.4665 logLik -48.1599 edf 10.5  eps 0.0069 iteration  91 qsel 5
#> AIC 116.7648 logLik -47.6714 edf 10.7  eps 0.1396 iteration  92 qsel 5
#> AIC 116.1922 logLik -47.3637 edf 10.7  eps 0.0261 iteration  93 qsel 5
#> AIC 115.6583 logLik -47.0968 edf 10.7  eps 0.0060 iteration  94 qsel 5
#> AIC 115.2587 logLik -46.8812 edf 10.7  eps 0.0362 iteration  95 qsel 5
#> AIC 114.9798 logLik -46.7417 edf 10.7  eps 0.0051 iteration  96 qsel 5
#> AIC 114.7844 logLik -46.6418 edf 10.7  eps 0.0444 iteration  97 qsel 5
#> AIC 114.6732 logLik -46.5862 edf 10.7  eps 0.0041 iteration  98 qsel 5
#> AIC 114.5974 logLik -46.5483 edf 10.7  eps 0.0115 iteration  99 qsel 5
#> AIC 114.5443 logLik -46.5217 edf 10.7  eps 0.0052 iteration 100 qsel 5
#> AIC 114.5066 logLik -46.5029 edf 10.7  eps 0.0034 iteration 101 qsel 5
#> AIC 114.4795 logLik -46.4894 edf 10.7  eps 0.0025 iteration 102 qsel 5
#> AIC 114.4578 logLik -46.4785 edf 10.7  eps 0.0032 iteration 103 qsel 5
#> AIC 114.4372 logLik -46.4682 edf 10.7  eps 0.0020 iteration 104 qsel 5
#> AIC 114.4222 logLik -46.4607 edf 10.7  eps 0.0016 iteration 105 qsel 5
#> AIC 114.4112 logLik -46.4552 edf 10.7  eps 0.0013 iteration 106 qsel 5
#> AIC 114.4032 logLik -46.4512 edf 10.7  eps 0.0011 iteration 107 qsel 5
#> AIC 114.3973 logLik -46.4482 edf 10.7  eps 0.0011 iteration 108 qsel 5
#> AIC 114.3930 logLik -46.4461 edf 10.7  eps 0.0013 iteration 109 qsel 5
#> AIC 114.3898 logLik -46.4445 edf 10.7  eps 0.0030 iteration 110 qsel 5
#> AIC 114.3874 logLik -46.4433 edf 10.7  eps 0.0021 iteration 111 qsel 5
#> AIC 114.3857 logLik -46.4424 edf 10.7  eps 0.0007 iteration 112 qsel 5
#> AIC 114.3844 logLik -46.4418 edf 10.7  eps 0.0004 iteration 113 qsel 5
#> AIC 114.3834 logLik -46.4413 edf 10.7  eps 0.0002 iteration 114 qsel 5
#> AIC 114.3827 logLik -46.4409 edf 10.7  eps 0.0001 iteration 115 qsel 5
#> AIC 114.3822 logLik -46.4407 edf 10.7  eps 0.0001 iteration 116 qsel 5
#> AIC 114.3818 logLik -46.4405 edf 10.7  eps 0.0000 iteration 117 qsel 5
#> AIC 114.3815 logLik -46.4403 edf 10.7  eps 0.0000 iteration 118 qsel 5
#> AIC 114.3813 logLik -46.4402 edf 10.7  eps 0.0000 iteration 119 qsel 5
#> AIC 114.3811 logLik -46.4402 edf 10.7  eps 0.0000 iteration 120 qsel 5
#> AIC 114.3810 logLik -46.4401 edf 10.7  eps 0.0000 iteration 121 qsel 5
#> AIC 114.3809 logLik -46.4401 edf 10.7  eps 0.0000 iteration 122 qsel 5
#> AIC 114.3809 logLik -46.4400 edf 10.7  eps 0.0000 iteration 123 qsel 5
#> AIC 114.3808 logLik -46.4400 edf 10.7  eps 0.0000 iteration 124 qsel 5
#> AIC 114.3808 logLik -46.4400 edf 10.7  eps 0.0000 iteration 125 qsel 5
#> AIC 114.3808 logLik -46.4400 edf 10.7  eps 0.0000 iteration 126 qsel 5
#> AIC 114.3808 logLik -46.4400 edf 10.7  eps 0.0000 iteration 127 qsel 5
#> AIC 114.3807 logLik -46.4400 edf 10.7  eps 0.0000 iteration 128 qsel 5
#> AIC 114.3807 logLik -46.4400 edf 10.7  eps 0.0000 iteration 129 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 130 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 131 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 132 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 133 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 134 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 135 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 136 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 137 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 138 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 139 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 140 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 141 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 142 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 143 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 144 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 145 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 146 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 147 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 148 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 149 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 150 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 151 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 152 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 153 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 154 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 155 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 156 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 157 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 158 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 159 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 160 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 161 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 162 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 163 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 164 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 165 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 166 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 167 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 168 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 169 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 170 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 171 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 172 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 173 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 174 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 175 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 176 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 177 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 178 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 179 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 180 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 181 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 182 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 183 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 184 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 185 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 186 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 187 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 188 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 189 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 190 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 191 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 192 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 193 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 194 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 195 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 196 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 197 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 198 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 199 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 200 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 201 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 202 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 203 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 204 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 205 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 206 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 207 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 208 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 209 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 210 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 211 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 212 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 213 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 214 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 215 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 216 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 217 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 218 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 219 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 220 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 221 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 222 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 223 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 224 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 225 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 226 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 227 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 228 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 229 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 230 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 231 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 232 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 233 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 234 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 235 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 236 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 237 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 238 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 239 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 240 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 241 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 242 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 243 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 244 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 245 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 246 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 247 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 248 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 249 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 250 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 251 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 252 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 253 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 254 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 255 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 256 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 257 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 258 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 259 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 260 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 261 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 262 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 263 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 264 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 265 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 266 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 267 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 268 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 269 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 270 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 271 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 272 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 273 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 274 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 275 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 276 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 277 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 278 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 279 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 280 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 281 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 282 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 283 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 284 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 285 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 286 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 287 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 288 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 289 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 290 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 291 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 292 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 293 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 294 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 295 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 296 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 297 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 298 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 299 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 300 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 301 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 302 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 303 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 304 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 305 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 306 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 307 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 308 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 309 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 310 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 311 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 312 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 313 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 314 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 315 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 316 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 317 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 318 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 319 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 320 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 321 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 322 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 323 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 324 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 325 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 326 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 327 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 328 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 329 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 330 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 331 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 332 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 333 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 334 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 335 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 336 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 337 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 338 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 339 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 340 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 341 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 342 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 343 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 344 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 345 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 346 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 347 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 348 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 349 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 350 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 351 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 352 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 353 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 354 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 355 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 356 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 357 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 358 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 359 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 360 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 361 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 362 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 363 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 364 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 365 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 366 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 367 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 368 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 369 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 370 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 371 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 372 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 373 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 374 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 375 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 376 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 377 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 378 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 379 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 380 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 381 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 382 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 383 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 384 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 385 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 386 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 387 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 388 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 389 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 390 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 391 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 392 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 393 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 394 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 395 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 396 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 397 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 398 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 399 qsel 5
#> AIC 114.3807 logLik -46.4399 edf 10.7  eps 0.0000 iteration 400 qsel 5
#> 
#>  elapsed time: 39.22sec


## Plot estimated effects.
## plot(b)

## Plot AIC and log-lik contributions.
## boost_plot(b, "AIC")
## boost_plot(b, "loglik.contrib")

## Out-of-sample selection of model terms.
set.seed(123)
d <- GAMart(n = 5000)

## Split data into training and testing
i <- sample(1:2, size = nrow(d), replace = TRUE)
dtest <- subset(d, i == 1)
dtrain <- subset(d, i == 2)

## Model formula
f <- list(
  num ~ s(x1) + s(x2) + s(x3),
  sigma ~ s(x1) + s(x2) + s(x3)
)

## Create model frame for out-of-sample selection.
sm <- boost_frame(f, train = dtrain, test = dtest, family = "gaussian")

## Out-of-sample selection function.
sfun <- function(parameters) {
  sm$parameters <- parameters
  p <- predict(sm, type = "parameter")
  -1 * sum(sm$family$d(dtest$num, p, log = TRUE))
}

## Start boosting with out-of-sample negative
## log-likelihood selection of model terms.
b <- bamlss(f, data = dtrain, sampler = FALSE, optimizer = opt_boost,
  selectfun = sfun, always = "best")
#> 
#> userIC 2413.338 logLik -598.545 eps 0.0531 iteration   2 qsel 1
#> userIC 2409.204 logLik -546.484 eps 0.0509 iteration   3 qsel 2
#> userIC 2408.146 logLik -501.945 eps 0.0545 iteration   4 qsel 2
#> userIC 2404.738 logLik -459.498 eps 0.0536 iteration   5 qsel 2
#> userIC 2403.927 logLik -423.168 eps 0.1060 iteration   6 qsel 2
#> userIC 2401.114 logLik -388.559 eps 0.1721 iteration   7 qsel 2
#> userIC 2399.574 logLik -356.674 eps 0.1006 iteration   8 qsel 3
#> userIC 2397.478 logLik -327.295 eps 0.3657 iteration   9 qsel 3
#> userIC 2395.112 logLik -299.251 eps 0.1390 iteration  10 qsel 3
#> userIC 2393.894 logLik -273.542 eps 0.1068 iteration  11 qsel 3
#> userIC 2392.253 logLik -249.784 eps 0.0981 iteration  12 qsel 3
#> userIC 2390.262 logLik -227.060 eps 0.2394 iteration  13 qsel 3
#> userIC 2389.281 logLik -206.297 eps 0.1422 iteration  14 qsel 3
#> userIC 2388.016 logLik -187.084 eps 0.1092 iteration  15 qsel 3
#> userIC 2386.337 logLik -168.671 eps 0.1048 iteration  16 qsel 3
#> userIC 2385.532 logLik -151.868 eps 0.1188 iteration  17 qsel 3
#> userIC 2384.575 logLik -136.330 eps 0.0938 iteration  18 qsel 3
#> userIC 2383.158 logLik -121.410 eps 0.1196 iteration  19 qsel 3
#> userIC 2382.483 logLik -107.780 eps 0.1400 iteration  20 qsel 3
#> userIC 2381.776 logLik -95.2136 eps 0.0898 iteration  21 qsel 3
#> userIC 2380.580 logLik -83.1257 eps 0.1321 iteration  22 qsel 3
#> userIC 2380.001 logLik -72.0359 eps 0.0879 iteration  23 qsel 3
#> userIC 2379.490 logLik -62.3061 eps 0.0778 iteration  24 qsel 3
#> userIC 2378.523 logLik -52.0815 eps 0.1159 iteration  25 qsel 3
#> userIC 2377.979 logLik -43.0269 eps 0.1050 iteration  26 qsel 3
#> userIC 2377.587 logLik -35.1386 eps 0.0663 iteration  27 qsel 3
#> userIC 2376.848 logLik -26.8749 eps 0.1179 iteration  28 qsel 3
#> userIC 2376.330 logLik -19.4503 eps 0.0809 iteration  29 qsel 3
#> userIC 2376.036 logLik -13.0555 eps 0.0874 iteration  30 qsel 3
#> userIC 2375.477 logLik  -6.3765 eps 0.0890 iteration  31 qsel 3
#> userIC 2374.982 logLik  -0.2570 eps 0.0854 iteration  32 qsel 3
#> userIC 2374.769 logLik   4.9264 eps 0.0690 iteration  33 qsel 3
#> userIC 2374.349 logLik  10.3245 eps 0.3250 iteration  34 qsel 3
#> userIC 2373.878 logLik  15.4001 eps 0.0592 iteration  35 qsel 3
#> userIC 2373.729 logLik  19.6011 eps 0.0460 iteration  36 qsel 3
#> userIC 2373.383 logLik  23.8714 eps 0.0804 iteration  37 qsel 3
#> userIC 2372.992 logLik  28.1978 eps 0.0558 iteration  38 qsel 3
#> userIC 2372.874 logLik  31.6019 eps 0.0490 iteration  39 qsel 3
#> userIC 2372.559 logLik  35.1998 eps 0.0411 iteration  40 qsel 3
#> userIC 2372.284 logLik  38.6987 eps 0.0569 iteration  41 qsel 3
#> userIC 2372.152 logLik  41.7811 eps 0.0983 iteration  42 qsel 3
#> userIC 2371.852 logLik  44.5202 eps 0.0783 iteration  43 qsel 3
#> userIC 2371.692 logLik  47.3496 eps 0.0628 iteration  44 qsel 3
#> userIC 2371.539 logLik  50.0094 eps 0.0344 iteration  45 qsel 3
#> userIC 2371.304 logLik  52.2288 eps 0.1278 iteration  46 qsel 3
#> userIC 2371.194 logLik  54.5166 eps 0.0387 iteration  47 qsel 3
#> userIC 2371.029 logLik  56.8489 eps 0.0669 iteration  48 qsel 3
#> userIC 2370.842 logLik  58.6465 eps 0.0471 iteration  49 qsel 3
#> userIC 2370.769 logLik  60.4962 eps 0.0331 iteration  50 qsel 3
#> userIC 2370.603 logLik  62.5801 eps 0.0328 iteration  51 qsel 3
#> userIC 2370.449 logLik  64.0356 eps 0.0312 iteration  52 qsel 3
#> userIC 2370.382 logLik  65.9383 eps 0.6300 iteration  53 qsel 3
#> userIC 2370.210 logLik  67.4205 eps 0.0266 iteration  54 qsel 3
#> userIC 2370.110 logLik  68.5985 eps 0.0424 iteration  55 qsel 3
#> userIC 2370.049 logLik  70.3715 eps 0.0226 iteration  56 qsel 3
#> userIC 2369.907 logLik  71.5704 eps 0.0226 iteration  57 qsel 3
#> userIC 2369.815 logLik  72.5233 eps 0.0275 iteration  58 qsel 3
#> userIC 2369.763 logLik  74.2077 eps 0.0460 iteration  59 qsel 3
#> userIC 2369.637 logLik  75.1778 eps 0.0180 iteration  60 qsel 3
#> userIC 2369.553 logLik  75.9481 eps 0.0152 iteration  61 qsel 3
#> userIC 2369.515 logLik  77.5642 eps 0.0875 iteration  62 qsel 3
#> userIC 2369.396 logLik  78.1818 eps 0.0530 iteration  63 qsel 3
#> userIC 2369.331 logLik  78.9716 eps 0.0240 iteration  64 qsel 3
#> userIC 2369.299 logLik  80.5106 eps 0.0273 iteration  65 qsel 3
#> userIC 2369.177 logLik  81.0100 eps 0.0129 iteration  66 qsel 3
#> userIC 2369.131 logLik  81.6499 eps 0.0393 iteration  67 qsel 3
#> userIC 2369.110 logLik  83.0792 eps 0.0206 iteration  68 qsel 3
#> userIC 2368.988 logLik  83.4832 eps 0.0116 iteration  69 qsel 3
#> userIC 2368.958 logLik  84.0026 eps 0.0171 iteration  70 qsel 3
#> userIC 2368.939 logLik  84.3331 eps 0.0104 iteration  71 qsel 3
#> userIC 2368.912 logLik  85.6200 eps 0.0154 iteration  72 qsel 3
#> userIC 2368.811 logLik  86.0426 eps 0.0104 iteration  73 qsel 3
#> userIC 2368.791 logLik  86.3101 eps 0.0350 iteration  74 qsel 3
#> userIC 2368.781 logLik  87.4468 eps 0.0227 iteration  75 qsel 3
#> userIC 2368.691 logLik  87.7915 eps 0.0118 iteration  76 qsel 3
#> userIC 2368.670 logLik  88.0082 eps 0.0117 iteration  77 qsel 3
#> userIC 2368.666 logLik  88.2900 eps 0.0087 iteration  78 qsel 3
#> userIC 2368.644 logLik 172.9726 eps 0.0244 iteration  79 qsel 4
#> userIC 2368.571 logLik 174.0924 eps 0.0188 iteration  80 qsel 4
#> userIC 2368.507 logLik 244.7876 eps 0.0214 iteration  81 qsel 4
#> userIC 2368.435 logLik 302.3209 eps 0.0186 iteration  82 qsel 4
#> userIC 2368.353 logLik 303.5135 eps 0.0124 iteration  83 qsel 4
#> userIC 2368.299 logLik 349.4761 eps 0.0161 iteration  84 qsel 4
#> userIC 2368.214 logLik 385.2245 eps 0.0138 iteration  85 qsel 4
#> userIC 2368.111 logLik 386.4442 eps 0.0161 iteration  86 qsel 4
#> userIC 2368.060 logLik 386.7164 eps 0.0133 iteration  87 qsel 4
#> userIC 2368.041 logLik 414.0993 eps 0.0118 iteration  88 qsel 4
#> userIC 2367.931 logLik 415.2377 eps 0.0097 iteration  89 qsel 4
#> userIC 2367.889 logLik 435.7674 eps 0.0100 iteration  90 qsel 4
#> userIC 2367.772 logLik 436.8224 eps 0.0081 iteration  91 qsel 4
#> userIC 2367.739 logLik 437.0671 eps 0.0070 iteration  92 qsel 4
#> userIC 2367.724 logLik 452.1933 eps 0.0085 iteration  93 qsel 4
#> userIC 2367.602 logLik 453.1626 eps 0.0091 iteration  94 qsel 4
#> userIC 2367.575 logLik 464.1367 eps 0.0071 iteration  95 qsel 5
#> userIC 2367.476 logLik 465.0218 eps 0.0185 iteration  96 qsel 5
#> userIC 2367.450 logLik 465.8041 eps 0.0086 iteration  97 qsel 5
#> userIC 2367.428 logLik 466.4965 eps 0.0114 iteration  98 qsel 5
#> userIC 2367.407 logLik 467.1109 eps 0.0099 iteration  99 qsel 5
#> userIC 2367.387 logLik 467.6570 eps 0.0289 iteration 100 qsel 5
#> userIC 2367.367 logLik 468.1428 eps 0.0119 iteration 101 qsel 5
#> userIC 2367.348 logLik 468.5742 eps 0.0086 iteration 102 qsel 5
#> userIC 2367.327 logLik 468.9557 eps 0.0120 iteration 103 qsel 5
#> userIC 2367.307 logLik 469.2909 eps 0.0072 iteration 104 qsel 5
#> userIC 2367.285 logLik 469.5800 eps 0.0056 iteration 105 qsel 5
#> userIC 2367.266 logLik 469.8162 eps 0.0048 iteration 106 qsel 5
#> userIC 2367.247 logLik 470.0114 eps 0.0046 iteration 107 qsel 5
#> userIC 2367.232 logLik 478.1570 eps 0.0060 iteration 108 qsel 5
#> userIC 2367.113 logLik 478.3300 eps 0.0316 iteration 109 qsel 5
#> userIC 2367.097 logLik 484.0327 eps 0.0050 iteration 110 qsel 5
#> userIC 2366.983 logLik 484.1843 eps 0.0473 iteration 111 qsel 5
#> userIC 2366.969 logLik 488.1785 eps 0.0041 iteration 112 qsel 5
#> userIC 2366.881 logLik 490.9084 eps 0.0033 iteration 113 qsel 5
#> userIC 2366.768 logLik 491.0446 eps 0.0092 iteration 114 qsel 5
#> userIC 2366.753 logLik 492.9112 eps 0.0027 iteration 115 qsel 5
#> userIC 2366.646 logLik 493.0268 eps 0.0043 iteration 116 qsel 5
#> userIC 2366.639 logLik 494.2997 eps 0.0022 iteration 117 qsel 5
#> userIC 2366.535 logLik 495.2359 eps 0.0018 iteration 118 qsel 5
#> userIC 2366.463 logLik 495.8433 eps 0.0014 iteration 119 qsel 5
#> userIC 2366.368 logLik 496.2647 eps 0.0011 iteration 120 qsel 5
#> userIC 2366.280 logLik 496.6283 eps 0.0009 iteration 121 qsel 5
#> userIC 2366.211 logLik 496.8505 eps 0.0007 iteration 122 qsel 5
#> userIC 2366.126 logLik 497.0148 eps 0.0006 iteration 123 qsel 5
#> userIC 2366.048 logLik 497.1397 eps 0.0005 iteration 124 qsel 5
#> userIC 2365.974 logLik 497.3073 eps 0.0005 iteration 125 qsel 5
#> userIC 2365.916 logLik 497.3918 eps 0.0003 iteration 126 qsel 5
#> userIC 2365.846 logLik 497.4619 eps 0.0003 iteration 127 qsel 5
#> userIC 2365.783 logLik 497.5215 eps 0.0002 iteration 128 qsel 5
#> userIC 2365.724 logLik 497.6429 eps 0.0004 iteration 129 qsel 5
#> userIC 2365.673 logLik 497.6923 eps 0.0002 iteration 130 qsel 5
#> userIC 2365.617 logLik 497.7365 eps 0.0001 iteration 131 qsel 5
#> userIC 2365.565 logLik 497.7763 eps 0.0001 iteration 132 qsel 5
#> userIC 2365.519 logLik 497.8126 eps 0.0001 iteration 133 qsel 5
#> userIC 2365.474 logLik 497.9173 eps 0.0003 iteration 134 qsel 5
#> userIC 2365.432 logLik 497.9503 eps 0.0001 iteration 135 qsel 5
#> userIC 2365.390 logLik 497.9755 eps 0.0106 iteration 136 qsel 5
#> userIC 2365.380 logLik 498.0168 eps 0.0233 iteration 137 qsel 5
#> 
#>  elapsed time:  3.34sec


## Plot curve of negative out-of-sample log-likelihood.
## boost_plot(b, which = "user")