从不同模型制作系数图

时间:2014-02-13 17:22:12

标签: r plot

我有一个制作系数图的问题。如何制作系数图,哪些系数来自不同类型的模型?例如

 > 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或任何推荐的系数绘图程序都对我有好处。

谢谢。

1 个答案:

答案 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"))