## Options

options(scipen=999)


# Chi Square Goodness of Fit Test

attach(Table15.1)

table(sex)

ftable(sex)

prop.table(table(sex))


chisq.test(table(sex), p = c(0.30,0.70))

chisq.test(table(sex), p = c(0.30,0.70), simulate.p.value=TRUE)




# Chi Square for Crosstabulationl

table(sex,ishopping)

chisq.test(table(sex,ishopping))

chisq.test(table(sex,ishopping), correct=TRUE)

chisq.test(table(sex,ishopping), simulate.p.value=TRUE)

# Package vcd

library(vcd)

fisher.test(table(sex,ishopping))

assocstats(table(sex,ishopping))




# Plots

barplot(table(sex,ishopping), names.arg=c("Male", "Female"), beside=TRUE)

mosaicplot(table(sex,ishopping),color=c("red", "blue"))