> #################### > file_path <- "https://www.sthda.com/sthda/RDoc/data/housetasks.txt" > housetasks <- read.delim(file_path, row.names = 1) > housetasks Wife Alternating Husband Jointly Laundry 156 14 2 4 Main_meal 124 20 5 4 Dinner 77 11 7 13 Breakfeast 82 36 15 7 Tidying 53 11 1 57 Dishes 32 24 4 53 Shopping 33 23 9 55 Official 12 46 23 15 Driving 10 51 75 3 Finances 13 13 21 66 Insurance 8 1 53 77 Repairs 0 3 160 2 Holidays 0 1 6 153 > > chisq <- chisq.test(housetasks) > chisq Pearson's Chi-squared test data: housetasks X-squared = 1944.5, df = 36, p-value < 2.2e-16 > > o <- chisq$observed > e <- chisq$expected > o Wife Alternating Husband Jointly Laundry 156 14 2 4 Main_meal 124 20 5 4 Dinner 77 11 7 13 Breakfeast 82 36 15 7 Tidying 53 11 1 57 Dishes 32 24 4 53 Shopping 33 23 9 55 Official 12 46 23 15 Driving 10 51 75 3 Finances 13 13 21 66 Insurance 8 1 53 77 Repairs 0 3 160 2 Holidays 0 1 6 153 > round(e,2) Wife Alternating Husband Jointly Laundry 60.55 25.63 38.45 51.37 Main_meal 52.64 22.28 33.42 44.65 Dinner 37.16 15.73 23.59 31.52 Breakfeast 48.17 20.39 30.58 40.86 Tidying 41.97 17.77 26.65 35.61 Dishes 38.88 16.46 24.69 32.98 Shopping 41.28 17.48 26.22 35.02 Official 33.03 13.98 20.97 28.02 Driving 47.82 20.24 30.37 40.57 Finances 38.88 16.46 24.69 32.98 Insurance 47.82 20.24 30.37 40.57 Repairs 56.77 24.03 36.05 48.16 Holidays 55.05 23.30 34.95 46.70 > cs.cal <- sum((e-o)^2/e) > cs.cal [1] 1944.456 > dim(housetasks) [1] 13 4 > dim(housetasks)-1 [1] 12 3 > pchisq(cs.cal, 12*3, lower.tail = F) [1] 0 > chisq Pearson's Chi-squared test data: housetasks X-squared = 1944.5, df = 36, p-value < 2.2e-16