Díky za sdílení úlohy. Ten nápad by šel vektorizovat, aby to bylo rychlejší. V tidyverse:
# rozebereme na 2 tabulky, prejmenujeme sloupecky, aby se to dalo zase slucovat
da <- d %>% filter(group == "A") %>% rename_with(~ paste0(.x, "_a"))
db <- d %>% filter(group == "B") %>% rename_with(~ paste0(.x, "_b"))
# dame obe tabulky vedle sebe do jedne, tu druhou (B) zamichame
dmix <- bind_cols(da, db %>% slice_sample(prop = 1))
# pocitame shrnuti, napriklad toto
dmix %>% summarise(b_larger_a_true_pc = mean(mq_true_b > mq_true_a))
# rozebereme na 2 tabulky, prejmenujeme sloupecky, aby se to dalo zase slucovat
da <- d %>% filter(group == "A") %>% rename_with(~ paste0(.x, "_a"))
db <- d %>% filter(group == "B") %>% rename_with(~ paste0(.x, "_b"))
# dame obe tabulky vedle sebe do jedne, tu druhou (B) zamichame
dmix <- bind_cols(da, db %>% slice_sample(prop = 1))
# pocitame shrnuti, napriklad toto
dmix %>% summarise(b_larger_a_true_pc = mean(mq_true_b > mq_true_a))