Another interesting feature of plot is that the content of the plot buffer is available in SIC variables: X_DATA, Y_DATA and W_DATA containing respectively the X values, the Y values and the weights associated to each data point. These values can be used in procedures to do your own computations.
Two options further enrich the possibilities:
Another useful option is /NOFRAME, that allow to not reset the plot before doing it, which in effect allows to overplot some another set of data on top of a previous plot (with another color if wanted).