I would like to do several scatterplots were I plot the correlation btw mortality and covid-cases and were I would like to identify the treatment and control by different colors, this will be done for each month. I don't know if I have done the best way but by using the following code I get a nice and neat scatterplot but were I have to do it separattly in four figures, what I now want to do is to merge all of them into one big figure but with several smaller plots. Any body knows how I can do this the best way?

twoway (scatter mortality_1000_treated new_cases_1000indiv ) (scatter mortality_1000_nontreated new_cases_1000indiv) if time== tm(2020m3)
twoway (scatter mortality_1000_treated new_cases_1000indiv ) (scatter mortality_1000_nontreated new_cases_1000indiv) if time== tm(2020m4)
twoway (scatter mortality_1000_treated new_cases_1000indiv ) (scatter mortality_1000_nontreated new_cases_1000indiv) if time== tm(2020m5)

/David