R语言编程啊

如果你可以把standard error算出来，就可以储存standard error。命令为

gen stderr=_se[bx]

其中bx是你的自变量名称，然后构造上限和下限，结合在一起绘制。

可以试试stata的twoway plot,用lpolyci这个选项:       
tw ( lpolyci coefficient        quantile,  ciplot(rline) pstyle(ci2))        (lpolyci coefficient quantile,  lp(shortdash))        ///
                 , legend(label(1 "95% CI") label(2 "2010") label(3 "95% CI") label(4 "2015"))                ///
                                        title("Coefficient estimate with 95% confidence interval", size(medium))