2016-01-13 22 views
-1

我有一组f4比率统计数据,针对不同人口/地区群体的标准误差。针对不同群体的误差线的绘图点估计

我想制作一个情节,其中X轴是不同的组,并在X轴上我有95%可信区间的误差线点估计。基本上,我想以图形方式展示我的兴趣群体的价值如何相对于彼此和其他群体出现。

我的数据基本是这样的:

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 

阿尔法将是点估计和SE将是标准的错误。 (Z和样本大小在这里不相关)。

任何人都可以提出一个这样做的好方法吗?谢谢!

(我所寻找的是一样的东西在这里图2b http://www.ncbi.nlm.nih.gov/pmc/articles/PMC3493647/figure/f2/除与误差条)。

+0

它更可能的是,我们将能够帮助你,如果你犯了一个最小的重复的例子,你的问题一起去。我们可以从中学习并使用它来向您展示如何解决您的问题。你可以看看[这个SO帖子](https://stackoverflow.com/questions/5963269/how-to-make-a-great-r-reproducible-example)关于如何在R中创建一个很好的重现例子此外,如果你概述你已经尝试过的东西,这也会很有帮助。 –

+0

我还没有尝试过任何东西,因为我不知道如何去做我需要做的事情(这就是我寻求帮助的原因)。 – MolecularAnthropologist

+0

您应该发布代码以进行数据输入。如果这是你不知道如何完成的部分,那么你应该先做一些关于如何进行基本R活动的搜索。这不是代码写作中心。研究“介绍R”作为第一步似乎是必要的。正如学习如何在R上搜索有很多答案的问题的答案一样。 –

回答

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

+0

谢谢!我有一些问题得到这个与我的数据表,其中有33个人口 我修改了代码来让它读取,但似乎是错误的。 'df < - structure(list(Group = structure(1:33,.Label = tbl $ Populations,class =“factor”),alpha = tbl $ alpha,SE = tbl $ SE,Z = tbl $ Z, Sample.size = tbl $ sample,.Names = c(“Group”,“alpha”,“SE”,“Z”,“Sample.size”),class =“data.frame”,row.names = c NA,-33L)))' '错误:ggplot2不知道如何处理班级列表中的数据 另外:警告消息: 在str.default(obj,...):'object'没有有效的级别()' – MolecularAnthropologist

+2

这就是为什么我们倾向于要求你提供[一个最小可重现的例子](http://stackoverflow.com/help/mcve)(是的,看看该页面)。更多的时候,人们在经历这个练习时回答自己的问题。两个指针,你很可能想重新排列你的数据从一个'list'结构到一个数据框,其次不要在你的'ggplot2'调用中使用'$'。 –

相关问题