Stuart-Maxwell-test for marginal homogeneity

Install required packages

coin, DescTools

wants <- c("coin", "DescTools")
has   <- wants %in% rownames(installed.packages())
if(any(!has)) install.packages(wants[!has])

MH-test

categ <- factor(1:3, labels=c("lo", "med", "hi"))
drug  <- rep(categ, c(30, 50, 20))
plac  <- rep(rep(categ, length(categ)), c(14,7,9, 5,26,19, 1,7,12))
cTab  <- table(drug, plac)
addmargins(cTab)

     plac
drug   lo med  hi Sum
  lo   14   7   9  30
  med   5  26  19  50
  hi    1   7  12  20
  Sum  20  40  40 100

library(coin)
mh_test(cTab, distribution=approximate(nresample=9999))

    Approximative Marginal Homogeneity Test

data:  response by
     conditions (drug, plac) 
     stratified by block
chi-squared = 12.139, p-value = 0.0011

Bhapkar-Test

library(DescTools)
BhapkarTest(cTab)

    Bhapkar test

data:  cTab
chi-squared = 13.816, df = 2, p-value = 0.0009999

Detach (automatically) loaded packages (if possible)

try(detach(package:coin))
try(detach(package:survival))
try(detach(package:DescTools))

Get the article source from GitHub

R markdown - markdown - R code - all posts