Calculate the c-value of the difference in PCCs produced by two hypotheses

Usage

compare_hypotheses(m1, m2, two_tailed)

Arguments

m1: an object of class "opafit" produced by a call to opa().
m2: an object of class "opafit" produced by a call to opa().
two_tailed: a boolean indicating whether the comparison is two-tailed.

Value

an object of class "opaHypothesisComparison".

Examples

dat <- data.frame(t1 = c(9, 4, 8, 10),
                  t2 = c(8, 8, 12, 10),
                  t3 = c(8, 5, 10, 11),
                  t4 = c(10, 5, 11, 12))
h1 <- hypothesis(c(1, 2, 3, 4))
h2 <- hypothesis(c(1, 4, 2, 3))
opamod1 <- opa(dat, h1)
opamod2 <- opa(dat, h2)
compare_hypotheses(opamod1, opamod2, two_tailed = TRUE)
#> ********* Hypothesis Comparison **********
#> H1: 1 2 3 4 
#> H2: 1 4 2 3 
#> H1 PCC: 62.5 
#> H2 PCC: 66.66667 
#> PCC difference: 4.166667 
#> cval: 0.895 
#> Comparison type: two-tailed