lt;!-- %%      main="Flow Duration Curve", xlab="\% Time flow equalled or exceeded",  --&gt;
&lt;!-- %%      ylab="Q, [m3/s]", ylim, yat=c(0.01, 0.1, 1), xat=c(0.01, 0.025, 0.05),  --&gt;
&lt;!-- %%      col=matrix(c(rep("black", NCOL(x)),  --&gt;
&lt;!-- %%      palette("default")[2:(NCOL(x)+1)]), byrow=FALSE, ncol=2), --&gt;
&lt;!-- %%      pch=matrix(rep(c(1, 15), NCOL(x)), byrow=TRUE, ncol=2), --&gt;
&lt;!-- %%      lwd=matrix(rep(c(1, 0.8), NCOL(x)), byrow=TRUE, ncol=2), --&gt;
&lt;!-- %%      lty=matr