tests/testthat/_snaps/Model-class.md

In Roche/crmPack: Object-Oriented Implementation of CRM Designs

MCMC computes correct values for LogisticNormal model

Code
  result@data
Output
  $alpha0
  [1] -1.955379 -1.955379 -1.955379 -2.325551

  $alpha1
  [1] 1.450219 1.450219 1.450219 1.059415

MCMC computes correct values for LogisticNormal model and empty data

Code
  result@data
Output
  $alpha0
  [1]  0.3948748  0.4142195  1.5222466 -1.5773076

  $alpha1
  [1]  0.9876408 -0.5152260 -0.6580012  1.4327070

MCMC computes correct values for LogisticLogNormal model

Code
  result@data
Output
  $alpha0
  [1] -1.296799 -1.296799 -1.296799 -1.680008

  $alpha1
  [1] 0.975694 0.975694 0.975694 0.651047

MCMC computes correct values for LogisticLogNormal model and empty data

Code
  result@data
Output
  $alpha0
  [1]  0.3948748  0.4142195  1.5222466 -1.5773076

  $alpha1
  [1] 2.6848927 0.5973656 0.5178854 4.1900261

MCMC computes correct values for LogisticLogNormalSub model

Code
  result@data
Output
  $alpha0
  [1] -2.965073 -2.965073 -2.965073 -2.965073

  $alpha1
  [1] 0.00974706 0.00974706 0.00974706 0.00974706

MCMC computes correct values for LogisticLogNormalSub model and empty data

Code
  result@data
Output
  $alpha0
  [1] -0.4517289 -2.2094272 -3.6537293  2.3394694

  $alpha1
  [1] 1745.91005   89.83612  619.23866   82.33520

MCMC computes correct values for ProbitLogNormal model

Code
  result@data
Output
  $alpha0
  [1] -1.709947 -1.709947 -1.709947 -1.988543

  $alpha1
  [1] 0.4101207 0.4101207 0.4101207 0.3056156

MCMC computes correct values for ProbitLogNormal model and empty data

Code
  result@data
Output
  $alpha0
  [1]  0.3948748  0.4142195  1.5222466 -1.5773076

  $alpha1
  [1] 2.6848927 0.5973656 0.5178854 4.1900261

MCMC computes correct values for ProbitLogNormalRel model

Code
  result@data
Output
  $alpha0
  [1] -0.8353588 -0.8353588 -0.8353588 -0.8353588

  $alpha1
  [1] 0.02201234 0.02201234 0.02201234 0.02201234

MCMC computes correct values for ProbitLogNormalRel model and empty data

Code
  result@data
Output
  $alpha0
  [1]  0.3948748  0.4142195  1.5222466 -1.5773076

  $alpha1
  [1] 2.6848927 0.5973656 0.5178854 4.1900261

MCMC computes correct values for LogisticLogNormalGrouped model

Code
  result@data
Output
  $alpha0
  [1] -1.848124 -1.848124 -2.195992 -2.195992

  $alpha1
  [1] 0.08165851 0.08165851 0.31496464 0.31496464

  $delta0
  [1] 0.1263061 0.1263061 0.1112928 0.1112928

  $delta1
  [1] 0.06280840 0.06280840 0.09082656 0.09082656

MCMC computes correct values for LogisticLogNormalGrouped model and empty data

Code
  result@data
Output
  $alpha0
  [1] -0.7258644 -2.3268647  1.2196258  0.9135351

  $alpha1
  [1] 1.3236188 0.2700508 1.8853673 0.2820018

  $delta0
  [1] 3.4298468 2.0426447 0.7105671 0.2355593

  $delta1
  [1] 0.3087583 0.2956680 3.1296169 1.7350122

MCMC computes correct values for LogisticKadane model

Code
  result@data
Output
  $gamma
  [1] 86.15670 87.65260 87.16488 77.19426

  $rho0
  [1] 0.15583729 0.04234820 0.01397354 0.05035614

MCMC computes correct values for LogisticKadane model and empty data

Code
  result@data
Output
  $gamma
  [1]  53.83622  41.13470 126.19369 137.71755

  $rho0
  [1] 0.12280089 0.29646860 0.02038947 0.12675423

MCMC computes correct values for LogisticKadaneBetaGamma model

Code
  result@data
Output
  $gamma
  [1] 3.430834 5.279510 5.967748 6.643819

  $rho0
  [1] 0.08783371 0.01565705 0.06100013 0.05312860

MCMC computes correct values for LogisticKadaneBetaGamma model and empty data

Code
  result@data
Output
  $gamma
  [1] 3.091929 5.915646 6.032858 3.804350

  $rho0
  [1] 0.0052913065 0.0004292867 0.0925659424 0.0113490098

MCMC computes correct values for LogisticNormalMixture model

Code
  result@data
Output
  $alpha0
  [1] -6.222009 -6.222009 -6.029068 -6.029068

  $alpha1
  [1] 1.189508 1.189508 1.245205 1.245205

  $w
  [1] 0.10161451 0.02275232 0.10965042 0.06311245

MCMC computes correct values for LogisticNormalMixture model and empty data

Code
  result@data
Output
  $alpha0
  [1] -2.6619787 -2.7232446 -0.1032151 -3.0446598

  $alpha1
  [1] 6.719249 3.381709 5.316616 7.827070

  $w
  [1] 0.35619738 0.01377870 0.04005730 0.06812647

MCMC computes correct values for LogisticNormalFixedMixture model

Code
  result@data
Output
  $alpha0
  [1] -1.482805 -1.482805 -0.716882 -1.132819

  $alpha1
  [1] 2.0964447 2.0964447 0.5026685 1.8211377

Code
  result_log_normal@data
Output
  $alpha0
  [1] -1.481202 -1.481202 -1.481202 -1.866221

  $alpha1
  [1] 0.9751036 0.9751036 0.9751036 3.3043857

MCMC computes correct values for LogisticNormalFixedMixture model and empty data

Code
  result@data
Output
  $alpha0
  [1] -0.6582401 -0.7511036 -0.9194915 -2.0027877

  $alpha1
  [1]  1.1678645  1.7176934 -0.2655944  2.1668228

Code
  result_log_normal@data
Output
  $alpha0
  [1] -0.6582401 -0.7511036 -0.9194915 -2.0027877

  $alpha1
  [1] 3.2151195 5.5716621 0.7667501 8.7305018

MCMC computes correct values for LogisticLogNormalMixture model

Code
  result@data
Output
  $alpha0
            [,1]       [,2]
  [1,] -0.939651 -0.6725854
  [2,] -0.939651 -0.6725854
  [3,] -2.371683 -1.2819273
  [4,] -2.371683 -1.0844256

  $alpha1
            [,1]      [,2]
  [1,] 0.4468436 0.7469892
  [2,] 0.4468436 0.7469892
  [3,] 0.3986522 1.1817339
  [4,] 0.3986522 0.6291422

  $comp
  [1] 1 1 1 1

MCMC computes correct values for LogisticLogNormalMixture model and empty data

Code
  result@data
Output
  $alpha0
              [,1]       [,2]
  [1,] -0.01379598 -1.1751964
  [2,] -0.24014099 -1.2185179
  [3,]  0.63412268  0.3440982
  [4,]  0.55101443 -0.4627752

  $alpha1
            [,1]      [,2]
  [1,] 3.8947301 0.5462308
  [2,] 0.7340741 5.3107587
  [3,] 8.5071807 2.2575712
  [4,] 6.2133610 2.8849464

  $comp
  [1] 1 2 1 1

MCMC computes correct values for DualEndpointRW model (fixed params)

Code
  result@data
