inst/doc/samplesizeCMH-introduction.R

In samplesizeCMH: Power and Sample Size Calculation for the Cochran-Mantel-Haenszel Test

## ----load-data-----------------------------------------------------------
data(Titanic, package = "datasets")
str(Titanic)

## ----partial-tables------------------------------------------------------
partial_tables <- margin.table(Titanic, c(2,4,1))
partial_tables

## ----marginal-table------------------------------------------------------
marginal_table <- margin.table(Titanic, c(2,4))
marginal_table

## ----prop-table----------------------------------------------------------
# Table proportions
prop.table(marginal_table)

# Row proportions
prop.table(marginal_table, 1)

# Column proportions
prop.table(marginal_table, 2)

## ----marginal-OR---------------------------------------------------------
library(samplesizeCMH)

odds.ratio(marginal_table)

## ----conditional-OR------------------------------------------------------
apply(partial_tables, 3, odds.ratio)

## ----mantelhaen----------------------------------------------------------
mantelhaen.test(partial_tables)

## ----BreslowDay----------------------------------------------------------
library(DescTools)

BreslowDayTest(x = partial_tables, OR = 1)