pairwiseMedianTest: Pairwise Mood's median tests

View source: R/pairwiseMedianTest.r

pairwiseMedianTestR Documentation

Pairwise Mood's median tests

Description

Conducts pairwise Mood's median tests across groups.

Usage

pairwiseMedianTest(
  formula = NULL,
  data = NULL,
  x = NULL,
  g = NULL,
  digits = 4,
  method = "fdr",
  ...
)

Arguments

formula

A formula indicating the measurement variable and the grouping variable. e.g. y ~ group.

data

The data frame to use.

x

The response variable as a vector.

g

The grouping variable as a vector.

digits

The number of significant digits to round output.

method

The p-value adjustment method to use for multiple tests. See stats::p.adjust.

...

Additional arguments passed to coin::median_test.

Details

The input should include either formula and data; or x, and g.

Mood's median test compares medians among two or more groups. See https://rcompanion.org/handbook/F_09.html for further discussion of this test.

The pairwiseMedianTest function can be used as a post-hoc method following an omnibus Mood's median test. It passes the data for pairwise groups to coin::median_test.

The output can be converted to a compact letter display, as in the example.

Value

A dataframe of the groups being compared, the p-values, and the adjusted p-values.

Note

The parsing of the formula is simplistic. The first variable on the left side is used as the measurement variable. The first variable on the right side is used for the grouping variable.

Author(s)

Salvatore Mangiafico, mangiafico@njaes.rutgers.edu

References

https://rcompanion.org/handbook/F_09.html

See Also

pairwiseMedianMatrix

Examples

data(PoohPiglet)
PoohPiglet$Speaker = factor(PoohPiglet$Speaker,
                     levels = c("Pooh", "Tigger", "Piglet"))
PT = pairwiseMedianTest(Likert ~ Speaker,
                        data   = PoohPiglet,
                        exact  = NULL,
                        method = "fdr")
PT                         
cldList(comparison = PT$Comparison,
        p.value    = PT$p.adjust,
        threshold  = 0.05)                         


rcompanion documentation built on Sept. 17, 2023, 5:07 p.m.