if (!require(fmsb)) {install.packages('fmsb', dep=T); require(fmsb)}
if (!require(vcd)) {install.packages('vcd', dep=T); require(vcd)}
data <- data.frame(
id=rep(c('MALE','FEMALE'),each=20),
resp=c("N", "S", "N", "N", "S", "S", "S", "S", "S", "N", "S", "S",
"S", "S", "SR", "SR", "N", "SR", "S", "S", "S", "N", "N", "S",
"S", "N", "S", "S", "S", "S", "S", "N", "S", "S", "S", "N", "S",
"SR", "S", "SR"))
table(data$id)
MALE <- data$resp[data$id=='MALE']
FEMALE <- data$resp[data$id=='FEMALE']
tab <- table(FEMALE,MALE); tab
tmp <- rbind(tab, FEMALE=margin.table(tab, 1))
tabtot <- cbind(tmp, MALE=margin.table(tmp, 1)); tabtot
Kappa(tab)
Kappa.test(tab)