* .do file for sbe15 insert by Wright and Royston

u dta/bpd, clear

xrigls bpd gawks, fp(m:2 2,s:1) cent(2.5 97.5) detail nogr
for M S Z C2.5 C97.5, any : ren @_gl @_sd
gr bpd C2.5_sd M_sd C97.5_sd gawks, c(.lll) s(oiii) sort xla /*
 */ yla(20,40,60,80,100) gap(5) saving(bpd1,replace)

xrigls bpd gawks, fp(m:2 2,s:df 0) cent(2.5 97.5) nogr
for M S Z C2.5 C97.5, any : ren @_gl @_co
gr bpd C2.5_co M_co C97.5_co gawks, c(.lll) s(oiii) sort xla /*
 */ yla(20,40,60,80,100) gap(5) saving(bpd2,replace)
gr Z_co gawks, xla yla gap(5) saving(bpd3,replace)

xrigls bpd gawks, fp(m:2 2,s:-2) cv cent(2.5 97.5) nogr
for M S Z C2.5 C97.5, any : ren @_gl @_cv
gen SD_cv=S_cv*M_cv
gr SD_cv S_sd gawks, c(ll) s(ii) so yla l1(Standard deviation) /*
 */ xla gap(5) saving(bpd4,replace)
gr C2.5_cv M_cv C97.5_cv C2.5_sd M_sd C97.5_sd gawks, c(llllll) /*
 */ s(iiiiii) sort xla yla(20,40,60,80,100) l1("Bi-parietal diameter (mm)") /*
 */ gap(5) saving(bpd5,replace)

xrigls bpd gawks, fp(m:2 2,s:1) cent(2.5 50 97.5) nogr se
for C2.5 C50 C97.5 ,any : gen l@=@_g-2*se@
for C2.5 C50 C97.5 ,any : gen u@=@_g+2*se@
gr  C2.5_gls  lC2.5 uC2.5 C50_gls  lC50 uC50 C97.5_gl  lC97.5 uC97.5 gawks /*
 */ if gawks>27, c(lllllllll) s(iiiiiiiii) xla yla /*
 */ l1("Bi-parietal diameter (mm)") gap(5) saving(bpd6,replace)