Learn R Programming

rcompanion (version 2.4.30)

pairwisePermutationMatrix: Pairwise two-sample independence tests with matrix output

Description

Conducts pairwise two-sample independence tests across groups.

Usage

pairwisePermutationMatrix(
  formula = NULL,
  data = NULL,
  x = NULL,
  g = NULL,
  method = "fdr",
  ...
)

Value

A list consisting of: A matrix of p-values; the p-value adjustment method; a matrix of adjusted p-values.

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.

method

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

...

Additional arguments passed to coin::independence_test.

Author

Salvatore Mangiafico, mangiafico@njaes.rutgers.edu

Details

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

This function is a wrapper for coin::independence_test, passing pairwise groups to the function. It's critical to read and understand the documentation for this function to understand its use and options.

For some options for common tests, see Horthorn et al., 2008.

References

http://rcompanion.org/handbook/K_02.html

Hothorn, T., K. Hornik, M.A. van de Wiel, and A. Zeileis. 2008. Implementing a Class of Permutation Tests: The coin Package. Journal of Statistical Software, 28(8), 1–23.

See Also

pairwisePermutationTest

Examples

Run this code
### Fisher-Pitman test

data(BrendonSmall)

library(coin)
                                 
independence_test(Sodium ~ Instructor, data = BrendonSmall, 
                  teststat = "quadratic") 
                                      
PT = pairwisePermutationMatrix(Sodium ~ Instructor,
                               data     = BrendonSmall,
                               teststat = "quadratic",
                               method   = "fdr")
PT

PA = PT$Adjusted
library(multcompView)
multcompLetters(PA,
                compare="<",
                threshold=0.05,
                Letters=letters)   

Run the code above in your browser using DataLab