Skip to contents

Extract posterior draws for specified parameters and gather them into a long format data frame with one row per draw per parameter.

Usage

gather_draws(object, ..., regex = FALSE, ndraws = NULL)

# S3 method for class 'ratiod_fit'
gather_draws(object, ..., regex = FALSE, ndraws = NULL)

Arguments

object

A ratiod_fit object

...

Parameter names to extract (unquoted).

regex

Logical; if TRUE, treat parameter names as regex patterns.

ndraws

Number of draws to return. If NULL, returns all.

Value

A data frame with columns .chain, .iteration, .draw, .variable, and .value.

Examples

# \donttest{
# Generate synthetic data
set.seed(333)
n <- 55
df <- data.frame(
  count = rpois(n, lambda = 24),
  effort = rgamma(n, shape = 12, rate = 1),
  depth = rnorm(n),
  site = sample(letters[1:4], n, replace = TRUE)
)

# Fit model
fit <- tratio(
  count | effort ~ depth + (1 | site),
  data = df,
  family = ratiod_poisson_gamma(),
  mode = "hmc",
  control = list(iter = 200, warmup = 100, chains = 1)
)
#> Inference: Exact (Tier 1)
#>   Backend: hmc
#> Fitting ratio model...
#>   Family: poisson_gamma
#>   Observations: 55
#>   Iterations: 200 (warmup: 100)
#> Running NUTS sampler...
#>   Parameters: 10
#>   Iterations: 200 (warmup: 100)
#>   Chains: 1 (cores: 1)

