Ooops, I didn't convert this one to text right for the list.
b1<- substitute( mu - d*sigma, list(d=*-round(dividers[1],2))* ) should be b1<- substitute( mu - d*sigma, list(d=-round(dividers[1],2)) ) and similarly for labels=*as.expression(c(b1,b2,b3,b4,b5))*, padj=-1) read labels=as.expression(c(b1,b2,b3,b4,b5)), padj=-1) Apologies Allan ______________________________________________ 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.