CIbars.R

cibars <- function(x, y, y.se, tinyfactor = 50, ...){
  
  lower <- y - 1.96*y.se
  upper <- y + 1.96*y.se
  
  segments(x0 = x, y0 = lower, x1 = x, y1 = upper, ...)
  tinybarwidth <- (max(x) - min(x))/tinyfactor
  
  segments(x0 = x - tinybarwidth, y0 = lower, x1 = x + tinybarwidth, y1 = lower, ...)
  segments(x0 = x - tinybarwidth, y0 = upper, x1 = x + tinybarwidth, y1 = upper, ...)
  }