View source: R/pairwise_comparisons.R
pairwise_comparisons | R Documentation |
Calculate parametric, non-parametric, robust, and Bayes Factor pairwise comparisons between group levels with corrections for multiple testing.
pairwise_comparisons( data, x, y, subject.id = NULL, type = "parametric", paired = FALSE, var.equal = FALSE, tr = 0.2, bf.prior = 0.707, p.adjust.method = "holm", k = 2L, ... )
data |
A dataframe (or a tibble) from which variables specified are to be taken. Other data types (e.g., matrix,table, array, etc.) will not be accepted. |
x |
The grouping (or independent) variable from the dataframe |
y |
The response (or outcome or dependent) variable from the
dataframe |
subject.id |
Relevant in case of a repeated measures or within-subjects
design ( |
type |
Type of statistic expected ( |
paired |
Logical that decides whether the experimental design is
repeated measures/within-subjects or between-subjects. The default is
|
var.equal |
a logical variable indicating whether to treat the
two variances as being equal. If |
tr |
Trim level for the mean when carrying out |
bf.prior |
A number between |
p.adjust.method |
Adjustment method for p-values for multiple
comparisons. Possible methods are: |
k |
Number of digits after decimal point (should be an integer)
(Default: |
... |
Additional arguments passed to other methods. |
A tibble dataframe containing two columns corresponding to group
levels being compared with each other (group1
and group2
) and p.value
column corresponding to this comparison. The dataframe will also contain a
p.value.label
column containing a label for this p-value, in case
this needs to be displayed in ggsignif::geom_ggsignif
. In addition to
these common columns across the different types of statistics, there will
be additional columns specific to the type
of test being run.
This function provides a unified syntax to carry out pairwise comparison tests and internally relies on other packages to carry out these tests. For more details about the included tests, see the documentation for the respective functions:
parametric : stats::pairwise.t.test()
(paired) and
PMCMRplus::gamesHowellTest()
(unpaired)
non-parametric :
PMCMRplus::durbinAllPairsTest()
(paired) and
PMCMRplus::kwAllPairsDunnTest()
(unpaired)
robust :
WRS2::rmmcp()
(paired) and WRS2::lincon()
(unpaired)
Bayes Factor : BayesFactor::ttestBF()
# for reproducibility set.seed(123) library(pairwiseComparisons) library(statsExpressions) # for data # show all columns and make the column titles bold # as a user, you don't need to do this; this is just for the package website options(tibble.width = Inf, pillar.bold = TRUE, pillar.neg = TRUE, pillar.subtle_num = TRUE) #------------------- between-subjects design ---------------------------- # parametric # if `var.equal = TRUE`, then Student's t-test will be run pairwise_comparisons( data = mtcars, x = cyl, y = wt, type = "parametric", var.equal = TRUE, paired = FALSE, p.adjust.method = "none" ) # if `var.equal = FALSE`, then Games-Howell test will be run pairwise_comparisons( data = mtcars, x = cyl, y = wt, type = "parametric", var.equal = FALSE, paired = FALSE, p.adjust.method = "bonferroni" ) # non-parametric (Dunn test) pairwise_comparisons( data = mtcars, x = cyl, y = wt, type = "nonparametric", paired = FALSE, p.adjust.method = "none" ) # robust (Yuen's trimmed means t-test) pairwise_comparisons( data = mtcars, x = cyl, y = wt, type = "robust", paired = FALSE, p.adjust.method = "fdr" ) # Bayes Factor (Student's t-test) pairwise_comparisons( data = mtcars, x = cyl, y = wt, type = "bayes", paired = FALSE ) #------------------- within-subjects design ---------------------------- # parametric (Student's t-test) pairwise_comparisons( data = bugs_long, x = condition, y = desire, subject.id = subject, type = "parametric", paired = TRUE, p.adjust.method = "BH" ) # non-parametric (Durbin-Conover test) pairwise_comparisons( data = bugs_long, x = condition, y = desire, subject.id = subject, type = "nonparametric", paired = TRUE, p.adjust.method = "BY" ) # robust (Yuen's trimmed means t-test) pairwise_comparisons( data = bugs_long, x = condition, y = desire, subject.id = subject, type = "robust", paired = TRUE, p.adjust.method = "hommel" ) # Bayes Factor (Student's t-test) pairwise_comparisons( data = bugs_long, x = condition, y = desire, subject.id = subject, type = "bayes", paired = TRUE )
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.