1

所以，我正在计算两个值之间的统计数据，一个是随机值，另一个是“真实”分数。因为它是成对的，所以我绘制了不同条件和状态的差异分布，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()

阴谋

4

0 回答 0