Output
  $betaW
           [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
  [1,] 37.62747 50.53208 54.27004 56.27200 58.47643 57.91182 58.26476 58.36793
  [2,] 37.55166 49.87282 54.34154 56.56615 58.78051 59.34977 59.99789 59.60488
  [3,] 37.03832 50.09082 54.26849 56.14897 58.40198 57.73258 57.91062 57.55694
  [4,] 37.71848 50.53490 54.83729 56.45795 58.44911 57.76732 57.30030 57.33680
           [,9]    [,10]    [,11]    [,12]    [,13]
  [1,] 58.38942 58.91251 59.30705 59.39655 58.68454
  [2,] 59.82375 60.57548 60.66172 61.18413 61.29934
  [3,] 57.93404 58.18794 58.54076 57.76713 57.42084
  [4,] 56.79581 55.70361 56.06342 56.13879 56.62076

  $betaZ
            [,1]       [,2]
  [1,] -2.698892 0.04732726
  [2,] -2.698892 0.04732726
  [3,] -2.698892 0.04732726
  [4,] -2.698892 0.04732726

  $delta
           [,1]     [,2]     [,3]     [,4]       [,5]       [,6]        [,7]
  [1,] 12.90460 3.737966 2.001957 2.204435 -0.5646161  0.3529426  0.10316457
  [2,] 12.32115 4.468726 2.224605 2.214358  0.5692620  0.6481230 -0.39301672
  [3,] 13.05249 4.177676 1.880479 2.253004 -0.6693999  0.1780457 -0.35368323
  [4,] 12.81642 4.302391 1.620659 1.991164 -0.6817942 -0.4670186  0.03649731
              [,8]       [,9]      [,10]       [,11]      [,12]
  [1,]  0.02149088  0.5230983 0.39453828  0.08950038 -0.7120159
  [2,]  0.21887260  0.7517329 0.08623658  0.52240780  0.1152141
  [3,]  0.37710644  0.2538971 0.35282056 -0.77363433 -0.3462910
  [4,] -0.54098389 -1.0922007 0.35981049  0.07536469  0.4819734

Code
  result_log_dose@data
Output
  $betaW
           [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
  [1,] 37.92798 49.81859 53.72168 55.92768 57.83541 57.72386 57.84462 57.97807
  [2,] 38.34152 50.28294 54.09302 55.53691 57.86823 58.31979 58.15104 57.96073
  [3,] 37.34457 50.82862 53.93501 56.76424 59.23104 58.80072 59.17973 59.05368
  [4,] 37.11602 50.46837 54.89170 56.69901 59.07928 59.73691 59.23760 58.83368
           [,9]    [,10]    [,11]    [,12]    [,13]
  [1,] 57.53570 57.29761 57.65197 58.74994 58.84660
  [2,] 57.95901 57.24585 57.26554 57.58131 56.62224
  [3,] 58.12606 58.19137 58.30478 58.61300 58.96139
  [4,] 58.56342 58.60440 58.48189 57.86828 58.11178

  $betaZ
            [,1]      [,2]
  [1,] -2.261079 0.4562101
  [2,] -2.261079 0.4562101
  [3,] -2.261079 0.4562101
  [4,] -2.261079 0.4562101

  $delta
           [,1]     [,2]     [,3]     [,4]       [,5]       [,6]       [,7]
  [1,] 11.89062 3.903092 2.205994 1.907736 -0.1115568  0.1207626  0.1334544
  [2,] 11.94142 3.810082 1.443886 2.331325  0.4515518 -0.1687490 -0.1903037
  [3,] 13.48406 3.106383 2.829239 2.466797 -0.4303261  0.3790168 -0.1260518
  [4,] 13.35236 4.423327 1.807311 2.380265  0.6576365 -0.4993145 -0.4039183
               [,8]        [,9]       [,10]      [,11]       [,12]
  [1,] -0.442370799 -0.23809251  0.35435512  1.0979752  0.09665914
  [2,] -0.001718259 -0.71316723  0.01969286  0.3157687 -0.95906729
  [3,] -0.927618955  0.06530940  0.11341143  0.3082144  0.34839812
  [4,] -0.270257937  0.04097757 -0.12251327 -0.6136042  0.24349751

MCMC computes correct values for DualEndpointRW model

Code
  result@data
Output
  $betaW
           [,1]     [,2]     [,3]     [,4]     [,5]     [,6]      [,7]      [,8]
  [1,] 45.89789 53.05191 66.44612 85.04402 66.94232 87.42737 102.96013 115.47159
  [2,] 43.82645 51.18066 68.66893 66.28827 75.02211 91.26635  94.27013  98.26580
  [3,] 40.50667 60.57052 60.16608 60.15285 60.98330 70.77219  77.56569  80.82604
  [4,] 35.80946 45.03871 44.58478 47.75190 63.52629 60.99178  69.46006  70.50504
            [,9]     [,10]     [,11]     [,12]     [,13]
  [1,] 120.82685 117.93494 118.17376 120.54576 120.47496
  [2,]  94.80926  99.03741  97.18131  99.02085 101.09911
  [3,]  83.35307  83.77029  93.51790 100.90913  93.32428
  [4,]  81.42350  88.23564  86.93943  70.28536  75.62961

  $betaZ
            [,1]       [,2]
  [1,] -2.089533 0.04811597
  [2,] -2.089533 0.04811597
  [3,] -1.359775 0.04435814
  [4,] -1.186205 0.03448847

  $delta
            [,1]       [,2]        [,3]        [,4]      [,5]      [,6]      [,7]
  [1,]  7.154014 13.3942071 18.59790278 -18.1017033 20.485050 15.532760 12.511462
  [2,]  7.354209 17.4882708 -2.38065917   8.7338402 16.244236  3.003776  3.995673
  [3,] 20.063851 -0.4044338 -0.01323763   0.8304503  9.788897  6.793498  3.260349
  [4,]  9.229245 -0.4539265  3.16711715  15.7743944 -2.534508  8.468280  1.044973
            [,8]       [,9]      [,10]      [,11]       [,12]
  [1,]  5.355268 -2.8919167  0.2388175   2.372007 -0.07080061
  [2,] -3.456543  4.2281569 -1.8561056   1.839540  2.07825773
  [3,]  2.527025  0.4172247  9.7476132   7.391231 -7.58485567
  [4,] 10.918465  6.8121410 -1.2962176 -16.654071  5.34425168

  $precBetaW
  [1] 0.2157918 0.1154706 0.2471180 0.4008307

  $precW
  [1] 0.001869948 0.001425776 0.001836746 0.002610231

  $rho
  [1]  0.5026811  0.7772394  0.3537694 -0.3502122

Code
  result_log_dose@data
Output
  $betaW
           [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
  [1,] 40.45741 69.56808 59.55233 79.13340 63.91507 46.05233 42.63871 28.26168
  [2,] 53.84549 62.20866 56.39624 64.75127 65.70342 68.02558 67.80433 57.85620
  [3,] 52.91517 52.72525 45.69151 41.54785 46.40927 54.91094 68.43927 74.22991
  [4,] 48.91754 59.63834 54.31261 54.45322 59.50827 62.02394 61.49123 56.93567
           [,9]    [,10]    [,11]    [,12]     [,13]
  [1,] 47.28609 62.44211 60.13589 52.29069  36.80784
  [2,] 57.93962 51.15583 42.19151 62.76152  74.34565
  [3,] 65.24259 71.93468 76.79675 91.77726 106.39135
  [4,] 55.56045 44.64745 48.05770 45.45153  44.58008

  $betaZ
            [,1]      [,2]
  [1,] -1.828305 0.4406148
  [2,] -1.828305 0.4406148
  [3,] -1.828305 0.4406148
  [4,] -1.828305 0.4406148

  $delta
             [,1]       [,2]       [,3]        [,4]       [,5]      [,6]
  [1,] 29.1106742 -10.015749 19.5810705 -15.2183359 -17.862738 -3.413620
  [2,]  8.3631685  -5.812421  8.3550353   0.9521496   2.322155 -0.221247
  [3,] -0.1899186  -7.033738 -4.1436554   4.8614117   8.501676 13.528324
  [4,] 10.7207979  -5.325733  0.1406142   5.0550447   2.515677 -0.532718
             [,7]       [,8]       [,9]     [,10]     [,11]       [,12]
  [1,] -14.377030 19.0244134  15.156014 -2.306212 -7.845210 -15.4828452
  [2,]  -9.948125  0.0834199  -6.783791 -8.964324 20.570011  11.5841288
  [3,]   5.790648 -8.9873284   6.692096  4.862073 14.980508  14.6140900
  [4,]  -4.555561 -1.3752176 -10.912995  3.410245 -2.606165  -0.8714534

  $precBetaW
  [1] 0.08906158 0.24302585 0.33020328 1.20661065

  $precW
  [1] 0.001839355 0.001311216 0.001120640 0.001377583

  $rho
  [1] 0.5169636 0.5470295 0.5800301 0.2189836

MCMC computes correct values for DualEndpointRW model with RW2

Code
  result@data
Output
  $betaW
           [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
  [1,] 30.32737 48.49422 58.12168 63.09706 64.87885 67.37127 70.76217 74.31984
  [2,] 30.68617 48.82946 58.57499 63.75225 66.06438 67.58474 69.86412 72.22545
  [3,] 29.51336 47.68636 57.57308 62.23465 64.54077 66.05562 67.43177 68.63151
  [4,] 28.90676 47.48893 57.46769 62.28280 64.73978 66.52814 67.11305 67.84362
           [,9]    [,10]    [,11]    [,12]    [,13]
  [1,] 77.56470 81.35809 84.81019 88.71579 92.43348
  [2,] 74.29452 75.99158 78.56877 80.40367 82.50585
  [3,] 69.19719 69.23396 69.27538 70.23891 71.68408
  [4,] 68.71528 69.49382 70.05987 71.12960 72.35257

  $betaZ
            [,1]       [,2]
  [1,] -2.242592 0.02871628
  [2,] -2.242592 0.02871628
  [3,] -1.757175 0.01148692
  [4,] -2.519110 0.02719582

  $delta
           [,1]     [,2]     [,3]     [,4]     [,5]      [,6]      [,7]      [,8]
  [1,] 18.16685 9.627453 4.975383 1.781788 2.492426 3.3909000 3.5576683 3.2448548
  [2,] 18.14329 9.745526 5.177257 2.312133 1.520355 2.2793882 2.3613280 2.0690683
  [3,] 18.17300 9.886723 4.661568 2.306116 1.514857 1.3761477 1.1997432 0.5656786
  [4,] 18.58217 9.978765 4.815104 2.456983 1.788359 0.5849097 0.7305757 0.8716594
             [,9]      [,10]     [,11]    [,12]
  [1,] 3.79339028 3.45210107 3.9056064 3.717690
  [2,] 1.69705765 2.57719225 1.8349016 2.102176
  [3,] 0.03677205 0.04141517 0.9635334 1.445165
  [4,] 0.77853552 0.56605477 1.0697225 1.222977

Code
  result_log_dose@data
Output
  $betaW
           [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
  [1,] 32.06799 47.95605 56.67823 61.76641 64.46203 66.23528 67.76664 69.82522
  [2,] 31.06149 47.87684 57.51957 62.69833 64.88478 67.29730 69.23368 71.25285
  [3,] 30.64417 47.81794 57.95373 62.49984 64.41398 66.77028 69.94141 72.50805
  [4,] 30.92944 47.91420 57.55503 61.44452 63.95053 66.26510 69.00486 71.21268
           [,9]    [,10]    [,11]    [,12]    [,13]
  [1,] 71.52412 73.23032 75.16740 76.92437 77.97551
  [2,] 73.02142 74.99042 77.17108 79.94857 81.70110
  [3,] 75.22159 78.31403 81.79406 84.68121 87.11906
  [4,] 73.67730 77.13788 80.78734 84.90042 88.46065

  $betaZ
            [,1]      [,2]
  [1,] -1.713626 0.3298758
  [2,] -1.713626 0.3298758
  [3,] -1.713626 0.3298758
  [4,] -1.713626 0.3298758

  $delta
           [,1]      [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
  [1,] 15.88807  8.722172 5.088188 2.695613 1.773252 1.531360 2.058584 1.698901
  [2,] 16.81535  9.642729 5.178764 2.186447 2.412528 1.936375 2.019171 1.768566
  [3,] 17.17377 10.135783 4.546111 1.914142 2.356298 3.171129 2.566648 2.713536
  [4,] 16.98476  9.640834 3.889486 2.506015 2.314566 2.739757 2.207820 2.464623
           [,9]    [,10]    [,11]    [,12]
  [1,] 1.706193 1.937086 1.756969 1.051138
  [2,] 1.969008 2.180659 2.777487 1.752525
  [3,] 3.092439 3.480032 2.887149 2.437850
  [4,] 3.460580 3.649458 4.113079 3.560226

MCMC computes correct values for DualEndpointRW model (fixed params) with RW2

Code
  result@data
Output
  $betaW
           [,1]     [,2]     [,3]     [,4]     [,5]     [,6]       [,7]
  [1,] 18.85216 39.67886 47.07057 48.07592 51.00794 58.72187  85.124404
  [2,] 27.10613 50.07723 70.37563 75.52978 76.99950 98.84815 144.464291
  [3,] 36.32607 57.36845 66.66144 66.69730 66.57335 78.36139  86.006482
  [4,] 41.38790 59.86394 57.79825 49.72352 37.86309 29.64807   9.392605
             [,8]      [,9]     [,10]     [,11]     [,12]     [,13]
  [1,] 118.105268 162.45184 219.97200 268.67755 309.41663 362.00415
  [2,] 194.117924 245.30292 306.35259 357.38392 409.96526 484.58834
  [3,] 105.196845 120.46583 123.20083 144.56411 168.25483 196.37063
  [4,]  -8.795638 -19.00679 -26.19848 -33.73996 -44.42018 -52.28546

  $betaZ
            [,1]       [,2]
  [1,] -1.746666 0.02987438
  [2,] -1.746666 0.02987438
  [3,] -2.027601 0.05389445
  [4,] -2.027601 0.05389445

  $delta
           [,1]      [,2]        [,3]        [,4]      [,5]       [,6]      [,7]
  [1,] 20.82670  7.391708  1.00535093   2.9320194  7.713929  26.402533  32.98086
  [2,] 22.97110 20.298400  5.15415276   1.4697179 21.848654  45.616141  49.65363
  [3,] 21.04238  9.292997  0.03585864  -0.1239471 11.788036   7.645091  19.19036
  [4,] 18.47604 -2.065687 -8.07473412 -11.8604310 -8.215020 -20.255464 -18.18824
            [,8]      [,9]     [,10]     [,11]     [,12]
  [1,]  44.34657 57.520164 48.705541  40.73909 52.587513
  [2,]  51.18500 61.049671 51.031328  52.58134 74.623076
  [3,]  15.26898  2.735001 21.363282  23.69072 28.115799
  [4,] -10.21115 -7.191692 -7.541483 -10.68022 -7.865276

  $precBetaW
  [1] 0.1855873 0.1786592 0.2768435 0.2804265

  $precW
  [1] 0.001458471 0.002268202 0.002229415 0.001081834

  $rho
  [1] -0.08830602 -0.10450848  0.02152350 -0.49952930

Code
  result_log_dose@data
Output
  $betaW
           [,1]     [,2]     [,3]     [,4]     [,5]      [,6]      [,7]      [,8]
  [1,] 32.65073 45.37149 53.17293 60.92134 64.89773  60.86839  65.83215  91.20007
  [2,] 41.71826 50.53644 56.49855 60.33543 72.70581  80.62209  89.22751 107.63254
  [3,] 52.51722 61.24755 66.81737 78.38312 88.49090 102.76567 120.54927 140.07753
  [4,] 38.99661 44.11117 49.75168 53.38641 44.08817  41.25106  40.85960  40.18481
            [,9]     [,10]    [,11]     [,12]     [,13]
  [1,] 123.75019 168.70935 228.0514 277.76967 318.79019
  [2,] 132.71762 156.72918 176.7843 204.34387 229.29783
  [3,] 162.92285 181.26264 196.4747 209.61837 226.96651
  [4,]  36.96778  36.03074  40.9033  40.23324  34.39948

  $betaZ
            [,1]      [,2]
  [1,] -1.244735 0.3580311
  [2,] -1.244735 0.3580311
  [3,] -1.244735 0.3580311
  [4,] -1.244735 0.3580311

  $delta
            [,1]     [,2]      [,3]      [,4]      [,5]       [,6]      [,7]
  [1,] 12.720760 7.801443  7.748407  3.976392 -4.029339  4.9637631 25.367919
  [2,]  8.818174 5.962111  3.836877 12.370382  7.916288  8.6054137 18.405036
  [3,]  8.730322 5.569827 11.565747 10.107784 14.274765 17.7835979 19.528261
  [4,]  5.114557 5.640509  3.634724 -9.298232 -2.837108 -0.3914679 -0.674784
            [,8]       [,9]     [,10]     [,11]     [,12]
  [1,] 32.550120 44.9591628 59.342045 49.718269 41.020523
  [2,] 25.085077 24.0115558 20.055094 27.559603 24.953952
  [3,] 22.845327 18.3397859 15.212062 13.143671 17.348133
  [4,] -3.217033 -0.9370433  4.872564 -0.670057 -5.833764

  $precBetaW
  [1] 0.1556914 0.9038538 0.4926834 0.4274567

  $precW
  [1] 0.0010851812 0.0009989832 0.0009079423 0.0013589435

  $rho
  [1] 0.03156483 0.59525044 0.76578093 0.39911145

MCMC computes correct values for DualEndpointRW model (fixed params, empty data)

Code
  result@data
Output
  $betaW
              [,1]        [,2]        [,3]        [,4]        [,5]       [,6]
  [1,]  -725.86444  -726.45202  -727.25438  -727.37988  -727.15350  -726.8426
  [2,]    59.50632    59.87105    60.19924    60.03736    59.91361    60.8442
  [3,] -1189.08277 -1188.35626 -1188.09319 -1187.45043 -1187.16567 -1187.2551
  [4,]   733.96049   733.92209   734.07999   733.54951   733.38720   733.0808
              [,7]        [,8]        [,9]       [,10]       [,11]       [,12]
  [1,]  -728.00604  -727.64892  -728.30349  -728.91275  -728.57788  -728.72518
  [2,]    61.07043    60.82423    60.88966    60.80774    60.16618    59.42379
  [3,] -1187.23146 -1187.57822 -1187.63407 -1186.97750 -1187.02293 -1187.00853
  [4,]   732.82127   732.83005   733.90648   733.49621   734.25865   733.87371
             [,13]
  [1,]  -728.53491
  [2,]    60.09169
  [3,] -1187.30303
  [4,]   733.33003

  $betaZ
             [,1]      [,2]
  [1,]  1.2325156 3.5979690
  [2,] -0.5088103 2.4388200
  [3,]  0.5511301 0.3890276
  [4,] -1.0176569 6.0064475

  $delta
              [,1]       [,2]       [,3]       [,4]        [,5]        [,6]
  [1,] -0.58758645 -0.8023568 -0.1255032  0.2263828  0.31088974 -1.16343233
  [2,]  0.36472627  0.3281960 -0.1618874 -0.1237429  0.93059106  0.22622599
  [3,]  0.72651451  0.2630678  0.6427632  0.2847582 -0.08942704  0.02363321
  [4,] -0.03840379  0.1579038 -0.5304832 -0.1623071 -0.30636545 -0.25956601
              [,7]        [,8]        [,9]       [,10]       [,11]      [,12]
  [1,]  0.35712269 -0.65457266 -0.60925895  0.33486735 -0.14730031  0.1902790
  [2,] -0.24620332  0.06543675 -0.08192282 -0.64155989 -0.74238665  0.6678930
  [3,] -0.34676091 -0.05584225  0.65656612 -0.04543084  0.01439569 -0.2944954
  [4,]  0.00877637  1.07643177 -0.41026823  0.76244412 -0.38494281 -0.5436842

Code
  result_log_dose@data
Output
  $betaW
              [,1]        [,2]        [,3]        [,4]        [,5]       [,6]
  [1,]  -725.86444  -726.45202  -727.25438  -727.37988  -727.15350  -726.8426
  [2,]    59.50632    59.87105    60.19924    60.03736    59.91361    60.8442
  [3,] -1189.08277 -1188.35626 -1188.09319 -1187.45043 -1187.16567 -1187.2551
  [4,]   733.96049   733.92209   734.07999   733.54951   733.38720   733.0808
              [,7]        [,8]        [,9]       [,10]       [,11]       [,12]
  [1,]  -728.00604  -727.64892  -728.30349  -728.91275  -728.57788  -728.72518
  [2,]    61.07043    60.82423    60.88966    60.80774    60.16618    59.42379
  [3,] -1187.23146 -1187.57822 -1187.63407 -1186.97750 -1187.02293 -1187.00853
  [4,]   732.82127   732.83005   733.90648   733.49621   734.25865   733.87371
             [,13]
  [1,]  -728.53491
  [2,]    60.09169
  [3,] -1187.30303
  [4,]   733.33003

  $betaZ
             [,1]      [,2]
  [1,]  1.2325156 3.5979690
  [2,] -0.5088103 2.4388200
  [3,]  0.5511301 0.3890276
  [4,] -1.0176569 6.0064475

  $delta
              [,1]       [,2]       [,3]       [,4]        [,5]        [,6]
  [1,] -0.58758645 -0.8023568 -0.1255032  0.2263828  0.31088974 -1.16343233
  [2,]  0.36472627  0.3281960 -0.1618874 -0.1237429  0.93059106  0.22622599
  [3,]  0.72651451  0.2630678  0.6427632  0.2847582 -0.08942704  0.02363321
  [4,] -0.03840379  0.1579038 -0.5304832 -0.1623071 -0.30636545 -0.25956601
              [,7]        [,8]        [,9]       [,10]       [,11]      [,12]
  [1,]  0.35712269 -0.65457266 -0.60925895  0.33486735 -0.14730031  0.1902790
  [2,] -0.24620332  0.06543675 -0.08192282 -0.64155989 -0.74238665  0.6678930
  [3,] -0.34676091 -0.05584225  0.65656612 -0.04543084  0.01439569 -0.2944954
  [4,]  0.00877637  1.07643177 -0.41026823  0.76244412 -0.38494281 -0.5436842

MCMC computes correct values for DualEndpointRW model (empty data)

Code
  result@data
Output
  $betaW
            [,1]      [,2]      [,3]      [,4]      [,5]      [,6]      [,7]
  [1,]  359.6244  353.5719  361.3905  363.1210  388.1782  381.1582  394.1862
  [2,] -166.0710 -170.1050 -174.3420 -178.2492 -181.1457 -185.8399 -178.1500
  [3,]  797.2342  798.1979  797.8669  806.9230  810.8529  805.4040  809.6491
  [4,] -759.2283 -759.7527 -761.9904 -763.6605 -774.3127 -779.4910 -785.5917
            [,8]      [,9]     [,10]     [,11]     [,12]     [,13]
  [1,]  417.6262  424.6957  420.8803  433.2697  445.1613  463.9299
  [2,] -177.6674 -178.1936 -177.9295 -183.6934 -181.0219 -190.4457
  [3,]  809.2378  810.9288  805.2480  803.5099  800.2291  797.4494
  [4,] -792.4441 -788.4116 -794.1396 -789.3519 -799.8292 -806.0798

  $betaZ
             [,1]      [,2]
  [1,] -1.2185179 5.3107587
  [2,]  0.6612293 2.3265189
  [3,]  1.3024961 0.9456107
  [4,] -0.1997936 3.4073663

  $delta
            [,1]       [,2]      [,3]       [,4]      [,5]      [,6]       [,7]
  [1,] -6.052450  7.8185661  1.730534  25.057217 -7.020062 13.028053 23.4400128
  [2,] -4.033993 -4.2369787 -3.907213  -2.896468 -4.694228  7.689898  0.4826465
  [3,]  0.963691 -0.3309774  9.056099   3.929909 -5.448930  4.245138 -0.4112659
  [4,] -0.524353 -2.2377166 -1.670103 -10.652162 -5.178337 -6.100732 -6.8523290
             [,8]       [,9]     [,10]      [,11]     [,12]
  [1,]  7.0694985 -3.8154217 12.389377  11.891589 18.768582
  [2,] -0.5262785  0.2640938 -5.763913   2.671526 -9.423777
  [3,]  1.6909561 -5.6808259 -1.738110  -3.280799 -2.779634
  [4,]  4.0324298 -5.7279604  4.787642 -10.477247 -6.250576

  $precBetaW
  [1] 0.05922806 1.06397010 0.87200774 0.72858362

  $precW
  [1] 0.09111477 0.32887065 0.75368176 0.60727117

  $rho
  [1] -0.05563211  0.28652884 -0.41335638 -0.58841977

Code
  result_log_dose@data
Output
  $betaW
            [,1]      [,2]      [,3]      [,4]      [,5]      [,6]      [,7]
  [1,]  359.6244  353.5719  361.3905  363.1210  388.1782  381.1582  394.1862
  [2,] -166.0710 -170.1050 -174.3420 -178.2492 -181.1457 -185.8399 -178.1500
  [3,]  797.2342  798.1979  797.8669  806.9230  810.8529  805.4040  809.6491
  [4,] -759.2283 -759.7527 -761.9904 -763.6605 -774.3127 -779.4910 -785.5917
            [,8]      [,9]     [,10]     [,11]     [,12]     [,13]
  [1,]  417.6262  424.6957  420.8803  433.2697  445.1613  463.9299
  [2,] -177.6674 -178.1936 -177.9295 -183.6934 -181.0219 -190.4457
  [3,]  809.2378  810.9288  805.2480  803.5099  800.2291  797.4494
  [4,] -792.4441 -788.4116 -794.1396 -789.3519 -799.8292 -806.0798

  $betaZ
             [,1]      [,2]
  [1,] -1.2185179 5.3107587
  [2,]  0.6612293 2.3265189
  [3,]  1.3024961 0.9456107
  [4,] -0.1997936 3.4073663

  $delta
            [,1]       [,2]      [,3]       [,4]      [,5]      [,6]       [,7]
  [1,] -6.052450  7.8185661  1.730534  25.057217 -7.020062 13.028053 23.4400128
  [2,] -4.033993 -4.2369787 -3.907213  -2.896468 -4.694228  7.689898  0.4826465
  [3,]  0.963691 -0.3309774  9.056099   3.929909 -5.448930  4.245138 -0.4112659
  [4,] -0.524353 -2.2377166 -1.670103 -10.652162 -5.178337 -6.100732 -6.8523290
             [,8]       [,9]     [,10]      [,11]     [,12]
  [1,]  7.0694985 -3.8154217 12.389377  11.891589 18.768582
  [2,] -0.5262785  0.2640938 -5.763913   2.671526 -9.423777
  [3,]  1.6909561 -5.6808259 -1.738110  -3.280799 -2.779634
  [4,]  4.0324298 -5.7279604  4.787642 -10.477247 -6.250576

  $precBetaW
  [1] 0.05922806 1.06397010 0.87200774 0.72858362

  $precW
  [1] 0.09111477 0.32887065 0.75368176 0.60727117

  $rho
  [1] -0.05563211  0.28652884 -0.41335638 -0.58841977

MCMC computes correct values for DualEndpointRW model with RW2 (empty data)

Code
  result@data
Output
  $betaW
              [,1]      [,2]        [,3]       [,4]      [,5]      [,6]      [,7]
  [1,]  -725.86444  506.6512  1738.36443  2969.9522  4201.766  5433.891  6664.853
  [2,]    59.50632 -449.3040  -957.78604 -1466.4300 -1975.198 -2483.035 -2990.646
  [3,] -1189.08277 -637.9527   -86.55955   465.4764  1017.797  1570.028  2122.283
  [4,]   733.96049 -283.6964 -1301.19537 -2319.2248 -3337.417 -4355.915 -5374.672
            [,8]      [,9]     [,10]     [,11]      [,12]      [,13]
  [1,]  7896.172  9126.836 10356.891 11587.280  12817.523  14047.956
  [2,] -3498.503 -4006.294 -4514.168 -5022.683  -5531.941  -6040.530
  [3,]  2674.191  3226.043  3778.552  4331.016   4883.494   5435.677
  [4,] -6393.421 -7411.094 -8429.177 -9446.497 -10464.202 -11482.451

  $betaZ
             [,1]       [,2]
  [1,]  0.2803695  0.8392922
  [2,] -0.1084857  5.6376489
  [3,] -1.9441049 11.6238390
  [4,]  0.7928335  2.5173100

  $delta
             [,1]       [,2]       [,3]       [,4]       [,5]       [,6]
  [1,]  1232.5156  1231.7133  1231.5878  1231.8141  1232.1250  1230.9616
  [2,]  -508.8103  -508.4821  -508.6440  -508.7677  -507.8371  -507.6109
  [3,]   551.1301   551.3931   552.0359   552.3207   552.2312   552.2549
  [4,] -1017.6569 -1017.4990 -1018.0295 -1018.1918 -1018.4981 -1018.7577
             [,7]       [,8]       [,9]      [,10]      [,11]      [,12]
  [1,]  1231.3187  1230.6641  1230.0549  1230.3898  1230.2425  1230.4327
  [2,]  -507.8571  -507.7917  -507.8736  -508.5151  -509.2575  -508.5896
  [3,]   551.9081   551.8523   552.5088   552.4634   552.4778   552.1833
  [4,] -1018.7489 -1017.6725 -1018.0828 -1017.3203 -1017.7053 -1018.2489

Code
  result_log_dose@data
Output
  $betaW
              [,1]      [,2]        [,3]       [,4]      [,5]      [,6]      [,7]
  [1,]  -725.86444  506.6512  1738.36443  2969.9522  4201.766  5433.891  6664.853
  [2,]    59.50632 -449.3040  -957.78604 -1466.4300 -1975.198 -2483.035 -2990.646
  [3,] -1189.08277 -637.9527   -86.55955   465.4764  1017.797  1570.028  2122.283
  [4,]   733.96049 -283.6964 -1301.19537 -2319.2248 -3337.417 -4355.915 -5374.672
            [,8]      [,9]     [,10]     [,11]      [,12]      [,13]
  [1,]  7896.172  9126.836 10356.891 11587.280  12817.523  14047.956
  [2,] -3498.503 -4006.294 -4514.168 -5022.683  -5531.941  -6040.530
  [3,]  2674.191  3226.043  3778.552  4331.016   4883.494   5435.677
  [4,] -6393.421 -7411.094 -8429.177 -9446.497 -10464.202 -11482.451

  $betaZ
             [,1]       [,2]
  [1,]  0.2803695  0.8392922
  [2,] -0.1084857  5.6376489
  [3,] -1.9441049 11.6238390
  [4,]  0.7928335  2.5173100

  $delta
             [,1]       [,2]       [,3]       [,4]       [,5]       [,6]
  [1,]  1232.5156  1231.7133  1231.5878  1231.8141  1232.1250  1230.9616
  [2,]  -508.8103  -508.4821  -508.6440  -508.7677  -507.8371  -507.6109
  [3,]   551.1301   551.3931   552.0359   552.3207   552.2312   552.2549
  [4,] -1017.6569 -1017.4990 -1018.0295 -1018.1918 -1018.4981 -1018.7577
             [,7]       [,8]       [,9]      [,10]      [,11]      [,12]
  [1,]  1231.3187  1230.6641  1230.0549  1230.3898  1230.2425  1230.4327
  [2,]  -507.8571  -507.7917  -507.8736  -508.5151  -509.2575  -508.5896
  [3,]   551.9081   551.8523   552.5088   552.4634   552.4778   552.1833
  [4,] -1018.7489 -1017.6725 -1018.0828 -1017.3203 -1017.7053 -1018.2489

MCMC computes correct values for DualEndpointRW model (fixed params, empty data) with RW2

Code
  result@data
Output
  $betaW
             [,1]        [,2]       [,3]      [,4]        [,5]       [,6]
  [1,]   359.6244  -815.57203 -1981.3298 -3132.413 -4310.39339 -5513.4078
  [2,]   452.4520   -39.95466  -524.9930 -1006.614 -1485.60685 -1961.3560
  [3,] -1312.6394  -966.39567  -604.1289  -254.852    96.63905   447.3697
  [4,]  -938.8200 -1676.01724 -2409.0232 -3143.089 -3875.95656 -4609.2988
            [,7]      [,8]      [,9]      [,10]      [,11]      [,12]      [,13]
  [1,] -6702.663 -7897.970 -9085.459 -10271.217 -11431.918 -12599.639 -13754.332
  [2,] -2440.265 -2931.387 -3426.103  -3924.593  -4426.563  -4931.113  -5439.844
  [3,]   818.907  1199.474  1567.521   1945.322   2322.177   2702.918   3070.607
  [4,] -5344.668 -6081.550 -6828.081  -7579.302  -8336.050  -9099.004  -9858.306

  $betaZ
             [,1]      [,2]
  [1,]  0.8165898 16.132993
  [2,] -1.3257977  1.426216
  [3,]  0.7220600  7.700468
  [4,] -0.7592283  1.664707

  $delta
             [,1]       [,2]       [,3]       [,4]       [,5]       [,6]
  [1,] -1175.1964 -1165.7578 -1151.0836 -1177.9800 -1203.0144 -1189.2547
  [2,]  -492.4066  -485.0384  -481.6212  -478.9926  -475.7491  -478.9088
  [3,]   346.2437   362.2668   349.2769   351.4911   350.7306   371.5374
  [4,]  -737.1972  -733.0060  -734.0661  -732.8673  -733.3422  -735.3692
             [,7]       [,8]       [,9]      [,10]      [,11]      [,12]
  [1,] -1195.3073 -1187.4887 -1185.7582 -1160.7010 -1167.7211 -1154.6930
  [2,]  -491.1224  -494.7157  -498.4897  -501.9700  -504.5500  -508.7312
  [3,]   380.5665   368.0474   377.8007   376.8558   380.7409   367.6889
  [4,]  -736.8820  -746.5308  -751.2214  -756.7475  -762.9544  -759.3018

  $precBetaW
  [1] 0.05922806 1.34102710 0.16519402 0.88798420

  $precW
  [1] 0.1691313 0.1954275 0.6113126 0.7285836

  $rho
  [1] -0.8574380  0.1725056 -0.0535309 -0.7464971

Code
  result_log_dose@data
Output
  $betaW
             [,1]        [,2]       [,3]      [,4]        [,5]       [,6]
  [1,]   359.6244  -815.57203 -1981.3298 -3132.413 -4310.39339 -5513.4078
  [2,]   452.4520   -39.95466  -524.9930 -1006.614 -1485.60685 -1961.3560
  [3,] -1312.6394  -966.39567  -604.1289  -254.852    96.63905   447.3697
  [4,]  -938.8200 -1676.01724 -2409.0232 -3143.089 -3875.95656 -4609.2988
            [,7]      [,8]      [,9]      [,10]      [,11]      [,12]      [,13]
  [1,] -6702.663 -7897.970 -9085.459 -10271.217 -11431.918 -12599.639 -13754.332
  [2,] -2440.265 -2931.387 -3426.103  -3924.593  -4426.563  -4931.113  -5439.844
  [3,]   818.907  1199.474  1567.521   1945.322   2322.177   2702.918   3070.607
  [4,] -5344.668 -6081.550 -6828.081  -7579.302  -8336.050  -9099.004  -9858.306

  $betaZ
             [,1]      [,2]
  [1,]  0.8165898 16.132993
  [2,] -1.3257977  1.426216
  [3,]  0.7220600  7.700468
  [4,] -0.7592283  1.664707

  $delta
             [,1]       [,2]       [,3]       [,4]       [,5]       [,6]
  [1,] -1175.1964 -1165.7578 -1151.0836 -1177.9800 -1203.0144 -1189.2547
  [2,]  -492.4066  -485.0384  -481.6212  -478.9926  -475.7491  -478.9088
  [3,]   346.2437   362.2668   349.2769   351.4911   350.7306   371.5374
  [4,]  -737.1972  -733.0060  -734.0661  -732.8673  -733.3422  -735.3692
             [,7]       [,8]       [,9]      [,10]      [,11]      [,12]
  [1,] -1195.3073 -1187.4887 -1185.7582 -1160.7010 -1167.7211 -1154.6930
  [2,]  -491.1224  -494.7157  -498.4897  -501.9700  -504.5500  -508.7312
  [3,]   380.5665   368.0474   377.8007   376.8558   380.7409   367.6889
  [4,]  -736.8820  -746.5308  -751.2214  -756.7475  -762.9544  -759.3018

  $precBetaW
  [1] 0.05922806 1.34102710 0.16519402 0.88798420

  $precW
  [1] 0.1691313 0.1954275 0.6113126 0.7285836

  $rho
  [1] -0.8574380  0.1725056 -0.0535309 -0.7464971

MCMC computes correct values for DualEndpointBeta model with fixed parameters

Code
  result@data
Output
  $betaW
       [,1]     [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13]
  [1,]   10 10.02243   10   10   10   10   10   10   10    10    10    10    10
  [2,]   10 10.02243   10   10   10   10   10   10   10    10    10    10    10
  [3,]   10 10.02243   10   10   10   10   10   10   10    10    10    10    10
  [4,]   10 10.02243   10   10   10   10   10   10   10    10    10    10    10

  $betaZ
            [,1]       [,2]
  [1,] -2.498969 0.03700095
  [2,] -2.498969 0.03700095
  [3,] -2.498969 0.03700095
  [4,] -2.498969 0.03700095

Code
  result_log_dose@data
Output
  $betaW
       [,1]     [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13]
  [1,]   10 10.02243   10   10   10   10   10   10   10    10    10    10    10
  [2,]   10 10.02243   10   10   10   10   10   10   10    10    10    10    10
  [3,]   10 10.02243   10   10   10   10   10   10   10    10    10    10    10
  [4,]   10 10.02243   10   10   10   10   10   10   10    10    10    10    10

  $betaZ
            [,1]      [,2]
  [1,] -1.791151 0.1417847
  [2,] -1.791151 0.1417847
  [3,] -1.791151 0.1417847
  [4,] -1.791151 0.1417847

MCMC computes correct values for DualEndpointBeta model

Code
  result@data
Output
  $E0
  [1] 2.959874 3.855516 2.877285 3.139652

  $Emax
  [1] 5.849261 4.195788 6.041242 6.927711

  $betaW
           [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
  [1,] 2.959874 2.961494 2.959874 2.959874 2.959874 2.959874 2.959874 2.959874
  [2,] 3.855516 3.855707 3.855516 3.855516 3.855516 3.855516 3.855516 3.855516
  [3,] 2.877285 2.879059 2.877285 2.877285 2.877285 2.877285 2.877285 2.877285
  [4,] 3.139652 3.141776 3.139652 3.139652 3.139652 3.139652 3.139652 3.139652
           [,9]    [,10]    [,11]    [,12]    [,13]
  [1,] 2.959874 2.959874 2.959874 2.959874 2.959874
  [2,] 3.855516 3.855516 3.855516 3.855516 3.855516
  [3,] 2.877285 2.877285 2.877285 2.877285 2.877285
  [4,] 3.139652 3.139652 3.139652 3.139652 3.139652

  $betaZ
            [,1]       [,2]
  [1,] -3.620806 0.07403348
  [2,] -3.931225 0.06756602
  [3,] -3.931225 0.06756602
  [4,] -3.931225 0.06756602

  $precW
  [1] 0.0004239878 0.0004673854 0.0005892920 0.0004534128

  $rho
  [1] -0.63946386 -0.42112478 -0.12308691  0.05374152

Code
  result_log_dose@data
Output
  $E0
  [1] 4.748688 5.187207 2.726835 3.550073

  $Emax
  [1] 8.254186 6.941639 7.529266 7.980235

  $betaW
           [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
  [1,] 4.748688 4.750654 4.748688 4.748688 4.748688 4.748688 4.748688 4.748688
  [2,] 5.187207 5.188191 5.187207 5.187207 5.187207 5.187207 5.187207 5.187207
  [3,] 2.726835 2.729528 2.726835 2.726835 2.726835 2.726835 2.726835 2.726835
  [4,] 3.550073 3.552557 3.550073 3.550073 3.550073 3.550073 3.550073 3.550073
           [,9]    [,10]    [,11]    [,12]    [,13]
  [1,] 4.748688 4.748688 4.748688 4.748688 4.748688
  [2,] 5.187207 5.187207 5.187207 5.187207 5.187207
  [3,] 2.726835 2.726835 2.726835 2.726835 2.726835
  [4,] 3.550073 3.550073 3.550073 3.550073 3.550073

  $betaZ
            [,1]      [,2]
  [1,] -2.117768 0.6925986
  [2,] -2.117768 0.6925986
  [3,] -2.117768 0.6925986
  [4,] -2.117768 0.6925986

  $precW
  [1] 0.0005884974 0.0002756555 0.0003533447 0.0004380398

  $rho
  [1] -0.5037564 -0.1568004 -0.7883783 -0.8992664

MCMC computes correct values for DualEndpointBeta model with fixed parameters (empty data)

Code
  result@data
Output
  $betaW
       [,1]     [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13]
  [1,]   10 10.02243   10   10   10   10   10   10   10    10    10    10    10
  [2,]   10 10.02243   10   10   10   10   10   10   10    10    10    10    10
  [3,]   10 10.02243   10   10   10   10   10   10   10    10    10    10    10
  [4,]   10 10.02243   10   10   10   10   10   10   10    10    10    10    10

  $betaZ
             [,1]     [,2]
  [1,] -0.7258644 9.323290
  [2,] -1.6047136 2.114870
  [3,] -2.3268647 5.552484
  [4,]  0.6697347 2.024655

Code
  result_log_dose@data
Output
  $betaW
       [,1]     [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13]
  [1,]   10 10.02243   10   10   10   10   10   10   10    10    10    10    10
  [2,]   10 10.02243   10   10   10   10   10   10   10    10    10    10    10
  [3,]   10 10.02243   10   10   10   10   10   10   10    10    10    10    10
  [4,]   10 10.02243   10   10   10   10   10   10   10    10    10    10    10

  $betaZ
             [,1]     [,2]
  [1,] -0.7258644 9.323290
  [2,] -1.6047136 2.114870
  [3,] -2.3268647 5.552484
  [4,]  0.6697347 2.024655

MCMC computes correct values for DualEndpointBeta model (empty data)

Code
  result@data
Output
  $E0
  [1] 2.860619 2.336103 4.946781 3.392726

  $Emax
  [1] 3.858561 2.878886 5.708037 3.806117

  $betaW
           [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
  [1,] 2.860619 2.861179 2.860619 2.860619 2.860619 2.860619 2.860619 2.860619
  [2,] 2.336103 2.336408 2.336103 2.336103 2.336103 2.336103 2.336103 2.336103
  [3,] 4.946781 4.947208 4.946781 4.946781 4.946781 4.946781 4.946781 4.946781
  [4,] 3.392726 3.392958 3.392726 3.392726 3.392726 3.392726 3.392726 3.392726
           [,9]    [,10]    [,11]    [,12]    [,13]
  [1,] 2.860619 2.860619 2.860619 2.860619 2.860619
  [2,] 2.336103 2.336103 2.336103 2.336103 2.336103
  [3,] 4.946781 4.946781 4.946781 4.946781 4.946781
  [4,] 3.392726 3.392726 3.392726 3.392726 3.392726

  $betaZ
             [,1]      [,2]
  [1,] -1.1751964 0.5462308
  [2,] -0.2946006 3.9771203
  [3,]  0.8896501 1.5289638
  [4,]  0.6563921 1.9664409

  $precW
  [1] 0.3589515 0.5428550 0.8249277 0.5744209

  $rho
  [1] -0.1434824 -0.4058777 -0.9311362  0.1910782

Code
  result_log_dose@data
Output
  $E0
  [1] 2.860619 2.336103 4.946781 3.392726

  $Emax
  [1] 3.858561 2.878886 5.708037 3.806117

  $betaW
           [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
  [1,] 2.860619 2.861179 2.860619 2.860619 2.860619 2.860619 2.860619 2.860619
  [2,] 2.336103 2.336408 2.336103 2.336103 2.336103 2.336103 2.336103 2.336103
  [3,] 4.946781 4.947208 4.946781 4.946781 4.946781 4.946781 4.946781 4.946781
  [4,] 3.392726 3.392958 3.392726 3.392726 3.392726 3.392726 3.392726 3.392726
           [,9]    [,10]    [,11]    [,12]    [,13]
  [1,] 2.860619 2.860619 2.860619 2.860619 2.860619
  [2,] 2.336103 2.336103 2.336103 2.336103 2.336103
  [3,] 4.946781 4.946781 4.946781 4.946781 4.946781
  [4,] 3.392726 3.392726 3.392726 3.392726 3.392726

  $betaZ
             [,1]      [,2]
  [1,] -1.1751964 0.5462308
  [2,] -0.2946006 3.9771203
  [3,]  0.8896501 1.5289638
  [4,]  0.6563921 1.9664409

  $precW
  [1] 0.3589515 0.5428550 0.8249277 0.5744209

  $rho
  [1] -0.1434824 -0.4058777 -0.9311362  0.1910782

MCMC computes correct values for DualEndpointEmax model with fixed parameters

Code
  result@data
Output
  $betaW
          [,1]     [,2] [,3]     [,4]     [,5]     [,6]     [,7]     [,8] [,9]
  [1,] 10.0002 14.44444   18 20.90909 23.33333 25.38462 27.14286 28.66667   30
  [2,] 10.0002 14.44444   18 20.90909 23.33333 25.38462 27.14286 28.66667   30
  [3,] 10.0002 14.44444   18 20.90909 23.33333 25.38462 27.14286 28.66667   30
  [4,] 10.0002 14.44444   18 20.90909 23.33333 25.38462 27.14286 28.66667   30
          [,10]    [,11]    [,12] [,13]
  [1,] 31.17647 32.22222 33.15789    34
  [2,] 31.17647 32.22222 33.15789    34
  [3,] 31.17647 32.22222 33.15789    34
  [4,] 31.17647 32.22222 33.15789    34

  $betaZ
            [,1]       [,2]
  [1,] -2.498969 0.03700095
  [2,] -2.498969 0.03700095
  [3,] -2.498969 0.03700095
  [4,] -2.498969 0.03700095

Code
  result_log_dose@data
Output
  $betaW
          [,1]     [,2] [,3]     [,4]     [,5]     [,6]     [,7]     [,8] [,9]
  [1,] 10.0002 14.44444   18 20.90909 23.33333 25.38462 27.14286 28.66667   30
  [2,] 10.0002 14.44444   18 20.90909 23.33333 25.38462 27.14286 28.66667   30
  [3,] 10.0002 14.44444   18 20.90909 23.33333 25.38462 27.14286 28.66667   30
  [4,] 10.0002 14.44444   18 20.90909 23.33333 25.38462 27.14286 28.66667   30
          [,10]    [,11]    [,12] [,13]
  [1,] 31.17647 32.22222 33.15789    34
  [2,] 31.17647 32.22222 33.15789    34
  [3,] 31.17647 32.22222 33.15789    34
  [4,] 31.17647 32.22222 33.15789    34

  $betaZ
            [,1]      [,2]
  [1,] -1.791151 0.1417847
  [2,] -1.791151 0.1417847
  [3,] -1.791151 0.1417847
  [4,] -1.791151 0.1417847

MCMC computes correct values for DualEndpointEmax model

Code
  result@data
Output
  $E0
  [1] 67.11412 63.00022 54.58665 57.35462

  $ED50
  [1] 260.2212 264.4737 258.8991 260.4557

  $Emax
  [1] 212.2862 208.7443 214.9901 225.6857

  $betaW
           [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
  [1,] 67.11417 68.49554 69.85093 71.18100 72.48646 73.76798 75.02623 76.26182
  [2,] 63.00027 64.36500 65.70446 67.01929 68.31017 69.57774 70.82264 72.04547
  [3,] 54.58671 56.12073 57.62575 59.10253 60.55184 61.97445 63.37109 64.74247
  [4,] 57.35468 58.95499 60.52522 62.06615 63.57860 65.06334 66.52112 67.95269
           [,9]    [,10]    [,11]    [,12]    [,13]
  [1,] 77.47536 78.66745 79.83863 80.98946 82.12047
  [2,] 73.24680 74.42720 75.58721 76.72736 77.84814
  [3,] 66.08926 67.41213 68.71169 69.98857 71.24334
  [4,] 69.35873 70.73993 72.09693 73.43037 74.74086

  $betaZ
            [,1]       [,2]
  [1,] -2.430269 0.02501872
  [2,] -2.430269 0.02501872
  [3,] -2.670751 0.04923801
  [4,] -2.670751 0.04923801

  $precW
  [1] 0.0007250610 0.0008827606 0.0009257190 0.0013743503

  $rho
  [1] -0.12229173  0.02420671  0.48289450  0.10104428

Code
  result_log_dose@data
Output
  $E0
  [1] 47.12090 39.47157 40.65994 42.60600

  $ED50
  [1] 202.9029 206.6280 200.5536 186.5548

  $Emax
  [1] 255.4720 256.1366 249.4947 248.5258

  $betaW
           [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
  [1,] 47.12100 49.65678 52.13168 54.54776 56.90710 59.21167 61.46337 63.66398
  [2,] 39.47167 42.06167 44.59057 47.06043 49.47328 51.83107 54.13567 56.38885
  [3,] 40.66005 43.23112 45.73976 48.18810 50.57831 52.91242 55.19239 57.42008
  [4,] 42.60611 45.32902 47.98096 50.56457 53.08246 55.53710 57.93085 60.26594
           [,9]    [,10]    [,11]    [,12]    [,13]
  [1,] 65.81523 67.91876 69.97614 71.98887 73.95840
  [2,] 58.59233 60.74772 62.85659 64.92041 66.94063
  [3,] 59.59727 61.72565 63.80686 65.84243 67.83386
  [4,] 62.54452 64.76859 66.94010 69.06088 71.13270

  $betaZ
            [,1]      [,2]
  [1,] -2.620468 0.7207425
  [2,] -2.620468 0.7207425
  [3,] -2.620468 0.7207425
  [4,] -2.620468 0.7207425

  $precW
  [1] 0.0009944088 0.0001887385 0.0006597842 0.0009435991

  $rho
  [1] -0.8478294 -0.9008383 -0.7634940 -0.7231665

MCMC computes correct values for DualEndpointEmax model with fixed parameters (empty data)

Code
  result@data
Output
  $betaW
          [,1]     [,2] [,3]     [,4]     [,5]     [,6]     [,7]     [,8] [,9]
  [1,] 10.0002 14.44444   18 20.90909 23.33333 25.38462 27.14286 28.66667   30
  [2,] 10.0002 14.44444   18 20.90909 23.33333 25.38462 27.14286 28.66667   30
  [3,] 10.0002 14.44444   18 20.90909 23.33333 25.38462 27.14286 28.66667   30
  [4,] 10.0002 14.44444   18 20.90909 23.33333 25.38462 27.14286 28.66667   30
          [,10]    [,11]    [,12] [,13]
  [1,] 31.17647 32.22222 33.15789    34
  [2,] 31.17647 32.22222 33.15789    34
  [3,] 31.17647 32.22222 33.15789    34
  [4,] 31.17647 32.22222 33.15789    34

  $betaZ
             [,1]     [,2]
  [1,] -0.7258644 9.323290
  [2,] -1.6047136 2.114870
  [3,] -2.3268647 5.552484
  [4,]  0.6697347 2.024655

Code
  result_log_dose@data
Output
  $betaW
          [,1]     [,2] [,3]     [,4]     [,5]     [,6]     [,7]     [,8] [,9]
  [1,] 10.0002 14.44444   18 20.90909 23.33333 25.38462 27.14286 28.66667   30
  [2,] 10.0002 14.44444   18 20.90909 23.33333 25.38462 27.14286 28.66667   30
  [3,] 10.0002 14.44444   18 20.90909 23.33333 25.38462 27.14286 28.66667   30
  [4,] 10.0002 14.44444   18 20.90909 23.33333 25.38462 27.14286 28.66667   30
          [,10]    [,11]    [,12] [,13]
  [1,] 31.17647 32.22222 33.15789    34
  [2,] 31.17647 32.22222 33.15789    34
  [3,] 31.17647 32.22222 33.15789    34
  [4,] 31.17647 32.22222 33.15789    34

  $betaZ
             [,1]     [,2]
  [1,] -0.7258644 9.323290
  [2,] -1.6047136 2.114870
  [3,] -2.3268647 5.552484
  [4,]  0.6697347 2.024655

MCMC computes correct values for DualEndpointEmax model (empty data)

Code
  result@data
Output
  $E0
  [1] 57.28534 34.03490 43.80971 77.73207

  $ED50
  [1] 132.7543 191.1940 238.8098 432.1697

  $Emax
  [1] 186.0619 434.8454 430.6047 194.9948

  $betaW
           [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
  [1,] 57.28543 59.66560 61.95948 64.17157 66.30619 68.36733 70.35872 72.28384
  [2,] 34.03511 39.20814 44.24955 49.16409 53.95650 58.63127 63.19268 67.64480
  [3,] 43.80987 47.81696 51.74202 55.58741 59.35552 63.04866 66.66905 70.21882
  [4,] 77.73210 78.40651 79.07323 79.73237 80.38405 81.02841 81.66557 82.29564
           [,9]    [,10]    [,11]    [,12]    [,13]
  [1,] 74.14595 75.94809 77.69311 79.38369 81.02233
  [2,] 71.99152 76.23654 80.38338 84.43542 88.39585
  [3,] 73.70000 77.11458 80.46443 83.75140 86.97722
  [4,] 82.91874 83.53499 84.14450 84.74738 85.34373

  $betaZ
             [,1]     [,2]
  [1,] -0.2510064 4.274921
  [2,]  1.2196258 1.931522
  [3,]  0.8042982 2.884946
  [4,]  0.2772333 2.381924

  $precW
  [1] 0.11792318 0.08988104 0.64591148 0.43738790

  $rho
  [1]  0.5447331 -0.2752328 -0.5923060  0.7406452

Code
  result_log_dose@data
Output
  $E0
  [1] 57.28534 34.03490 43.80971 77.73207

  $ED50
  [1] 132.7543 191.1940 238.8098 432.1697

  $Emax
  [1] 186.0619 434.8454 430.6047 194.9948

  $betaW
           [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
  [1,] 57.28543 59.66560 61.95948 64.17157 66.30619 68.36733 70.35872 72.28384
  [2,] 34.03511 39.20814 44.24955 49.16409 53.95650 58.63127 63.19268 67.64480
  [3,] 43.80987 47.81696 51.74202 55.58741 59.35552 63.04866 66.66905 70.21882
  [4,] 77.73210 78.40651 79.07323 79.73237 80.38405 81.02841 81.66557 82.29564
           [,9]    [,10]    [,11]    [,12]    [,13]
  [1,] 74.14595 75.94809 77.69311 79.38369 81.02233
  [2,] 71.99152 76.23654 80.38338 84.43542 88.39585
  [3,] 73.70000 77.11458 80.46443 83.75140 86.97722
  [4,] 82.91874 83.53499 84.14450 84.74738 85.34373

  $betaZ
             [,1]     [,2]
  [1,] -0.2510064 4.274921
  [2,]  1.2196258 1.931522
  [3,]  0.8042982 2.884946
  [4,]  0.2772333 2.381924

  $precW
  [1] 0.11792318 0.08988104 0.64591148 0.43738790

  $rho
  [1]  0.5447331 -0.2752328 -0.5923060  0.7406452

MCMC computes correct values for LogisticIndepBeta model

Code
  result@data
Output
  $phi1
  [1] -0.2140691 -0.9963186 -0.9963186 -1.7178745

  $phi2
  [1] 0.27003542 0.01697991 0.01697991 0.37729886

MCMC computes correct values for LogisticIndepBeta model (empty data)

Code
  result@data
Output
  $phi1
  [1] -1.775733 -5.848248 -1.266877 -5.483583

  $phi2
  [1] 0.3449701 1.0647430 0.3513896 1.0368889

MCMC computes correct values for Effloglog model

Code
  result@data
Output
  $theta1
  [1] -1.248083 -2.088724 -4.963567 -1.699768

  $theta2
  [1] 1.423735 1.870617 4.347765 1.936372

  $nu
  [1] 2.472057 2.242145 1.656834 2.802397

MCMC computes correct values for Effloglog model (empty data)

Code
  result@data
Output
  $theta1
  [1] -0.8217308 -1.2444975 -2.4870098 -0.9905060

  $theta2
  [1] 1.742472 1.935503 3.022090 1.968669

  $nu
  [1] 20.533736 15.128041  4.359347 29.198563

MCMC computes correct values for EffFlexi model

Code
  result@data
Output
  $ExpEff
            [,1]      [,2]      [,3]        [,4]       [,5]     [,6]       [,7]
  [1,] 0.7672193 0.4879123 0.4928340  2.06456500  1.7002369 2.016392  1.9229294
  [2,] 0.7562453 0.5330042 0.4214018 -0.15655774  2.3823607 3.507772  4.6037886
  [3,] 0.7576968 0.5111914 0.4668699 -0.04876909  0.6188389 2.042441  2.6289971
  [4,] 0.7568636 0.5091830 0.4689064  0.06682629 -0.3135319 1.489857 -0.3968666
            [,8]       [,9]     [,10]    [,11]    [,12]
  [1,]  3.058090  1.1998861 2.1821223 2.693873 2.502903
  [2,]  5.363741  6.3654165 4.2922960 2.907385 2.506099
  [3,]  1.958357  0.6000557 1.1226388 2.599905 2.501354
  [4,] -2.172702 -0.3058473 0.5452364 2.069557 2.505000

  $sigma2W
             [,1]
  [1,] 0.10913016
  [2,] 0.07964666
  [3,] 0.11801705
  [4,] 0.13254077

  $sigma2betaW
           [,1]
  [1,] 2.291492
  [2,] 2.608590
  [3,] 2.252753
  [4,] 1.962545

MCMC computes correct values for EffFlexi model (RW2)

Code
  result@data
Output
  $ExpEff
            [,1]      [,2]      [,3]       [,4]       [,5]      [,6]      [,7]
  [1,] 0.7630782 0.5095909 0.4481623 -0.5470574 -0.6979189 -2.864326 -3.946862
  [2,] 0.7607543 0.5166401 0.4381781  0.2468354 -2.7113638 -3.252579 -3.619454
  [3,] 0.7664345 0.5001368 0.4609672  2.5330036  5.1413807  7.359390  8.452245
  [4,] 0.7697145 0.4941193 0.4663217  0.1844835  3.0573439  5.446762  6.386212
            [,8]       [,9]      [,10]     [,11]    [,12]
  [1,] -3.652219 -2.1922017 -0.4877764 0.8099165 2.512499
  [2,] -2.773986 -0.3137721  0.1227855 0.6939263 2.513033
  [3,] 10.099099  6.5574661  3.1359366 1.8899065 2.511890
  [4,]  6.906622  6.3523570  7.1189684 5.5187513 2.512011

  $sigma2W
             [,1]
  [1,] 0.21716195
  [2,] 0.06111584
  [3,] 0.10309552
  [4,] 0.16816591

  $sigma2betaW
           [,1]
  [1,] 2.559479
  [2,] 3.328569
  [3,] 2.210422
  [4,] 2.056716

MCMC computes correct values for EffFlexi model (empty data)

Code
  result@data
Output
  $ExpEff
           [,1]       [,2]       [,3]       [,4]        [,5]        [,6]
  [1,] 1.225490  3.0839844 2.12928349  0.2700333  0.41273635 -0.34108462
  [2,] 1.227582  2.6616783 3.62952402  4.1859148  1.88733382 -1.00489873
  [3,] 1.232744  2.0292269 1.08522425  1.0795338 -0.09933442 -0.07668353
  [4,] 1.228831 -0.9402123 0.05784494 -0.8959898 -0.57887080 -0.83072784
             [,7]      [,8]        [,9]     [,10]    [,11]    [,12]
  [1,]  1.2629508 0.7880890  0.73875705 0.8797659 1.855001 2.510510
  [2,]  1.7970138 2.6919593  1.21092403 0.5493481 1.882342 2.508418
  [3,] -0.7743345 0.4562472 -0.22742142 0.3832043 2.386059 2.503256
  [4,]  1.5963565 1.8285942  0.03311922 2.3840162 2.930087 2.507169

  $sigma2W
             [,1]
  [1,] 0.69508707
  [2,] 0.09938619
  [3,] 1.02791665
  [4,] 0.14289676

  $sigma2betaW
           [,1]
  [1,] 2.398108
  [2,] 2.139745
  [3,] 2.465330
  [4,] 2.049230

MCMC computes correct values for EffFlexi model (empty data, RW2)

Code
  result@data
Output
  $ExpEff
        [,1]      [,2]       [,3]       [,4]      [,5]      [,6]      [,7]
  [1,] 1.223  6.803911  11.788163  13.949884  14.21939  12.00250  10.10042
  [2,] 1.223 -2.267596  -5.612351  -8.690355 -10.12993 -11.24248 -10.74001
  [3,] 1.223  7.085129  10.158879  11.207407  11.93213  12.92086  12.90089
  [4,] 1.223 -8.357629 -15.530921 -18.221162 -20.41940 -21.15290 -20.27631
            [,8]       [,9]      [,10]     [,11] [,12]
  [1,]   9.04445   7.301129   5.409379  4.147922 2.513
  [2,] -11.41255 -10.705392  -8.707142 -3.117653 2.513
  [3,]  10.58026   8.926475   5.930708  2.546946 2.513
  [4,] -21.44332 -19.769558 -14.047895 -7.380435 2.513

  $sigma2W
             [,1]
  [1,] 0.95061320
  [2,] 0.12591440
  [3,] 0.06719078
  [4,] 0.02936185

  $sigma2betaW
           [,1]
  [1,] 2.051256
  [2,] 2.233393
  [3,] 3.469060
  [4,] 2.977588

MCMC computes correct values for DALogisticLogNormal model

Code
  result@data
Output
  $alpha0
  [1] -0.8866238  0.2638399 -1.5478437 -0.5337333

  $alpha1
  [1] 0.4124879 0.4451248 0.1405851 1.1509924

  $lambda
               [,1]         [,2]        [,3]
  [1,] 1.865608e-11 0.0003936759 0.007446998
  [2,] 5.690263e-10 0.0554967361 0.015927912
  [3,] 5.510949e-10 0.1019245232 0.287915936
  [4,] 5.504115e-10 0.0007306373 0.025253135

MCMC computes correct values for DALogisticLogNormal model and empty data

Code
  result@data
Output
  $alpha0
  [1] -0.7258644  0.7181814  0.3449151 -1.3669881

  $alpha1
  [1] 9.3232903 3.7472384 6.0757060 0.7219512

  $lambda
               [,1]       [,2]       [,3]
  [1,] 4.704793e-03 0.44544093 0.06189658
  [2,] 1.687200e-01 1.13156451 0.22390135
  [3,] 3.471273e-02 0.01260272 1.66933173
  [4,] 3.430421e-05 1.14915573 0.38080491

MCMC computes correct values for TITELogisticLogNormal model (linear)

Code
  result@data
Output
  $alpha0
  [1] -1.379114 -1.379114 -1.379114 -1.780157

  $alpha1
  [1] 0.5018051 0.5018051 0.5018051 0.3285917

MCMC computes correct values for TITELogisticLogNormal model (linear, data 2)

Code
  result@data
Output
  $alpha0
  [1] -0.1830693 -0.1830693 -0.1830693 -0.5301295

  $alpha1
  [1] 0.5047497 0.5047497 0.5047497 0.3499039

MCMC computes correct values for TITELogisticLogNormal model (linear, empty data)

Code
  result@data
Output
  $alpha0
  [1] -0.7258644 -1.6047136 -2.3268647  0.6697347

  $alpha1
  [1] 9.323290 2.114870 5.552484 2.024655

MCMC computes correct values for TITELogisticLogNormal model (adaptive)

Code
  result@data
Output
  $alpha0
  [1] -1.389973 -1.389973 -1.389973 -1.797080

  $alpha1
  [1] 0.5051006 0.5051006 0.5051006 0.3286390

MCMC computes correct values for TITELogisticLogNormal model (adaptive, data 2)

Code
  result@data
Output
  $alpha0
  [1] -0.1830693 -0.1830693 -0.1830693 -0.5301295

  $alpha1
  [1] 0.5047497 0.5047497 0.5047497 0.3499039

MCMC computes correct values for TITELogisticLogNormal model (adaptive, empty data)

Code
  result@data
Output
  $alpha0
  [1] -0.7258644 -1.6047136 -2.3268647  0.6697347

  $alpha1
  [1] 9.323290 2.114870 5.552484 2.024655

MCMC computes correct values for OneParLogNormalPrior model

Code
  result@data
Output
  $alpha
  [1] 0.2149170 0.3987768 1.0718321 0.8496711

MCMC computes correct values for OneParLogNormalPrior model and empty data

Code
  result@data
Output
  $alpha
  [1] -1.0265273  0.3965024 -2.2694077  0.6403072

MCMC computes correct values for OneParExpPrior model

Code
  result@data
Output
  $theta
  [1] 0.7980177 1.7746948 0.9487709 0.9468932

MCMC computes correct values for OneParExpPrior model and empty data

Code
  result@data
Output
  $theta
  [1] 0.37759092 0.07285336 0.21803431 0.61478103

MCMC computes correct values for FractionalCRM model

Code
  result@data
Output
  $alpha
  [1] -0.8525427 -0.5374850  0.1646873 -0.5913058

MCMC computes correct values for FractionalCRM model and empty data

Code
  result@data
Output
  $alpha
  [1] -1.0265273  0.3965024 -2.2694077  0.6403072

Roche/crmPack documentation built on June 30, 2024, 1:31 a.m.