按照陈传波老师在《基尼系数的测算与分解--excel算法与stata程序》一文,在stata中输入如下程序
#delimit;
use"D:\BaiduNetdiskDownload\Stata17\nhb97.dta",clear;
egen tinc=sum(inc);
gen mi=inc/hhs;
sort mi;
egen thh=sum(hhs);
gen pi=hhs/thh;
save temp,replace;
program drop_all;
program define gini;
egen zsr=sum($sr);
gen u$sr=zsr/thh;
gen wi=$sr/zsr;
gen w$sr=zsr/tinc;
gen qi=sum(wi);
egen gini1=pi*(2*Qi-wi);
egen gini2=sum(gini1);
gen C$sr=1-gini2;
gen C$srW$sr=C$sr*w$sr;
drop zsr wi qi gini1 gini2;
save temp,replace;
end;
use temp,clear;
global sr"bc";
gini;
global sr"jy";
gini;
global sr"zy";
gini;
global sr"cc";
gini;
gen G=CbcWbc+CjyWjy+CzyWzy+CccWcc;
gen shGbc=CbcWbc/G;
gen shGjy=CjyWjy/G;
gen shGzy=CzyWzy/G;
gen shGcc=CccWcc/G;
list u *C * W* sh* in l;
执行后,stata一直报错(如下图所示)

想请大神帮我看看,是哪里出了问题,或者有没有其他更简便的程序可以计算按收入来源分组的基尼系数分解,非常感谢。