library(plot2)
par(pch = 19)

mod = lm(mpg ~ hp + factor(cyl), mtcars)
coefs = data.frame(names(coef(mod)), coef(mod), confint(mod))
coefs = setNames(coefs, c("x", "y", "ymin", "ymax"))

with(
  coefs,
  plot2(
    x = x,
    y = y,
    ymin = ymin,
    ymax = ymax,
    type = "pointrange"
  )
)

with(
  coefs,
  plot2(
    x = x,
    y = y,
    ymin = ymin,
    ymax = ymax,
    type = "errorbar"
  )
)