Блог

Итак, я вычисляю статистику между двумя значениями, одним случайным, а другим — «реальными» баллами. Поскольку она является своего рода парной, я строю распределения различий для разных условий и состояний, facet_grid но меня также интересует маргинальное (глобальное распределение, т.Е. margin = T ). Теперь я добавляю уровень достоверности для t-теста (имеет ли это смысл, кстати?), Поскольку они нормально распределены для каждого условия. Тем не менее, я не хочу, чтобы горизонтальные полосы были забиты в маргинальном распределении (all) . Возможно ли это? На самом деле было бы неплохо добавить только одну горизонтальную полосу, чтобы указать доверительный предел для глобального распределения.

Код для создания проблемы приведен ниже, вместе с выходным графиком:

 library(tidyverse)
library(rstatix)
main_tb <- tibble("scores" = rnorm(300, 0.6, 0.1), 
                  "rnd_scores" = rnorm(300, 0.4, 0.1), 
                  "condition" = sample(c("A", "B", "C"), 300, replace = T), 
                  "state" = sample(letters[5:10], 300, replace = T)) %>% 
  mutate(diff = scores - rnd_scores) %>% 
  print(n = 5)

conf_tb <- main_tb %>% 
  select(-diff) %>% 
  pivot_longer(names_to = "type", values_to = "fin_scores", c(-condition, -state)) %>% 
  group_by(condition) %>% 
  rstatix::t_test(fin_scores ~ type, detailed = T) %>% 
  rstatix::add_significance() %>% print(n = 5)

ggplot(main_tb, aes(y = diff, x = state))   
  geom_boxplot()   
  geom_hline(data = conf_tb, aes(yintercept = -conf.low))  
  facet_grid(~condition, margins = T)  
  theme_classic()