Nothing
## ----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)
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.