library("plm")
> data("Grunfeld",package="plm")
> grun.fe<-plm(inv~value+capital,data=Grunfeld,model="within"
+ )
> summary(grun.fe)
Oneway (individual) effect Within Model
Call:
plm(formula = inv ~ value + capital, data = Grunfeld, model = "within")
Balanced Panel: n=10, T=20, N=200
Residuals :
    Min.  1st Qu.   Median  3rd Qu.     Max.
-184.000  -17.600    0.563   19.200  251.000
Coefficients :
        Estimate Std. Error t-value  Pr(>|t|)   
value   0.110124   0.011857  9.2879 < 2.2e-16 ***
capital 0.310065   0.017355 17.8666 < 2.2e-16 ***

1、应该怎么用公式计算出value和capital对应的置信区间？

2、假如我改变数据运行这样的模型5次，设为m=c(1,2,3,4,5),那么得到对应value和capital值各有5个，如设得到value值5个，为value=c(1.9,1.8,1.5,1,2)
且对应的置信区间也各有5个。
我现在想画包含置信区间的图形（以m为x轴，value为y轴），这样图形应该怎么实现？类似如下的图形。

简单的说，就是已经有了m和value，及value对应的置信区间，应该如何画如下的图形。

谢谢各位帮忙