В дополнение к ответу @Rui Barradas, вот код, чтобы сделать это с базовой графикой:
df <- read.table(text = "variable lowerci beta upperci
a 0.86471152 0.87615515 0.88759878
b 0.966626566 0.970159119 0.973691671
c 0.94946211 0.955502661 0.961543211
", header = TRUE, as.is = TRUE)
plot(df[, 3], pch = 19, cex = 2, ylim = range(df[, -1]))
for (i in c(1:nrow(df))) {
lines(c(i, i), c(df[i, 2], df[i, 4]))
}