####################
file_path <- "https://www.sthda.com/sthda/RDoc/data/housetasks.txt"
housetasks <- read.delim(file_path, row.names = 1)
housetasks

chisq <- chisq.test(housetasks)
chisq

o <- chisq$observed
e <- chisq$expected
o
round(e,2)
cs.cal <- sum((e-o)^2/e)
cs.cal
dim(housetasks)
dim(housetasks)-1
pchisq(cs.cal, 12*3, lower.tail = F)
chisq