Hi,
I want to create two different plot color for "inc" and "mort" by "cancer_type". Using the codes below, I only get one (red) color. How do I get this done?


ID outcome lnor_m lnlci_m lnuci_m cancer_type _ES _seES _LCI _UCI
1 inc .2390169 .0392207 .4317824 Lung .2390169 .10014514 .0392207 .4317824
2 mort .3784365 .2311117 .5247285 Lung .3784365 .07490363 .2311117 .5247285
3 inc .7275486 .5988365 .8586616 Lung .7275486 .06628312 .5988365 .8586616
4 mort .3001046 .10436 .5007753 Lung .3001046 .1011282 .10436 .5007753
5 inc .1739534 .0099503 .3364722 Lung .1739534 .08329793 .0099503 .3364722
6 mort .285179 .1397619 .4317824 Lung .285179 .07449639 .1397619 .4317824
7 inc .5187938 .1988509 .8458682 Bladder .5187938 .16505848 .1988509 .8458682
8 mort .1310282 -.1508229 .4187103 Bladder .1310282 .14529175 -.1508229 .4187103
9 inc .4317824 .0953102 .6830968 Bladder .4317824 .14994832 .0953102 .6830968
10 mort .2926696 -.0618754 .6471033 Bladder .2926696 .18086523 -.0618754 .6471033
11 inc .4187103 .14842 .6881347 Bladder .4187103 .13768485 .14842 .6881347
12 mort .6523252 .4054651 .8960881 Bladder .6523252 .12516122 .4054651 .8960881
13 inc -.356675 -.9416085 .2070142 Lip-mouth-pharynx -.356675 .29302139 -.9416085 .2070142
14 mort .14842 -.2357223 .5187938 Lip-mouth-pharynx .14842 .19248212 -.2357223 .5187938
15 inc .285179 -.0512933 .6151857 Lip-mouth-pharynx .285179 .17002327 -.0512933 .6151857
16 mort .2390169 .10436 .3784365 Lip-mouth-pharynx .2390169 .06991875 .10436 .3784365
17 inc or no .1222176 .0099503 .2311117 Lip-mouth-pharynx .1222176 .05641976 .0099503 .2311117
18 mort .1397619 .0392207 .2390169 Lip-mouth-pharynx .1397619 .05096936 .0392207 .2390169

metan lnor_m lnlci_m lnuci_m if men!=., ///
by(cancer) nosubgroup noov nohet nostat null(1) eform nobox double graphregion(fcolor(white)) ///
label(namevar=outcome) t1(inc-mort) ///
boxopt(msymbol(o) mcolor(red if _education=="inc") ///
xlabel(0.5 1 2 3 4) xsize(4) ysize(20)

Thanks a lot!

Deependra