how to create a scatter plot of two categorical variables by assigning one to shape and other to color. in this example
-
Login or Register
- Log in with
clear input str3 location str5 fraction str3 rep float total_af_ppb "BMT" "heavy" "one" .634899 "BMT" "heavy" "one" .9944572 "BMT" "heavy" "one" .7497677 "BMT" "heavy" "one" .1736788 "BMT" "heavy" "one" .6107705 "BMT" "heavy" "two" 2.5754216 "BMT" "heavy" "two" 2.3678162 "BMT" "heavy" "two" 2.3005245 "BMT" "heavy" "two" 2.007538 "BMT" "heavy" "two" 2.670137 "BMT" "light" "one" 1.4241406 "BMT" "light" "one" 1.953762 "BMT" "light" "one" 1.0867478 "BMT" "light" "one" 1.8949648 "BMT" "light" "one" 1.5890286 "BMT" "light" "two" 3.400583 "BMT" "light" "two" 3.66549 "BMT" "light" "two" 3.419839 "BMT" "light" "two" 3.747205 "BMT" "light" "two" 3.7190144 "BSA" "heavy" "one" 104.08141 "BSA" "heavy" "one" 106.1555 "BSA" "heavy" "one" 101.74577 "BSA" "heavy" "one" 103.61765 "BSA" "heavy" "one" 101.339 "BSA" "heavy" "two" 300.01364 "BSA" "heavy" "two" 302.571 "BSA" "heavy" "two" 306.5174 "BSA" "heavy" "two" 309.2521 "BSA" "heavy" "two" 308.2334 "BSA" "light" "one" 209.2294 "BSA" "light" "one" 207.48042 "BSA" "light" "one" 205.21414 "BSA" "light" "one" 204.02216 "BSA" "light" "one" 208.682 "BSA" "light" "two" 402.7266 "BSA" "light" "two" 407.2395 "BSA" "light" "two" 407.9555 "BSA" "light" "two" 408.9251 "BSA" "light" "two" 407.0788 "HBY" "heavy" "one" 53.65269 "HBY" "heavy" "one" 59.3105 "HBY" "heavy" "one" 56.21681 "HBY" "heavy" "one" 58.00435 "HBY" "heavy" "one" 54.79837 "HBY" "heavy" "two" 151.42947 "HBY" "heavy" "two" 158.34344 "HBY" "heavy" "two" 153.43124 "HBY" "heavy" "two" 155.90686 "HBY" "heavy" "two" 153.2953 "HBY" "light" "one" 106.9963 "HBY" "light" "one" 108.14297 "HBY" "light" "one" 103.42973 "HBY" "light" "one" 101.07968 "HBY" "light" "one" 109.67175 "HBY" "light" "two" 201.2855 "HBY" "light" "two" 207.57854 "HBY" "light" "two" 202.5003 "HBY" "light" "two" 209.26913 "HBY" "light" "two" 207.11842 "KKA" "heavy" "one" 203.59343 "KKA" "heavy" "one" 209.01303 "KKA" "heavy" "one" 203.89214 "KKA" "heavy" "one" 214.2715 "KKA" "heavy" "one" 204.90623 "KKA" "heavy" "two" 415.3449 "KKA" "heavy" "two" 407.3071 "KKA" "heavy" "two" 405.4138 "KKA" "heavy" "two" 419.8226 "KKA" "heavy" "two" 413.7026 "KKA" "light" "one" 310.0553 "KKA" "light" "one" 313.80765 "KKA" "light" "one" 317.272 "KKA" "light" "one" 300.80927 "KKA" "light" "one" 303.68445 "KKA" "light" "two" 508.3976 "KKA" "light" "two" 512.951 "KKA" "light" "two" 518.2063 "KKA" "light" "two" 513.61847 "KKA" "light" "two" 517.13763 "KSM" "heavy" "one" 150.96303 "KSM" "heavy" "one" 162.586 "KSM" "heavy" "one" 159.3123 "KSM" "heavy" "one" 156.06264 "KSM" "heavy" "one" 164.67955 "KSM" "heavy" "two" 355.4415 "KSM" "heavy" "two" 354.7074 "KSM" "heavy" "two" 359.5711 "KSM" "heavy" "two" 352.8049 "KSM" "heavy" "two" 357.5802 "KSM" "light" "one" 257.91446 "KSM" "light" "one" 261.78012 "KSM" "light" "one" 257.07602 "KSM" "light" "one" 253.44977 "KSM" "light" "one" 261.96524 "KSM" "light" "two" 452.4741 "KSM" "light" "two" 463.9942 "KSM" "light" "two" 455.999 "KSM" "light" "two" 464.823 "KSM" "light" "two" 463.9318 "SYA" "heavy" "one" 786.4038 "SYA" "heavy" "one" 723.6804 "SYA" "heavy" "one" 753.365 "SYA" "heavy" "one" 798.243 "SYA" "heavy" "one" 721.7005 "SYA" "heavy" "two" 2342.782 "SYA" "heavy" "two" 2377.8247 "SYA" "heavy" "two" 2395.336 "SYA" "heavy" "two" 2405.7007 "SYA" "heavy" "two" 2413.2275 "SYA" "light" "one" 1552.4447 "SYA" "light" "one" 1536.377 "SYA" "light" "one" 1615.82 "SYA" "light" "one" 1580.086 "SYA" "light" "one" 1617.3347 "SYA" "light" "two" 3143.638 "SYA" "light" "two" 3187.648 "SYA" "light" "two" 3134.591 "SYA" "light" "two" 3149.664 "SYA" "light" "two" 3202.1956 end
stripplot total_af_ppb, over(location) separate(fraction) by(rep) jitter(2 2) msymbol (o o) vertical
separate total, by(rep) veryshortlabel encode location, gen(n_location) twoway (scatter total_af_ppb? n_location if fraction == "light", xlabel(1 2 3 4 5 6, val) msymbol(o d) mcolor(red red) jitter(5)) /// (scatter total_af_ppb? n_location if fraction == "heavy", xlabel(1 2 3 4 5 6, val) msymbol(o d) mcolor(blue blue) jitter(5)), /// legend(label(1 "One (Light)") label(2 "Two (Light)") label(3 "One (Heavy)") label(4 "Two (Heavy)"))
clear input str3 location str5 fraction str3 rep float total_af_ppb "BMT" "heavy" "one" .634899 "BMT" "heavy" "one" .9944572 "BMT" "heavy" "one" .7497677 "BMT" "heavy" "one" .1736788 "BMT" "heavy" "one" .6107705 "BMT" "heavy" "two" 2.5754216 "BMT" "heavy" "two" 2.3678162 "BMT" "heavy" "two" 2.3005245 "BMT" "heavy" "two" 2.007538 "BMT" "heavy" "two" 2.670137 "BMT" "light" "one" 1.4241406 "BMT" "light" "one" 1.953762 "BMT" "light" "one" 1.0867478 "BMT" "light" "one" 1.8949648 "BMT" "light" "one" 1.5890286 "BMT" "light" "two" 3.400583 "BMT" "light" "two" 3.66549 "BMT" "light" "two" 3.419839 "BMT" "light" "two" 3.747205 "BMT" "light" "two" 3.7190144 "BSA" "heavy" "one" 104.08141 "BSA" "heavy" "one" 106.1555 "BSA" "heavy" "one" 101.74577 "BSA" "heavy" "one" 103.61765 "BSA" "heavy" "one" 101.339 "BSA" "heavy" "two" 300.01364 "BSA" "heavy" "two" 302.571 "BSA" "heavy" "two" 306.5174 "BSA" "heavy" "two" 309.2521 "BSA" "heavy" "two" 308.2334 "BSA" "light" "one" 209.2294 "BSA" "light" "one" 207.48042 "BSA" "light" "one" 205.21414 "BSA" "light" "one" 204.02216 "BSA" "light" "one" 208.682 "BSA" "light" "two" 402.7266 "BSA" "light" "two" 407.2395 "BSA" "light" "two" 407.9555 "BSA" "light" "two" 408.9251 "BSA" "light" "two" 407.0788 "HBY" "heavy" "one" 53.65269 "HBY" "heavy" "one" 59.3105 "HBY" "heavy" "one" 56.21681 "HBY" "heavy" "one" 58.00435 "HBY" "heavy" "one" 54.79837 "HBY" "heavy" "two" 151.42947 "HBY" "heavy" "two" 158.34344 "HBY" "heavy" "two" 153.43124 "HBY" "heavy" "two" 155.90686 "HBY" "heavy" "two" 153.2953 "HBY" "light" "one" 106.9963 "HBY" "light" "one" 108.14297 "HBY" "light" "one" 103.42973 "HBY" "light" "one" 101.07968 "HBY" "light" "one" 109.67175 "HBY" "light" "two" 201.2855 "HBY" "light" "two" 207.57854 "HBY" "light" "two" 202.5003 "HBY" "light" "two" 209.26913 "HBY" "light" "two" 207.11842 "KKA" "heavy" "one" 203.59343 "KKA" "heavy" "one" 209.01303 "KKA" "heavy" "one" 203.89214 "KKA" "heavy" "one" 214.2715 "KKA" "heavy" "one" 204.90623 "KKA" "heavy" "two" 415.3449 "KKA" "heavy" "two" 407.3071 "KKA" "heavy" "two" 405.4138 "KKA" "heavy" "two" 419.8226 "KKA" "heavy" "two" 413.7026 "KKA" "light" "one" 310.0553 "KKA" "light" "one" 313.80765 "KKA" "light" "one" 317.272 "KKA" "light" "one" 300.80927 "KKA" "light" "one" 303.68445 "KKA" "light" "two" 508.3976 "KKA" "light" "two" 512.951 "KKA" "light" "two" 518.2063 "KKA" "light" "two" 513.61847 "KKA" "light" "two" 517.13763 "KSM" "heavy" "one" 150.96303 "KSM" "heavy" "one" 162.586 "KSM" "heavy" "one" 159.3123 "KSM" "heavy" "one" 156.06264 "KSM" "heavy" "one" 164.67955 "KSM" "heavy" "two" 355.4415 "KSM" "heavy" "two" 354.7074 "KSM" "heavy" "two" 359.5711 "KSM" "heavy" "two" 352.8049 "KSM" "heavy" "two" 357.5802 "KSM" "light" "one" 257.91446 "KSM" "light" "one" 261.78012 "KSM" "light" "one" 257.07602 "KSM" "light" "one" 253.44977 "KSM" "light" "one" 261.96524 "KSM" "light" "two" 452.4741 "KSM" "light" "two" 463.9942 "KSM" "light" "two" 455.999 "KSM" "light" "two" 464.823 "KSM" "light" "two" 463.9318 "SYA" "heavy" "one" 786.4038 "SYA" "heavy" "one" 723.6804 "SYA" "heavy" "one" 753.365 "SYA" "heavy" "one" 798.243 "SYA" "heavy" "one" 721.7005 "SYA" "heavy" "two" 2342.782 "SYA" "heavy" "two" 2377.8247 "SYA" "heavy" "two" 2395.336 "SYA" "heavy" "two" 2405.7007 "SYA" "heavy" "two" 2413.2275 "SYA" "light" "one" 1552.4447 "SYA" "light" "one" 1536.377 "SYA" "light" "one" 1615.82 "SYA" "light" "one" 1580.086 "SYA" "light" "one" 1617.3347 "SYA" "light" "two" 3143.638 "SYA" "light" "two" 3187.648 "SYA" "light" "two" 3134.591 "SYA" "light" "two" 3149.664 "SYA" "light" "two" 3202.1956 end egen group = group(rep fraction), label niceloglabels total_af_ppb , style(13) local(show) egen median = median(total), by(location) egen order = group(median) labmask order, values(location) set scheme s1color stripplot total_af_ppb, over(order) vertical box cumul cumprob centre separate(group) ysc(log) ms(Oh + Oh +) mc(red red blue blue) yla(`show', ang(h)) xla(, noticks) xtitle("") legend(col(1) order(7 6 5 4) pos(3)) refline reflevel(gmean) name(NJC, replace) stripplot total_af_ppb, over(location) separate(fraction) by(rep) jitter(2 2) msymbol (o o) vertical name(IP, replace) separate total, by(rep) veryshortlabel encode location, gen(n_location) twoway (scatter total_af_ppb? n_location if fraction == "light", xlabel(1 2 3 4 5 6, val) msymbol(o d) mcolor(red red) jitter(5)) /// (scatter total_af_ppb? n_location if fraction == "heavy", xlabel(1 2 3 4 5 6, val) msymbol(o d) mcolor(blue blue) jitter(5)), /// legend(label(1 "One (Light)") label(2 "Two (Light)") label(3 "One (Heavy)") label(4 "Two (Heavy)")) name(WW, replace)
Comment