我有一组针对不同人群/地区群体的标准误差的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/,除了错误条。)
答案 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)