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.

x

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).

y

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).

k

Number of digits after decimal point (should be an integer) (Default: k = 2L).

conf.level

Scalar between 0 and 1. If unspecified, the defaults return 95% lower and upper confidence intervals (0.95).

beta

bending constant (Default: 0.1). For more, see WRS2::pbcor().

type

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.

bf.prior

A number between 0.5 and 2 (default 0.707), the prior width to use in calculating Bayes factors.

output

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)