From 13c12f87b82badcf2b4c9964b06fc88ae1ddd3d7 Mon Sep 17 00:00:00 2001 From: Indrajeet Patil Date: Sun, 25 Oct 2020 13:12:55 +0100 Subject: [PATCH 1/2] initialize refactoring --- API | 3 +- DESCRIPTION | 4 +- NAMESPACE | 1 + R/bf_contingency_tab.R | 67 ++++--- R/bf_expr_template.R | 71 +++++++ R/bf_meta_random.R | 47 ++--- R/helpers_bf_tests.R | 52 ++--- README.md | 2 +- _pkgdown.yml | 1 + codemeta.json | 8 +- docs/index.html | 2 +- docs/pkgdown.yml | 2 +- docs/reference/bf_contingency_tab.html | 22 ++- docs/reference/bf_corr_test.html | 14 +- docs/reference/bf_expr.html | 5 +- docs/reference/bf_expr_template.html | 235 +++++++++++++++++++++++ docs/reference/bf_extractor.html | 44 ++--- docs/reference/bf_meta_random.html | 30 +-- docs/reference/bf_oneway_anova.html | 48 +++-- docs/reference/bf_ttest.html | 32 ++- docs/reference/bugs_long.html | 28 +-- docs/reference/index.html | 6 + docs/sitemap.xml | 3 + man/bf_contingency_tab.Rd | 6 +- man/bf_expr_template.Rd | 57 ++++++ tests/testthat/test-bf_contingency_tab.R | 77 ++++++-- tests/testthat/test-bf_oneway_anova.R | 1 + 27 files changed, 622 insertions(+), 246 deletions(-) create mode 100644 R/bf_expr_template.R create mode 100644 docs/reference/bf_expr_template.html create mode 100644 man/bf_expr_template.Rd diff --git a/API b/API index 5cd052b..4bfb8e7 100644 --- a/API +++ b/API @@ -2,9 +2,10 @@ ## Exported functions -bf_contingency_tab(data, x, y = NULL, counts = NULL, ratio = NULL, sampling.plan = "indepMulti", fixed.margin = "rows", prior.concentration = 1, top.text = NULL, output = "dataframe", k = 2L, ...) +bf_contingency_tab(data, x, y = NULL, counts = NULL, ratio = NULL, sampling.plan = "indepMulti", fixed.margin = "rows", prior.concentration = 1, top.text = NULL, output = "dataframe", conf.level = 0.95, k = 2L, ...) bf_corr_test(data, x, y, bf.prior = 0.707, top.text = NULL, output = "dataframe", k = 2L, ...) bf_expr(bf.object, k = 2L, conf.level = 0.95, conf.method = "hdi", centrality = "median", top.text = NULL, anova.design = FALSE, ...) +bf_expr_template(top.text, bf.value, bf.prior = 0.707, prior.type = quote(italic("r")["Cauchy"]^"JZS"), estimate.type = quote(delta), estimate, estimate.LB, estimate.UB, centrality = "median", conf.level = 0.95, conf.method = "HDI", k = 2L, ...) bf_extractor(bf.object, conf.method = "hdi", centrality = "median", conf.level = 0.95, ...) bf_meta(data, d = prior("norm", c(mean = 0, sd = 0.3)), tau = prior("invgamma", c(shape = 1, scale = 0.15)), k = 2L, conf.level = 0.95, output = "dataframe", top.text = NULL, ...) bf_meta_random(data, d = prior("norm", c(mean = 0, sd = 0.3)), tau = prior("invgamma", c(shape = 1, scale = 0.15)), k = 2L, conf.level = 0.95, output = "dataframe", top.text = NULL, ...) diff --git a/DESCRIPTION b/DESCRIPTION index 27fd4fe..3b0d3ec 100644 --- a/DESCRIPTION +++ b/DESCRIPTION @@ -21,8 +21,6 @@ URL: https://indrajeetpatil.github.io/tidyBF/, https://github.com/IndrajeetPatil/tidyBF BugReports: https://github.com/IndrajeetPatil/tidyBF/issues Remotes: - easystats/insight, - easystats/bayestestR, easystats/parameters, easystats/performance, easystats/effectsize, @@ -41,7 +39,7 @@ Imports: bayestestR (>= 0.7.5), dplyr, effectsize (>= 0.3.9.009), - insight (>= 0.10.1), + insight (>= 0.10.0), ipmisc (>= 4.1.0), lme4, metaBMA, diff --git a/NAMESPACE b/NAMESPACE index 7a40330..02e079f 100644 --- a/NAMESPACE +++ b/NAMESPACE @@ -11,6 +11,7 @@ export(as_tibble) export(bf_contingency_tab) export(bf_corr_test) export(bf_expr) +export(bf_expr_template) export(bf_extractor) export(bf_meta) export(bf_meta_random) diff --git a/R/bf_contingency_tab.R b/R/bf_contingency_tab.R index ef15bb9..24ab2f3 100644 --- a/R/bf_contingency_tab.R +++ b/R/bf_contingency_tab.R @@ -24,7 +24,8 @@ #' This means if there are two levels this will be `ratio = c(0.5,0.5)` or if #' there are four levels this will be `ratio = c(0.25,0.25,0.25,0.25)`, etc. #' @param counts A string naming a variable in data containing counts, or `NULL` -#' if each row represents a single observation (Default). +#' if each row represents a single observation. +#' @inheritParams bf_expr #' #' @importFrom BayesFactor contingencyTableBF logMeanExpLogs #' @importFrom dplyr pull select rename mutate tibble @@ -70,6 +71,7 @@ bf_contingency_tab <- function(data, prior.concentration = 1, top.text = NULL, output = "dataframe", + conf.level = 0.95, k = 2L, ...) { @@ -122,6 +124,21 @@ bf_contingency_tab <- function(data, fixed.margin = fixed.margin, prior.concentration = prior.concentration ) + + # Bayes Factor expression + bf01_expr <- + bf_expr_template( + top.text = top.text, + bf.value = -log(df$bf10[[1]]), + bf.prior = prior.concentration, + prior.type = quote(italic("a")["Gunel-Dickey"]), + estimate.type = quote(widehat(italic("V"))["Cramer"]), + estimate = df$estimate[[1]], + estimate.LB = df$conf.low[[1]], + estimate.UB = df$conf.high[[1]], + conf.level = conf.level, + k = k + ) } # ---------------------------- goodness of fit ---------------------------- @@ -175,35 +192,33 @@ bf_contingency_tab <- function(data, df <- tibble(bf10 = exp(pr_y_h1 - pr_y_h0)) %>% dplyr::mutate(log_e_bf10 = log(bf10), prior.concentration = prior.concentration) - } - # ========================= top.text preparation ============================= - - # final expression - bf01_expr <- - substitute( - atop( - displaystyle(top.text), - expr = paste( - "log"["e"], - "(BF"["01"], - ") = ", - bf, - ", ", - italic("a")["Gunel-Dickey"], - " = ", - a + # final expression + bf01_expr <- + substitute( + atop( + displaystyle(top.text), + expr = paste( + "log"["e"], + "(BF"["01"], + ") = ", + bf, + ", ", + italic("a")["Gunel-Dickey"], + " = ", + a + ) + ), + env = list( + top.text = top.text, + bf = specify_decimal_p(x = -log(df$bf10[[1]]), k = k), + a = specify_decimal_p(x = df$prior.concentration[[1]], k = k) ) - ), - env = list( - top.text = top.text, - bf = specify_decimal_p(x = -log(df$bf10[[1]]), k = k), - a = specify_decimal_p(x = df$prior.concentration[[1]], k = k) ) - ) - # the final expression - if (is.null(top.text)) bf01_expr <- bf01_expr$expr + # the final expression + if (is.null(top.text)) bf01_expr <- bf01_expr$expr + } # return the expression or the dataframe switch( diff --git a/R/bf_expr_template.R b/R/bf_expr_template.R new file mode 100644 index 0000000..2dd41a6 --- /dev/null +++ b/R/bf_expr_template.R @@ -0,0 +1,71 @@ +#' @title Expression template for Bayes Factor results +#' @name bf_expr_template +#' +#' @param bf.value Numeric value corresponding to Bayes Factor. +#' @param prior.type A character that specifies the prior type. +#' @param estimate.type A character that specifies the relevant effect size. +#' @param estimate,estimate.LB,estimate.UB Values of posterior estimates and +#' their credible intervals. +#' @param ... Currently ignored. +#' @inheritParams bf_expr +#' @inheritParams bf_ttest +#' +#' @export + +bf_expr_template <- function(top.text, + bf.value, + bf.prior = 0.707, + prior.type = quote(italic("r")["Cauchy"]^"JZS"), + estimate.type = quote(delta), + estimate, + estimate.LB, + estimate.UB, + centrality = "median", + conf.level = 0.95, + conf.method = "HDI", + k = 2L, + ...) { + # prepare the Bayes Factor message + bf01_expr <- + substitute( + atop( + displaystyle(top.text), + expr = paste( + "log"["e"], + "(BF"["01"], + ") = ", + bf, + ", ", + widehat(italic(estimate.type))[centrality]^"posterior", + " = ", + estimate, + ", CI"[conf.level]^conf.method, + " [", + estimate.LB, + ", ", + estimate.UB, + "]", + ", ", + prior.type, + " = ", + bf.prior + ) + ), + env = list( + top.text = top.text, + estimate.type = estimate.type, + centrality = centrality, + conf.level = paste0(conf.level * 100, "%"), + conf.method = toupper(conf.method), + bf = specify_decimal_p(x = bf.value, k = k), + estimate = specify_decimal_p(x = estimate, k = k), + estimate.LB = specify_decimal_p(x = estimate.LB, k = k), + estimate.UB = specify_decimal_p(x = estimate.UB, k = k), + prior.type = prior.type, + bf.prior = specify_decimal_p(x = bf.prior, k = k) + ) + ) + + # return the final expression + if (is.null(top.text)) bf01_expr$expr else bf01_expr +} diff --git a/R/bf_meta_random.R b/R/bf_meta_random.R index 34446b1..0232ff6 100644 --- a/R/bf_meta_random.R +++ b/R/bf_meta_random.R @@ -89,45 +89,20 @@ bf_meta_random <- function(data, as_tibble(mod$estimates, rownames = "term") %>% dplyr::mutate(.data = ., bf10 = mod$BF["random_H1", "random_H0"]) - # prepare the Bayes factor message + # Bayes Factor expression bf01_expr <- - substitute( - atop(displaystyle(top.text), - expr = paste( - "log"["e"], - "(BF"["01"], - ") = ", - bf, - ", ", - widehat(italic(delta))["mean"]^"posterior", - " = ", - estimate, - ", CI"[conf.level]^"HDI", - " [", - estimate.LB, - ", ", - estimate.UB, - "]", - ", ", - italic("r")["Cauchy"]^"JZS", - " = ", - bf.prior - ) - ), - env = list( - top.text = top.text, - bf = specify_decimal_p(x = -log(df$bf10[[1]]), k = k), - conf.level = paste0(conf.level * 100, "%"), - estimate = specify_decimal_p(x = df$mean[[1]], k = k), - estimate.LB = specify_decimal_p(x = df$hpd95_lower[[1]], k = k), - estimate.UB = specify_decimal_p(x = df$hpd95_upper[[1]], k = k), - bf.prior = specify_decimal_p(x = mod$jzs$rscale_discrete[[1]], k = k) - ) + bf_expr_template( + top.text = top.text, + bf.value = -log(df$bf10[[1]]), + bf.prior = mod$jzs$rscale_discrete[[1]], + estimate = df$mean[[1]], + estimate.LB = df$hpd95_lower[[1]], + estimate.UB = df$hpd95_upper[[1]], + centrality = "mean", + conf.level = conf.level, + k = k ) - # the final expression - if (is.null(top.text)) bf01_expr <- bf01_expr$expr - # return the text results or the dataframe with results switch( EXPR = output, diff --git a/R/helpers_bf_tests.R b/R/helpers_bf_tests.R index a5596de..2aa6fb3 100644 --- a/R/helpers_bf_tests.R +++ b/R/helpers_bf_tests.R @@ -182,47 +182,21 @@ bf_expr <- function(bf.object, c(df$r2[[1]], df$r2.conf.low[[1]], df$r2.conf.high[[1]], df_prior$prior.scale[[1]]) } - # prepare the Bayes Factor message + # Bayes Factor expression bf01_expr <- - substitute( - atop(displaystyle(top.text), - expr = paste( - "log"["e"], - "(BF"["01"], - ") = ", - bf, - ", ", - widehat(italic(estimate.type))[centrality]^"posterior", - " = ", - estimate, - ", CI"[conf.level]^conf.method, - " [", - estimate.LB, - ", ", - estimate.UB, - "]", - ", ", - italic("r")["Cauchy"]^"JZS", - " = ", - bf.prior - ) - ), - env = list( - top.text = top.text, - estimate.type = estimate.type, - centrality = centrality, - conf.level = paste0(conf.level * 100, "%"), - conf.method = toupper(conf.method), - bf = specify_decimal_p(x = -log(df$bf10[[1]]), k = k), - estimate = specify_decimal_p(x = estimate, k = k), - estimate.LB = specify_decimal_p(x = estimate.LB, k = k), - estimate.UB = specify_decimal_p(x = estimate.UB, k = k), - bf.prior = specify_decimal_p(x = bf.prior, k = k) - ) + bf_expr_template( + top.text = top.text, + bf.value = -log(df$bf10[[1]]), + bf.prior = bf.prior, + estimate.type = estimate.type, + estimate = estimate, + estimate.LB = estimate.LB, + estimate.UB = estimate.UB, + centrality = centrality, + conf.level = conf.level, + conf.method = conf.method, + k = k ) - - # return the final expression - if (is.null(top.text)) bf01_expr$expr else bf01_expr } #' @name meta_data_check diff --git a/README.md b/README.md index db437ec..637783b 100644 --- a/README.md +++ b/README.md @@ -21,7 +21,7 @@ Status](https://ci.appveyor.com/api/projects/status/github/IndrajeetPatil/tidyBF [![Project Status: Active - The project has reached a stable, usable state and is being actively developed.](http://www.repostatus.org/badges/latest/active.svg)](https://www.repostatus.org/#active) -[![Last-changedate](https://img.shields.io/badge/last%20change-2020--10--24-yellowgreen.svg)](https://github.com/IndrajeetPatil/tidyBF/commits/master) +[![Last-changedate](https://img.shields.io/badge/last%20change-2020--10--25-yellowgreen.svg)](https://github.com/IndrajeetPatil/tidyBF/commits/master) [![minimal R version](https://img.shields.io/badge/R%3E%3D-3.6.0-6666ff.svg)](https://cran.r-project.org/) [![Coverage diff --git a/_pkgdown.yml b/_pkgdown.yml index 033c030..583bf7b 100644 --- a/_pkgdown.yml +++ b/_pkgdown.yml @@ -26,6 +26,7 @@ reference: contents: - bf_expr - bf_extractor + - bf_expr_template - meta_data_check - title: Reexported functions diff --git a/codemeta.json b/codemeta.json index 00fb2e9..90f89df 100644 --- a/codemeta.json +++ b/codemeta.json @@ -140,7 +140,7 @@ "name": "Comprehensive R Archive Network (CRAN)", "url": "https://cran.r-project.org" }, - "sameAs": "https://github.com/easystats/bayestestR" + "sameAs": "https://CRAN.R-project.org/package=bayestestR" }, { "@type": "SoftwareApplication", @@ -171,14 +171,14 @@ "@type": "SoftwareApplication", "identifier": "insight", "name": "insight", - "version": ">= 0.10.1", + "version": ">= 0.10.0", "provider": { "@id": "https://cran.r-project.org", "@type": "Organization", "name": "Comprehensive R Archive Network (CRAN)", "url": "https://cran.r-project.org" }, - "sameAs": "https://github.com/easystats/insight" + "sameAs": "https://CRAN.R-project.org/package=insight" }, { "@type": "SoftwareApplication", @@ -283,5 +283,5 @@ "name": "Comprehensive R Archive Network (CRAN)", "url": "https://cran.r-project.org" }, - "fileSize": "1761.6KB" + "fileSize": "1787.638KB" } diff --git a/docs/index.html b/docs/index.html index 6dcc207..553fd04 100644 --- a/docs/index.html +++ b/docs/index.html @@ -507,7 +507,7 @@

