# Examine correlation and regression coefficients
 
pwcorr x y z // x&y: r=.70; x&z:  r=.61
reg y x // b=.70
reg z x // b=4.04
 
# Make graph
scatter y x, msize(large) mcolor(stgreen%80) mlcolor(lime) mlwidth(medium) jitter(10) ///
xlab( , glpattern(solid) glcolor(gs14) glwidth(thin)) ///
ylab( , glpattern(solid) glcolor(gs14) glwidth(thin)) /// 
|| lfit y x, lwidth(medthick) lcolor(red%80) || ///
scatter z x, msize(large) mcolor(stblue%80) mlcolor(cyan) mlwidth(medium) jitter(10) || ///
lfit z x, lcolor(stc10%90) lwidth(medthick) ///
legend(order (1 "Y Variable" 2 "{&beta}=.70" 3 "Z Variable" 4 "{&beta}=4.04")) ///
legend(region(fcolor(gs15) lcolor(gs10))) ///
xsize(6.5) ysize(4.5) graphregion(margin(vsmall)) ///
title(Correlation vs. Regression, box fcolor(black) color(white) span bexpand) ///
scheme(white_jet) /// note: requires installing "schemepack"
text(75 48.5 "{stSans: {it:r}} = .70", placement(c) justification(left) size(medium) box fcolor(edkblue) color(white)) ///
text(50 48.5 "{stSans: {it:r}} = .61", placement(c) justification(left) size(medium) box fcolor(edkblue) color(white)) ///
xtitle(X Variable)

IyBFeGFtaW5lIGNvcnJlbGF0aW9uIGFuZCByZWdyZXNzaW9uIGNvZWZmaWNpZW50cwoKcHdjb3JyIHggeSB6IC8vIHgmeTogcj0uNzA7IHgmejogIHI9LjYxCnJlZyB5IHggLy8gYj0uNzAKcmVnIHogeCAvLyBiPTQuMDQKCiMgTWFrZSBncmFwaApzY2F0dGVyIHkgeCwgbXNpemUobGFyZ2UpIG1jb2xvcihzdGdyZWVuJTgwKSBtbGNvbG9yKGxpbWUpIG1sd2lkdGgobWVkaXVtKSBqaXR0ZXIoMTApIC8vLwp4bGFiKCAsIGdscGF0dGVybihzb2xpZCkgZ2xjb2xvcihnczE0KSBnbHdpZHRoKHRoaW4pKSAvLy8KeWxhYiggLCBnbHBhdHRlcm4oc29saWQpIGdsY29sb3IoZ3MxNCkgZ2x3aWR0aCh0aGluKSkgLy8vIAp8fCBsZml0IHkgeCwgbHdpZHRoKG1lZHRoaWNrKSBsY29sb3IocmVkJTgwKSB8fCAvLy8Kc2NhdHRlciB6IHgsIG1zaXplKGxhcmdlKSBtY29sb3Ioc3RibHVlJTgwKSBtbGNvbG9yKGN5YW4pIG1sd2lkdGgobWVkaXVtKSBqaXR0ZXIoMTApIHx8IC8vLwpsZml0IHogeCwgbGNvbG9yKHN0YzEwJTkwKSBsd2lkdGgobWVkdGhpY2spIC8vLwpsZWdlbmQob3JkZXIgKDEgIlkgVmFyaWFibGUiIDIgInsmYmV0YX09LjcwIiAzICJaIFZhcmlhYmxlIiA0ICJ7JmJldGF9PTQuMDQiKSkgLy8vCmxlZ2VuZChyZWdpb24oZmNvbG9yKGdzMTUpIGxjb2xvcihnczEwKSkpIC8vLwp4c2l6ZSg2LjUpIHlzaXplKDQuNSkgZ3JhcGhyZWdpb24obWFyZ2luKHZzbWFsbCkpIC8vLwp0aXRsZShDb3JyZWxhdGlvbiB2cy4gUmVncmVzc2lvbiwgYm94IGZjb2xvcihibGFjaykgY29sb3Iod2hpdGUpIHNwYW4gYmV4cGFuZCkgLy8vCnNjaGVtZSh3aGl0ZV9qZXQpIC8vLyBub3RlOiByZXF1aXJlcyBpbnN0YWxsaW5nICJzY2hlbWVwYWNrIgp0ZXh0KDc1IDQ4LjUgIntzdFNhbnM6IHtpdDpyfX0gPSAuNzAiLCBwbGFjZW1lbnQoYykganVzdGlmaWNhdGlvbihsZWZ0KSBzaXplKG1lZGl1bSkgYm94IGZjb2xvcihlZGtibHVlKSBjb2xvcih3aGl0ZSkpIC8vLwp0ZXh0KDUwIDQ4LjUgIntzdFNhbnM6IHtpdDpyfX0gPSAuNjEiLCBwbGFjZW1lbnQoYykganVzdGlmaWNhdGlvbihsZWZ0KSBzaXplKG1lZGl1bSkgYm94IGZjb2xvcihlZGtibHVlKSBjb2xvcih3aGl0ZSkpIC8vLwp4dGl0bGUoWCBWYXJpYWJsZSk=