cap log close
clear
log using StataLog, replace 

include "http://staff.pubhealth.ku.dk/~linearpredictors/datafiles/readFever.do"

/* Tab 7.4.4 */

gen newweight=weight if weight>0
gen ga=(gwbirth>31)+(gwbirth>32)+(gwbirth>33)+(gwbirth>34)+(gwbirth>35)+(gwbirth>36)+(gwbirth>37)+(gwbirth>38)+(gwbirth>39)+(gwbirth>40)
gen parity2=(parity>0) if parity<.

gen bw1=newweight if ga==1
summarize bw1, detail
gen sga1=(newweight<r(p5)) if ga<2  & newweight<.
recode sga1 (0=0) (1=1) (.=0)

gen bw2=newweight if ga==2
summarize bw2, detail
gen sga2=(newweight<r(p5)) if ga==2  & newweight<.
recode sga2 (0=0) (1=1) (.=0)

gen bw3=newweight if ga==3
summarize bw3, detail
gen sga3=(newweight<r(p5)) if ga==3  & newweight<.
recode sga3 (0=0) (1=1) (.=0)

gen bw4=newweight if ga==4
summarize bw4, detail
gen sga4=(newweight<r(p5)) if ga==4  & newweight<.
recode sga4 (0=0) (1=1) (.=0)

gen bw5=newweight if ga==5
summarize bw5, detail
gen sga5=(newweight<r(p5)) if ga==5  & newweight<.
recode sga5 (0=0) (1=1) (.=0)

gen bw6=newweight if ga==6
summarize bw6, detail
gen sga6=(newweight<r(p5)) if ga==6  & newweight<.
recode sga6 (0=0) (1=1) (.=0)

gen bw7=newweight if ga==7
summarize bw7, detail
gen sga7=(newweight<r(p5)) if ga==7  & newweight<.
recode sga7 (0=0) (1=1) (.=0)

gen bw8=newweight if ga==8
summarize bw8, detail
gen sga8=(newweight<r(p5)) if ga==8  & newweight<.
recode sga8 (0=0) (1=1) (.=0)

gen bw9=newweight if ga==9
summarize bw9, detail
gen sga9=(newweight<r(p5)) if ga==9  & newweight<.
recode sga9 (0=0) (1=1) (.=0)

gen bw10=newweight if ga==10
summarize bw10, detail
gen sga10=(newweight<r(p5)) if ga==10  & newweight<.
recode sga10 (0=0) (1=1) (.=0)

gen sga=sga1*(ga<2)+sga2*(ga==2)+sga3*(ga==3)+sga4*(ga==4)+sga5*(ga==5)+sga6*(ga==6)+sga7*(ga==7)+sga8*(ga==8)+sga9*(ga==9)+sga10*(ga==10) if newweight<.

table sga smoke , col row
table sga parity2 , col row
table sga parity2, by(smoke) col row
table smoke, c(mean sga)
table parity2, c(mean sga)
table parity smoke, c(mean sga)

log close