Skip to contents

Obtain posterior draws of meta-d' model parameters

Source: R/linpred_draws_metad.R
linpred_draws_metad.Rd

Given a data frame and a meta-d' model, adds estimates of all model parameters. For linpred_draws_metad and add_linpred_draws_metad, parameters are returned in a tidy tibble with one row per posterior draw. For linpred_rvars_metad and add_linpred_rvars_metad, parameters are returned as posterior::rvars, with one row per row in newdata.

Usage

linpred_draws_metad(object, newdata, ..., pivot_longer = FALSE)

add_linpred_draws_metad(newdata, object, ..., pivot_longer = FALSE)

linpred_rvars_metad(object, newdata, ..., pivot_longer = FALSE)

add_linpred_rvars_metad(newdata, object, ..., pivot_longer = FALSE)

Arguments

object

The brms model with the metad family

newdata

A data frame from which to generate posterior predictions

...

Additional arguments passed to tidybayes::add_linpred_draws or tidybayes::add_linpred_rvars

pivot_longer

Return the draws in long format?

  • if TRUE, resulting data frame has one row per posterior draw per model parameter

  • if FALSE (default), resulting data frame has one row per posterior draw

Value

a tibble containing posterior draws of model parameters with the following columns:

  • .row: the row of newdata

  • .chain, .iteration, .draw: for linpred_draws_metad, identifiers for the posterior sample

  • .variable, .value: if pivot_longer=TRUE, .variable identifies different meta-d' model parameters and .value stores posterior samples

  • M, dprime, c, meta_dprime, meta_c, meta_c2_0_<k>, meta_c2_1_<k>: if pivot_longer=FALSE, posterior samples of all meta-d' model parameters

See also

tidybayes::linpred_draws(), tidybayes::linpred_rvars()

Examples

# \donttest{
newdata <- tidyr::tibble(.row = 1)

# obtain model parameters (wide format)
# equivalent to `add_linpred_draws_metad(newdata, example_model())`
linpred_draws_metad(example_model(), newdata)
#> # A tibble: 1,000 × 15
#> # Groups:   .row [1]
#>     .row .chain .iteration .draw     M dprime        c meta_dprime   meta_c
#>    <int>  <int>      <int> <int> <dbl>  <dbl>    <dbl>       <dbl>    <dbl>
#>  1     1     NA         NA     1 1.07   0.977  0.0754        1.05   0.0754 
#>  2     1     NA         NA     2 0.923  1.02  -0.0103        0.938 -0.0103 
#>  3     1     NA         NA     3 1.07   0.953  0.0707        1.02   0.0707 
#>  4     1     NA         NA     4 1.33   0.972  0.0192        1.29   0.0192 
#>  5     1     NA         NA     5 1.19   0.894  0.0137        1.06   0.0137 
#>  6     1     NA         NA     6 1.01   1.04   0.0365        1.05   0.0365 
#>  7     1     NA         NA     7 1.20   0.959 -0.00256       1.15  -0.00256
#>  8     1     NA         NA     8 1.03   1.02   0.0226        1.05   0.0226 
#>  9     1     NA         NA     9 1.15   1.06   0.0646        1.21   0.0646 
#> 10     1     NA         NA    10 0.974  1.01   0.0836        0.984  0.0836 
#> # ℹ 990 more rows
#> # ℹ 6 more variables: meta_c2_0_1 <dbl>, meta_c2_0_2 <dbl>, meta_c2_0_3 <dbl>,
#> #   meta_c2_1_1 <dbl>, meta_c2_1_2 <dbl>, meta_c2_1_3 <dbl>

# obtain model parameters (long format)
# equivalent to `add_linpred_draws_metad(newdata, example_model(), pivot_longer = TRUE)`
linpred_draws_metad(example_model(), newdata, pivot_longer = TRUE)
#> # A tibble: 11,000 × 6
#> # Groups:   .row, .variable [11]
#>     .row .chain .iteration .draw .variable    .value
#>    <int>  <int>      <int> <int> <chr>         <dbl>
#>  1     1     NA         NA     1 M            1.07  
#>  2     1     NA         NA     1 dprime       0.977 
#>  3     1     NA         NA     1 c            0.0754
#>  4     1     NA         NA     1 meta_dprime  1.05  
#>  5     1     NA         NA     1 meta_c       0.0754
#>  6     1     NA         NA     1 meta_c2_0_1 -0.429 
#>  7     1     NA         NA     1 meta_c2_0_2 -0.967 
#>  8     1     NA         NA     1 meta_c2_0_3 -1.65  
#>  9     1     NA         NA     1 meta_c2_1_1  0.561 
#> 10     1     NA         NA     1 meta_c2_1_2  1.06  
#> # ℹ 10,990 more rows

# obtain model parameters (wide format, posterior::rvar)
# equivalent to `add_linpred_rvars_metad(newdata, example_model())`
linpred_rvars_metad(example_model(), newdata)
#> # A tibble: 1 × 12
#> # Groups:   .row [1]
#>    .row           M     dprime              c meta_dprime         meta_c
#>   <dbl>  <rvar[1d]> <rvar[1d]>     <rvar[1d]>  <rvar[1d]>     <rvar[1d]>
#> 1     1  1.1 ± 0.15  1 ± 0.083  0.025 ± 0.042  1.1 ± 0.12  0.025 ± 0.042
#> # ℹ 6 more variables: meta_c2_0_1 <rvar[1d]>, meta_c2_0_2 <rvar[1d]>,
#> #   meta_c2_0_3 <rvar[1d]>, meta_c2_1_1 <rvar[1d]>, meta_c2_1_2 <rvar[1d]>,
#> #   meta_c2_1_3 <rvar[1d]>

# obtain model parameters (long format, posterior::rvar)
# equivalent to `add_linpred_rvars_metad(newdata, example_model(), pivot_longer = TRUE)`
linpred_rvars_metad(example_model(), newdata, pivot_longer = TRUE)
#> # A tibble: 11 × 3
#> # Groups:   .row, .variable [11]
#>     .row .variable            .value
#>    <dbl> <chr>            <rvar[1d]>
#>  1     1 M             1.094 ± 0.149
#>  2     1 dprime        1.025 ± 0.083
#>  3     1 c             0.025 ± 0.042
#>  4     1 meta_dprime   1.114 ± 0.121
#>  5     1 meta_c        0.025 ± 0.042
#>  6     1 meta_c2_0_1  -0.453 ± 0.043
#>  7     1 meta_c2_0_2  -0.967 ± 0.050
#>  8     1 meta_c2_0_3  -1.601 ± 0.065
#>  9     1 meta_c2_1_1   0.504 ± 0.041
#> 10     1 meta_c2_1_2   1.042 ± 0.048
#> 11     1 meta_c2_1_3   1.585 ± 0.062
# }