Hello, I have another question. I would like to plot something within a self written function (plotdensity) called by tapply t <- tapply(mat, classes, plotdensity) Now I would like to add each plot the class/level as title. How can I do this? Antje