首页
学习
活动
专区
圈层
工具
发布
首页
学习
活动
专区
圈层
工具
MCP广场
社区首页 >问答首页 >如何从ggplot2 R中的自举模型中绘制具有估计平均值和置信区间的原始数据点?

如何从ggplot2 R中的自举模型中绘制具有估计平均值和置信区间的原始数据点?
EN

Stack Overflow用户
提问于 2019-07-23 03:31:21
回答 1查看 143关注 0票数 0

我有几个因素的数据。我运行了一个带有均值和置信区间的自举模型。我想绘制两个数据点以及模型结果均值和置信区间。我怎样才能做到这一点呢?

代码语言:javascript
运行
复制
yval <- c(0.10,0.12,0.50,0.13,0.20,0.15,0.68,0.71,0.90,0.10,0.11,0.12,0.14,0.15,0.16,0.20,0.14,0.15,0.66,0.222)
type <- c("A","S","A","S","A","S","A","S","A","S","A","S","A","S","A","S","A","S","A","S")
status <- c("L","H","L","H","L","H","L","H","L","H","L","H","L","H","L","H","L","H","L","H")


data <- data.frame(yval, type, status)

然后我尝试了一个很长很复杂的模型,它给出了这个矩阵

代码语言:javascript
运行
复制
                      est        2.5%     97.5%
    (Intercept)  0.14593333 -0.10347804 0.3817681
    typeA        0.26295217 -0.01156827 0.5503199
    statusH      0.09181279 -0.18906489 0.3748384

我希望将原始数据绘制为模型中具有适当平均值和置信区间的点。请不要建议在ggplot2中使用内置引导置信区间。我想特别使用上面的几个。

代码语言:javascript
运行
复制
ggplot(data, aes(x=status, y= yval, color = status)) + theme_bw()  +  geom_point(position = pd, size = 1.5) + facet_grid(.~type)
EN

回答 1

Stack Overflow用户

回答已采纳

发布于 2019-07-23 04:03:41

你在找geom_errorbar吗?

代码语言:javascript
运行
复制
data2<-data.frame(yval= c(0.10,0.12,0.50,0.13,0.20,0.15,0.68,0.71,0.90,0.10,0.11,0.12,0.14,0.15,0.16,0.20,0.14,0.15,0.66,0.222),
                  typeA =c("A","S","A","S","A","S","A","S","A","S","A","S","A","S","A","S","A","S","A","S"),
                  statusH=c("L","H","L","H","L","H","L","H","L","H","L","H","L","H","L","H","L","H","L","H")

)
data2<-data2%>%
  mutate(.,type=mapvalues(data2$type, from = c("A", "S"), to = c("typeA", "statusH")))

data<-read.table(text="
type est        l2.5     u97.5
  (Intercept)  0.14593333 -0.10347804 0.3817681
typeA        0.26295217 -0.01156827 0.5503199
statusH      0.09181279 -0.18906489 0.3748384
",header=T)
pd=position_dodge(10)
ggplot(data, aes(x=type, y= est, color = type)) + 
  theme_bw()  +
  geom_point(size = 3.5,position=pd) + 
geom_errorbar(aes(ymin=l2.5, ymax=est+u97.5), width=0.1,position = pd)+
    geom_point(data=data2, aes(x=type,y=yval),size=1)

我可能遗漏了一些东西,但听起来你可以直接使用第二个矩阵。

票数 1
EN
页面原文内容由Stack Overflow提供。腾讯云小微IT领域专用引擎提供翻译支持
原文链接:

https://stackoverflow.com/questions/57152612

复制
相关文章

相似问题

领券
问题归档专栏文章快讯文章归档关键词归档开发者手册归档开发者手册 Section 归档