O código abaixo usa dados novos, já que os dados originais não foram fornecidos na pergunta.
time <- 3:10
tq1_mean <- 3:10
sd <- sd(tq1_mean)
plot (time, tq1_mean, type ="p", col = "red", ylim=c(0, 15))
# desenha o traço vertical
segments(time,tq1_mean - sd, time, tq1_mean + sd)
# Define o tamanho dos traços horizontais
tamanho <- 0.1
# desenha o traço horizontal superior
segments(time - tamanho, tq1_mean - sd, time + tamanho, tq1_mean - sd)
# desenha o traço horizontal inferior
segments(time - tamanho, tq1_mean + sd, time + tamanho, tq1_mean + sd)