LBASF

clear all
import delimited gunssuicide.csv

gen logsuicide = log(suicides)
gen logplans = log(plans)
gen logattempts = log(attempts)

// list in 1/10

reg ownership logsuicide

graph twoway (scatter logsuicide ownership) (lfit logsuicide ownership), xtitle("Percent of Household Gun Ownership") ytitle("Logged Suicide Rates per 100,000") title("Suicides vs Household Gun Ownership") plotregion(color(white)) graphregion(color(white)) leg(off)