I am creating an albatross plot (basically a scatter of p-values vs sample size) of compiled literature. I have two dimensions which I would like to differentiate in the plot. One is whether the paper is rated as "weak" or "strong" and the other is what the broad topic of the paper is, 1 of 5 options. Using the Albatross command with the by() option you can tell STATA to differentiate the points visually on the plot by one variable. It also has a color option, but color just seems to use the variable specified in the by() option.

I would ideally like to specify my rating variable in the by() option, and then have different colors differentiate the topic of the paper. I initially just tried this:

albatross SampleSize p_value direction, type(correlation) by(Rating) color(Topic)

However, the color option does not allow for a varlist or varname to be used with it.

I am hoping to learn either a way to work around this issue or if there is a module I can install that gives more power when creating plots and editing their appearance. Thank you!