根据这个公式写代码,有没有大佬会啊,我的数据是这样的
capture program drop marshall
program define marshall
version 15
gen s=0
syntax varlist (min=1 max=1)
local p: word 1 of `varlist'
sum `p'
local t=r(N)
local u=r(sum)
forvalue i=1/`t'{
forvalue j=1/`t'{
qui replace s=s+abs(`p'[`i']-`p'[`j'])
}
}
local mgini=s/(2*`u'*(`t'-1))
drop s
dis _n in g _dup(20) "=" in y "Marshall GINI is " in g _dup(20) "="
dis `mgini'
end