u humer, clear

* SD constant
xriml humer gest, fp(m:0.5 3) di(n) cent(5 95) nogr
gr humer C5_ml M_ml C95_ml gest, c(.lll) s(oiii) so xla yla bord /*
 */ gap(5) t1([a]) saving(fig1a,replace)
ren Z_ml Z
lab var Z "Z-scores"
gr Z gest, xla yla bord gap(5) t1([b]) saving(fig1b,replace)
xriqtest Z gest, par(m:3,s:1) gen
gr _sd _group, c(J) s(i)  yli(1) l1("_sd") bord gap(5) xla yla saving(fig2,replace)
permband Z gest, mom(2) seed(125) bord gap(5) xla yla saving(fig3,replace)
 
* SD linear
xriml humer gest, fp(m:0.5 3,s:1) di(n) cent(5 95) nogr
ren Z_ml Znew
lab var Znew "new Z-scores"
xriqtest Znew gest, par(m:3,s:2)
permband Znew gest, mom(3) seed(129) yline(0) xla yla bord gap(5) saving(fig4,replace)