Hi Giovanni, The basic idea is:
classiclimits <- aes(x=x[1:100],ymax = classiccis[1:100,e,p, ymin=classiccis[1:100,e,p,2], colour = "classic") ownlimits <- aes(x=x[1:100]+0.4,ymax = owncis[1:100,e,p,1], ymin=owncis[1:100,e,p,2], colour = "own") rbootlimits <- aes(x=x[1:100]+0.8,ymax = rbootcis[1:100,e,p,1], ymin=rbootcis[1:100,e,p,2], colour = "rboot") And then remove the explicit colours elsewhere in your code. Hadley -- Assistant Professor / Dobelman Family Junior Chair Department of Statistics / Rice University http://had.co.nz/ ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.