# 練習問題
# 1.
Ctable <- xtabs(Freq ~ management + own, data = JobSatisfaction)
# (a)
chisq.test(Ctable)
# (b)
fisher.test(Ctable)

# 2.
Ctable <- xtabs(Freq ~ Sex + Survived, data = Titanic[,,"Child",])
chisq.test(Ctable)

# 3.
Ctable <- xtabs(Freq ~ Class + Survived, data = Titanic)[-4,]
result <- chisq.test(Ctable)
v <- result$expected * (1 - rowSums(Ctable)/sum(Ctable)) %o%
  (1 - colSums(Ctable)/sum(Ctable))
(Ctable - result$expected) / sqrt(v)

# 4.
lbl_test(Hospital)