我有一个制作系数图的问题。如何制作系数图,哪些系数来自不同类型的模型?例如
> med.out1$d.avg
[1] 0.03856374
> med.out1$d.avg.ci
2.5% 97.5%
0.02406029 0.05213770
> med.out2$d.avg
[1] 0.03703407
> med.out2$d.avg.ci
2.5% 97.5%
0.02320818 0.05294310
> med.out3$d.avg
[1] 0.03631961
> med.out3$d.avg.ci
2.5% 97.5%
0.02158095 0.05110621
通过这些估算,如何将所有这三个系数组合在一起?基本绘图程序,ggplot或任何推荐的系数绘图程序都对我有好处。
谢谢。
答案 0 :(得分:1)
有一个更易于重现的问题会很好,但有点像:
mList <- list(med.out1,med.out2,med.out3)
est <- sapply(mList,"[[","d.avg")
lwr <- sapply(mList,function(x) x$d.avg.ci[1])
upr <- sapply(mList,function(x) x$d.avg.ci[2])
library(plotrix)
plotCI(1:3,est,ui=upr,li=lwr,axes=FALSE)
axis(side=2)
axis(side=1,at=1:3,labels=c("med.out1","med.out2","med.out3"))