df <- structure(list(ID = structure(c(1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L,
                                      9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 17L, 18L, 19L, 1L, 2L,
                                      3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L,
                                      17L, 18L, 19L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L,
                                      12L, 13L, 14L, 15L, 16L, 17L, 18L, 19L, 1L, 2L, 3L, 4L, 5L, 6L,
                                      7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 17L, 18L, 19L,
                                      1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L,
                                      15L, 16L, 17L, 18L, 19L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L,
                                      10L, 11L, 12L, 13L, 14L, 15L, 16L, 17L, 18L, 19L),
                                    .Label = c("1", "2", "3", "4", "5", "6", "7", "8", "9", "10", "11", "12", "13",
                                               "14", "15", "16", "17", "18", "19"),
                                    class = "factor"),
                     cong = structure(c(1L,
                                        1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
                                        1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,
                                        2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L,
                                        3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
                                        1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L,
                                        2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L,
                                        3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L,
                                        3L),
                                      .Label = c("c", "ic", "n"),
                                      class = "factor"),
                     goNogo = structure(c(1L,
                                          1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
                                          1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
                                          1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
                                          1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,
                                          2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,
                                          2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,
                                          2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,
                                          2L),
                                        .Label = c("g", "ng"),
                                        class = "factor"),
                     rt = c(632, 416,
                            533, 419, 431, 387, 429, 394, 358, 405, 373, 389, 361, 477, 454,
                            574, 455, 492, 389, 722, 431, 545, 468, 450, 405, 439, 454, 374,
                            403, 388, 395, 381, 482, 474, 658, 470, 560, 423, 579, 442, 532,
                            411, 481, 399, 458, 439, 369, 395, 386, 374, 360, 474, 465, 615,
                            469, 494, 403, 533, 611, 502, 519, 391, 470, 371, 453, 381, 438,
                            385, 508, 657, 488, 667, 395, 538, 359, 417, 598, 701, 535, 534,
                            413, 523, 421, 620, 396, 446, 389, 511, 788, 554, 793, 410, 581,
                            395, 449, 531, 609, 468, 530, 414, 480, 387, 502, 382, 473, 386,
                            524, 682, 502, 602, 395, 571, 405, 413)),
                class = "data.frame", row.names = c(NA,-114L))


# Analysis ----------------------------------------------------------------

library(afex)
library(BayesFactor)
library(BFEffects)

rmanova <- aov_ez('ID','rt',df,within = c('cong','goNogo'))
rmanova


BF_norandomeffects <- generalTestBF(rt ~ cong*goNogo + ID,df,
                                    whichRandom = c('ID'),
                                    whichModels = 'withmain')
inclusionBF(BF_norandomeffects, match_models = TRUE)


BF_withrandomeffects <- generalTestBF(rt ~ cong*goNogo*ID,df,
                                      whichRandom = c('ID','cong:ID','goNogo:ID','cong:goNogo:ID'),
                                      whichModels = 'withmain')
inclusionBF(BF_withrandomeffects, match_models = TRUE)