* !hillp.ado version 1.0.1 June/2000    (STB-56: sg147)
* Manuel G Scotto
* Hill Plot
program define hillp
version 6.0
local varlist "req ex min(1) max(1)"
local if "opt"
local in "opt"
parse "`*'"
parse "`varlist'", parse (" ")
local X `1'
preserve
if ("`if'"!="") {
        qui keep `if'
}
if ("`in'"!="") {
        qui keep `in'
}
postfile hill value y using points, replace
sort `X'
gen xl=log(`X')
local k=2
local num=xl[_N]
local y=1
while (`k'<=_N){
local num=`num'+xl[_N-`k'+1]
local value=`num'/`k'
local value=`value'-xl[_N-`k'+1]
local value=1/`value'
post hill `value' `y'
local y=`y'+1
local k=`k'+1
}
postclose hill
clear
use points
label var value "hill estimate of alpha"
label var y "k"
graph value y, s(.) border
end