# Gather all fixed effects
gather_draws(fit, beta_num, beta_denom)
#>     .chain .iteration .draw     .variable        .value
#> 1        1          1     1   beta_num.1.  3.1329319284
#> 2        1          2     2   beta_num.1.  3.1486022812
#> 3        1          3     3   beta_num.1.  3.1432601972
#> 4        1          4     4   beta_num.1.  3.1263841539
#> 5        1          5     5   beta_num.1.  3.1682768215
#> 6        1          6     6   beta_num.1.  3.1663165196
#> 7        1          7     7   beta_num.1.  3.1336968503
#> 8        1          8     8   beta_num.1.  3.2129158470
#> 9        1          9     9   beta_num.1.  3.1120753300
#> 10       1         10    10   beta_num.1.  3.1787526672
#> 11       1         11    11   beta_num.1.  3.3003977714
#> 12       1         12    12   beta_num.1.  3.3201137589
#> 13       1         13    13   beta_num.1.  3.2777393719
#> 14       1         14    14   beta_num.1.  3.2975941346
#> 15       1         15    15   beta_num.1.  3.2230293594
#> 16       1         16    16   beta_num.1.  3.2319967480
#> 17       1         17    17   beta_num.1.  3.1702185950
#> 18       1         18    18   beta_num.1.  3.1696019654
#> 19       1         19    19   beta_num.1.  3.1607235562
#> 20       1         20    20   beta_num.1.  3.1850110346
#> 21       1         21    21   beta_num.1.  3.1559430042
#> 22       1         22    22   beta_num.1.  3.1827484818
#> 23       1         23    23   beta_num.1.  3.1781010591
#> 24       1         24    24   beta_num.1.  3.1725244552
#> 25       1         25    25   beta_num.1.  3.1529799397
#> 26       1         26    26   beta_num.1.  3.2205042983
#> 27       1         27    27   beta_num.1.  3.1987548628
#> 28       1         28    28   beta_num.1.  3.1761891347
#> 29       1         29    29   beta_num.1.  3.1870166356
#> 30       1         30    30   beta_num.1.  3.1859588638
#> 31       1         31    31   beta_num.1.  3.1802715169
#> 32       1         32    32   beta_num.1.  3.1414064820
#> 33       1         33    33   beta_num.1.  3.1414064820
#> 34       1         34    34   beta_num.1.  3.2148169334
#> 35       1         35    35   beta_num.1.  3.1885033431
#> 36       1         36    36   beta_num.1.  3.1306451269
#> 37       1         37    37   beta_num.1.  3.1074462798
#> 38       1         38    38   beta_num.1.  3.1435647790
#> 39       1         39    39   beta_num.1.  3.1071444969
#> 40       1         40    40   beta_num.1.  3.1070384333
#> 41       1         41    41   beta_num.1.  3.1330372260
#> 42       1         42    42   beta_num.1.  3.2087636943
#> 43       1         43    43   beta_num.1.  3.2045126075
#> 44       1         44    44   beta_num.1.  3.1594953017
#> 45       1         45    45   beta_num.1.  3.1296419829
#> 46       1         46    46   beta_num.1.  3.1681939747
#> 47       1         47    47   beta_num.1.  3.1995224548
#> 48       1         48    48   beta_num.1.  3.2167888005
#> 49       1         49    49   beta_num.1.  3.1598885136
#> 50       1         50    50   beta_num.1.  3.1728546933
#> 51       1         51    51   beta_num.1.  3.1626301884
#> 52       1         52    52   beta_num.1.  3.1829000407
#> 53       1         53    53   beta_num.1.  3.1955486823
#> 54       1         54    54   beta_num.1.  3.2138276585
#> 55       1         55    55   beta_num.1.  3.2802181443
#> 56       1         56    56   beta_num.1.  3.2878678621
#> 57       1         57    57   beta_num.1.  3.3209957142
#> 58       1         58    58   beta_num.1.  3.2799537345
#> 59       1         59    59   beta_num.1.  3.3825116534
#> 60       1         60    60   beta_num.1.  3.4010331027
#> 61       1         61    61   beta_num.1.  3.2434963884
#> 62       1         62    62   beta_num.1.  3.3361308081
#> 63       1         63    63   beta_num.1.  3.2285337034
#> 64       1         64    64   beta_num.1.  3.2073180683
#> 65       1         65    65   beta_num.1.  3.2462124337
#> 66       1         66    66   beta_num.1.  3.2009534293
#> 67       1         67    67   beta_num.1.  3.1659457203
#> 68       1         68    68   beta_num.1.  3.1880253588
#> 69       1         69    69   beta_num.1.  3.2477955593
#> 70       1         70    70   beta_num.1.  3.2120426380
#> 71       1         71    71   beta_num.1.  3.2032986179
#> 72       1         72    72   beta_num.1.  3.1739842941
#> 73       1         73    73   beta_num.1.  3.1400553637
#> 74       1         74    74   beta_num.1.  3.2130915335
#> 75       1         75    75   beta_num.1.  3.1559416798
#> 76       1         76    76   beta_num.1.  3.1430040639
#> 77       1         77    77   beta_num.1.  3.1613578533
#> 78       1         78    78   beta_num.1.  3.1858799402
#> 79       1         79    79   beta_num.1.  3.1460915537
#> 80       1         80    80   beta_num.1.  3.2011863320
#> 81       1         81    81   beta_num.1.  3.2070306316
#> 82       1         82    82   beta_num.1.  3.1619971001
#> 83       1         83    83   beta_num.1.  3.2169392019
#> 84       1         84    84   beta_num.1.  3.2647441185
#> 85       1         85    85   beta_num.1.  3.1251802711
#> 86       1         86    86   beta_num.1.  3.2072167761
#> 87       1         87    87   beta_num.1.  3.1836717853
#> 88       1         88    88   beta_num.1.  3.1608026866
#> 89       1         89    89   beta_num.1.  3.0988291320
#> 90       1         90    90   beta_num.1.  3.1678125366
#> 91       1         91    91   beta_num.1.  3.1139521044
#> 92       1         92    92   beta_num.1.  3.1371039919
#> 93       1         93    93   beta_num.1.  3.1327253292
#> 94       1         94    94   beta_num.1.  3.1298187402
#> 95       1         95    95   beta_num.1.  3.1211346415
#> 96       1         96    96   beta_num.1.  3.1353905825
#> 97       1         97    97   beta_num.1.  3.1013882122
#> 98       1         98    98   beta_num.1.  3.1953241971
#> 99       1         99    99   beta_num.1.  3.2101621980
#> 100      1        100   100   beta_num.1.  3.1635337587
#> 101      1          1     1   beta_num.2. -0.0719734664
#> 102      1          2     2   beta_num.2. -0.0274504634
#> 103      1          3     3   beta_num.2. -0.0551527420
#> 104      1          4     4   beta_num.2. -0.0348767846
#> 105      1          5     5   beta_num.2. -0.0318397203
#> 106      1          6     6   beta_num.2. -0.0226692426
#> 107      1          7     7   beta_num.2. -0.0335885572
#> 108      1          8     8   beta_num.2. -0.0311749251
#> 109      1          9     9   beta_num.2. -0.0642089233
#> 110      1         10    10   beta_num.2. -0.0632724752
#> 111      1         11    11   beta_num.2. -0.0051942410
#> 112      1         12    12   beta_num.2. -0.0188643860
#> 113      1         13    13   beta_num.2. -0.0392167605
#> 114      1         14    14   beta_num.2. -0.0629560807
#> 115      1         15    15   beta_num.2. -0.0032710570
#> 116      1         16    16   beta_num.2. -0.0465701808
#> 117      1         17    17   beta_num.2. -0.0295395663
#> 118      1         18    18   beta_num.2. -0.0321007146
#> 119      1         19    19   beta_num.2. -0.0146370954
#> 120      1         20    20   beta_num.2. -0.0003511561
#> 121      1         21    21   beta_num.2. -0.0364385144
#> 122      1         22    22   beta_num.2. -0.0399106141
#> 123      1         23    23   beta_num.2. -0.0422339147
#> 124      1         24    24   beta_num.2. -0.0597879946
#> 125      1         25    25   beta_num.2. -0.0172114185
#> 126      1         26    26   beta_num.2. -0.0277618784
#> 127      1         27    27   beta_num.2. -0.0419166706
#> 128      1         28    28   beta_num.2. -0.0356458526
#> 129      1         29    29   beta_num.2. -0.0215422126
#> 130      1         30    30   beta_num.2. -0.0740030287
#> 131      1         31    31   beta_num.2. -0.0153935563
#> 132      1         32    32   beta_num.2. -0.0352860987
#> 133      1         33    33   beta_num.2. -0.0352860987
#> 134      1         34    34   beta_num.2. -0.0345262607
#> 135      1         35    35   beta_num.2. -0.0139040137
#> 136      1         36    36   beta_num.2. -0.0787028239
#> 137      1         37    37   beta_num.2. -0.0390023051
#> 138      1         38    38   beta_num.2. -0.0689853150
#> 139      1         39    39   beta_num.2. -0.0277386260
#> 140      1         40    40   beta_num.2. -0.0629440146
#> 141      1         41    41   beta_num.2. -0.0500089563
#> 142      1         42    42   beta_num.2. -0.0512104730
#> 143      1         43    43   beta_num.2. -0.0470934107
#> 144      1         44    44   beta_num.2.  0.0001734759
#> 145      1         45    45   beta_num.2. -0.0157422843
#> 146      1         46    46   beta_num.2.  0.0148165911
#> 147      1         47    47   beta_num.2.  0.0195736435
#> 148      1         48    48   beta_num.2. -0.0547735001
#> 149      1         49    49   beta_num.2. -0.0797705929
#> 150      1         50    50   beta_num.2. -0.0101676327
#> 151      1         51    51   beta_num.2. -0.0530723674
#> 152      1         52    52   beta_num.2. -0.0082046341
#> 153      1         53    53   beta_num.2. -0.0186806477
#> 154      1         54    54   beta_num.2. -0.0118660502
#> 155      1         55    55   beta_num.2. -0.0592467678
#> 156      1         56    56   beta_num.2. -0.0408653793
#> 157      1         57    57   beta_num.2. -0.0683121089
#> 158      1         58    58   beta_num.2. -0.0979553993
#> 159      1         59    59   beta_num.2. -0.0664845913
#> 160      1         60    60   beta_num.2. -0.0415678217
#> 161      1         61    61   beta_num.2. -0.0742948610
#> 162      1         62    62   beta_num.2. -0.0680035064
#> 163      1         63    63   beta_num.2. -0.0704096854
#> 164      1         64    64   beta_num.2.  0.0694829653
#> 165      1         65    65   beta_num.2. -0.0810106638
#> 166      1         66    66   beta_num.2. -0.0836937181
#> 167      1         67    67   beta_num.2. -0.0602832799
#> 168      1         68    68   beta_num.2. -0.0093674358
#> 169      1         69    69   beta_num.2. -0.0095856564
#> 170      1         70    70   beta_num.2.  0.0067348840
#> 171      1         71    71   beta_num.2. -0.0288903320
#> 172      1         72    72   beta_num.2. -0.0395126957
#> 173      1         73    73   beta_num.2. -0.0015830578
#> 174      1         74    74   beta_num.2. -0.0490907547
#> 175      1         75    75   beta_num.2. -0.0723835385
#> 176      1         76    76   beta_num.2. -0.0367440526
#> 177      1         77    77   beta_num.2. -0.0150309835
#> 178      1         78    78   beta_num.2. -0.0548615796
#> 179      1         79    79   beta_num.2. -0.0481993019
#> 180      1         80    80   beta_num.2. -0.0260312854
#> 181      1         81    81   beta_num.2. -0.0183444920
#> 182      1         82    82   beta_num.2.  0.0218034429
#> 183      1         83    83   beta_num.2. -0.0575033301
#> 184      1         84    84   beta_num.2. -0.0460246838
#> 185      1         85    85   beta_num.2. -0.0585800986
#> 186      1         86    86   beta_num.2.  0.0489736846
#> 187      1         87    87   beta_num.2.  0.0369806653
#> 188      1         88    88   beta_num.2.  0.0064916672
#> 189      1         89    89   beta_num.2. -0.0212301549
#> 190      1         90    90   beta_num.2. -0.0404134683
#> 191      1         91    91   beta_num.2. -0.0253568040
#> 192      1         92    92   beta_num.2. -0.0245437736
#> 193      1         93    93   beta_num.2. -0.0463734306
#> 194      1         94    94   beta_num.2. -0.0164130621
#> 195      1         95    95   beta_num.2. -0.0194300240
#> 196      1         96    96   beta_num.2. -0.0311281631
#> 197      1         97    97   beta_num.2. -0.0104185345
#> 198      1         98    98   beta_num.2. -0.0303695149
#> 199      1         99    99   beta_num.2. -0.0004444491
#> 200      1        100   100   beta_num.2. -0.0797194824
#> 201      1          1     1 beta_denom.1.  2.4423945257
#> 202      1          2     2 beta_denom.1.  2.4017887937
#> 203      1          3     3 beta_denom.1.  2.3845685900
#> 204      1          4     4 beta_denom.1.  2.4224352020
#> 205      1          5     5 beta_denom.1.  2.3620063267
#> 206      1          6     6 beta_denom.1.  2.4102247183
#> 207      1          7     7 beta_denom.1.  2.4116079082
#> 208      1          8     8 beta_denom.1.  2.5172420554
#> 209      1          9     9 beta_denom.1.  2.5514928034
#> 210      1         10    10 beta_denom.1.  2.5007751174
#> 211      1         11    11 beta_denom.1.  2.5616979418
#> 212      1         12    12 beta_denom.1.  2.5779980603
#> 213      1         13    13 beta_denom.1.  2.5451976656
#> 214      1         14    14 beta_denom.1.  2.6612478771
#> 215      1         15    15 beta_denom.1.  2.5594234412
#> 216      1         16    16 beta_denom.1.  2.4646914692
#> 217      1         17    17 beta_denom.1.  2.4132133721
#> 218      1         18    18 beta_denom.1.  2.4894037015
#> 219      1         19    19 beta_denom.1.  2.4780369088
#> 220      1         20    20 beta_denom.1.  2.4330357337
#> 221      1         21    21 beta_denom.1.  2.4474314456
#> 222      1         22    22 beta_denom.1.  2.4463225967
#> 223      1         23    23 beta_denom.1.  2.3903376439
#> 224      1         24    24 beta_denom.1.  2.3639679482
#> 225      1         25    25 beta_denom.1.  2.5506490202
#> 226      1         26    26 beta_denom.1.  2.4363578739
#> 227      1         27    27 beta_denom.1.  2.4363788064
#> 228      1         28    28 beta_denom.1.  2.5009457784
#> 229      1         29    29 beta_denom.1.  2.5144832239
#> 230      1         30    30 beta_denom.1.  2.4440206161
#> 231      1         31    31 beta_denom.1.  2.4372789836
#> 232      1         32    32 beta_denom.1.  2.4080874898
#> 233      1         33    33 beta_denom.1.  2.4080874898
#> 234      1         34    34 beta_denom.1.  2.5007585672
#> 235      1         35    35 beta_denom.1.  2.3753192593
#> 236      1         36    36 beta_denom.1.  2.4520936812
#> 237      1         37    37 beta_denom.1.  2.3673166629
#> 238      1         38    38 beta_denom.1.  2.4374144532
#> 239      1         39    39 beta_denom.1.  2.3459980211
#> 240      1         40    40 beta_denom.1.  2.4244366825
#> 241      1         41    41 beta_denom.1.  2.4038555333
#> 242      1         42    42 beta_denom.1.  2.4459879203
#> 243      1         43    43 beta_denom.1.  2.4606704227
#> 244      1         44    44 beta_denom.1.  2.4448945232
#> 245      1         45    45 beta_denom.1.  2.3933473138
#> 246      1         46    46 beta_denom.1.  2.3986469051
#> 247      1         47    47 beta_denom.1.  2.3883947551
#> 248      1         48    48 beta_denom.1.  2.4213842436
#> 249      1         49    49 beta_denom.1.  2.4124777881
#> 250      1         50    50 beta_denom.1.  2.4757883452
#> 251      1         51    51 beta_denom.1.  2.3860292681
#> 252      1         52    52 beta_denom.1.  2.4686815368
#> 253      1         53    53 beta_denom.1.  2.4372498859
#> 254      1         54    54 beta_denom.1.  2.4671709241
#> 255      1         55    55 beta_denom.1.  2.5619799389
#> 256      1         56    56 beta_denom.1.  2.5615449286
#> 257      1         57    57 beta_denom.1.  2.6453175308
#> 258      1         58    58 beta_denom.1.  2.6437364975
#> 259      1         59    59 beta_denom.1.  2.6236233529
#> 260      1         60    60 beta_denom.1.  2.5730754403
#> 261      1         61    61 beta_denom.1.  2.5933318929
#> 262      1         62    62 beta_denom.1.  2.5867597803
#> 263      1         63    63 beta_denom.1.  2.4846755378
#> 264      1         64    64 beta_denom.1.  2.5304373342
#> 265      1         65    65 beta_denom.1.  2.4901887837
#> 266      1         66    66 beta_denom.1.  2.4949068625
#> 267      1         67    67 beta_denom.1.  2.4522881025
#> 268      1         68    68 beta_denom.1.  2.4562751385
#> 269      1         69    69 beta_denom.1.  2.4816351727
#> 270      1         70    70 beta_denom.1.  2.4510962320
#> 271      1         71    71 beta_denom.1.  2.4482555775
#> 272      1         72    72 beta_denom.1.  2.4637681327
#> 273      1         73    73 beta_denom.1.  2.4169860742
#> 274      1         74    74 beta_denom.1.  2.4752698695
#> 275      1         75    75 beta_denom.1.  2.4816256741
#> 276      1         76    76 beta_denom.1.  2.5297839130
#> 277      1         77    77 beta_denom.1.  2.3991331420
#> 278      1         78    78 beta_denom.1.  2.4353384392
#> 279      1         79    79 beta_denom.1.  2.4054633261
#> 280      1         80    80 beta_denom.1.  2.4618265336
#> 281      1         81    81 beta_denom.1.  2.4488197111
#> 282      1         82    82 beta_denom.1.  2.4508712574
#> 283      1         83    83 beta_denom.1.  2.4465001559
#> 284      1         84    84 beta_denom.1.  2.4827260508
#> 285      1         85    85 beta_denom.1.  2.4833421363
#> 286      1         86    86 beta_denom.1.  2.3911593110
#> 287      1         87    87 beta_denom.1.  2.4710390608
#> 288      1         88    88 beta_denom.1.  2.3608701273
#> 289      1         89    89 beta_denom.1.  2.4048026127
#> 290      1         90    90 beta_denom.1.  2.4050086879
#> 291      1         91    91 beta_denom.1.  2.5198570382
#> 292      1         92    92 beta_denom.1.  2.4525702186
#> 293      1         93    93 beta_denom.1.  2.3931855136
#> 294      1         94    94 beta_denom.1.  2.3956185226
#> 295      1         95    95 beta_denom.1.  2.3860444812
#> 296      1         96    96 beta_denom.1.  2.4261511402
#> 297      1         97    97 beta_denom.1.  2.3978429691
#> 298      1         98    98 beta_denom.1.  2.3560647391
#> 299      1         99    99 beta_denom.1.  2.4828987776
#> 300      1        100   100 beta_denom.1.  2.4780947674
#> 301      1          1     1 beta_denom.2. -0.0078161313
#> 302      1          2     2 beta_denom.2. -0.0022649586
#> 303      1          3     3 beta_denom.2.  0.0637693571
#> 304      1          4     4 beta_denom.2. -0.0494107701
#> 305      1          5     5 beta_denom.2.  0.0475237827
#> 306      1          6     6 beta_denom.2.  0.0977717853
#> 307      1          7     7 beta_denom.2.  0.0749902872
#> 308      1          8     8 beta_denom.2.  0.0965453818
#> 309      1          9     9 beta_denom.2.  0.0619326450
#> 310      1         10    10 beta_denom.2.  0.0126633258
#> 311      1         11    11 beta_denom.2.  0.0436041830
#> 312      1         12    12 beta_denom.2. -0.0260422249
#> 313      1         13    13 beta_denom.2. -0.0077166669
#> 314      1         14    14 beta_denom.2.  0.0075389908
#> 315      1         15    15 beta_denom.2.  0.0107760366
#> 316      1         16    16 beta_denom.2. -0.0440742928
#> 317      1         17    17 beta_denom.2.  0.0043565722
#> 318      1         18    18 beta_denom.2.  0.0509830151
#> 319      1         19    19 beta_denom.2.  0.0624905940
#> 320      1         20    20 beta_denom.2.  0.0065296670
#> 321      1         21    21 beta_denom.2.  0.0269811108
#> 322      1         22    22 beta_denom.2.  0.0011985183
#> 323      1         23    23 beta_denom.2.  0.0872154150
#> 324      1         24    24 beta_denom.2.  0.0855483126
#> 325      1         25    25 beta_denom.2. -0.0369418581
#> 326      1         26    26 beta_denom.2.  0.0621101536
#> 327      1         27    27 beta_denom.2.  0.0662859686
#> 328      1         28    28 beta_denom.2. -0.0349723874
#> 329      1         29    29 beta_denom.2. -0.0254451847
#> 330      1         30    30 beta_denom.2. -0.0738034601
#> 331      1         31    31 beta_denom.2. -0.0303943368
#> 332      1         32    32 beta_denom.2.  0.0076201839
#> 333      1         33    33 beta_denom.2.  0.0076201839
#> 334      1         34    34 beta_denom.2.  0.0303996396
#> 335      1         35    35 beta_denom.2.  0.0105455389
#> 336      1         36    36 beta_denom.2.  0.0691509500
#> 337      1         37    37 beta_denom.2.  0.0058163296
#> 338      1         38    38 beta_denom.2.  0.0122271280
#> 339      1         39    39 beta_denom.2.  0.0679062966
#> 340      1         40    40 beta_denom.2. -0.0082741763
#> 341      1         41    41 beta_denom.2.  0.0920257711
#> 342      1         42    42 beta_denom.2. -0.0163635141
#> 343      1         43    43 beta_denom.2.  0.0070829625
#> 344      1         44    44 beta_denom.2. -0.0011968698
#> 345      1         45    45 beta_denom.2. -0.0174506017
#> 346      1         46    46 beta_denom.2.  0.0164162476
#> 347      1         47    47 beta_denom.2. -0.0412710801
#> 348      1         48    48 beta_denom.2.  0.0016245865
#> 349      1         49    49 beta_denom.2.  0.0536890585
#> 350      1         50    50 beta_denom.2. -0.0355274678
#> 351      1         51    51 beta_denom.2.  0.0837253932
#> 352      1         52    52 beta_denom.2. -0.0143805553
#> 353      1         53    53 beta_denom.2.  0.0655042147
#> 354      1         54    54 beta_denom.2.  0.0463540168
#> 355      1         55    55 beta_denom.2. -0.0353056305
#> 356      1         56    56 beta_denom.2.  0.0030013391
#> 357      1         57    57 beta_denom.2. -0.0272641109
#> 358      1         58    58 beta_denom.2. -0.0554655836
#> 359      1         59    59 beta_denom.2. -0.1358273258
#> 360      1         60    60 beta_denom.2.  0.1267842824
#> 361      1         61    61 beta_denom.2.  0.0826108375
#> 362      1         62    62 beta_denom.2.  0.0749311463
#> 363      1         63    63 beta_denom.2. -0.0515768724
#> 364      1         64    64 beta_denom.2.  0.1130934729
#> 365      1         65    65 beta_denom.2.  0.0177795623
#> 366      1         66    66 beta_denom.2.  0.0245424909
#> 367      1         67    67 beta_denom.2.  0.0674614795
#> 368      1         68    68 beta_denom.2.  0.0343110002
#> 369      1         69    69 beta_denom.2.  0.0658671116
#> 370      1         70    70 beta_denom.2. -0.0061076515
#> 371      1         71    71 beta_denom.2.  0.0488995726
#> 372      1         72    72 beta_denom.2.  0.0583834179
#> 373      1         73    73 beta_denom.2.  0.0634563792
#> 374      1         74    74 beta_denom.2. -0.0329305445
#> 375      1         75    75 beta_denom.2. -0.0204427146
#> 376      1         76    76 beta_denom.2. -0.0076179133
#> 377      1         77    77 beta_denom.2.  0.0546765971
#> 378      1         78    78 beta_denom.2.  0.0026852685
#> 379      1         79    79 beta_denom.2.  0.0166023938
#> 380      1         80    80 beta_denom.2.  0.0130129996
#> 381      1         81    81 beta_denom.2.  0.0415093058
#> 382      1         82    82 beta_denom.2.  0.0406851352
#> 383      1         83    83 beta_denom.2. -0.0469143779
#> 384      1         84    84 beta_denom.2.  0.0313174487
#> 385      1         85    85 beta_denom.2.  0.0527196289
#> 386      1         86    86 beta_denom.2.  0.0412722008
#> 387      1         87    87 beta_denom.2.  0.0081892865
#> 388      1         88    88 beta_denom.2.  0.0378669048
#> 389      1         89    89 beta_denom.2.  0.0399440573
#> 390      1         90    90 beta_denom.2.  0.0017369397
#> 391      1         91    91 beta_denom.2. -0.0390511919
#> 392      1         92    92 beta_denom.2. -0.0759772464
#> 393      1         93    93 beta_denom.2.  0.0470265355
#> 394      1         94    94 beta_denom.2.  0.0305813763
#> 395      1         95    95 beta_denom.2.  0.0375877227
#> 396      1         96    96 beta_denom.2.  0.0222978750
#> 397      1         97    97 beta_denom.2.  0.0503010588
#> 398      1         98    98 beta_denom.2.  0.0190012811
#> 399      1         99    99 beta_denom.2.  0.0235121415
#> 400      1        100   100 beta_denom.2.  0.0369568942

# Use for plotting
if (requireNamespace("ggplot2", quietly = TRUE)) {
  library(ggplot2)
  gather_draws(fit, beta_num) |>
    ggplot(aes(x = .value)) +
    geom_histogram() +
    facet_wrap(~ .variable)
}
#> `stat_bin()` using `bins = 30`. Pick better value `binwidth`.

# }