2

我心里有一个情节想要创作，但不知道如何成功实现这个目标。

我有 2 个数据框，一个包含每个因子级别的平均值，另一个包含这些级别之间的成对差异。

contrasts <- data.frame(Level1 = c("setosa", "setosa", "versicolor"),
                        Level2 = c("versicolor", "virginica", "virginica"),
                        Diff = c(0.65, 0.46, -0.20),
                        CI_low = c(0.53, 0.35, -0.32),
                        CI_high = c(0.75, 0.56, -0.09))

means <- data.frame(Species = c("setosa", "versicolor", "virginica"),
                    Mean = c(3.42, .77, 2.97))

我的目标是使用该方法作为三角形的起点，该三角形将“投影”到相应对比度的水平上，其高度将等于 CI (CI_low和CI_high)。所以它看起来像这样（请原谅我的油漆）：

在此处输入图像描述

使用以下内容，我轻松地添加了初始点：

library(tidyverse)

means %>%
  ggplot() + 
  geom_point(aes(x = Species, y= Mean)) + 
  geom_ribbon(data=contrasts, aes(x=Level1, ymin=CI_low, ymax=CI_high))

但是我在添加三角形时遇到了麻烦。有任何想法吗？非常感谢！

编辑

感谢 Yuriy Barvinchenko，他提供了获取此代码的代码：

contrasts %>% 
  bind_cols(id=1:3) %>% 
  inner_join(means, by=c('Level1' = 'Species')) %>% 
  select(id, x=Level1, y=Mean) %>% 
  bind_rows( (contrasts %>% 
                bind_cols(id=1:3) %>% 
                select(id, x=Level2, y=CI_low)),
             (contrasts %>% 
                bind_cols(id=1:3) %>% 
                select(id, x=Level2, y=CI_high))) %>% 
  ggplot(aes(x = x, y= y, group=id)) + 
  geom_polygon()

然而，基于手段，我预计中间层（杂色）是“最低的”，而在那个图中，它是弗吉尼亚州，它是最低值。

4

1 回答 1

3

如果我正确理解你的问题，你需要这样的代码：

contrasts <- tibble(Level1 = c("setosa", "setosa", "versicolor"),
                        Level2 = c("versicolor", "virginica", "virginica"),
                        Diff = c(0.65, 0.46, -0.20),
                        CI_low = c(0.53, 0.35, -0.32),
                        CI_high = c(0.75, 0.56, -0.09))

means <- tibble(Species = c("setosa", "versicolor", "virginica"),
                                            Mean = c(3.42, .77, 2.97))

library(tidyverse)

contrasts %>% 
  bind_cols(id=1:3) %>% 
  inner_join(means, by=c('Level1' = 'Species')) %>% 
  select(id, x=Level1, y=Mean) %>% 
  bind_rows( (contrasts %>% 
                bind_cols(id=1:3) %>% 
                select(id, x=Level2, y=CI_low)),
             (contrasts %>% 
                bind_cols(id=1:3) %>% 
                select(id, x=Level2, y=CI_high))) %>% 
  ggplot(aes(x = x, y= y, group=id)) + 
  geom_polygon()

请注意，我使用tibble()而不是data.frame()为了避免因素，以便更容易加入这些表。

于 2019-04-25T07:37:01.187 回答