program define lsfit
    version 4.0

    local x     = "`1'"
    local y     = "`2'"

    tempvar z

    qui summarize `x'
    local xbar=_result(3)
    local sxx =_result(4)*(_N-1)
    di `xbar'
    qui summarize `y'
    local ybar=_result(3)
    local syy =_result(4)*(_N-1)

    gen `z'=(`x'-`xbar')*(`y'-`ybar')
    replace `z'=sum(`z')
    local sxy=`z'[_N]

    global slope = `sxy'/`sxx'
    global int   = `ybar'-$slope*`xbar'
    global corr  = $slope*sqrt(`sxx'/`syy')

end
exit