# NOT RUN {
a = rnorm(100)
b = 0.5*a + rnorm(100)
c = -0.5*b + rnorm(100)
est = summary(lm(a ~ c + b))
errbar(est$coefficients, x.shift = -.2)
errbar(est$coefficients, x.shift = .2, add = TRUE, col = 2, bar.lty = 2, pch=15)
# }
Run the code above in your browser using DataLab