Thank you all for your answers. I eventually use the points command plot(MATRIX[,1]) points(MATRIX[,2]) points(MATRIX[,3]) ... with matplot I had numbers instead of points with type="p" and par(new=TRUE) makes complications with the axis.