Hi there

Does anyone know how to rename the titles of each graph when using the by() option?

For example, using the code below, I would like each graph to be labelled with its group name + mean read + mean write.

Code:
use https://stats.idre.ucla.edu/stat/stata/notes/hsb2, clear
bysort ses female: egen meanread = mean(read)
bysort ses female: egen meanwrite = mean(write)
replace meanread = round(meanread)
replace meanwrite = round(meanwrite)
twoway (scatter read write), by(ses female, cols(2))
So, for instance, the first graph would be entitled as follows:

low, male
avg read: 49
avg write: 47

on three separate lines as above.

This seems like an impossible challenge!

I would be extremely grateful for any advice.

With thanks