Random Effects P-Spline

This smooth constructor implements the random effects representation of a P-spline.

# S3 method for sr.smooth.spec
smooth.construct(object, data, knots, ...)

Arguments

object, data, knots

See smooth.construct.

...

Currently not used.

Value

See smooth.construct

See also

bamlss, predict.bamlss, opt_bfit, opt_boost

Examples

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

## Estimate model.
f <- num ~ x1 + x2 + x3 + s2(x1,bs="sr") + s2(x2,bs="sr") + s2(x3,bs="sr")

b <- bamlss(f, data = d, optimizer = boost, sampler = FALSE)
#> 
#> logLik -241.891 eps 0.0867 iteration   2 qsel 1
#> logLik -231.740 eps 0.1021 iteration   3 qsel 1
#> logLik -223.509 eps 0.7382 iteration   4 qsel 1
#> logLik -216.833 eps 0.3125 iteration   5 qsel 1
#> logLik -213.573 eps 0.0977 iteration   6 qsel 2
#> logLik -208.175 eps 0.2492 iteration   7 qsel 2
#> logLik -202.847 eps 0.2142 iteration   8 qsel 3
#> logLik -200.247 eps 0.1136 iteration   9 qsel 3
#> logLik -195.921 eps 0.3445 iteration  10 qsel 3
#> logLik -191.615 eps 0.1065 iteration  11 qsel 3
#> logLik -184.375 eps 0.1418 iteration  12 qsel 4
#> logLik -180.920 eps 0.8809 iteration  13 qsel 4
#> logLik -177.444 eps 0.0765 iteration  14 qsel 4
#> logLik -171.569 eps 0.1322 iteration  15 qsel 4
#> logLik -168.371 eps 0.0108 iteration  16 qsel 4
#> logLik -165.499 eps 0.0888 iteration  17 qsel 4
#> logLik -162.571 eps 0.1337 iteration  18 qsel 4
#> logLik -157.600 eps 0.3013 iteration  19 qsel 4
#> logLik -154.054 eps 0.0111 iteration  20 qsel 4
#> logLik -151.126 eps 0.0098 iteration  21 qsel 4
#> logLik -148.638 eps 0.0583 iteration  22 qsel 4
#> logLik -146.040 eps 0.0090 iteration  23 qsel 4
#> logLik -143.371 eps 0.0728 iteration  24 qsel 4
#> logLik -141.047 eps 0.0084 iteration  25 qsel 4
#> logLik -139.135 eps 0.0074 iteration  26 qsel 4
#> logLik -137.566 eps 0.0066 iteration  27 qsel 4
#> logLik -136.017 eps 0.0576 iteration  28 qsel 4
#> logLik -131.264 eps 0.1594 iteration  29 qsel 4
#> logLik -129.012 eps 0.0559 iteration  30 qsel 4
#> logLik -126.644 eps 0.0527 iteration  31 qsel 4
#> logLik -122.765 eps 0.0857 iteration  32 qsel 4
#> logLik -120.602 eps 0.0077 iteration  33 qsel 4
#> logLik -118.742 eps 0.0342 iteration  34 qsel 4
#> logLik -116.744 eps 0.0494 iteration  35 qsel 4
#> logLik -113.464 eps 0.0978 iteration  36 qsel 4
#> logLik -111.221 eps 0.0077 iteration  37 qsel 4
#> logLik -109.377 eps 0.0068 iteration  38 qsel 4
#> logLik -107.791 eps 0.0886 iteration  39 qsel 4
#> logLik -106.046 eps 0.0385 iteration  40 qsel 4
#> logLik -103.175 eps 0.0389 iteration  41 qsel 4
#> logLik -101.295 eps 0.0068 iteration  42 qsel 4
#> logLik -99.7507 eps 0.0061 iteration  43 qsel 4
#> logLik -98.4061 eps 0.0231 iteration  44 qsel 4
#> logLik -96.8818 eps 0.0362 iteration  45 qsel 4
#> logLik -94.3702 eps 0.0415 iteration  46 qsel 4
#> logLik -92.8183 eps 0.0060 iteration  47 qsel 4
#> logLik -91.5456 eps 0.0053 iteration  48 qsel 4
#> logLik -90.4105 eps 0.0206 iteration  49 qsel 4
#> logLik -89.0771 eps 0.0304 iteration  50 qsel 4
#> logLik -86.8777 eps 0.0758 iteration  51 qsel 4
#> logLik -85.6143 eps 0.0053 iteration  52 qsel 4
#> logLik -84.5795 eps 0.0047 iteration  53 qsel 4
#> logLik -83.6250 eps 0.0188 iteration  54 qsel 4
#> logLik -82.4547 eps 0.0298 iteration  55 qsel 4
#> logLik -81.5210 eps 0.0044 iteration  56 qsel 4
#> logLik -80.7574 eps 0.0039 iteration  57 qsel 4
#> logLik -80.1334 eps 0.0035 iteration  58 qsel 4
#> logLik -79.5990 eps 0.0144 iteration  59 qsel 4
#> logLik -77.6340 eps 0.0810 iteration  60 qsel 4
#> logLik -76.8249 eps 0.0228 iteration  61 qsel 4
#> logLik -75.7784 eps 0.0334 iteration  62 qsel 4
#> logLik -74.1096 eps 0.0647 iteration  63 qsel 4
#> logLik -73.4008 eps 0.0037 iteration  64 qsel 4
#> logLik -72.7370 eps 0.0568 iteration  65 qsel 4
#> logLik -71.8209 eps 0.0263 iteration  66 qsel 4
#> logLik -71.1897 eps 0.0035 iteration  67 qsel 4
#> logLik -70.6744 eps 0.0031 iteration  68 qsel 4
#> logLik -70.2540 eps 0.0028 iteration  69 qsel 4
#> logLik -69.9019 eps 0.0181 iteration  70 qsel 4
#> logLik -69.5489 eps 0.0134 iteration  71 qsel 4
#> logLik -69.1888 eps 0.0025 iteration  72 qsel 4
#> logLik -68.8422 eps 0.0253 iteration  73 qsel 4
#> logLik -68.4957 eps 0.0131 iteration  74 qsel 4
#> logLik -68.1589 eps 0.0336 iteration  75 qsel 4
#> logLik -67.8230 eps 0.0153 iteration  76 qsel 4
#> logLik -67.4957 eps 0.0233 iteration  77 qsel 4
#> logLik -67.1630 eps 0.0024 iteration  78 qsel 4
#> logLik -66.8333 eps 0.0156 iteration  79 qsel 4
#> logLik -66.5115 eps 0.0194 iteration  80 qsel 4
#> logLik -66.1918 eps 0.0157 iteration  81 qsel 4
#> logLik -65.6199 eps 0.3729 iteration  82 qsel 4
#> logLik -65.3073 eps 0.0157 iteration  83 qsel 4
#> logLik -64.9948 eps 0.0023 iteration  84 qsel 4
#> logLik -64.6828 eps 0.0103 iteration  85 qsel 4
#> logLik -63.8129 eps 0.0261 iteration  86 qsel 4
#> logLik -63.5059 eps 0.0242 iteration  87 qsel 4
#> logLik -63.2062 eps 0.0101 iteration  88 qsel 4
#> logLik -62.9085 eps 0.0150 iteration  89 qsel 4
#> logLik -62.6086 eps 0.0023 iteration  90 qsel 4
#> logLik -62.3141 eps 0.0128 iteration  91 qsel 4
#> logLik -62.0220 eps 0.0139 iteration  92 qsel 4
#> logLik -61.7360 eps 0.0116 iteration  93 qsel 4
#> logLik -61.4526 eps 0.0143 iteration  94 qsel 4
#> logLik -61.1749 eps 0.0104 iteration  95 qsel 4
#> logLik -60.9000 eps 0.0781 iteration  96 qsel 4
#> logLik -60.6197 eps 0.0022 iteration  97 qsel 4
#> logLik -60.3470 eps 0.0118 iteration  98 qsel 4
#> logLik -60.0773 eps 0.0158 iteration  99 qsel 4
#> logLik -59.8126 eps 0.0276 iteration 100 qsel 4
#> logLik -59.5509 eps 0.0204 iteration 101 qsel 4
#> logLik -59.2939 eps 0.0133 iteration 102 qsel 4
#> logLik -58.8219 eps 0.0152 iteration 103 qsel 4
#> logLik -58.5569 eps 0.0021 iteration 104 qsel 4
#> logLik -58.3002 eps 0.0258 iteration 105 qsel 4
#> logLik -57.4907 eps 0.0226 iteration 106 qsel 4
#> logLik -57.2426 eps 0.0090 iteration 107 qsel 4
#> logLik -56.9938 eps 0.0143 iteration 108 qsel 4
#> logLik -56.7472 eps 0.0020 iteration 109 qsel 4
#> logLik -56.5035 eps 0.0088 iteration 110 qsel 4
#> logLik -56.2598 eps 0.0156 iteration 111 qsel 4
#> logLik -56.0230 eps 0.0135 iteration 112 qsel 4
#> logLik -55.7868 eps 0.0122 iteration 113 qsel 4
#> logLik -55.5566 eps 0.0105 iteration 114 qsel 4
#> logLik -55.3276 eps 0.0119 iteration 115 qsel 4
#> logLik -55.0999 eps 0.0019 iteration 116 qsel 4
#> logLik -54.8741 eps 0.0699 iteration 117 qsel 4
#> logLik -54.6497 eps 0.0124 iteration 118 qsel 4
#> logLik -54.4304 eps 0.0105 iteration 119 qsel 4
#> logLik -54.2128 eps 0.0142 iteration 120 qsel 4
#> logLik -53.9998 eps 0.0086 iteration 121 qsel 4
#> logLik -53.7887 eps 0.0233 iteration 122 qsel 4
#> logLik -53.5796 eps 0.0018 iteration 123 qsel 4
#> logLik -53.3707 eps 0.0083 iteration 124 qsel 4
#> logLik -52.5901 eps 0.0206 iteration 125 qsel 4
#> logLik -52.3834 eps 0.0205 iteration 126 qsel 4
#> logLik -51.9919 eps 0.0173 iteration 127 qsel 4
#> logLik -51.7925 eps 0.0119 iteration 128 qsel 4
#> logLik -51.5908 eps 0.0018 iteration 129 qsel 4
#> logLik -51.3886 eps 0.0128 iteration 130 qsel 4
#> logLik -51.1929 eps 0.0208 iteration 131 qsel 4
#> logLik -50.9970 eps 0.0122 iteration 132 qsel 4
#> logLik -50.8068 eps 0.0125 iteration 133 qsel 4
#> logLik -50.6168 eps 0.0130 iteration 134 qsel 4
#> logLik -50.4320 eps 0.0087 iteration 135 qsel 4
#> logLik -50.2478 eps 0.0178 iteration 136 qsel 4
#> logLik -50.0608 eps 0.0017 iteration 137 qsel 4
#> logLik -49.8796 eps 0.0080 iteration 138 qsel 4
#> logLik -49.6994 eps 0.0454 iteration 139 qsel 4
#> logLik -49.5234 eps 0.0125 iteration 140 qsel 4
#> logLik -49.3486 eps 0.0260 iteration 141 qsel 4
#> logLik -49.1776 eps 0.0075 iteration 142 qsel 4
#> logLik -48.3963 eps 0.0279 iteration 143 qsel 4
#> logLik -48.2165 eps 0.0017 iteration 144 qsel 4
#> logLik -48.0455 eps 0.0156 iteration 145 qsel 4
#> logLik -47.7188 eps 0.1629 iteration 146 qsel 4
#> logLik -47.5541 eps 0.0103 iteration 147 qsel 4
#> logLik -47.3884 eps 0.0154 iteration 148 qsel 4
#> logLik -47.2282 eps 0.0081 iteration 149 qsel 4
#> logLik -47.0662 eps 0.0016 iteration 150 qsel 4
#> logLik -46.9043 eps 0.0205 iteration 151 qsel 4
#> logLik -46.7473 eps 0.0076 iteration 152 qsel 4
#> logLik -46.5904 eps 0.0305 iteration 153 qsel 4
#> logLik -46.4379 eps 0.0190 iteration 154 qsel 4
#> logLik -46.2857 eps 0.0144 iteration 155 qsel 4
#> logLik -46.1375 eps 0.0070 iteration 156 qsel 4
#> logLik -45.9900 eps 0.0144 iteration 157 qsel 4
#> logLik -45.8419 eps 0.0015 iteration 158 qsel 4
#> logLik -45.6967 eps 0.0066 iteration 159 qsel 4
#> logLik -45.5525 eps 0.0146 iteration 160 qsel 4
#> logLik -44.7181 eps 0.0276 iteration 161 qsel 4
#> logLik -44.5785 eps 0.0240 iteration 162 qsel 4
#> logLik -44.4388 eps 0.0125 iteration 163 qsel 4
#> logLik -44.2978 eps 0.0015 iteration 164 qsel 4
#> logLik -44.1610 eps 0.0098 iteration 165 qsel 4
#> logLik -44.0245 eps 0.0127 iteration 166 qsel 4
#> logLik -43.7507 eps 0.0144 iteration 167 qsel 4
#> logLik -43.6187 eps 0.0084 iteration 168 qsel 4
#> logLik -43.4864 eps 0.0150 iteration 169 qsel 4
#> logLik -43.3581 eps 0.0180 iteration 170 qsel 4
#> logLik -43.2299 eps 0.0132 iteration 171 qsel 4
#> logLik -43.1002 eps 0.0014 iteration 172 qsel 4
#> logLik -42.9745 eps 0.0099 iteration 173 qsel 4
#> logLik -42.8493 eps 0.0125 iteration 174 qsel 4
#> logLik -42.7272 eps 0.0071 iteration 175 qsel 4
#> logLik -42.6058 eps 0.0147 iteration 176 qsel 4
#> logLik -42.4872 eps 0.0062 iteration 177 qsel 4
#> logLik -42.3694 eps 0.0223 iteration 178 qsel 4
#> logLik -41.4220 eps 0.0504 iteration 179 qsel 4
#> logLik -41.2933 eps 0.0014 iteration 180 qsel 4
#> logLik -41.1783 eps 0.0050 iteration 181 qsel 4
#> logLik -41.0632 eps 0.0086 iteration 182 qsel 4
#> logLik -40.9514 eps 0.0047 iteration 183 qsel 4
#> logLik -40.8399 eps 0.0080 iteration 184 qsel 4
#> logLik -40.7293 eps 0.0013 iteration 185 qsel 4
#> logLik -40.6199 eps 0.0046 iteration 186 qsel 4
#> logLik -40.5110 eps 0.0082 iteration 187 qsel 4
#> logLik -40.2788 eps 0.0092 iteration 188 qsel 4
#> logLik -40.1733 eps 0.0046 iteration 189 qsel 4
#> logLik -40.0678 eps 0.0086 iteration 190 qsel 4
#> logLik -39.9652 eps 0.0047 iteration 191 qsel 4
#> logLik -39.8629 eps 0.0157 iteration 192 qsel 4
#> logLik -39.7623 eps 0.0012 iteration 193 qsel 4
#> logLik -39.6620 eps 0.0065 iteration 194 qsel 4
#> logLik -39.5621 eps 0.0089 iteration 195 qsel 4
#> logLik -38.4116 eps 0.0254 iteration 196 qsel 4
#> logLik -38.3135 eps 0.0012 iteration 197 qsel 4
#> logLik -38.2159 eps 0.0073 iteration 198 qsel 4
#> logLik -38.1185 eps 0.0246 iteration 199 qsel 4
#> logLik -38.0237 eps 0.0053 iteration 200 qsel 4
#> logLik -37.9293 eps 0.0071 iteration 201 qsel 4
#> logLik -37.8372 eps 0.0041 iteration 202 qsel 4
#> logLik -37.7456 eps 0.0062 iteration 203 qsel 4
#> logLik -37.6562 eps 0.0038 iteration 204 qsel 4
#> logLik -37.5673 eps 0.0058 iteration 205 qsel 4
#> logLik -37.4789 eps 0.0011 iteration 206 qsel 4
#> logLik -37.3915 eps 0.0036 iteration 207 qsel 4
#> logLik -37.3048 eps 0.0055 iteration 208 qsel 4
#> logLik -37.1043 eps 0.0065 iteration 209 qsel 4
#> logLik -37.0201 eps 0.0036 iteration 210 qsel 4
#> logLik -36.9360 eps 0.0054 iteration 211 qsel 4
#> logLik -35.8320 eps 0.0154 iteration 212 qsel 4
#> logLik -35.7413 eps 0.0011 iteration 213 qsel 4
#> logLik -35.6589 eps 0.0036 iteration 214 qsel 4
#> logLik -35.5769 eps 0.0054 iteration 215 qsel 4
#> logLik -35.4969 eps 0.0035 iteration 216 qsel 4
#> logLik -35.4174 eps 0.0053 iteration 217 qsel 4
#> logLik -35.3397 eps 0.0034 iteration 218 qsel 4
#> logLik -35.2616 eps 0.0010 iteration 219 qsel 4
#> logLik -35.1840 eps 0.0051 iteration 220 qsel 4
#> logLik -35.1081 eps 0.0034 iteration 221 qsel 4
#> logLik -35.0329 eps 0.0051 iteration 222 qsel 4
#> logLik -34.9591 eps 0.0033 iteration 223 qsel 4
#> logLik -34.8861 eps 0.0050 iteration 224 qsel 4
#> logLik -33.9006 eps 0.0153 iteration 225 qsel 4
#> logLik -33.8227 eps 0.0010 iteration 226 qsel 4
#> logLik -33.7506 eps 0.0046 iteration 227 qsel 4
#> logLik -33.6729 eps 0.0043 iteration 228 qsel 5
#> logLik -32.8726 eps 0.0121 iteration 229 qsel 5
#> logLik -32.7988 eps 0.0096 iteration 230 qsel 5
#> logLik -32.7248 eps 0.0010 iteration 231 qsel 5
#> logLik -32.6511 eps 0.0064 iteration 232 qsel 5
#> logLik -32.5808 eps 0.0042 iteration 233 qsel 5
#> logLik -32.5101 eps 0.0075 iteration 234 qsel 5
#> logLik -32.4417 eps 0.0049 iteration 235 qsel 5
#> logLik -32.3733 eps 0.0062 iteration 236 qsel 5
#> logLik -32.1974 eps 0.0077 iteration 237 qsel 5
#> logLik -32.1315 eps 0.0036 iteration 238 qsel 5
#> logLik -32.0652 eps 0.0075 iteration 239 qsel 5
#> logLik -31.1547 eps 0.0156 iteration 240 qsel 5
#> logLik -31.0781 eps 0.0010 iteration 241 qsel 5
#> logLik -31.0135 eps 0.0037 iteration 242 qsel 5
#> logLik -30.9398 eps 0.0042 iteration 243 qsel 5
#> logLik -30.0843 eps 0.0289 iteration 244 qsel 5
#> logLik -30.0117 eps 0.0010 iteration 245 qsel 5
#> logLik -29.9438 eps 0.0060 iteration 246 qsel 5
#> logLik -29.8744 eps 0.0052 iteration 247 qsel 5
#> logLik -29.8098 eps 0.0042 iteration 248 qsel 5
#> logLik -29.7443 eps 0.0041 iteration 249 qsel 5
#> logLik -29.6828 eps 0.0009 iteration 250 qsel 5
#> logLik -29.6211 eps 0.0036 iteration 251 qsel 5
#> logLik -29.5575 eps 0.0074 iteration 252 qsel 5
#> logLik -28.6986 eps 0.0241 iteration 253 qsel 5
#> logLik -28.6379 eps 0.0047 iteration 254 qsel 5
#> logLik -28.5683 eps 0.0040 iteration 255 qsel 5
#> logLik -28.5079 eps 0.0009 iteration 256 qsel 5
#> logLik -28.4494 eps 0.0042 iteration 257 qsel 5
#> logLik -28.3838 eps 0.0038 iteration 258 qsel 5
#> logLik -28.3278 eps 0.0040 iteration 259 qsel 5
#> logLik -28.2663 eps 0.0038 iteration 260 qsel 5
#> logLik -27.4755 eps 0.0151 iteration 261 qsel 5
#> logLik -27.4131 eps 0.0038 iteration 262 qsel 5
#> logLik -27.3539 eps 0.0009 iteration 263 qsel 5
#> logLik -27.2951 eps 0.0043 iteration 264 qsel 5
#> logLik -27.2359 eps 0.0037 iteration 265 qsel 5
#> logLik -27.1801 eps 0.0095 iteration 266 qsel 5
#> logLik -27.1241 eps 0.0045 iteration 267 qsel 5
#> logLik -26.3871 eps 0.0181 iteration 268 qsel 5
#> logLik -26.3271 eps 0.0039 iteration 269 qsel 5
#> logLik -26.1659 eps 0.0110 iteration 270 qsel 5
#> logLik -26.1068 eps 0.0009 iteration 271 qsel 5
#> logLik -26.0536 eps 0.0057 iteration 272 qsel 5
#> logLik -25.9974 eps 0.0039 iteration 273 qsel 5
#> logLik -25.9444 eps 0.0120 iteration 274 qsel 5
#> logLik -25.8900 eps 0.0042 iteration 275 qsel 5
#> logLik -25.1474 eps 0.0197 iteration 276 qsel 5
#> logLik -25.0883 eps 0.0043 iteration 277 qsel 5
#> logLik -25.0310 eps 0.0009 iteration 278 qsel 5
#> logLik -24.9787 eps 0.0033 iteration 279 qsel 5
#> logLik -24.9228 eps 0.0040 iteration 280 qsel 5
#> logLik -24.8729 eps 0.0032 iteration 281 qsel 5
#> logLik -24.8203 eps 0.0038 iteration 282 qsel 5
#> logLik -24.1428 eps 0.0386 iteration 283 qsel 5
#> logLik -24.0883 eps 0.0008 iteration 284 qsel 5
#> logLik -24.0348 eps 0.0069 iteration 285 qsel 5
#> logLik -23.9841 eps 0.0088 iteration 286 qsel 5
#> logLik -23.9323 eps 0.0048 iteration 287 qsel 5
#> logLik -23.8830 eps 0.0152 iteration 288 qsel 5
#> logLik -23.8329 eps 0.0044 iteration 289 qsel 5
#> logLik -23.7849 eps 0.0123 iteration 290 qsel 5
#> logLik -23.7364 eps 0.0046 iteration 291 qsel 5
#> logLik -23.0485 eps 0.0178 iteration 292 qsel 5
#> logLik -22.9935 eps 0.0045 iteration 293 qsel 5
#> logLik -22.9396 eps 0.0008 iteration 294 qsel 5
#> logLik -22.8927 eps 0.0098 iteration 295 qsel 5
#> logLik -22.8408 eps 0.0056 iteration 296 qsel 5
#> logLik -22.7960 eps 0.0047 iteration 297 qsel 5
#> logLik -22.7473 eps 0.0184 iteration 298 qsel 5
#> logLik -22.5918 eps 0.0146 iteration 299 qsel 5
#> logLik -21.9609 eps 0.0115 iteration 300 qsel 5
#> logLik -21.9085 eps 0.0008 iteration 301 qsel 5
#> logLik -21.8612 eps 0.0054 iteration 302 qsel 5
#> logLik -21.8117 eps 0.0035 iteration 303 qsel 5
#> logLik -21.7666 eps 0.0054 iteration 304 qsel 5
#> logLik -21.7199 eps 0.0035 iteration 305 qsel 5
#> logLik -21.6758 eps 0.0008 iteration 306 qsel 5
#> logLik -21.6328 eps 0.0056 iteration 307 qsel 5
#> logLik -21.5885 eps 0.0035 iteration 308 qsel 5
#> logLik -21.0911 eps 0.0175 iteration 309 qsel 5
#> logLik -21.0495 eps 0.0068 iteration 310 qsel 5
#> logLik -21.0025 eps 0.0053 iteration 311 qsel 5
#> logLik -20.4105 eps 0.0131 iteration 312 qsel 5
#> logLik -20.3644 eps 0.0050 iteration 313 qsel 5
#> logLik -20.3179 eps 0.0052 iteration 314 qsel 5
#> logLik -20.2712 eps 0.0008 iteration 315 qsel 5
#> logLik -20.2272 eps 0.0051 iteration 316 qsel 5
#> logLik -20.1831 eps 0.0054 iteration 317 qsel 5
#> logLik -20.1414 eps 0.0054 iteration 318 qsel 5
#> logLik -20.0996 eps 0.0058 iteration 319 qsel 5
#> logLik -20.0600 eps 0.0129 iteration 320 qsel 5
#> logLik -20.0204 eps 0.0007 iteration 321 qsel 5
#> logLik -19.9808 eps 0.0065 iteration 322 qsel 5
#> logLik -19.9431 eps 0.0055 iteration 323 qsel 5
#> logLik -19.3823 eps 0.1204 iteration 324 qsel 5
#> logLik -19.3439 eps 0.0130 iteration 325 qsel 5
#> logLik -19.3006 eps 0.0053 iteration 326 qsel 5
#> logLik -18.8673 eps 0.0130 iteration 327 qsel 5
#> logLik -18.8265 eps 0.0007 iteration 328 qsel 5
#> logLik -18.7859 eps 0.0053 iteration 329 qsel 5
#> logLik -18.7433 eps 0.0044 iteration 330 qsel 5
#> logLik -18.7046 eps 0.0046 iteration 331 qsel 5
#> logLik -18.6644 eps 0.0040 iteration 332 qsel 5
#> logLik -18.6276 eps 0.0042 iteration 333 qsel 5
#> logLik -18.5897 eps 0.0039 iteration 334 qsel 5
#> logLik -18.4271 eps 0.0082 iteration 335 qsel 5
#> logLik -17.9098 eps 0.0120 iteration 336 qsel 5
#> logLik -17.8696 eps 0.0007 iteration 337 qsel 5
#> logLik -17.8313 eps 0.0036 iteration 338 qsel 5
#> logLik -17.7908 eps 0.0029 iteration 339 qsel 5
#> logLik -17.7543 eps 0.0034 iteration 340 qsel 5
#> logLik -17.7162 eps 0.0027 iteration 341 qsel 5
#> logLik -17.6814 eps 0.0032 iteration 342 qsel 5
#> logLik -17.6454 eps 0.0027 iteration 343 qsel 5
#> logLik -17.2704 eps 0.0170 iteration 344 qsel 5
#> logLik -17.2334 eps 0.0007 iteration 345 qsel 5
#> logLik -17.1995 eps 0.0034 iteration 346 qsel 5
#> logLik -17.1613 eps 0.0039 iteration 347 qsel 5
#> logLik -17.1289 eps 0.0036 iteration 348 qsel 5
#> logLik -17.0924 eps 0.0048 iteration 349 qsel 5
#> logLik -16.6101 eps 0.0098 iteration 350 qsel 5
#> logLik -16.5731 eps 0.0037 iteration 351 qsel 5
#> logLik -16.5359 eps 0.0030 iteration 352 qsel 5
#> logLik -16.5007 eps 0.0035 iteration 353 qsel 5
#> logLik -16.4654 eps 0.0007 iteration 354 qsel 5
#> logLik -16.4300 eps 0.0029 iteration 355 qsel 5
#> logLik -16.3965 eps 0.0034 iteration 356 qsel 5
#> logLik -16.3631 eps 0.0027 iteration 357 qsel 5
#> logLik -16.3313 eps 0.0035 iteration 358 qsel 5
#> logLik -16.2997 eps 0.0026 iteration 359 qsel 5
#> logLik -16.2695 eps 0.0040 iteration 360 qsel 5
#> logLik -15.8175 eps 0.0095 iteration 361 qsel 5
#> logLik -15.7843 eps 0.0006 iteration 362 qsel 5
#> logLik -15.7531 eps 0.0033 iteration 363 qsel 5
#> logLik -15.7183 eps 0.0027 iteration 364 qsel 5
#> logLik -15.4067 eps 0.0873 iteration 365 qsel 5
#> logLik -15.3740 eps 0.0048 iteration 366 qsel 5
#> logLik -15.3395 eps 0.0073 iteration 367 qsel 5
#> logLik -15.3084 eps 0.0049 iteration 368 qsel 5
#> logLik -15.2758 eps 0.0063 iteration 369 qsel 5
#> logLik -15.2455 eps 0.0006 iteration 370 qsel 5
#> logLik -15.2157 eps 0.0054 iteration 371 qsel 5
#> logLik -15.1850 eps 0.0119 iteration 372 qsel 5
#> logLik -14.9898 eps 0.0350 iteration 373 qsel 5
#> logLik -14.5695 eps 0.0104 iteration 374 qsel 5
#> logLik -14.5382 eps 0.0042 iteration 375 qsel 5
#> logLik -14.5058 eps 0.0034 iteration 376 qsel 5
#> logLik -14.4760 eps 0.0190 iteration 377 qsel 5
#> logLik -14.4454 eps 0.0037 iteration 378 qsel 5
#> logLik -14.4155 eps 0.0006 iteration 379 qsel 5
#> logLik -14.3871 eps 0.0045 iteration 380 qsel 5
#> logLik -14.3581 eps 0.0032 iteration 381 qsel 5
#> logLik -14.0881 eps 0.0096 iteration 382 qsel 5
#> logLik -14.0607 eps 0.0025 iteration 383 qsel 5
#> logLik -14.0297 eps 0.0044 iteration 384 qsel 5
#> logLik -14.0030 eps 0.0006 iteration 385 qsel 5
#> logLik -13.9766 eps 0.0025 iteration 386 qsel 5
#> logLik -13.9474 eps 0.0121 iteration 387 qsel 5
#> logLik -13.5501 eps 0.0108 iteration 388 qsel 5
#> logLik -13.5199 eps 0.0036 iteration 389 qsel 5
#> logLik -13.4903 eps 0.0025 iteration 390 qsel 5
#> logLik -13.4617 eps 0.0040 iteration 391 qsel 5
#> logLik -13.4337 eps 0.0025 iteration 392 qsel 5
#> logLik -13.4066 eps 0.0338 iteration 393 qsel 5
#> logLik -13.3801 eps 0.0041 iteration 394 qsel 5
#> logLik -13.3544 eps 0.0037 iteration 395 qsel 5
#> logLik -13.3288 eps 0.0006 iteration 396 qsel 5
#> logLik -13.3036 eps 0.0024 iteration 397 qsel 5
#> logLik -13.0673 eps 0.0087 iteration 398 qsel 5
#> logLik -12.7780 eps 0.0093 iteration 399 qsel 5
#> logLik -12.4116 eps 0.0109 iteration 400 qsel 5
#> 
#>  elapsed time:  0.88sec


plot(b)