## ----include = FALSE---------------------------------------------------------- knitr::opts_chunk$set( collapse = TRUE, comment = "#>" ) ## ----setup-------------------------------------------------------------------- library(DFBA) ## ----eval = FALSE------------------------------------------------------------- # G1 <- c(96.49, 96.78, 97.26, 98.85, 99.75, 100.14, 101.15, 101.39, 102.58, # 107.22, 107.70, 113.26) # # G2 <- c(101.16, 102.09, 103.14, 104.70, 105.27, 108.22, 108.32, 108.51, 109.88, # 110.32, 110.55, 113.42) # # set.seed(1) # # dfba_mann_whitney(E = G1, # C = G2, # hide_progress = TRUE) ## ----echo = FALSE------------------------------------------------------------- G1 <- c(96.49, 96.78, 97.26, 98.85, 99.75, 100.14, 101.15, 101.39, 102.58, 107.22, 107.70, 113.26) G2 <- c(101.16, 102.09, 103.14, 104.70, 105.27, 108.22, 108.32, 108.51, 109.88, 110.32, 110.55, 113.42) load("mw_ex1") #mw_ex1 ## ----eval = FALSE------------------------------------------------------------- # set.seed(2) # dfba_mann_whitney(E = G1, # C = G2, # hide_progress = TRUE) # # mw_ex2 <- dfba_mann_whitney(E = G1, # C = G2, # hide_progress = TRUE) # # ## ----echo = FALSE------------------------------------------------------------- load("mw_ex2") mw_ex2 ## ----eval = FALSE------------------------------------------------------------- # plot(dfba_mann_whitney(E = G1, # C = G2, # hide_progress = TRUE)) ## ----echo = FALSE, fig.width = 7---------------------------------------------- plot(mw_ex2) ## ----------------------------------------------------------------------------- set.seed(1) dfba_mann_whitney(E = G1, C = G2, method = "large") ## ----------------------------------------------------------------------------- A1B1 <- c(11.541, 11.854, 11.313, 14.201, 11.333, 11.583, 11.223) A2B1 <- c(11.210, 11.117, 12.967, 12.514, 11.232, 13.585, 11.023) A3B1 <- c(4.762, 2.323, 5.890, 2.722, 2.499, 2.534, 2.016) A1B2 <- c(1.500, 1.562, 1.444, 1.822, 1.802, 1.075, 1.464) A2B2 <- c(2.663, 1.503, 1.086, 1.459, 1.296, 1.009, 3.316) A3B2 <- c(11.067, 11.117, 10.180, 10.060, 10.664, 10.074, 10.355) ## ----echo=FALSE, fig.width=7-------------------------------------------------- yB1 <- c(mean(A1B1), mean(A2B1), mean(A3B1)) yB2 <- c(mean(A1B2), mean(A2B2), mean(A3B2)) xA <- c(1, 2, 3) plot(xA, yB1, type = "o", xlab = "A levels", ylab = "group means", main = "B1 black; B2 red", ylim = c(0, 12)) lines(xA, yB2, type = "o", col = "firebrick") ## ----------------------------------------------------------------------------- dfba_mann_whitney(E = c(A1B1, A2B1, A3B1), C = c(A1B2, A2B2, A3B2)) ## ----------------------------------------------------------------------------- dfba_mann_whitney(E = c(A1B1, A2B1, A3B2), C = c(A1B2, A2B2, A3B1))