http://www.commres.net/ 2026-05-01T07:33:21+00:00 http://www.commres.net/ http://www.commres.net/_media/wiki/logo.png text/html 2025-12-02T22:42:26+00:00 Anonymous (anonymous@undisclosed.example.com) http://www.commres.net/r/chi-square_test/code01?rev=1764715346&do=diff # install.packages("MASS") library(MASS) print(str(survey)) ??MASS::survevy wh.clap = data.frame(survey$Clap,survey$W.Hnd) wh.clap.tbl = table(survey$Clap,survey$W.Hnd) wh.clap.tbl chisq <- chisq.test(wh.clap.tbl) chisq o <- chisq$observed e <- chisq$expected o round(e, 2) res <- (e-o)/sqrt(e) round(res,2) chisq.cal <- sum(res^2) p.val <- pchisq(chisq.cal, df=2, lower.tail = F) chisq.cal p.val chisq barplot(wh.clap.tbl, beside = TRUE, col = c("red", "lightgreen", "blue"), … text/html 2025-12-02T22:48:10+00:00 Anonymous (anonymous@undisclosed.example.com) http://www.commres.net/r/chi-square_test/code02?rev=1764715690&do=diff #################### 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 text/html 2025-12-02T22:49:35+00:00 Anonymous (anonymous@undisclosed.example.com) http://www.commres.net/r/chi-square_test/code03?rev=1764715775&do=diff # or the below works also chisq <- chisq.test(df$treatment, df$improvement, correct=F) chisq o <- chisq$observed e <- chisq$expected o round(e,2) r <- (e-o)/sqrt(e) # this is called residuals round(r, 3) round(chisq$residuals, 3) cs.tmp <- sum(chisq$residuals^2) cs.tmp cs.cal <- sum((e-o)^2/e) pchisq(cs.cal, 1,lower.tail = F) cs.cal chisq library(corrplot) corrplot(chisq$residuals, is.cor = FALSE) contrib <- 100*chisq$residuals^2/chisq$statistic round(contrib, 3) # Visualize the contributio… text/html 2025-12-02T22:41:55+00:00 Anonymous (anonymous@undisclosed.example.com) http://www.commres.net/r/chi-square_test/output01?rev=1764715315&do=diff > # install.packages("MASS") > library(MASS) > > print(str(survey)) 'data.frame': 237 obs. of 12 variables: $ Sex : Factor w/ 2 levels "Female","Male": 1 2 2 2 2 1 2 1 2 2 ... $ Wr.Hnd: num 18.5 19.5 18 18.8 20 18 17.7 17 20 18.5 ... $ NW.Hnd: num 18 20.5 13.3 18.9 20 17.7 17.7 17.3 19.5 18.5 ... $ W.Hnd : Factor w/ 2 levels "Left","Right": 2 1 2 2 2 2 2 2 2 2 ... $ Fold : Factor w/ 3 levels "L on R","Neither",..: 3 3 1 3 2 1 1 3 3 3 ... $ Pulse : int 92 104 87 NA 35 64 83 74 … text/html 2025-12-02T22:47:56+00:00 Anonymous (anonymous@undisclosed.example.com) http://www.commres.net/r/chi-square_test/output02?rev=1764715676&do=diff > #################### > 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 … text/html 2025-12-05T20:23:33+00:00 Anonymous (anonymous@undisclosed.example.com) http://www.commres.net/r/chi-square_test/output03?rev=1764966213&do=diff > > #Author DataFlair > df <- read.csv("https://goo.gl/j6lRXD") #Reading CSV > head(df) id treatment improvement 1 1 treated improved 2 2 treated improved 3 3 not-treated improved 4 4 treated improved 5 5 treated not-improved 6 6 treated not-improved > tail(df) id treatment improvement 100 100 not-treated improved 101 101 treated improved 102 102 treated improved 103 103 not-treated not-improved 104 104 treated …