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

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 -169.168 eps 0.0479 iteration   2 qsel 1
#> logLik -162.007 eps 0.0457 iteration   3 qsel 1
#> logLik -156.202 eps 0.0450 iteration   4 qsel 1
#> logLik -151.496 eps 0.0461 iteration   5 qsel 1
#> logLik -144.078 eps 0.0537 iteration   6 qsel 2
#> logLik -140.232 eps 0.0555 iteration   7 qsel 2
#> logLik -137.908 eps 0.0284 iteration   8 qsel 3
#> logLik -131.939 eps 0.3467 iteration   9 qsel 3
#> logLik -128.795 eps 0.0972 iteration  10 qsel 3
#> logLik -125.232 eps 0.0728 iteration  11 qsel 4
#> logLik -120.197 eps 0.1062 iteration  12 qsel 4
#> logLik -117.624 eps 0.2880 iteration  13 qsel 4
#> logLik -114.664 eps 0.0678 iteration  14 qsel 4
#> logLik -112.894 eps 0.0067 iteration  15 qsel 4
#> logLik -111.426 eps 0.0425 iteration  16 qsel 4
#> logLik -107.242 eps 0.1577 iteration  17 qsel 4
#> logLik -105.074 eps 0.0919 iteration  18 qsel 4
#> logLik -102.532 eps 0.0945 iteration  19 qsel 4
#> logLik -99.0027 eps 0.1529 iteration  20 qsel 4
#> logLik -96.6826 eps 0.0077 iteration  21 qsel 4
#> logLik -94.7744 eps 0.0068 iteration  22 qsel 4
#> logLik -92.8825 eps 0.0672 iteration  23 qsel 4
#> logLik -90.6119 eps 0.0909 iteration  24 qsel 4
#> logLik -87.4351 eps 0.1025 iteration  25 qsel 4
#> logLik -85.4212 eps 0.0069 iteration  26 qsel 4
#> logLik -83.7666 eps 0.0061 iteration  27 qsel 4
#> logLik -82.1227 eps 0.1426 iteration  28 qsel 4
#> logLik -80.0954 eps 0.0486 iteration  29 qsel 4
#> logLik -78.5392 eps 0.0059 iteration  30 qsel 4
#> logLik -77.2630 eps 0.0052 iteration  31 qsel 4
#> logLik -76.2175 eps 0.0046 iteration  32 qsel 4
#> logLik -75.2996 eps 0.0249 iteration  33 qsel 4
#> logLik -73.8612 eps 0.0409 iteration  34 qsel 4
#> logLik -70.9064 eps 0.0603 iteration  35 qsel 4
#> logLik -69.0600 eps 0.0574 iteration  36 qsel 4
#> logLik -67.8807 eps 0.2180 iteration  37 qsel 4
#> logLik -65.3759 eps 0.0770 iteration  38 qsel 4
#> logLik -63.8059 eps 0.0564 iteration  39 qsel 4
#> logLik -62.3700 eps 0.0054 iteration  40 qsel 4
#> logLik -61.1929 eps 0.0048 iteration  41 qsel 4
#> logLik -60.1763 eps 0.0640 iteration  42 qsel 4
#> logLik -57.9354 eps 0.0892 iteration  43 qsel 4
#> logLik -56.5237 eps 0.0730 iteration  44 qsel 4
#> logLik -55.3482 eps 0.0048 iteration  45 qsel 4
#> logLik -54.3857 eps 0.0043 iteration  46 qsel 4
#> logLik -53.5131 eps 0.0612 iteration  47 qsel 4
#> logLik -51.5059 eps 0.0810 iteration  48 qsel 4
#> logLik -50.2338 eps 0.1963 iteration  49 qsel 4
#> logLik -49.2788 eps 0.0042 iteration  50 qsel 4
#> logLik -48.4976 eps 0.0038 iteration  51 qsel 4
#> logLik -47.7515 eps 0.0555 iteration  52 qsel 4
#> logLik -45.9462 eps 0.0618 iteration  53 qsel 4
#> logLik -44.7961 eps 0.0412 iteration  54 qsel 4
#> logLik -44.0247 eps 0.0037 iteration  55 qsel 4
#> logLik -43.3945 eps 0.0033 iteration  56 qsel 4
#> logLik -42.8799 eps 0.0030 iteration  57 qsel 4
#> logLik -42.4602 eps 0.0026 iteration  58 qsel 4
#> logLik -42.0736 eps 0.0177 iteration  59 qsel 4
#> logLik -41.4188 eps 0.1432 iteration  60 qsel 4
#> logLik -39.7811 eps 0.0404 iteration  61 qsel 4
#> logLik -38.7051 eps 0.0327 iteration  62 qsel 4
#> logLik -38.1660 eps 0.0285 iteration  63 qsel 4
#> logLik -36.7064 eps 0.0306 iteration  64 qsel 4
#> logLik -35.7535 eps 0.0735 iteration  65 qsel 4
#> logLik -35.2269 eps 0.0030 iteration  66 qsel 4
#> logLik -34.7755 eps 0.0186 iteration  67 qsel 4
#> logLik -33.4233 eps 0.0288 iteration  68 qsel 4
#> logLik -32.5555 eps 0.0270 iteration  69 qsel 4
#> logLik -32.0491 eps 0.0029 iteration  70 qsel 4
#> logLik -31.6359 eps 0.0026 iteration  71 qsel 4
#> logLik -31.2520 eps 0.0160 iteration  72 qsel 4
#> logLik -30.9059 eps 0.0023 iteration  73 qsel 4
#> logLik -30.6239 eps 0.0021 iteration  74 qsel 4
#> logLik -30.3943 eps 0.0019 iteration  75 qsel 4
#> logLik -30.1779 eps 0.0167 iteration  76 qsel 4
#> logLik -29.9611 eps 0.0098 iteration  77 qsel 4
#> logLik -29.7495 eps 0.0315 iteration  78 qsel 4
#> logLik -29.5381 eps 0.0099 iteration  79 qsel 4
#> logLik -29.3312 eps 0.1370 iteration  80 qsel 4
#> logLik -29.1249 eps 0.0105 iteration  81 qsel 4
#> logLik -28.9149 eps 0.0018 iteration  82 qsel 4
#> logLik -28.7108 eps 0.0187 iteration  83 qsel 4
#> logLik -27.8605 eps 0.0325 iteration  84 qsel 4
#> logLik -27.6571 eps 0.0148 iteration  85 qsel 4
#> logLik -27.4576 eps 0.0175 iteration  86 qsel 4
#> logLik -27.2591 eps 0.0784 iteration  87 qsel 4
#> logLik -27.0591 eps 0.0017 iteration  88 qsel 4
#> logLik -26.8624 eps 0.0195 iteration  89 qsel 4
#> logLik -26.6668 eps 0.0130 iteration  90 qsel 4
#> logLik -26.4746 eps 0.0341 iteration  91 qsel 4
#> logLik -26.2838 eps 0.0117 iteration  92 qsel 4
#> logLik -26.0960 eps 0.0232 iteration  93 qsel 4
#> logLik -25.7633 eps 0.0161 iteration  94 qsel 4
#> logLik -25.5769 eps 0.0154 iteration  95 qsel 4
#> logLik -25.3889 eps 0.0017 iteration  96 qsel 4
#> logLik -25.2038 eps 0.0137 iteration  97 qsel 4
#> logLik -25.0201 eps 0.0127 iteration  98 qsel 4
#> logLik -24.8393 eps 0.0131 iteration  99 qsel 4
#> logLik -24.6600 eps 0.0125 iteration 100 qsel 4
#> logLik -24.4834 eps 0.0128 iteration 101 qsel 4
#> logLik -24.3084 eps 0.0130 iteration 102 qsel 4
#> logLik -24.1358 eps 0.0126 iteration 103 qsel 4
#> logLik -23.9622 eps 0.0016 iteration 104 qsel 4
#> logLik -23.7899 eps 0.0146 iteration 105 qsel 4
#> logLik -22.9805 eps 0.0333 iteration 106 qsel 4
#> logLik -22.8105 eps 0.0200 iteration 107 qsel 4
#> logLik -22.6420 eps 0.0139 iteration 108 qsel 4
#> logLik -22.4760 eps 0.0147 iteration 109 qsel 4
#> logLik -22.3112 eps 0.0015 iteration 110 qsel 4
#> logLik -22.1453 eps 0.0249 iteration 111 qsel 4
#> logLik -21.9817 eps 0.0135 iteration 112 qsel 4
#> logLik -21.8198 eps 0.0113 iteration 113 qsel 4
#> logLik -21.6600 eps 0.0134 iteration 114 qsel 4
#> logLik -21.5019 eps 0.0098 iteration 115 qsel 4
#> logLik -21.2202 eps 0.0228 iteration 116 qsel 4
#> logLik -21.0642 eps 0.0128 iteration 117 qsel 4
#> logLik -20.9097 eps 0.0087 iteration 118 qsel 4
#> logLik -20.7540 eps 0.0015 iteration 119 qsel 4
#> logLik -20.6004 eps 0.0128 iteration 120 qsel 4
#> logLik -20.4483 eps 0.0082 iteration 121 qsel 4
#> logLik -20.2982 eps 0.0138 iteration 122 qsel 4
#> logLik -20.1498 eps 0.0082 iteration 123 qsel 4
#> logLik -20.0032 eps 0.0184 iteration 124 qsel 4
#> logLik -19.8582 eps 0.0086 iteration 125 qsel 4
#> logLik -19.7151 eps 0.0304 iteration 126 qsel 4
#> logLik -18.9369 eps 0.0624 iteration 127 qsel 4
#> logLik -18.7826 eps 0.0015 iteration 128 qsel 4
#> logLik -18.6398 eps 0.0098 iteration 129 qsel 4
#> logLik -18.4988 eps 0.0632 iteration 130 qsel 4
#> logLik -18.3594 eps 0.0123 iteration 131 qsel 4
#> logLik -18.2217 eps 0.0172 iteration 132 qsel 4
#> logLik -18.0855 eps 0.0103 iteration 133 qsel 4
#> logLik -17.9499 eps 0.0014 iteration 134 qsel 4
#> logLik -17.8144 eps 0.0194 iteration 135 qsel 4
#> logLik -17.6804 eps 0.0089 iteration 136 qsel 4
#> logLik -17.5481 eps 0.0290 iteration 137 qsel 4
#> logLik -17.4172 eps 0.0094 iteration 138 qsel 4
#> logLik -17.1783 eps 0.0779 iteration 139 qsel 4
#> logLik -17.0491 eps 0.0165 iteration 140 qsel 4
#> logLik -16.9213 eps 0.0084 iteration 141 qsel 4
#> logLik -16.7951 eps 0.0223 iteration 142 qsel 4
#> logLik -16.6687 eps 0.0013 iteration 143 qsel 4
#> logLik -16.5430 eps 0.0087 iteration 144 qsel 4
#> logLik -16.4188 eps 0.0302 iteration 145 qsel 4
#> logLik -16.2960 eps 0.0140 iteration 146 qsel 4
#> logLik -16.1747 eps 0.0134 iteration 147 qsel 4
#> logLik -16.0548 eps 0.0111 iteration 148 qsel 4
#> logLik -15.2820 eps 0.0452 iteration 149 qsel 4
#> logLik -15.1605 eps 0.0013 iteration 150 qsel 4
#> logLik -15.0412 eps 0.0149 iteration 151 qsel 4
#> logLik -14.9233 eps 0.0096 iteration 152 qsel 4
#> logLik -14.8068 eps 0.0145 iteration 153 qsel 4
#> logLik -14.6916 eps 0.0102 iteration 154 qsel 4
#> logLik -14.5778 eps 0.0556 iteration 155 qsel 4
#> logLik -14.4653 eps 0.0153 iteration 156 qsel 4
#> logLik -14.3542 eps 0.0129 iteration 157 qsel 4
#> logLik -14.2444 eps 0.0190 iteration 158 qsel 4
#> logLik -14.1337 eps 0.0012 iteration 159 qsel 4
#> logLik -14.0244 eps 0.0117 iteration 160 qsel 4
#> logLik -13.9164 eps 0.0121 iteration 161 qsel 4
#> logLik -13.7134 eps 0.0149 iteration 162 qsel 4
#> logLik -13.6067 eps 0.0120 iteration 163 qsel 4
#> logLik -13.5011 eps 0.0931 iteration 164 qsel 4
#> logLik -13.3969 eps 0.0113 iteration 165 qsel 4
#> logLik -13.2939 eps 0.0110 iteration 166 qsel 4
#> logLik -13.1921 eps 0.0114 iteration 167 qsel 4
#> logLik -12.0647 eps 0.0567 iteration 168 qsel 4
#> logLik -11.2884 eps 0.0340 iteration 169 qsel 4
#> logLik -11.1606 eps 0.0013 iteration 170 qsel 4
#> logLik -11.0568 eps 0.0012 iteration 171 qsel 4
#> logLik -10.9641 eps 0.0075 iteration 172 qsel 4
#> logLik -10.8641 eps 0.0328 iteration 173 qsel 4
#> logLik -10.7728 eps 0.0069 iteration 174 qsel 4
#> logLik -10.6759 eps 0.0133 iteration 175 qsel 4
#> logLik -10.5862 eps 0.0068 iteration 176 qsel 4
#> logLik -10.4921 eps 0.0269 iteration 177 qsel 4
#> logLik -10.4009 eps 0.0011 iteration 178 qsel 4
#> logLik -10.3123 eps 0.0074 iteration 179 qsel 4
#> logLik -10.2202 eps 0.0462 iteration 180 qsel 4
#> logLik -10.0472 eps 0.0350 iteration 181 qsel 4
#> logLik  -9.9602 eps 0.0086 iteration 182 qsel 4
#> logLik  -9.8706 eps 0.0174 iteration 183 qsel 4
#> logLik  -9.7854 eps 0.0076 iteration 184 qsel 4
#> logLik  -9.6982 eps 0.0423 iteration 185 qsel 4
#> logLik  -8.4346 eps 0.0924 iteration 186 qsel 4
#> logLik  -7.6021 eps 0.0337 iteration 187 qsel 4
#> logLik  -7.4942 eps 0.0012 iteration 188 qsel 4
#> logLik  -7.4066 eps 0.0011 iteration 189 qsel 4
#> logLik  -7.3276 eps 0.0076 iteration 190 qsel 4
#> logLik  -7.2421 eps 0.0113 iteration 191 qsel 4
#> logLik  -7.1644 eps 0.0338 iteration 192 qsel 4
#> logLik  -7.0816 eps 0.0109 iteration 193 qsel 4
#> logLik  -7.0051 eps 0.0073 iteration 194 qsel 4
#> logLik  -6.9247 eps 0.0091 iteration 195 qsel 4
#> logLik  -6.8482 eps 0.0010 iteration 196 qsel 4
#> logLik  -6.7727 eps 0.0064 iteration 197 qsel 4
#> logLik  -6.6941 eps 0.0088 iteration 198 qsel 4
#> logLik  -5.5430 eps 0.0209 iteration 199 qsel 4
#> logLik  -5.3953 eps 0.0102 iteration 200 qsel 4
#> logLik  -5.3178 eps 0.0010 iteration 201 qsel 4
#> logLik  -5.2466 eps 0.0060 iteration 202 qsel 4
#> logLik  -5.1701 eps 0.0124 iteration 203 qsel 4
#> logLik  -5.1000 eps 0.0062 iteration 204 qsel 4
#> logLik  -5.0257 eps 0.0159 iteration 205 qsel 4
#> logLik  -4.0493 eps 0.0239 iteration 206 qsel 4
#> logLik  -3.9729 eps 0.0010 iteration 207 qsel 4
#> logLik  -3.9044 eps 0.0051 iteration 208 qsel 4
#> logLik  -3.8319 eps 0.0074 iteration 209 qsel 4
#> logLik  -3.7646 eps 0.0050 iteration 210 qsel 4
#> logLik  -3.6942 eps 0.0073 iteration 211 qsel 4
#> logLik  -2.6533 eps 0.0193 iteration 212 qsel 4
#> logLik  -2.5774 eps 0.0010 iteration 213 qsel 4
#> logLik  -2.5125 eps 0.0050 iteration 214 qsel 4
#> logLik  -2.4438 eps 0.0077 iteration 215 qsel 4
#> logLik  -2.3801 eps 0.0050 iteration 216 qsel 4
#> logLik  -2.3134 eps 0.0079 iteration 217 qsel 4
#> logLik  -2.2491 eps 0.0009 iteration 218 qsel 4
#> logLik  -2.1861 eps 0.0050 iteration 219 qsel 4
#> logLik  -2.1209 eps 0.0088 iteration 220 qsel 4
#> logLik  -2.0592 eps 0.0050 iteration 221 qsel 4
#> logLik  -1.9958 eps 0.0129 iteration 222 qsel 4
#> logLik  -1.9353 eps 0.0052 iteration 223 qsel 4
#> logLik  -1.8736 eps 0.0125 iteration 224 qsel 4
#> logLik  -0.9330 eps 0.0192 iteration 225 qsel 4
#> logLik  -0.8680 eps 0.0009 iteration 226 qsel 4
#> logLik  -0.8082 eps 0.0067 iteration 227 qsel 4
#> logLik  -0.7477 eps 0.0073 iteration 228 qsel 4
#> logLik  -0.6196 eps 0.0085 iteration 229 qsel 4
#> logLik  -0.5610 eps 0.0050 iteration 230 qsel 4
#> logLik  -0.5021 eps 0.0069 iteration 231 qsel 4
#> logLik   0.3969 eps 0.0347 iteration 232 qsel 4
#> logLik   0.4619 eps 0.0009 iteration 233 qsel 4
#> logLik   0.5185 eps 0.0056 iteration 234 qsel 4
#> logLik   0.5762 eps 0.0078 iteration 235 qsel 4
#> logLik   0.6317 eps 0.0054 iteration 236 qsel 4
#> logLik   0.6879 eps 0.0073 iteration 237 qsel 4
#> logLik   0.7428 eps 0.0008 iteration 238 qsel 4
#> logLik   0.7973 eps 0.0052 iteration 239 qsel 4
#> logLik   0.8524 eps 0.0070 iteration 240 qsel 4
#> logLik   1.7027 eps 0.0248 iteration 241 qsel 4
#> logLik   1.7573 eps 0.0061 iteration 242 qsel 4
#> logLik   1.8111 eps 0.0386 iteration 243 qsel 4
#> logLik   1.8651 eps 0.0008 iteration 244 qsel 4
#> logLik   1.9187 eps 0.0058 iteration 245 qsel 4
#> logLik   1.9714 eps 0.0097 iteration 246 qsel 4
#> logLik   2.0237 eps 0.0055 iteration 247 qsel 4
#> logLik   2.0753 eps 0.0084 iteration 248 qsel 4
#> logLik   2.1263 eps 0.0054 iteration 249 qsel 4
#> logLik   2.1767 eps 0.0087 iteration 250 qsel 4
#> logLik   2.2266 eps 0.0054 iteration 251 qsel 4
#> logLik   3.0237 eps 0.0241 iteration 252 qsel 4
#> logLik   3.0777 eps 0.0008 iteration 253 qsel 4
#> logLik   3.1272 eps 0.0095 iteration 254 qsel 4
#> logLik   3.2381 eps 0.0111 iteration 255 qsel 4
#> logLik   3.2859 eps 0.0088 iteration 256 qsel 4
#> logLik   3.3341 eps 0.0088 iteration 257 qsel 4
#> logLik   4.0913 eps 0.0242 iteration 258 qsel 4
#> logLik   4.1440 eps 0.0008 iteration 259 qsel 4
#> logLik   4.1936 eps 0.0056 iteration 260 qsel 4
#> logLik   4.2410 eps 0.0086 iteration 261 qsel 4
#> logLik   4.2892 eps 0.0055 iteration 262 qsel 4
#> logLik   4.3356 eps 0.0131 iteration 263 qsel 4
#> logLik   4.3826 eps 0.0056 iteration 264 qsel 4
#> logLik   4.4280 eps 0.0090 iteration 265 qsel 4
#> logLik   4.4738 eps 0.0054 iteration 266 qsel 4
#> logLik   4.5194 eps 0.0007 iteration 267 qsel 4
#> logLik   4.5641 eps 0.0124 iteration 268 qsel 4
#> logLik   4.6089 eps 0.0056 iteration 269 qsel 4
#> logLik   5.3309 eps 0.5096 iteration 270 qsel 4
#> logLik   5.3746 eps 0.0214 iteration 271 qsel 4
#> logLik   5.4184 eps 0.0007 iteration 272 qsel 4
#> logLik   5.4612 eps 0.0098 iteration 273 qsel 4
#> logLik   5.5040 eps 0.0235 iteration 274 qsel 4
#> logLik   6.1807 eps 0.0238 iteration 275 qsel 4
#> logLik   6.2255 eps 0.0127 iteration 276 qsel 4
#> logLik   6.2675 eps 0.0106 iteration 277 qsel 4
#> logLik   6.3111 eps 0.0432 iteration 278 qsel 4
#> logLik   6.3538 eps 0.0007 iteration 279 qsel 4
#> logLik   6.3951 eps 0.0116 iteration 280 qsel 4
#> logLik   6.4377 eps 0.0165 iteration 281 qsel 4
#> logLik   6.4783 eps 0.0187 iteration 282 qsel 4
#> logLik   6.5197 eps 0.0184 iteration 283 qsel 4
#> logLik   6.5594 eps 0.0131 iteration 284 qsel 4
#> logLik   6.5998 eps 0.0330 iteration 285 qsel 4
#> logLik   7.2524 eps 0.0205 iteration 286 qsel 4
#> logLik   7.3503 eps 0.0073 iteration 287 qsel 4
#> logLik   7.3927 eps 0.0007 iteration 288 qsel 4
#> logLik   7.4317 eps 0.0073 iteration 289 qsel 4
#> logLik   7.4702 eps 0.0054 iteration 290 qsel 4
#> logLik   7.5083 eps 0.0060 iteration 291 qsel 4
#> logLik   7.5459 eps 0.0042 iteration 292 qsel 4
#> logLik   7.5831 eps 0.0058 iteration 293 qsel 4
#> logLik   8.1807 eps 0.0197 iteration 294 qsel 4
#> logLik   8.2212 eps 0.0007 iteration 295 qsel 4
#> logLik   8.2617 eps 0.0041 iteration 296 qsel 4
#> logLik   8.2985 eps 0.2928 iteration 297 qsel 4
#> logLik   8.3378 eps 0.0052 iteration 298 qsel 4
#> logLik   8.3740 eps 0.0081 iteration 299 qsel 4
#> logLik   8.4122 eps 0.0058 iteration 300 qsel 4
#> logLik   9.0051 eps 0.0138 iteration 301 qsel 4
#> logLik   9.0438 eps 0.0007 iteration 302 qsel 4
#> logLik   9.0795 eps 0.0068 iteration 303 qsel 4
#> logLik   9.1157 eps 0.0042 iteration 304 qsel 4
#> logLik   9.1506 eps 0.0080 iteration 305 qsel 4
#> logLik   9.1859 eps 0.0039 iteration 306 qsel 4
#> logLik   9.2201 eps 0.0076 iteration 307 qsel 4
#> logLik   9.2544 eps 0.0041 iteration 308 qsel 4
#> logLik   9.2879 eps 0.0143 iteration 309 qsel 4
#> logLik   9.3214 eps 0.0039 iteration 310 qsel 4
#> logLik   9.3548 eps 0.0006 iteration 311 qsel 4
#> logLik   9.3876 eps 0.0059 iteration 312 qsel 4
#> logLik   9.4205 eps 0.0035 iteration 313 qsel 4
#> logLik   9.9377 eps 0.0384 iteration 314 qsel 4
#> logLik  10.4780 eps 0.0235 iteration 315 qsel 4
#> logLik  10.5130 eps 0.0006 iteration 316 qsel 4
#> logLik  10.5452 eps 0.0079 iteration 317 qsel 4
#> logLik  10.5806 eps 0.0064 iteration 318 qsel 4
#> logLik  10.6695 eps 0.0140 iteration 319 qsel 4
#> logLik  10.7013 eps 0.0068 iteration 320 qsel 4
#> logLik  10.7356 eps 0.0042 iteration 321 qsel 4
#> logLik  10.7669 eps 0.0053 iteration 322 qsel 4
#> logLik  10.8002 eps 0.0035 iteration 323 qsel 4
#> logLik  10.8310 eps 0.0048 iteration 324 qsel 4
#> logLik  10.8633 eps 0.0032 iteration 325 qsel 4
#> logLik  10.8941 eps 0.0006 iteration 326 qsel 4
#> logLik  10.9244 eps 0.0045 iteration 327 qsel 4
#> logLik  10.9559 eps 0.0031 iteration 328 qsel 4
#> logLik  10.9857 eps 0.0044 iteration 329 qsel 4
#> logLik  11.0163 eps 0.0030 iteration 330 qsel 4
#> logLik  11.0455 eps 0.0043 iteration 331 qsel 4
#> logLik  11.0753 eps 0.0030 iteration 332 qsel 4
#> logLik  11.5658 eps 0.0122 iteration 333 qsel 4
#> logLik  11.5952 eps 0.0006 iteration 334 qsel 4
#> logLik  11.6238 eps 0.0045 iteration 335 qsel 4
#> logLik  12.0622 eps 0.0289 iteration 336 qsel 4
#> logLik  12.0945 eps 0.0039 iteration 337 qsel 4
#> logLik  12.1228 eps 0.0065 iteration 338 qsel 4
#> logLik  12.1540 eps 0.0034 iteration 339 qsel 4
#> logLik  12.1819 eps 0.0052 iteration 340 qsel 4
#> logLik  12.2121 eps 0.0032 iteration 341 qsel 4
#> logLik  12.2399 eps 0.0006 iteration 342 qsel 4
#> logLik  12.2675 eps 0.0049 iteration 343 qsel 4
#> logLik  12.2968 eps 0.0031 iteration 344 qsel 4
#> logLik  12.3240 eps 0.0051 iteration 345 qsel 4
#> logLik  12.3525 eps 0.0030 iteration 346 qsel 4
#> logLik  12.3791 eps 0.0096 iteration 347 qsel 4
#> logLik  12.4068 eps 0.0033 iteration 348 qsel 4
#> logLik  12.4917 eps 0.0062 iteration 349 qsel 4
#> logLik  12.9389 eps 0.0117 iteration 350 qsel 4
#> logLik  12.9658 eps 0.0005 iteration 351 qsel 4
#> logLik  12.9920 eps 0.0045 iteration 352 qsel 4
#> logLik  13.0183 eps 0.0031 iteration 353 qsel 4
#> logLik  13.0439 eps 0.0045 iteration 354 qsel 4
#> logLik  13.0696 eps 0.0030 iteration 355 qsel 4
#> logLik  13.4299 eps 0.0198 iteration 356 qsel 4
#> logLik  13.4550 eps 0.0078 iteration 357 qsel 4
#> logLik  13.4838 eps 0.0040 iteration 358 qsel 4
#> logLik  13.5087 eps 0.0005 iteration 359 qsel 4
#> logLik  13.5336 eps 0.0049 iteration 360 qsel 4
#> logLik  13.5615 eps 0.0032 iteration 361 qsel 4
#> logLik  13.5861 eps 0.0044 iteration 362 qsel 4
#> logLik  13.6131 eps 0.0030 iteration 363 qsel 4
#> logLik  13.6373 eps 0.0041 iteration 364 qsel 4
#> logLik  13.6634 eps 0.0029 iteration 365 qsel 4
#> logLik  13.6873 eps 0.0039 iteration 366 qsel 4
#> logLik  13.7126 eps 0.0029 iteration 367 qsel 4
#> logLik  14.1206 eps 0.0251 iteration 368 qsel 4
#> logLik  14.1444 eps 0.0005 iteration 369 qsel 4
#> logLik  14.1679 eps 0.0118 iteration 370 qsel 4
#> logLik  14.1919 eps 0.0041 iteration 371 qsel 4
#> logLik  14.2150 eps 0.0065 iteration 372 qsel 4
#> logLik  14.2383 eps 0.0049 iteration 373 qsel 4
#> logLik  14.2609 eps 0.0052 iteration 374 qsel 4
#> logLik  14.2837 eps 0.0132 iteration 375 qsel 4
#> logLik  14.5815 eps 0.0140 iteration 376 qsel 4
#> logLik  14.6035 eps 0.0044 iteration 377 qsel 4
#> logLik  14.6293 eps 0.0033 iteration 378 qsel 4
#> logLik  14.6514 eps 0.0005 iteration 379 qsel 4
#> logLik  14.6733 eps 0.0044 iteration 380 qsel 4
#> logLik  14.6982 eps 0.0032 iteration 381 qsel 4
#> logLik  14.7199 eps 0.0044 iteration 382 qsel 4
#> logLik  14.7439 eps 0.0032 iteration 383 qsel 4
#> logLik  15.1176 eps 0.0123 iteration 384 qsel 4
#> logLik  15.2010 eps 0.0081 iteration 385 qsel 4
#> logLik  15.2224 eps 0.0048 iteration 386 qsel 4
#> logLik  15.2451 eps 0.0038 iteration 387 qsel 4
#> logLik  15.2663 eps 0.0005 iteration 388 qsel 4
#> logLik  15.2874 eps 0.0045 iteration 389 qsel 4
#> logLik  15.3095 eps 0.0036 iteration 390 qsel 4
#> logLik  15.3302 eps 0.0042 iteration 391 qsel 4
#> logLik  15.3517 eps 0.0034 iteration 392 qsel 4
#> logLik  15.3721 eps 0.0040 iteration 393 qsel 4
#> logLik  15.3929 eps 0.0032 iteration 394 qsel 4
#> logLik  15.6578 eps 0.0228 iteration 395 qsel 4
#> logLik  15.6777 eps 0.0040 iteration 396 qsel 4
#> logLik  15.7013 eps 0.0358 iteration 397 qsel 4
#> logLik  15.7211 eps 0.0041 iteration 398 qsel 4
#> logLik  15.7439 eps 0.0071 iteration 399 qsel 4
#> logLik  16.0847 eps 0.0202 iteration 400 qsel 4
#> 
#>  elapsed time:  1.52sec


plot(b)