Making expression for correlation analysis

expr_corr_test(
data,
x,
y,
k = 2L,
conf.level = 0.95,
beta = 0.1,
type = "parametric",
bf.prior = 0.707,
output = "expression",
...
)

Arguments

data A dataframe (or a tibble) from which variables specified are to be taken. A matrix or tables will not be accepted. The column in data containing the explanatory variable to be plotted on the x-axis. Can be entered either as a character string (e.g., "x") or as a bare expression (e.g, x). The column in data containing the response (outcome) variable to be plotted on the y-axis. Can be entered either as a character string (e.g., "y") or as a bare expression (e.g, y). Number of digits after decimal point (should be an integer) (Default: k = 2L). Scalar between 0 and 1. If unspecified, the defaults return 95% lower and upper confidence intervals (0.95). bending constant (Default: 0.1). For more, see WRS2::pbcor(). Type of association between paired samples required (""parametric": Pearson's product moment correlation coefficient" or ""nonparametric": Spearman's rho" or ""robust": percentage bend correlation coefficient" or ""bayes": Bayes Factor for Pearson's r"). Corresponding abbreviations are also accepted: "p" (for parametric/pearson), "np" (nonparametric/spearman), "r" (robust), "bf" (for bayes factor), resp. A number between 0.5 and 2 (default 0.707), the prior width to use in calculating Bayes factors. If "expression", will return expression with statistical details, while "dataframe" will return a dataframe containing the results. Additional arguments passed to parameters::model_parameters.BFBayesFactor().

Value

Expression containing results from correlation test with confidence intervals for the correlation coefficient estimate. Results are extracted via correlation::correlation.

References

For more details, see- https://indrajeetpatil.github.io/statsExpressions/articles/stats_details.html

Examples

# for reproducibility
set.seed(123)
library(statsExpressions)

# without changing defaults
expr_corr_test(
data = ggplot2::midwest,
x = area,
y = percblack,
type = "parametric"
)
#> paste(italic("t")["Student"], "(", "435", ") = ", "-0.97", ", ",
#>     italic("p"), " = ", "0.332", ", ", widehat(italic("r"))["Pearson"],
#>     " = ", "-0.05", ", CI"["95%"], " [", "-0.14", ", ", "0.05",
#>     "]", ", ", italic("n")["pairs"], " = ", 437L)
# changing defaults
expr_corr_test(
data = ggplot2::midwest,
x = area,
y = percblack,
beta = 0.2,
type = "robust"
)
#> paste(italic("t")["Student"], "(", "435", ") = ", "0.05", ", ",
#>     italic("p"), " = ", "0.960", ", ", widehat(italic(rho))["pb"],
#>     " = ", "0.00", ", CI"["95%"], " [", "-0.09", ", ", "0.10",
#>     "]", ", ", italic("n")["pairs"], " = ", 437L)