使用不同组的误差条绘制点估计值

时间:2016-01-13 22:29:23

标签: r

我有一组针对不同人群/地区群体的标准误差的f4比率统计数据。

我想制作一个图,其中X轴是不同的组,在X轴上,我有点估计值,误差条为95%CI。基本上,我想以图形方式显示我感兴趣的人群的价值相对于彼此和其他群体的价值。

我的数据基本上是这样的:

Group   alpha   SE  Z   Sample size
Pop1    0.029000    0.003589    8.116   9
Pop2    0.031868    0.003498    8.231   9
Pop3    0.028969    0.003765    7.942   8
Pop4    0.030651    0.003479    8.792   10

alpha是点估计,SE是标准误差。 (Z和样本大小在这里无关。)

有人能建议这样做的好方法吗?谢谢!

(我正在寻找的是类似于图2b的http://www.ncbi.nlm.nih.gov/pmc/articles/PMC3493647/figure/f2/,除了错误条。)

1 个答案:

答案 0 :(得分:1)

这有用吗?

# > dput(df)
df <- structure(list(Group = structure(1:4, .Label = c("Pop1", "Pop2", 
"Pop3", "Pop4"), class = "factor"), alpha = c(0.029, 0.031868, 
0.028969, 0.030651), SE = c(0.003589, 0.003498, 0.003765, 0.003479
), Z = c(8.116, 8.231, 7.942, 8.792), Sample.size = c(9L, 9L, 
8L, 10L)), .Names = c("Group", "alpha", "SE", "Z", "Sample.size"
), class = "data.frame", row.names = c(NA, -4L))
# > df
#   Group    alpha       SE     Z Sample.size
# 1  Pop1 0.029000 0.003589 8.116           9
# 2  Pop2 0.031868 0.003498 8.231           9
# 3  Pop3 0.028969 0.003765 7.942           8
# 4  Pop4 0.030651 0.003479 8.792          10

# install.packages("ggplot2", dependencies = TRUE)
require(ggplot2)

ggplot(df, aes(y = Group, x = alpha, xmin = alpha - SE,  
               xmax = alpha + SE, label = Group, 
               colour = as.factor(Sample.size))) + 
      geom_point(colour = "black") + geom_text(hjust = 1.2) + theme_classic() +
      theme(axis.title = element_blank(), axis.ticks = element_blank(), 
      axis.text.y = element_blank(), legend.position = "none") # +
      geom_errorbarh(height = .1)

sigh