>I would like to plot datapoints according to the sex of the person, e.g. >circle for a girl and square for a boy, like this: > > > Have you tried using stripplot (a lattice plot) instead? You can then use the group argument to display different groups differently. eg. stripplot(age ~ family, data=famdata, group=sex) Hadley