matrix define A=J(5,80,0) 
 
 
 

 
local fifi  " de10 dk10 es10 gr10"  
 
 
 
foreach toto in `fifi' {  
local perso "$`toto'p" 
local house "$`toto'h" 
qui use hid ppopwgt age sex relation educ  nchildren immigr educ_c using `perso'  , clear 
qui joinby hid using `house'  
keep hid ppopwgt  dname age sex relation educ year iso2 hpopwgt  dhi hxit hits hicvip pension nchildren immigr educ_c npers  regio 
local save   "t`toto'"  
qui save  `save' , replace 
quietly : gen dpi=dhi 
quietly : gen d3=. 
quietly : keep if dpi>0 & dpi!=. 
 
 
 
quietly : gen country=iso2 
capture :gen hw= hpopwgt 
quietly :gen d4= npers   

keep if age>=25 & age<=54 

quietly :gen wgt=int(hw*d4) 
quietly :gen rev=exp(1+.01*rnormal())  
quietly :replace rev=((dpi)/sqrt(d4)) 
replace hxit=0 if hxit==. 
replace hits=0 if hits==. 
 
 
quietly :gen red=(dhi+hxit-hits)/dhi  
quietly :gen bri=(dhi+hxit-hits)  
  
local ii = `ii' + 1 
 
qui { 
quietly :sort country year rev 
quietly :gen x=. 
quietly :gen v=. 
quietly :gen pop=. 
quietly :gen COUNTRY = dname 
gen pw=int(wgt)  
qui su rev [w=pw]  
qui replace rev=rev/r(mean) 
qui replace v=rev 
 
local ctry "`toto'" 
noi di `ii' "    " "`ctry'" 
 
qui {  
su v  [fw=wgt]  if COUNTRY == "`ctry'" 
qui ineqdeco v  [fw=wgt]  if COUNTRY == "`ctry'" 
matrix A[`ii',1]= r(a2)  
matrix A[`ii',2]= r(a1)  
matrix A[`ii',3]= r(ahalf)  
matrix A[`ii',4]= r(ge2)  
matrix A[`ii',5]= r(ge1)  
matrix A[`ii',6]= r(ge0)  
matrix A[`ii',7]= r(gem1)  
matrix A[`ii',8]= r(gini)  
matrix A[`ii',13]= r(p90)/r(p50)  
matrix A[`ii',14]= r(p50)/r(p10)  
matrix A[`ii',15]= r(p90)/r(p75)/(r(p75)/r(p50))  
matrix A[`ii',16]= r(p25)/r(p10)/(r(p50)/r(p25)) 
matrix A[`ii',17]= r(p75p25) 
 
gen lv=ln(v) 
su lv [fw=wgt]  if COUNTRY == "`ctry'" , d 
matrix A[`ii',10]= r(skewness)  
matrix A[`ii',11]= r(kurtosis)  
matrix A[`ii',12]= r(sd)  
 
su lv if COUNTRY == "`ctry'" , d 
 
matrix A[`ii',19]=r(N) 
noi di r(N) 
drop lv 
 
quietly : svyset _n [pweight=wgt], vce(linearized) singleunit(missing) 
capture :  ipolfw v   
matrix A[`ii',18]=e(est) 
svyset, clear 
 
povdeco v   [fw=wgt]  if  COUNTRY == "`ctry'"  , pline(.5) 
matrix A[`ii',21]=r(fgt2) 
matrix A[`ii',22]=r(fgt1) 
matrix A[`ii',23]=r(fgt0) 
matrix A[`ii',24]= r(meangappoor) 
matrix A[`ii',25]= r(meanpoor)  
 
qui su v  [fw=pw]   , d 
qui replace rev=v/r(p50) 
qui gen mc=rev>.75 & rev<1.25 
qui su mc  [fw=pw]   
matrix A[`ii',26]=r(mean) 
qui gen mcl=rev>.50 & rev<.75 
qui su mcl  [fw=pw]   
matrix A[`ii',27]=r(mean) 
qui gen mcu=rev>1.25 & rev<2 
qui su mcu  [fw=pw]   
matrix A[`ii',28]=r(mean) 
qui gen poo=rev<.5  
qui su poo  [fw=pw]   
matrix A[`ii',29]=r(mean) 
qui gen ric=rev>2 
qui su ric  [fw=pw]   
matrix A[`ii',30]=r(mean) 
 
   
capture graph drop ga 
 
isograph  v  [w=wgt]  if  COUNTRY == "`ctry'"
mat def is=r(isograph)
 
forvalue h=41/52 {
local g=`h'-40
matrix A[`ii',`h']=is[`g',2]
}
 
 
  
isograph rev [pweight=wgt] if COUNTRY == "`ctry'"
mat def is=r(isograph)
 
forvalue h=61/72 {
local g=`h'-60
matrix A[`ii',`h']=is[`g',2]
}

} 
 
} 
} 
 
noi mat li A