Dev status

  • Licence
  • DOI
  • Project Status: Active - The project has reached a stable, usable state and is being actively developed.
    • -
  • Last-changedate
    • +
  • Last-changedate
  • minimal R version
  • Coverage Status
  • Coverage Status
    • diff --git a/docs/pkgdown.yml b/docs/pkgdown.yml index 600ac1a..63aacf3 100644 --- a/docs/pkgdown.yml +++ b/docs/pkgdown.yml @@ -2,7 +2,7 @@ pandoc: 2.11.0.2 pkgdown: 1.6.1 pkgdown_sha: ~ articles: {} -last_built: 2020-10-24T14:51Z +last_built: 2020-10-25T12:04Z urls: reference: https://indrajeetpatil.github.io/tidyBF//reference article: https://indrajeetpatil.github.io/tidyBF//articles diff --git a/docs/reference/bf_contingency_tab.html b/docs/reference/bf_contingency_tab.html index d9e4912..4b68c3e 100644 --- a/docs/reference/bf_contingency_tab.html +++ b/docs/reference/bf_contingency_tab.html @@ -141,6 +141,7 @@

    Bayes Factor for contingency table analysis

    prior.concentration = 1, top.text = NULL, output = "dataframe", + conf.level = 0.95, k = 2L, ... ) @@ -167,7 +168,7 @@

    Arg counts

    A string naming a variable in data containing counts, or NULL -if each row represents a single observation (Default).

    +if each row represents a single observation.

    ratio @@ -206,6 +207,11 @@

    Arg

    If "expression", will return expression with statistical details, while "dataframe" will return a dataframe containing the results.

    + + + conf.level +

    Confidence/Credible Interval (CI) level. Default to 0.95 +(95%).

    k @@ -237,16 +243,16 @@

    Examp y = cyl, output = "dataframe" ) -
    #> # A tibble: 1 x 15 +
    #> # A tibble: 1 x 15 #> bf10 log_e_bf10 term estimate ci.width conf.low conf.high pd rope.ci -#> <dbl> <dbl> <chr> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> -#> 1 16.8 2.82 Cramers_v 0.476 95 0.207 0.701 1 89 +#> <dbl> <dbl> <chr> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> +#> 1 16.8 2.82 Cramers_v 0.476 95 0.207 0.701 1 89 #> rope.low rope.high rope.percentage sampling.plan fixed.margin -#> <dbl> <dbl> <dbl> <chr> <chr> -#> 1 -0.1 0.1 0 indepMulti rows +#> <dbl> <dbl> <dbl> <chr> <chr> +#> 1 -0.1 0.1 0 indepMulti rows #> prior.concentration -#> <dbl> -#> 1 1
    +#> <dbl> +#> 1 1
    # ------------------ goodness of fit tests -------------------------------- # to get expression diff --git a/docs/reference/bf_corr_test.html b/docs/reference/bf_corr_test.html index 9c3c2f1..b29fc70 100644 --- a/docs/reference/bf_corr_test.html +++ b/docs/reference/bf_corr_test.html @@ -207,13 +207,13 @@

    Examp bf.prior = 0.3, output = "dataframe" ) -

    #> # A tibble: 1 x 13 +
    #> # A tibble: 1 x 13 #> term estimate conf.low conf.high pd rope.percentage prior.distribution -#> <chr> <dbl> <dbl> <dbl> <dbl> <dbl> <chr> -#> 1 rho -0.182 -0.614 0.296 0.773 0.266 cauchy +#> <chr> <dbl> <dbl> <dbl> <dbl> <dbl> <chr> +#> 1 rho -0.182 -0.614 0.296 0.773 0.266 cauchy #> prior.location prior.scale effect component bf10 log_e_bf10 -#> <dbl> <dbl> <chr> <chr> <dbl> <dbl> -#> 1 0 0.3 fixed conditional 0.828 -0.189
    +#> <dbl> <dbl> <chr> <chr> <dbl> <dbl> +#> 1 0 0.3 fixed conditional 0.828 -0.189
    # to get expression bf_corr_test( data = anscombe, @@ -222,9 +222,7 @@

    Examp bf.prior = 0.8, output = "expression" ) -

    #> paste("log"["e"], "(BF"["01"], ") = ", "0.50", ", ", widehat(italic(rho))["median"]^"posterior", -#> " = ", "-0.23", ", CI"["95%"]^"HDI", " [", "-0.72", ", ", -#> "0.30", "]", ", ", italic("r")["Cauchy"]^"JZS", " = ", "0.80")
    +
    #> atop(displaystyle("Note: Iris dataset"), expr = paste("log"["e"], -#> "(BF"["01"], ") = ", "-99.770", ", ", widehat(italic(rho))["median"]^"posterior", -#> " = ", "0.862", ", CI"["95%"]^"HDI", " [", "0.818", ", ", -#> "0.900", "]", ", ", italic("r")["Cauchy"]^"JZS", " = ", "0.333"))
    +
    #> # A tibble: 6 x 18 +
    #> # A tibble: 6 x 18 #> term estimate conf.low conf.high pd rope.percentage -#> <chr> <dbl> <dbl> <dbl> <dbl> <dbl> -#> 1 mu 5.84 5.76 5.93 1 0 -#> 2 Species-setosa -0.828 -0.941 -0.706 1 0 -#> 3 Species-versicolor 0.0916 -0.0275 0.208 0.936 0.434 -#> 4 Species-virginica 0.737 0.617 0.850 1 0 -#> 5 sig2 0.268 0.211 0.337 1 0 -#> 6 g_Species 2.04 0.179 14.2 1 0 +#> <chr> <dbl> <dbl> <dbl> <dbl> <dbl> +#> 1 mu 5.84 5.76 5.93 1 0 +#> 2 Species-setosa -0.828 -0.941 -0.706 1 0 +#> 3 Species-versicolor 0.0916 -0.0275 0.208 0.936 0.434 +#> 4 Species-virginica 0.737 0.617 0.850 1 0 +#> 5 sig2 0.268 0.211 0.337 1 0 +#> 6 g_Species 2.04 0.179 14.2 1 0 #> prior.distribution prior.location prior.scale effect component bf10 -#> <chr> <dbl> <dbl> <chr> <chr> <dbl> -#> 1 NA NA NA fixed extra 1.43e28 -#> 2 NA NA NA fixed conditional 1.43e28 -#> 3 NA NA NA fixed conditional 1.43e28 -#> 4 NA NA NA fixed conditional 1.43e28 -#> 5 NA NA NA fixed extra 1.43e28 -#> 6 NA NA NA fixed extra 1.43e28 +#> <chr> <dbl> <dbl> <chr> <chr> <dbl> +#> 1 NA NA NA fixed extra 1.43e28 +#> 2 NA NA NA fixed conditional 1.43e28 +#> 3 NA NA NA fixed conditional 1.43e28 +#> 4 NA NA NA fixed conditional 1.43e28 +#> 5 NA NA NA fixed extra 1.43e28 +#> 6 NA NA NA fixed extra 1.43e28 #> log_e_bf10 r2 std.dev ci.width r2.conf.low r2.conf.high -#> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> -#> 1 64.8 0.611 0.0311 95 0.544 0.667 -#> 2 64.8 0.611 0.0311 95 0.544 0.667 -#> 3 64.8 0.611 0.0311 95 0.544 0.667 -#> 4 64.8 0.611 0.0311 95 0.544 0.667 -#> 5 64.8 0.611 0.0311 95 0.544 0.667 -#> 6 64.8 0.611 0.0311 95 0.544 0.667
    +#> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> +#> 1 64.8 0.611 0.0311 95 0.544 0.667 +#> 2 64.8 0.611 0.0311 95 0.544 0.667 +#> 3 64.8 0.611 0.0311 95 0.544 0.667 +#> 4 64.8 0.611 0.0311 95 0.544 0.667 +#> 5 64.8 0.611 0.0311 95 0.544 0.667 +#> 6 64.8 0.611 0.0311 95 0.544 0.667
    #> Loading required package: Rcpp
    # creating a dataframe (df <- structure( @@ -260,14 +260,14 @@

    Examp row.names = c(NA, -5L), class = c("tbl_df", "tbl", "data.frame") )) -

    #> # A tibble: 5 x 3 +
    #> # A tibble: 5 x 3 #> study estimate std.error -#> <chr> <dbl> <dbl> -#> 1 1 0.382 0.0466 -#> 2 2 0.781 0.0330 -#> 3 3 0.426 0.0363 -#> 4 4 0.558 0.0481 -#> 5 5 0.956 0.0622
    +#> <chr> <dbl> <dbl> +#> 1 1 0.382 0.0466 +#> 2 2 0.781 0.0330 +#> 3 3 0.426 0.0363 +#> 4 4 0.558 0.0481 +#> 5 5 0.956 0.0622
    # to get dataframe bf_meta_random( data = df, @@ -279,15 +279,15 @@

    Examp )

    #> Warning: There were 2 divergent transitions after warmup. See #> http://mc-stan.org/misc/warnings.html#divergent-transitions-after-warmup -#> to find out why this is a problem and how to eliminate them.
    #> Warning: Examine the pairs() plot to diagnose sampling problems
    #> # A tibble: 2 x 11 +#> to find out why this is a problem and how to eliminate them.
    #> Warning: Examine the pairs() plot to diagnose sampling problems
    #> # A tibble: 2 x 11 #> term mean sd `2.5%` `50%` `97.5%` hpd95_lower hpd95_upper n_eff Rhat -#> <chr> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> -#> 1 d 0.535 0.119 0.273 0.544 0.751 0.287 0.762 682. 1.00 -#> 2 tau 0.277 0.118 0.133 0.249 0.602 0.119 0.518 431. 1.00 +#> <chr> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> +#> 1 d 0.535 0.119 0.273 0.544 0.751 0.287 0.762 682. 1.00 +#> 2 tau 0.277 0.118 0.133 0.249 0.602 0.119 0.518 431. 1.00 #> bf10 -#> <dbl> -#> 1 28.2 -#> 2 28.2
    # } +#> <dbl> +#> 1 28.2 +#> 2 28.2
    # }
    diff --git a/docs/reference/bf_oneway_anova.html b/docs/reference/bf_oneway_anova.html index 96097ad..5094dd3 100644 --- a/docs/reference/bf_oneway_anova.html +++ b/docs/reference/bf_oneway_anova.html @@ -221,31 +221,31 @@

    Examp bf.prior = 0.8, output = "dataframe" ) -
    #> # A tibble: 6 x 18 +
    #> # A tibble: 6 x 18 #> term estimate conf.low conf.high pd rope.percentage -#> <chr> <dbl> <dbl> <dbl> <dbl> <dbl> -#> 1 mu 5.84 5.76 5.93 1 0 -#> 2 Species-setosa -0.829 -0.943 -0.708 1 0 -#> 3 Species-versicolor 0.0917 -0.0277 0.208 0.937 0.434 -#> 4 Species-virginica 0.738 0.619 0.852 1 0 -#> 5 sig2 0.267 0.210 0.335 1 0 -#> 6 g_Species 2.21 0.290 15.5 1 0 +#> <chr> <dbl> <dbl> <dbl> <dbl> <dbl> +#> 1 mu 5.84 5.76 5.93 1 0 +#> 2 Species-setosa -0.829 -0.943 -0.708 1 0 +#> 3 Species-versicolor 0.0917 -0.0277 0.208 0.937 0.434 +#> 4 Species-virginica 0.738 0.619 0.852 1 0 +#> 5 sig2 0.267 0.210 0.335 1 0 +#> 6 g_Species 2.21 0.290 15.5 1 0 #> prior.distribution prior.location prior.scale effect component bf10 -#> <chr> <dbl> <dbl> <chr> <chr> <dbl> -#> 1 NA NA NA fixed extra 2.03e28 -#> 2 NA NA NA fixed conditional 2.03e28 -#> 3 NA NA NA fixed conditional 2.03e28 -#> 4 NA NA NA fixed conditional 2.03e28 -#> 5 NA NA NA fixed extra 2.03e28 -#> 6 NA NA NA fixed extra 2.03e28 +#> <chr> <dbl> <dbl> <chr> <chr> <dbl> +#> 1 NA NA NA fixed extra 2.03e28 +#> 2 NA NA NA fixed conditional 2.03e28 +#> 3 NA NA NA fixed conditional 2.03e28 +#> 4 NA NA NA fixed conditional 2.03e28 +#> 5 NA NA NA fixed extra 2.03e28 +#> 6 NA NA NA fixed extra 2.03e28 #> log_e_bf10 r2 std.dev ci.width r2.conf.low r2.conf.high -#> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> -#> 1 65.2 0.612 0.0310 95 0.545 0.668 -#> 2 65.2 0.612 0.0310 95 0.545 0.668 -#> 3 65.2 0.612 0.0310 95 0.545 0.668 -#> 4 65.2 0.612 0.0310 95 0.545 0.668 -#> 5 65.2 0.612 0.0310 95 0.545 0.668 -#> 6 65.2 0.612 0.0310 95 0.545 0.668
    +#> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> +#> 1 65.2 0.612 0.0310 95 0.545 0.668 +#> 2 65.2 0.612 0.0310 95 0.545 0.668 +#> 3 65.2 0.612 0.0310 95 0.545 0.668 +#> 4 65.2 0.612 0.0310 95 0.545 0.668 +#> 5 65.2 0.612 0.0310 95 0.545 0.668 +#> 6 65.2 0.612 0.0310 95 0.545 0.668
    # within-subjects ------------------------------------------------------- # to get expression @@ -257,9 +257,7 @@

    Examp paired = TRUE, output = "expression" ) -

    #> paste("log"["e"], "(BF"["01"], ") = ", "-21.04", ", ", widehat(italic(R^"2"))["median"]^"posterior", -#> " = ", "0.53", ", CI"["95%"]^"HDI", " [", "0.46", ", ", "0.59", -#> "]", ", ", italic("r")["Cauchy"]^"JZS", " = ", "0.71")
    +
    #> # A tibble: 1 x 13 +
    #> # A tibble: 1 x 13 #> term estimate conf.low conf.high pd rope.percentage -#> <chr> <dbl> <dbl> <dbl> <dbl> <dbl> -#> 1 Difference -1.27 -1.79 -0.728 1 0 +#> <chr> <dbl> <dbl> <dbl> <dbl> <dbl> +#> 1 Difference -1.27 -1.79 -0.728 1 0 #> prior.distribution prior.location prior.scale effect component bf10 -#> <chr> <dbl> <dbl> <chr> <chr> <dbl> -#> 1 cauchy 0 0.88 fixed conditional 1579. +#> <chr> <dbl> <dbl> <chr> <chr> <dbl> +#> 1 cauchy 0 0.88 fixed conditional 1579. #> log_e_bf10 -#> <dbl> -#> 1 7.36
    +#> <dbl> +#> 1 7.36
    # to get expression (within-subjects) bf_ttest( data = dplyr::filter(bugs_long, condition %in% c("LDLF", "LDHF")), @@ -256,16 +256,16 @@

    Examp bf.prior = 0.880, output = "dataframe" ) -

    #> # A tibble: 1 x 13 +
    #> # A tibble: 1 x 13 #> term estimate conf.low conf.high pd rope.percentage -#> <chr> <dbl> <dbl> <dbl> <dbl> <dbl> -#> 1 Difference -1.64 -2.12 -1.14 1 0 +#> <chr> <dbl> <dbl> <dbl> <dbl> <dbl> +#> 1 Difference -1.64 -2.12 -1.14 1 0 #> prior.distribution prior.location prior.scale effect component bf10 -#> <chr> <dbl> <dbl> <chr> <chr> <dbl> -#> 1 cauchy 0 0.88 fixed conditional 4605067. +#> <chr> <dbl> <dbl> <chr> <chr> <dbl> +#> 1 cauchy 0 0.88 fixed conditional 4605067. #> log_e_bf10 -#> <dbl> -#> 1 15.3
    +#> <dbl> +#> 1 15.3
    # ------------------- one-samples test ----------------------------------- # to get expression @@ -275,9 +275,7 @@

    Examp test.value = 5.85, output = "expression" ) -

    #> paste("log"["e"], "(BF"["01"], ") = ", "2.39", ", ", widehat(italic(delta))["median"]^"posterior", -#> " = ", "0.01", ", CI"["95%"]^"HDI", " [", "-0.12", ", ", -#> "0.15", "]", ", ", italic("r")["Cauchy"]^"JZS", " = ", "0.71")
    + -

    + theme(text = element_text(size = 12)) +#> Error: passing unknown arguments: conf.level.
    diff --git a/docs/pkgdown.yml b/docs/pkgdown.yml index 63aacf3..3d9bbfe 100644 --- a/docs/pkgdown.yml +++ b/docs/pkgdown.yml @@ -2,7 +2,7 @@ pandoc: 2.11.0.2 pkgdown: 1.6.1 pkgdown_sha: ~ articles: {} -last_built: 2020-10-25T12:04Z +last_built: 2020-10-25T15:01Z urls: reference: https://indrajeetpatil.github.io/tidyBF//reference article: https://indrajeetpatil.github.io/tidyBF//articles diff --git a/docs/reference/bf_contingency_tab.html b/docs/reference/bf_contingency_tab.html index 4b68c3e..90ba323 100644 --- a/docs/reference/bf_contingency_tab.html +++ b/docs/reference/bf_contingency_tab.html @@ -141,7 +141,6 @@

    Bayes Factor for contingency table analysis

    prior.concentration = 1, top.text = NULL, output = "dataframe", - conf.level = 0.95, k = 2L, ... ) @@ -207,11 +206,6 @@

    Arg

    If "expression", will return expression with statistical details, while "dataframe" will return a dataframe containing the results.

    - - - conf.level -

    Confidence/Credible Interval (CI) level. Default to 0.95 -(95%).

    k @@ -243,16 +237,13 @@

    Examp y = cyl, output = "dataframe" ) -

    #> # A tibble: 1 x 15 +
    #> # A tibble: 1 x 13 #> bf10 log_e_bf10 term estimate ci.width conf.low conf.high pd rope.ci -#> <dbl> <dbl> <chr> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> -#> 1 16.8 2.82 Cramers_v 0.476 95 0.207 0.701 1 89 -#> rope.low rope.high rope.percentage sampling.plan fixed.margin -#> <dbl> <dbl> <dbl> <chr> <chr> -#> 1 -0.1 0.1 0 indepMulti rows -#> prior.concentration -#> <dbl> -#> 1 1
    +#> <dbl> <dbl> <chr> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> +#> 1 16.8 2.82 Cramers_v 0.476 95 0.207 0.701 1 89 +#> rope.low rope.high rope.percentage prior.scale +#> <dbl> <dbl> <dbl> <dbl> +#> 1 -0.1 0.1 0 1
    # ------------------ goodness of fit tests -------------------------------- # to get expression diff --git a/docs/reference/bf_corr_test.html b/docs/reference/bf_corr_test.html index b29fc70..9c3c2f1 100644 --- a/docs/reference/bf_corr_test.html +++ b/docs/reference/bf_corr_test.html @@ -207,13 +207,13 @@

    Examp bf.prior = 0.3, output = "dataframe" ) -

    #> # A tibble: 1 x 13 +
    #> # A tibble: 1 x 13 #> term estimate conf.low conf.high pd rope.percentage prior.distribution -#> <chr> <dbl> <dbl> <dbl> <dbl> <dbl> <chr> -#> 1 rho -0.182 -0.614 0.296 0.773 0.266 cauchy +#> <chr> <dbl> <dbl> <dbl> <dbl> <dbl> <chr> +#> 1 rho -0.182 -0.614 0.296 0.773 0.266 cauchy #> prior.location prior.scale effect component bf10 log_e_bf10 -#> <dbl> <dbl> <chr> <chr> <dbl> <dbl> -#> 1 0 0.3 fixed conditional 0.828 -0.189
    +#> <dbl> <dbl> <chr> <chr> <dbl> <dbl> +#> 1 0 0.3 fixed conditional 0.828 -0.189
    # to get expression bf_corr_test( data = anscombe, @@ -222,7 +222,9 @@

    Examp bf.prior = 0.8, output = "expression" ) -

    +
    #> paste("log"["e"], "(BF"["01"], ") = ", "0.50", ", ", widehat(italic(rho))["median"]^"posterior", +#> " = ", "-0.23", ", CI"["95%"]^"HDI", " [", "-0.72", ", ", +#> "0.30", "]", ", ", italic("r")["Cauchy"]^"JZS", " = ", "0.80")
    +
    #> atop(displaystyle("Note: Iris dataset"), expr = paste("log"["e"], +#> "(BF"["01"], ") = ", "-99.770", ", ", widehat(italic(rho))["median"]^"posterior", +#> " = ", "0.862", ", CI"["95%"]^"HDI", " [", "0.818", ", ", +#> "0.900", "]", ", ", italic("r")["Cauchy"]^"JZS", " = ", "0.333"))
    #> # A tibble: 6 x 18 +
    #> # A tibble: 6 x 18 #> term estimate conf.low conf.high pd rope.percentage -#> <chr> <dbl> <dbl> <dbl> <dbl> <dbl> -#> 1 mu 5.84 5.76 5.93 1 0 -#> 2 Species-setosa -0.828 -0.941 -0.706 1 0 -#> 3 Species-versicolor 0.0916 -0.0275 0.208 0.936 0.434 -#> 4 Species-virginica 0.737 0.617 0.850 1 0 -#> 5 sig2 0.268 0.211 0.337 1 0 -#> 6 g_Species 2.04 0.179 14.2 1 0 +#> <chr> <dbl> <dbl> <dbl> <dbl> <dbl> +#> 1 mu 5.84 5.76 5.93 1 0 +#> 2 Species-setosa -0.828 -0.941 -0.706 1 0 +#> 3 Species-versicolor 0.0916 -0.0275 0.208 0.936 0.434 +#> 4 Species-virginica 0.737 0.617 0.850 1 0 +#> 5 sig2 0.268 0.211 0.337 1 0 +#> 6 g_Species 2.04 0.179 14.2 1 0 #> prior.distribution prior.location prior.scale effect component bf10 -#> <chr> <dbl> <dbl> <chr> <chr> <dbl> -#> 1 NA NA NA fixed extra 1.43e28 -#> 2 NA NA NA fixed conditional 1.43e28 -#> 3 NA NA NA fixed conditional 1.43e28 -#> 4 NA NA NA fixed conditional 1.43e28 -#> 5 NA NA NA fixed extra 1.43e28 -#> 6 NA NA NA fixed extra 1.43e28 +#> <chr> <dbl> <dbl> <chr> <chr> <dbl> +#> 1 NA NA NA fixed extra 1.43e28 +#> 2 NA NA NA fixed conditional 1.43e28 +#> 3 NA NA NA fixed conditional 1.43e28 +#> 4 NA NA NA fixed conditional 1.43e28 +#> 5 NA NA NA fixed extra 1.43e28 +#> 6 NA NA NA fixed extra 1.43e28 #> log_e_bf10 r2 std.dev ci.width r2.conf.low r2.conf.high -#> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> -#> 1 64.8 0.611 0.0311 95 0.544 0.667 -#> 2 64.8 0.611 0.0311 95 0.544 0.667 -#> 3 64.8 0.611 0.0311 95 0.544 0.667 -#> 4 64.8 0.611 0.0311 95 0.544 0.667 -#> 5 64.8 0.611 0.0311 95 0.544 0.667 -#> 6 64.8 0.611 0.0311 95 0.544 0.667
    +#> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> +#> 1 64.8 0.611 0.0311 95 0.544 0.667 +#> 2 64.8 0.611 0.0311 95 0.544 0.667 +#> 3 64.8 0.611 0.0311 95 0.544 0.667 +#> 4 64.8 0.611 0.0311 95 0.544 0.667 +#> 5 64.8 0.611 0.0311 95 0.544 0.667 +#> 6 64.8 0.611 0.0311 95 0.544 0.667
    #> Loading required package: Rcpp
    + # creating a dataframe (df <- structure( @@ -260,14 +253,14 @@

    Examp row.names = c(NA, -5L), class = c("tbl_df", "tbl", "data.frame") )) -

    #> # A tibble: 5 x 3 +
    #> # A tibble: 5 x 3 #> study estimate std.error -#> <chr> <dbl> <dbl> -#> 1 1 0.382 0.0466 -#> 2 2 0.781 0.0330 -#> 3 3 0.426 0.0363 -#> 4 4 0.558 0.0481 -#> 5 5 0.956 0.0622
    +#> <chr> <dbl> <dbl> +#> 1 1 0.382 0.0466 +#> 2 2 0.781 0.0330 +#> 3 3 0.426 0.0363 +#> 4 4 0.558 0.0481 +#> 5 5 0.956 0.0622
    # to get dataframe bf_meta_random( data = df, @@ -279,15 +272,15 @@

    Examp )

    #> Warning: There were 2 divergent transitions after warmup. See #> http://mc-stan.org/misc/warnings.html#divergent-transitions-after-warmup -#> to find out why this is a problem and how to eliminate them.
    #> Warning: Examine the pairs() plot to diagnose sampling problems
    #> # A tibble: 2 x 11 -#> term mean sd `2.5%` `50%` `97.5%` hpd95_lower hpd95_upper n_eff Rhat -#> <chr> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> -#> 1 d 0.535 0.119 0.273 0.544 0.751 0.287 0.762 682. 1.00 -#> 2 tau 0.277 0.118 0.133 0.249 0.602 0.119 0.518 431. 1.00 -#> bf10 -#> <dbl> -#> 1 28.2 -#> 2 28.2
    # } +#> to find out why this is a problem and how to eliminate them.
    #> Warning: Examine the pairs() plot to diagnose sampling problems
    #> # A tibble: 2 x 13 +#> term estimate sd `2.5%` `50%` `97.5%` conf.low conf.high n_eff Rhat bf10 +#> <chr> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> +#> 1 d 0.535 0.119 0.273 0.544 0.751 0.287 0.762 682. 1.00 28.2 +#> 2 tau 0.277 0.118 0.133 0.249 0.602 0.119 0.518 431. 1.00 28.2 +#> prior.scale log_e_bf10 +#> <dbl> <dbl> +#> 1 0.707 3.34 +#> 2 0.707 3.34
    # }
    diff --git a/docs/reference/bf_oneway_anova.html b/docs/reference/bf_oneway_anova.html index 5094dd3..96097ad 100644 --- a/docs/reference/bf_oneway_anova.html +++ b/docs/reference/bf_oneway_anova.html @@ -221,31 +221,31 @@

    Examp bf.prior = 0.8, output = "dataframe" ) -
    #> # A tibble: 6 x 18 +
    #> # A tibble: 6 x 18 #> term estimate conf.low conf.high pd rope.percentage -#> <chr> <dbl> <dbl> <dbl> <dbl> <dbl> -#> 1 mu 5.84 5.76 5.93 1 0 -#> 2 Species-setosa -0.829 -0.943 -0.708 1 0 -#> 3 Species-versicolor 0.0917 -0.0277 0.208 0.937 0.434 -#> 4 Species-virginica 0.738 0.619 0.852 1 0 -#> 5 sig2 0.267 0.210 0.335 1 0 -#> 6 g_Species 2.21 0.290 15.5 1 0 +#> <chr> <dbl> <dbl> <dbl> <dbl> <dbl> +#> 1 mu 5.84 5.76 5.93 1 0 +#> 2 Species-setosa -0.829 -0.943 -0.708 1 0 +#> 3 Species-versicolor 0.0917 -0.0277 0.208 0.937 0.434 +#> 4 Species-virginica 0.738 0.619 0.852 1 0 +#> 5 sig2 0.267 0.210 0.335 1 0 +#> 6 g_Species 2.21 0.290 15.5 1 0 #> prior.distribution prior.location prior.scale effect component bf10 -#> <chr> <dbl> <dbl> <chr> <chr> <dbl> -#> 1 NA NA NA fixed extra 2.03e28 -#> 2 NA NA NA fixed conditional 2.03e28 -#> 3 NA NA NA fixed conditional 2.03e28 -#> 4 NA NA NA fixed conditional 2.03e28 -#> 5 NA NA NA fixed extra 2.03e28 -#> 6 NA NA NA fixed extra 2.03e28 +#> <chr> <dbl> <dbl> <chr> <chr> <dbl> +#> 1 NA NA NA fixed extra 2.03e28 +#> 2 NA NA NA fixed conditional 2.03e28 +#> 3 NA NA NA fixed conditional 2.03e28 +#> 4 NA NA NA fixed conditional 2.03e28 +#> 5 NA NA NA fixed extra 2.03e28 +#> 6 NA NA NA fixed extra 2.03e28 #> log_e_bf10 r2 std.dev ci.width r2.conf.low r2.conf.high -#> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> -#> 1 65.2 0.612 0.0310 95 0.545 0.668 -#> 2 65.2 0.612 0.0310 95 0.545 0.668 -#> 3 65.2 0.612 0.0310 95 0.545 0.668 -#> 4 65.2 0.612 0.0310 95 0.545 0.668 -#> 5 65.2 0.612 0.0310 95 0.545 0.668 -#> 6 65.2 0.612 0.0310 95 0.545 0.668
    +#> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> +#> 1 65.2 0.612 0.0310 95 0.545 0.668 +#> 2 65.2 0.612 0.0310 95 0.545 0.668 +#> 3 65.2 0.612 0.0310 95 0.545 0.668 +#> 4 65.2 0.612 0.0310 95 0.545 0.668 +#> 5 65.2 0.612 0.0310 95 0.545 0.668 +#> 6 65.2 0.612 0.0310 95 0.545 0.668
    # within-subjects ------------------------------------------------------- # to get expression @@ -257,7 +257,9 @@

    Examp paired = TRUE, output = "expression" ) -

    +
    #> paste("log"["e"], "(BF"["01"], ") = ", "-21.04", ", ", widehat(italic(R^"2"))["median"]^"posterior", +#> " = ", "0.53", ", CI"["95%"]^"HDI", " [", "0.46", ", ", "0.59", +#> "]", ", ", italic("r")["Cauchy"]^"JZS", " = ", "0.71")
    #> # A tibble: 1 x 13 +
    #> # A tibble: 1 x 13 #> term estimate conf.low conf.high pd rope.percentage -#> <chr> <dbl> <dbl> <dbl> <dbl> <dbl> -#> 1 Difference -1.27 -1.79 -0.728 1 0 +#> <chr> <dbl> <dbl> <dbl> <dbl> <dbl> +#> 1 Difference -1.27 -1.79 -0.728 1 0 #> prior.distribution prior.location prior.scale effect component bf10 -#> <chr> <dbl> <dbl> <chr> <chr> <dbl> -#> 1 cauchy 0 0.88 fixed conditional 1579. +#> <chr> <dbl> <dbl> <chr> <chr> <dbl> +#> 1 cauchy 0 0.88 fixed conditional 1579. #> log_e_bf10 -#> <dbl> -#> 1 7.36
    +#> <dbl> +#> 1 7.36
    # to get expression (within-subjects) bf_ttest( data = dplyr::filter(bugs_long, condition %in% c("LDLF", "LDHF")), @@ -256,16 +256,16 @@

    Examp bf.prior = 0.880, output = "dataframe" ) -

    #> # A tibble: 1 x 13 +
    #> # A tibble: 1 x 13 #> term estimate conf.low conf.high pd rope.percentage -#> <chr> <dbl> <dbl> <dbl> <dbl> <dbl> -#> 1 Difference -1.64 -2.12 -1.14 1 0 +#> <chr> <dbl> <dbl> <dbl> <dbl> <dbl> +#> 1 Difference -1.64 -2.12 -1.14 1 0 #> prior.distribution prior.location prior.scale effect component bf10 -#> <chr> <dbl> <dbl> <chr> <chr> <dbl> -#> 1 cauchy 0 0.88 fixed conditional 4605067. +#> <chr> <dbl> <dbl> <chr> <chr> <dbl> +#> 1 cauchy 0 0.88 fixed conditional 4605067. #> log_e_bf10 -#> <dbl> -#> 1 15.3
    +#> <dbl> +#> 1 15.3
    # ------------------- one-samples test ----------------------------------- # to get expression @@ -275,7 +275,9 @@

    Examp test.value = 5.85, output = "expression" ) -

    +
    #> paste("log"["e"], "(BF"["01"], ") = ", "2.39", ", ", widehat(italic(delta))["median"]^"posterior", +#> " = ", "0.01", ", CI"["95%"]^"HDI", " [", "-0.12", ", ", +#> "0.15", "]", ", ", italic("r")["Cauchy"]^"JZS", " = ", "0.71")