R Under development (unstable) (2025-09-25 r88874 ucrt) -- "Unsuffered Consequences"
Copyright (C) 2025 The R Foundation for Statistical Computing
Platform: x86_64-w64-mingw32/x64

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.

> # Clean environment and load libraries
> rm(list = ls())  # Clear environment to avoid conflicts
> 
> 
> # Source files
> library(MixStable)

Attaching package: 'MixStable'

The following object is masked from 'package:graphics':

    clip

> 
> # Safe mean extractor with NA-check
> safe_mean <- function(x) {
+   vals <- unlist(x)
+   if (!is.numeric(vals) || all(is.na(vals)) || length(vals) == 0) return(NA_real_)
+   return(mean(vals, na.rm = TRUE))
+ }
> # === Run Metropolis-Hastings and plot ===
> set.seed(123)
> xx <- seq(-8, 11, length.out = 200)
> xx_true <- seq(-8, 11, length.out = 200)
> yy_true <- r_stable_pdf(xx_true, 1.5, -0.5, 4, 5)
> 
> # Generate synthetic data
> lok <- rnorm(1000, mean = 0, sd = 1)
> 
> # ⚠️ Increase iterations to ensure sufficient samples
> mh_result <- metropolis_hastings(CDF, iterations = 1000, lok = lok, aa = c(1, 1), proposal_std = 0.1)
Iteration 1 / 1000 completed.
Iteration 2 / 1000 completed.
Iteration 3 / 1000 completed.
Iteration 4 / 1000 completed.
Iteration 5 / 1000 completed.
Iteration 6 / 1000 completed.
Iteration 7 / 1000 completed.
Iteration 8 / 1000 completed.
Iteration 9 / 1000 completed.
Iteration 10 / 1000 completed.
Iteration 11 / 1000 completed.
Iteration 12 / 1000 completed.
Iteration 13 / 1000 completed.
Iteration 14 / 1000 completed.
Iteration 15 / 1000 completed.
Iteration 16 / 1000 completed.
Iteration 17 / 1000 completed.
Iteration 18 / 1000 completed.
Iteration 19 / 1000 completed.
Iteration 20 / 1000 completed.
Iteration 21 / 1000 completed.
Iteration 22 / 1000 completed.
Iteration 23 / 1000 completed.
Iteration 24 / 1000 completed.
Iteration 25 / 1000 completed.
Iteration 26 / 1000 completed.
Iteration 27 / 1000 completed.
Iteration 28 / 1000 completed.
Iteration 29 / 1000 completed.
Iteration 30 / 1000 completed.
Iteration 31 / 1000 completed.
Iteration 32 / 1000 completed.
Iteration 33 / 1000 completed.
Iteration 34 / 1000 completed.
Iteration 35 / 1000 completed.
Iteration 36 / 1000 completed.
Iteration 37 / 1000 completed.
Iteration 38 / 1000 completed.
Iteration 39 / 1000 completed.
Iteration 40 / 1000 completed.
Iteration 41 / 1000 completed.
Iteration 42 / 1000 completed.
Iteration 43 / 1000 completed.
Iteration 44 / 1000 completed.
Iteration 45 / 1000 completed.
Iteration 46 / 1000 completed.
Iteration 47 / 1000 completed.
Iteration 48 / 1000 completed.
Iteration 49 / 1000 completed.
Iteration 50 / 1000 completed.
Iteration 51 / 1000 completed.
Iteration 52 / 1000 completed.
Iteration 53 / 1000 completed.
Iteration 54 / 1000 completed.
Iteration 55 / 1000 completed.
Iteration 56 / 1000 completed.
Iteration 57 / 1000 completed.
Iteration 58 / 1000 completed.
Iteration 59 / 1000 completed.
Iteration 60 / 1000 completed.
Iteration 61 / 1000 completed.
Iteration 62 / 1000 completed.
Iteration 63 / 1000 completed.
Iteration 64 / 1000 completed.
Iteration 65 / 1000 completed.
Iteration 66 / 1000 completed.
Iteration 67 / 1000 completed.
Iteration 68 / 1000 completed.
Iteration 69 / 1000 completed.
Iteration 70 / 1000 completed.
Iteration 71 / 1000 completed.
Iteration 72 / 1000 completed.
Iteration 73 / 1000 completed.
Iteration 74 / 1000 completed.
Iteration 75 / 1000 completed.
Iteration 76 / 1000 completed.
Iteration 77 / 1000 completed.
Iteration 78 / 1000 completed.
Iteration 79 / 1000 completed.
Iteration 80 / 1000 completed.
Iteration 81 / 1000 completed.
Iteration 82 / 1000 completed.
Iteration 83 / 1000 completed.
Iteration 84 / 1000 completed.
Iteration 85 / 1000 completed.
Iteration 86 / 1000 completed.
Iteration 87 / 1000 completed.
Iteration 88 / 1000 completed.
Iteration 89 / 1000 completed.
Iteration 90 / 1000 completed.
Iteration 91 / 1000 completed.
Iteration 92 / 1000 completed.
Iteration 93 / 1000 completed.
Iteration 94 / 1000 completed.
Iteration 95 / 1000 completed.
Iteration 96 / 1000 completed.
Iteration 97 / 1000 completed.
Iteration 98 / 1000 completed.
Iteration 99 / 1000 completed.
Iteration 100 / 1000 completed.
Iteration 101 / 1000 completed.
Iteration 102 / 1000 completed.
Iteration 103 / 1000 completed.
Iteration 104 / 1000 completed.
Iteration 105 / 1000 completed.
Iteration 106 / 1000 completed.
Iteration 107 / 1000 completed.
Iteration 108 / 1000 completed.
Iteration 109 / 1000 completed.
Iteration 110 / 1000 completed.
Iteration 111 / 1000 completed.
Iteration 112 / 1000 completed.
Iteration 113 / 1000 completed.
Iteration 114 / 1000 completed.
Iteration 115 / 1000 completed.
Iteration 116 / 1000 completed.
Iteration 117 / 1000 completed.
Iteration 118 / 1000 completed.
Iteration 119 / 1000 completed.
Iteration 120 / 1000 completed.
Iteration 121 / 1000 completed.
Iteration 122 / 1000 completed.
Iteration 123 / 1000 completed.
Iteration 124 / 1000 completed.
Iteration 125 / 1000 completed.
Iteration 126 / 1000 completed.
Iteration 127 / 1000 completed.
Iteration 128 / 1000 completed.
Iteration 129 / 1000 completed.
Iteration 130 / 1000 completed.
Iteration 131 / 1000 completed.
Iteration 132 / 1000 completed.
Iteration 133 / 1000 completed.
Iteration 134 / 1000 completed.
Iteration 135 / 1000 completed.
Iteration 136 / 1000 completed.
Iteration 137 / 1000 completed.
Iteration 138 / 1000 completed.
Iteration 139 / 1000 completed.
Iteration 140 / 1000 completed.
Iteration 141 / 1000 completed.
Iteration 142 / 1000 completed.
Iteration 143 / 1000 completed.
Iteration 144 / 1000 completed.
Iteration 145 / 1000 completed.
Iteration 146 / 1000 completed.
Iteration 147 / 1000 completed.
Iteration 148 / 1000 completed.
Iteration 149 / 1000 completed.
Iteration 150 / 1000 completed.
Iteration 151 / 1000 completed.
Iteration 152 / 1000 completed.
Iteration 153 / 1000 completed.
Iteration 154 / 1000 completed.
Iteration 155 / 1000 completed.
Iteration 156 / 1000 completed.
Iteration 157 / 1000 completed.
Iteration 158 / 1000 completed.
Iteration 159 / 1000 completed.
Iteration 160 / 1000 completed.
Iteration 161 / 1000 completed.
Iteration 162 / 1000 completed.
Iteration 163 / 1000 completed.
Iteration 164 / 1000 completed.
Iteration 165 / 1000 completed.
Iteration 166 / 1000 completed.
Iteration 167 / 1000 completed.
Iteration 168 / 1000 completed.
Iteration 169 / 1000 completed.
Iteration 170 / 1000 completed.
Iteration 171 / 1000 completed.
Iteration 172 / 1000 completed.
Iteration 173 / 1000 completed.
Iteration 174 / 1000 completed.
Iteration 175 / 1000 completed.
Iteration 176 / 1000 completed.
Iteration 177 / 1000 completed.
Iteration 178 / 1000 completed.
Iteration 179 / 1000 completed.
Iteration 180 / 1000 completed.
Iteration 181 / 1000 completed.
Iteration 182 / 1000 completed.
Iteration 183 / 1000 completed.
Iteration 184 / 1000 completed.
Iteration 185 / 1000 completed.
Iteration 186 / 1000 completed.
Iteration 187 / 1000 completed.
Iteration 188 / 1000 completed.
Iteration 189 / 1000 completed.
Iteration 190 / 1000 completed.
Iteration 191 / 1000 completed.
Iteration 192 / 1000 completed.
Iteration 193 / 1000 completed.
Iteration 194 / 1000 completed.
Iteration 195 / 1000 completed.
Iteration 196 / 1000 completed.
Iteration 197 / 1000 completed.
Iteration 198 / 1000 completed.
Iteration 199 / 1000 completed.
Iteration 200 / 1000 completed.
Iteration 201 / 1000 completed.
Iteration 202 / 1000 completed.
Iteration 203 / 1000 completed.
Iteration 204 / 1000 completed.
Iteration 205 / 1000 completed.
Iteration 206 / 1000 completed.
Iteration 207 / 1000 completed.
Iteration 208 / 1000 completed.
Iteration 209 / 1000 completed.
Iteration 210 / 1000 completed.
Iteration 211 / 1000 completed.
Iteration 212 / 1000 completed.
Iteration 213 / 1000 completed.
Iteration 214 / 1000 completed.
Iteration 215 / 1000 completed.
Iteration 216 / 1000 completed.
Iteration 217 / 1000 completed.
Iteration 218 / 1000 completed.
Iteration 219 / 1000 completed.
Iteration 220 / 1000 completed.
Iteration 221 / 1000 completed.
Iteration 222 / 1000 completed.
Iteration 223 / 1000 completed.
Iteration 224 / 1000 completed.
Iteration 225 / 1000 completed.
Iteration 226 / 1000 completed.
Iteration 227 / 1000 completed.
Iteration 228 / 1000 completed.
Iteration 229 / 1000 completed.
Iteration 230 / 1000 completed.
Iteration 231 / 1000 completed.
Iteration 232 / 1000 completed.
Iteration 233 / 1000 completed.
Iteration 234 / 1000 completed.
Iteration 235 / 1000 completed.
Iteration 236 / 1000 completed.
Iteration 237 / 1000 completed.
Iteration 238 / 1000 completed.
Iteration 239 / 1000 completed.
Iteration 240 / 1000 completed.
Iteration 241 / 1000 completed.
Iteration 242 / 1000 completed.
Iteration 243 / 1000 completed.
Iteration 244 / 1000 completed.
Iteration 245 / 1000 completed.
Iteration 246 / 1000 completed.
Iteration 247 / 1000 completed.
Iteration 248 / 1000 completed.
Iteration 249 / 1000 completed.
Iteration 250 / 1000 completed.
Iteration 251 / 1000 completed.
Iteration 252 / 1000 completed.
Iteration 253 / 1000 completed.
Iteration 254 / 1000 completed.
Iteration 255 / 1000 completed.
Iteration 256 / 1000 completed.
Iteration 257 / 1000 completed.
Iteration 258 / 1000 completed.
Iteration 259 / 1000 completed.
Iteration 260 / 1000 completed.
Iteration 261 / 1000 completed.
Iteration 262 / 1000 completed.
Iteration 263 / 1000 completed.
Iteration 264 / 1000 completed.
Iteration 265 / 1000 completed.
Iteration 266 / 1000 completed.
Iteration 267 / 1000 completed.
Iteration 268 / 1000 completed.
Iteration 269 / 1000 completed.
Iteration 270 / 1000 completed.
Iteration 271 / 1000 completed.
Iteration 272 / 1000 completed.
Iteration 273 / 1000 completed.
Iteration 274 / 1000 completed.
Iteration 275 / 1000 completed.
Iteration 276 / 1000 completed.
Iteration 277 / 1000 completed.
Iteration 278 / 1000 completed.
Iteration 279 / 1000 completed.
Iteration 280 / 1000 completed.
Iteration 281 / 1000 completed.
Iteration 282 / 1000 completed.
Iteration 283 / 1000 completed.
Iteration 284 / 1000 completed.
Iteration 285 / 1000 completed.
Iteration 286 / 1000 completed.
Iteration 287 / 1000 completed.
Iteration 288 / 1000 completed.
Iteration 289 / 1000 completed.
Iteration 290 / 1000 completed.
Iteration 291 / 1000 completed.
Iteration 292 / 1000 completed.
Iteration 293 / 1000 completed.
Iteration 294 / 1000 completed.
Iteration 295 / 1000 completed.
Iteration 296 / 1000 completed.
Iteration 297 / 1000 completed.
Iteration 298 / 1000 completed.
Iteration 299 / 1000 completed.
Iteration 300 / 1000 completed.
Iteration 301 / 1000 completed.
Iteration 302 / 1000 completed.
Iteration 303 / 1000 completed.
Iteration 304 / 1000 completed.
Iteration 305 / 1000 completed.
Iteration 306 / 1000 completed.
Iteration 307 / 1000 completed.
Iteration 308 / 1000 completed.
Iteration 309 / 1000 completed.
Iteration 310 / 1000 completed.
Iteration 311 / 1000 completed.
Iteration 312 / 1000 completed.
Iteration 313 / 1000 completed.
Iteration 314 / 1000 completed.
Iteration 315 / 1000 completed.
Iteration 316 / 1000 completed.
Iteration 317 / 1000 completed.
Iteration 318 / 1000 completed.
Iteration 319 / 1000 completed.
Iteration 320 / 1000 completed.
Iteration 321 / 1000 completed.
Iteration 322 / 1000 completed.
Iteration 323 / 1000 completed.
Iteration 324 / 1000 completed.
Iteration 325 / 1000 completed.
Iteration 326 / 1000 completed.
Iteration 327 / 1000 completed.
Iteration 328 / 1000 completed.
Iteration 329 / 1000 completed.
Iteration 330 / 1000 completed.
Iteration 331 / 1000 completed.
Iteration 332 / 1000 completed.
Iteration 333 / 1000 completed.
Iteration 334 / 1000 completed.
Iteration 335 / 1000 completed.
Iteration 336 / 1000 completed.
Iteration 337 / 1000 completed.
Iteration 338 / 1000 completed.
Iteration 339 / 1000 completed.
Iteration 340 / 1000 completed.
Iteration 341 / 1000 completed.
Iteration 342 / 1000 completed.
Iteration 343 / 1000 completed.
Iteration 344 / 1000 completed.
Iteration 345 / 1000 completed.
Iteration 346 / 1000 completed.
Iteration 347 / 1000 completed.
Iteration 348 / 1000 completed.
Iteration 349 / 1000 completed.
Iteration 350 / 1000 completed.
Iteration 351 / 1000 completed.
Iteration 352 / 1000 completed.
Iteration 353 / 1000 completed.
Iteration 354 / 1000 completed.
Iteration 355 / 1000 completed.
Iteration 356 / 1000 completed.
Iteration 357 / 1000 completed.
Iteration 358 / 1000 completed.
Iteration 359 / 1000 completed.
Iteration 360 / 1000 completed.
Iteration 361 / 1000 completed.
Iteration 362 / 1000 completed.
Iteration 363 / 1000 completed.
Iteration 364 / 1000 completed.
Iteration 365 / 1000 completed.
Iteration 366 / 1000 completed.
Iteration 367 / 1000 completed.
Iteration 368 / 1000 completed.
Iteration 369 / 1000 completed.
Iteration 370 / 1000 completed.
Iteration 371 / 1000 completed.
Iteration 372 / 1000 completed.
Iteration 373 / 1000 completed.
Iteration 374 / 1000 completed.
Iteration 375 / 1000 completed.
Iteration 376 / 1000 completed.
Iteration 377 / 1000 completed.
Iteration 378 / 1000 completed.
Iteration 379 / 1000 completed.
Iteration 380 / 1000 completed.
Iteration 381 / 1000 completed.
Iteration 382 / 1000 completed.
Iteration 383 / 1000 completed.
Iteration 384 / 1000 completed.
Iteration 385 / 1000 completed.
Iteration 386 / 1000 completed.
Iteration 387 / 1000 completed.
Iteration 388 / 1000 completed.
Iteration 389 / 1000 completed.
Iteration 390 / 1000 completed.
Iteration 391 / 1000 completed.
Iteration 392 / 1000 completed.
Iteration 393 / 1000 completed.
Iteration 394 / 1000 completed.
Iteration 395 / 1000 completed.
Iteration 396 / 1000 completed.
Iteration 397 / 1000 completed.
Iteration 398 / 1000 completed.
Iteration 399 / 1000 completed.
Iteration 400 / 1000 completed.
Iteration 401 / 1000 completed.
Iteration 402 / 1000 completed.
Iteration 403 / 1000 completed.
Iteration 404 / 1000 completed.
Iteration 405 / 1000 completed.
Iteration 406 / 1000 completed.
Iteration 407 / 1000 completed.
Iteration 408 / 1000 completed.
Iteration 409 / 1000 completed.
Iteration 410 / 1000 completed.
Iteration 411 / 1000 completed.
Iteration 412 / 1000 completed.
Iteration 413 / 1000 completed.
Iteration 414 / 1000 completed.
Iteration 415 / 1000 completed.
Iteration 416 / 1000 completed.
Iteration 417 / 1000 completed.
Iteration 418 / 1000 completed.
Iteration 419 / 1000 completed.
Iteration 420 / 1000 completed.
Iteration 421 / 1000 completed.
Iteration 422 / 1000 completed.
Iteration 423 / 1000 completed.
Iteration 424 / 1000 completed.
Iteration 425 / 1000 completed.
Iteration 426 / 1000 completed.
Iteration 427 / 1000 completed.
Iteration 428 / 1000 completed.
Iteration 429 / 1000 completed.
Iteration 430 / 1000 completed.
Iteration 431 / 1000 completed.
Iteration 432 / 1000 completed.
Iteration 433 / 1000 completed.
Iteration 434 / 1000 completed.
Iteration 435 / 1000 completed.
Iteration 436 / 1000 completed.
Iteration 437 / 1000 completed.
Iteration 438 / 1000 completed.
Iteration 439 / 1000 completed.
Iteration 440 / 1000 completed.
Iteration 441 / 1000 completed.
Iteration 442 / 1000 completed.
Iteration 443 / 1000 completed.
Iteration 444 / 1000 completed.
Iteration 445 / 1000 completed.
Iteration 446 / 1000 completed.
Iteration 447 / 1000 completed.
Iteration 448 / 1000 completed.
Iteration 449 / 1000 completed.
Iteration 450 / 1000 completed.
Iteration 451 / 1000 completed.
Iteration 452 / 1000 completed.
Iteration 453 / 1000 completed.
Iteration 454 / 1000 completed.
Iteration 455 / 1000 completed.
Iteration 456 / 1000 completed.
Iteration 457 / 1000 completed.
Iteration 458 / 1000 completed.
Iteration 459 / 1000 completed.
Iteration 460 / 1000 completed.
Iteration 461 / 1000 completed.
Iteration 462 / 1000 completed.
Iteration 463 / 1000 completed.
Iteration 464 / 1000 completed.
Iteration 465 / 1000 completed.
Iteration 466 / 1000 completed.
Iteration 467 / 1000 completed.
Iteration 468 / 1000 completed.
Iteration 469 / 1000 completed.
Iteration 470 / 1000 completed.
Iteration 471 / 1000 completed.
Iteration 472 / 1000 completed.
Iteration 473 / 1000 completed.
Iteration 474 / 1000 completed.
Iteration 475 / 1000 completed.
Iteration 476 / 1000 completed.
Iteration 477 / 1000 completed.
Iteration 478 / 1000 completed.
Iteration 479 / 1000 completed.
Iteration 480 / 1000 completed.
Iteration 481 / 1000 completed.
Iteration 482 / 1000 completed.
Iteration 483 / 1000 completed.
Iteration 484 / 1000 completed.
Iteration 485 / 1000 completed.
Iteration 486 / 1000 completed.
Iteration 487 / 1000 completed.
Iteration 488 / 1000 completed.
Iteration 489 / 1000 completed.
Iteration 490 / 1000 completed.
Iteration 491 / 1000 completed.
Iteration 492 / 1000 completed.
Iteration 493 / 1000 completed.
Iteration 494 / 1000 completed.
Iteration 495 / 1000 completed.
Iteration 496 / 1000 completed.
Iteration 497 / 1000 completed.
Iteration 498 / 1000 completed.
Iteration 499 / 1000 completed.
Iteration 500 / 1000 completed.
Iteration 501 / 1000 completed.
Iteration 502 / 1000 completed.
Iteration 503 / 1000 completed.
Iteration 504 / 1000 completed.
Iteration 505 / 1000 completed.
Iteration 506 / 1000 completed.
Iteration 507 / 1000 completed.
Iteration 508 / 1000 completed.
Iteration 509 / 1000 completed.
Iteration 510 / 1000 completed.
Iteration 511 / 1000 completed.
Iteration 512 / 1000 completed.
Iteration 513 / 1000 completed.
Iteration 514 / 1000 completed.
Iteration 515 / 1000 completed.
Iteration 516 / 1000 completed.
Iteration 517 / 1000 completed.
Iteration 518 / 1000 completed.
Iteration 519 / 1000 completed.
Iteration 520 / 1000 completed.
Iteration 521 / 1000 completed.
Iteration 522 / 1000 completed.
Iteration 523 / 1000 completed.
Iteration 524 / 1000 completed.
Iteration 525 / 1000 completed.
Iteration 526 / 1000 completed.
Iteration 527 / 1000 completed.
Iteration 528 / 1000 completed.
Iteration 529 / 1000 completed.
Iteration 530 / 1000 completed.
Iteration 531 / 1000 completed.
Iteration 532 / 1000 completed.
Iteration 533 / 1000 completed.
Iteration 534 / 1000 completed.
Iteration 535 / 1000 completed.
Iteration 536 / 1000 completed.
Iteration 537 / 1000 completed.
Iteration 538 / 1000 completed.
Iteration 539 / 1000 completed.
Iteration 540 / 1000 completed.
Iteration 541 / 1000 completed.
Iteration 542 / 1000 completed.
Iteration 543 / 1000 completed.
Iteration 544 / 1000 completed.
Iteration 545 / 1000 completed.
Iteration 546 / 1000 completed.
Iteration 547 / 1000 completed.
Iteration 548 / 1000 completed.
Iteration 549 / 1000 completed.
Iteration 550 / 1000 completed.
Iteration 551 / 1000 completed.
Iteration 552 / 1000 completed.
Iteration 553 / 1000 completed.
Iteration 554 / 1000 completed.
Iteration 555 / 1000 completed.
Iteration 556 / 1000 completed.
Iteration 557 / 1000 completed.
Iteration 558 / 1000 completed.
Iteration 559 / 1000 completed.
Iteration 560 / 1000 completed.
Iteration 561 / 1000 completed.
Iteration 562 / 1000 completed.
Iteration 563 / 1000 completed.
Iteration 564 / 1000 completed.
Iteration 565 / 1000 completed.
Iteration 566 / 1000 completed.
Iteration 567 / 1000 completed.
Iteration 568 / 1000 completed.
Iteration 569 / 1000 completed.
Iteration 570 / 1000 completed.
Iteration 571 / 1000 completed.
Iteration 572 / 1000 completed.
Iteration 573 / 1000 completed.
Iteration 574 / 1000 completed.
Iteration 575 / 1000 completed.
Iteration 576 / 1000 completed.
Iteration 577 / 1000 completed.
Iteration 578 / 1000 completed.
Iteration 579 / 1000 completed.
Iteration 580 / 1000 completed.
Iteration 581 / 1000 completed.
Iteration 582 / 1000 completed.
Iteration 583 / 1000 completed.
Iteration 584 / 1000 completed.
Iteration 585 / 1000 completed.
Iteration 586 / 1000 completed.
Iteration 587 / 1000 completed.
Iteration 588 / 1000 completed.
Iteration 589 / 1000 completed.
Iteration 590 / 1000 completed.
Iteration 591 / 1000 completed.
Iteration 592 / 1000 completed.
Iteration 593 / 1000 completed.
Iteration 594 / 1000 completed.
Iteration 595 / 1000 completed.
Iteration 596 / 1000 completed.
Iteration 597 / 1000 completed.
Iteration 598 / 1000 completed.
Iteration 599 / 1000 completed.
Iteration 600 / 1000 completed.
Iteration 601 / 1000 completed.
Iteration 602 / 1000 completed.
Iteration 603 / 1000 completed.
Iteration 604 / 1000 completed.
Iteration 605 / 1000 completed.
Iteration 606 / 1000 completed.
Iteration 607 / 1000 completed.
Iteration 608 / 1000 completed.
Iteration 609 / 1000 completed.
Iteration 610 / 1000 completed.
Iteration 611 / 1000 completed.
Iteration 612 / 1000 completed.
Iteration 613 / 1000 completed.
Iteration 614 / 1000 completed.
Iteration 615 / 1000 completed.
Iteration 616 / 1000 completed.
Iteration 617 / 1000 completed.
Iteration 618 / 1000 completed.
Iteration 619 / 1000 completed.
Iteration 620 / 1000 completed.
Iteration 621 / 1000 completed.
Iteration 622 / 1000 completed.
Iteration 623 / 1000 completed.
Iteration 624 / 1000 completed.
Iteration 625 / 1000 completed.
Iteration 626 / 1000 completed.
Iteration 627 / 1000 completed.
Iteration 628 / 1000 completed.
Iteration 629 / 1000 completed.
Iteration 630 / 1000 completed.
Iteration 631 / 1000 completed.
Iteration 632 / 1000 completed.
Iteration 633 / 1000 completed.
Iteration 634 / 1000 completed.
Iteration 635 / 1000 completed.
Iteration 636 / 1000 completed.
Iteration 637 / 1000 completed.
Iteration 638 / 1000 completed.
Iteration 639 / 1000 completed.
Iteration 640 / 1000 completed.
Iteration 641 / 1000 completed.
Iteration 642 / 1000 completed.
Iteration 643 / 1000 completed.
Iteration 644 / 1000 completed.
Iteration 645 / 1000 completed.
Iteration 646 / 1000 completed.
Iteration 647 / 1000 completed.
Iteration 648 / 1000 completed.
Iteration 649 / 1000 completed.
Iteration 650 / 1000 completed.
Iteration 651 / 1000 completed.
Iteration 652 / 1000 completed.
Iteration 653 / 1000 completed.
Iteration 654 / 1000 completed.
Iteration 655 / 1000 completed.
Iteration 656 / 1000 completed.
Iteration 657 / 1000 completed.
Iteration 658 / 1000 completed.
Iteration 659 / 1000 completed.
Iteration 660 / 1000 completed.
Iteration 661 / 1000 completed.
Iteration 662 / 1000 completed.
Iteration 663 / 1000 completed.
Iteration 664 / 1000 completed.
Iteration 665 / 1000 completed.
Iteration 666 / 1000 completed.
Iteration 667 / 1000 completed.
Iteration 668 / 1000 completed.
Iteration 669 / 1000 completed.
Iteration 670 / 1000 completed.
Iteration 671 / 1000 completed.
Iteration 672 / 1000 completed.
Iteration 673 / 1000 completed.
Iteration 674 / 1000 completed.
Iteration 675 / 1000 completed.
Iteration 676 / 1000 completed.
Iteration 677 / 1000 completed.
Iteration 678 / 1000 completed.
Iteration 679 / 1000 completed.
Iteration 680 / 1000 completed.
Iteration 681 / 1000 completed.
Iteration 682 / 1000 completed.
Iteration 683 / 1000 completed.
Iteration 684 / 1000 completed.
Iteration 685 / 1000 completed.
Iteration 686 / 1000 completed.
Iteration 687 / 1000 completed.
Iteration 688 / 1000 completed.
Iteration 689 / 1000 completed.
Iteration 690 / 1000 completed.
Iteration 691 / 1000 completed.
Iteration 692 / 1000 completed.
Iteration 693 / 1000 completed.
Iteration 694 / 1000 completed.
Iteration 695 / 1000 completed.
Iteration 696 / 1000 completed.
Iteration 697 / 1000 completed.
Iteration 698 / 1000 completed.
Iteration 699 / 1000 completed.
Iteration 700 / 1000 completed.
Iteration 701 / 1000 completed.
Iteration 702 / 1000 completed.
Iteration 703 / 1000 completed.
Iteration 704 / 1000 completed.
Iteration 705 / 1000 completed.
Iteration 706 / 1000 completed.
Iteration 707 / 1000 completed.
Iteration 708 / 1000 completed.
Iteration 709 / 1000 completed.
Iteration 710 / 1000 completed.
Iteration 711 / 1000 completed.
Iteration 712 / 1000 completed.
Iteration 713 / 1000 completed.
Iteration 714 / 1000 completed.
Iteration 715 / 1000 completed.
Iteration 716 / 1000 completed.
Iteration 717 / 1000 completed.
Iteration 718 / 1000 completed.
Iteration 719 / 1000 completed.
Iteration 720 / 1000 completed.
Iteration 721 / 1000 completed.
Iteration 722 / 1000 completed.
Iteration 723 / 1000 completed.
Iteration 724 / 1000 completed.
Iteration 725 / 1000 completed.
Iteration 726 / 1000 completed.
Iteration 727 / 1000 completed.
Iteration 728 / 1000 completed.
Iteration 729 / 1000 completed.
Iteration 730 / 1000 completed.
Iteration 731 / 1000 completed.
Iteration 732 / 1000 completed.
Iteration 733 / 1000 completed.
Iteration 734 / 1000 completed.
Iteration 735 / 1000 completed.
Iteration 736 / 1000 completed.
Iteration 737 / 1000 completed.
Iteration 738 / 1000 completed.
Iteration 739 / 1000 completed.
Iteration 740 / 1000 completed.
Iteration 741 / 1000 completed.
Iteration 742 / 1000 completed.
Iteration 743 / 1000 completed.
Iteration 744 / 1000 completed.
Iteration 745 / 1000 completed.
Iteration 746 / 1000 completed.
Iteration 747 / 1000 completed.
Iteration 748 / 1000 completed.
Iteration 749 / 1000 completed.
Iteration 750 / 1000 completed.
Iteration 751 / 1000 completed.
Iteration 752 / 1000 completed.
Iteration 753 / 1000 completed.
Iteration 754 / 1000 completed.
Iteration 755 / 1000 completed.
Iteration 756 / 1000 completed.
Iteration 757 / 1000 completed.
Iteration 758 / 1000 completed.
Iteration 759 / 1000 completed.
Iteration 760 / 1000 completed.
Iteration 761 / 1000 completed.
Iteration 762 / 1000 completed.
Iteration 763 / 1000 completed.
Iteration 764 / 1000 completed.
Iteration 765 / 1000 completed.
Iteration 766 / 1000 completed.
Iteration 767 / 1000 completed.
Iteration 768 / 1000 completed.
Iteration 769 / 1000 completed.
Iteration 770 / 1000 completed.
Iteration 771 / 1000 completed.
Iteration 772 / 1000 completed.
Iteration 773 / 1000 completed.
Iteration 774 / 1000 completed.
Iteration 775 / 1000 completed.
Iteration 776 / 1000 completed.
Iteration 777 / 1000 completed.
Iteration 778 / 1000 completed.
Iteration 779 / 1000 completed.
Iteration 780 / 1000 completed.
Iteration 781 / 1000 completed.
Iteration 782 / 1000 completed.
Iteration 783 / 1000 completed.
Iteration 784 / 1000 completed.
Iteration 785 / 1000 completed.
Iteration 786 / 1000 completed.
Iteration 787 / 1000 completed.
Iteration 788 / 1000 completed.
Iteration 789 / 1000 completed.
Iteration 790 / 1000 completed.
Iteration 791 / 1000 completed.
Iteration 792 / 1000 completed.
Iteration 793 / 1000 completed.
Iteration 794 / 1000 completed.
Iteration 795 / 1000 completed.
Iteration 796 / 1000 completed.
Iteration 797 / 1000 completed.
Iteration 798 / 1000 completed.
Iteration 799 / 1000 completed.
Iteration 800 / 1000 completed.
Iteration 801 / 1000 completed.
Iteration 802 / 1000 completed.
Iteration 803 / 1000 completed.
Iteration 804 / 1000 completed.
Iteration 805 / 1000 completed.
Iteration 806 / 1000 completed.
Iteration 807 / 1000 completed.
Iteration 808 / 1000 completed.
Iteration 809 / 1000 completed.
Iteration 810 / 1000 completed.
Iteration 811 / 1000 completed.
Iteration 812 / 1000 completed.
Iteration 813 / 1000 completed.
Iteration 814 / 1000 completed.
Iteration 815 / 1000 completed.
Iteration 816 / 1000 completed.
Iteration 817 / 1000 completed.
Iteration 818 / 1000 completed.
Iteration 819 / 1000 completed.
Iteration 820 / 1000 completed.
Iteration 821 / 1000 completed.
Iteration 822 / 1000 completed.
Iteration 823 / 1000 completed.
Iteration 824 / 1000 completed.
Iteration 825 / 1000 completed.
Iteration 826 / 1000 completed.
Iteration 827 / 1000 completed.
Iteration 828 / 1000 completed.
Iteration 829 / 1000 completed.
Iteration 830 / 1000 completed.
Iteration 831 / 1000 completed.
Iteration 832 / 1000 completed.
Iteration 833 / 1000 completed.
Iteration 834 / 1000 completed.
Iteration 835 / 1000 completed.
Iteration 836 / 1000 completed.
Iteration 837 / 1000 completed.
Iteration 838 / 1000 completed.
Iteration 839 / 1000 completed.
Iteration 840 / 1000 completed.
Iteration 841 / 1000 completed.
Iteration 842 / 1000 completed.
Iteration 843 / 1000 completed.
Iteration 844 / 1000 completed.
Iteration 845 / 1000 completed.
Iteration 846 / 1000 completed.
Iteration 847 / 1000 completed.
Iteration 848 / 1000 completed.
Iteration 849 / 1000 completed.
Iteration 850 / 1000 completed.
Iteration 851 / 1000 completed.
Iteration 852 / 1000 completed.
Iteration 853 / 1000 completed.
Iteration 854 / 1000 completed.
Iteration 855 / 1000 completed.
Iteration 856 / 1000 completed.
Iteration 857 / 1000 completed.
Iteration 858 / 1000 completed.
Iteration 859 / 1000 completed.
Iteration 860 / 1000 completed.
Iteration 861 / 1000 completed.
Iteration 862 / 1000 completed.
Iteration 863 / 1000 completed.
Iteration 864 / 1000 completed.
Iteration 865 / 1000 completed.
Iteration 866 / 1000 completed.
Iteration 867 / 1000 completed.
Iteration 868 / 1000 completed.
Iteration 869 / 1000 completed.
Iteration 870 / 1000 completed.
Iteration 871 / 1000 completed.
Iteration 872 / 1000 completed.
Iteration 873 / 1000 completed.
Iteration 874 / 1000 completed.
Iteration 875 / 1000 completed.
Iteration 876 / 1000 completed.
Iteration 877 / 1000 completed.
Iteration 878 / 1000 completed.
Iteration 879 / 1000 completed.
Iteration 880 / 1000 completed.
Iteration 881 / 1000 completed.
Iteration 882 / 1000 completed.
Iteration 883 / 1000 completed.
Iteration 884 / 1000 completed.
Iteration 885 / 1000 completed.
Iteration 886 / 1000 completed.
Iteration 887 / 1000 completed.
Iteration 888 / 1000 completed.
Iteration 889 / 1000 completed.
Iteration 890 / 1000 completed.
Iteration 891 / 1000 completed.
Iteration 892 / 1000 completed.
Iteration 893 / 1000 completed.
Iteration 894 / 1000 completed.
Iteration 895 / 1000 completed.
Iteration 896 / 1000 completed.
Iteration 897 / 1000 completed.
Iteration 898 / 1000 completed.
Iteration 899 / 1000 completed.
Iteration 900 / 1000 completed.
Iteration 901 / 1000 completed.
Iteration 902 / 1000 completed.
Iteration 903 / 1000 completed.
Iteration 904 / 1000 completed.
Iteration 905 / 1000 completed.
Iteration 906 / 1000 completed.
Iteration 907 / 1000 completed.
Iteration 908 / 1000 completed.
Iteration 909 / 1000 completed.
Iteration 910 / 1000 completed.
Iteration 911 / 1000 completed.
Iteration 912 / 1000 completed.
Iteration 913 / 1000 completed.
Iteration 914 / 1000 completed.
Iteration 915 / 1000 completed.
Iteration 916 / 1000 completed.
Iteration 917 / 1000 completed.
Iteration 918 / 1000 completed.
Iteration 919 / 1000 completed.
Iteration 920 / 1000 completed.
Iteration 921 / 1000 completed.
Iteration 922 / 1000 completed.
Iteration 923 / 1000 completed.
Iteration 924 / 1000 completed.
Iteration 925 / 1000 completed.
Iteration 926 / 1000 completed.
Iteration 927 / 1000 completed.
Iteration 928 / 1000 completed.
Iteration 929 / 1000 completed.
Iteration 930 / 1000 completed.
Iteration 931 / 1000 completed.
Iteration 932 / 1000 completed.
Iteration 933 / 1000 completed.
Iteration 934 / 1000 completed.
Iteration 935 / 1000 completed.
Iteration 936 / 1000 completed.
Iteration 937 / 1000 completed.
Iteration 938 / 1000 completed.
Iteration 939 / 1000 completed.
Iteration 940 / 1000 completed.
Iteration 941 / 1000 completed.
Iteration 942 / 1000 completed.
Iteration 943 / 1000 completed.
Iteration 944 / 1000 completed.
Iteration 945 / 1000 completed.
Iteration 946 / 1000 completed.
Iteration 947 / 1000 completed.
Iteration 948 / 1000 completed.
Iteration 949 / 1000 completed.
Iteration 950 / 1000 completed.
Iteration 951 / 1000 completed.
Iteration 952 / 1000 completed.
Iteration 953 / 1000 completed.
Iteration 954 / 1000 completed.
Iteration 955 / 1000 completed.
Iteration 956 / 1000 completed.
Iteration 957 / 1000 completed.
Iteration 958 / 1000 completed.
Iteration 959 / 1000 completed.
Iteration 960 / 1000 completed.
Iteration 961 / 1000 completed.
Iteration 962 / 1000 completed.
Iteration 963 / 1000 completed.
Iteration 964 / 1000 completed.
Iteration 965 / 1000 completed.
Iteration 966 / 1000 completed.
Iteration 967 / 1000 completed.
Iteration 968 / 1000 completed.
Iteration 969 / 1000 completed.
Iteration 970 / 1000 completed.
Iteration 971 / 1000 completed.
Iteration 972 / 1000 completed.
Iteration 973 / 1000 completed.
Iteration 974 / 1000 completed.
Iteration 975 / 1000 completed.
Iteration 976 / 1000 completed.
Iteration 977 / 1000 completed.
Iteration 978 / 1000 completed.
Iteration 979 / 1000 completed.
Iteration 980 / 1000 completed.
Iteration 981 / 1000 completed.
Iteration 982 / 1000 completed.
Iteration 983 / 1000 completed.
Iteration 984 / 1000 completed.
Iteration 985 / 1000 completed.
Iteration 986 / 1000 completed.
Iteration 987 / 1000 completed.
Iteration 988 / 1000 completed.
Iteration 989 / 1000 completed.
Iteration 990 / 1000 completed.
Iteration 991 / 1000 completed.
Iteration 992 / 1000 completed.
Iteration 993 / 1000 completed.
Iteration 994 / 1000 completed.
Iteration 995 / 1000 completed.
Iteration 996 / 1000 completed.
Iteration 997 / 1000 completed.
Iteration 998 / 1000 completed.
Iteration 999 / 1000 completed.
Iteration 1000 / 1000 completed.
> 
> 
> # Plot results using final MCMC samples
> plot_results(mh_result$M2_w1, mh_result$M2_alpha1, mh_result$M2_beta1,
+              mh_result$M2_delta1, mh_result$M2_omega1,
+              mh_result$M2_w2, mh_result$M2_alpha2, mh_result$M2_beta2,
+              mh_result$M2_delta2, mh_result$M2_omega2,
+              xx, xx_true, yy_true)
null device 
          1 
Warning message:
In plot_results(mh_result$M2_w1, mh_result$M2_alpha1, mh_result$M2_beta1,  :
  Invalid scale parameter for component 1: -0.469318010196124
> 
> proc.time()
   user  system elapsed 
 213.43    6.25   28.09