Расчет коэффициента Джини по данным LIS (в Stata)

Мне нужно рассчитать коэффициент Джини на основе данных о располагаемом личном доходе в LIS. Согласно учебному документу LIS, код Stata для этого:

    di "** INCOME DISTRIBUTION II – Exercise 13 **"
    program define bottop
    qui sum ey [w=hweight*d4]
    replace ey = .01*r(mean) if ey<.01*r(mean)
    qui sum dpi [w=hweight*d4], de
    replace ey = (10*r(p50)/(d4^.5)) if dpi>10*r(p50)
    end
    foreach file in $us00h $fi00h {
    display "`file'"
    use hweight d4 dpi if (!mi(dpi) & !(dpi==0)) using "`file'", clear
    gen ey=dpi/(d4^0.5)
    bottop
    ineqdeco ey [w=hweight*d4]
    }

Я просто скопировал и вставил этот код из учебного документа. Фрагменты

    qui sum ey [w=hweight*d4]
    replace ey=0.01*r(mean) if ey<0.01*r(mean)

и

    qui sum dpi [w=hweight*d4], de
    replace ey=(10*r(p50)/(d4^0.5)) if dpi>10*r(p50)

соответственно нижнее и верхнее кодирование.

Когда я попытался запустить этот код, переменная hweight не была найдена. Кто-нибудь знает, какое новое имя hweight в LIS? Или может кто-нибудь предложить, как я мог бы иначе преодолеть этот тупик?

Я знаком со stata, но сложность этого кода мне непонятна.

Очень признателен.