数据文件中：fid为家庭编号，countyid为所在区县的编号，fincomp为家庭人均收入。我想基于家庭的人均收入计算出每个区县的基尼系数，并将生成的基尼系数保存到表中，即增加一个变量为每个家庭所在区县对应的基尼系数。数据为截面数据。我的程序如下：gen inequality=0
    qui ineqdeco fincomp , by(countyid)
qui levelsof countyid, local(countyid)
                         foreach j of local countyid {
             replace inequality=r(gini_j) if  countyid==`j'
｝

但循环的结果不对，请各位老师帮忙指导下，程序应该怎么写。