This do-file creates Table 1 (p. 60) in the SSR paper.

// Task: descriptive analysis; 2005 cgss; rural sample
// Project: Chinese Famine & Self-rated Health

**********************************************************************
*** 0. Program setup
**********************************************************************

version 13.0
clear all
macro drop _all
set linesize 80
set more off

**********************************************************************
*** 1. Identify sample
**********************************************************************
*********************************
*** 1.1. open data
*********************************

use "2014-02-24_clean10a_SocScienceRes_58notinfamine", clear

*********************************
*** 1.2. sample
*********************************

keep if rural_newfm == 1 & fmparty < .
drop if rural_b1f  == 0 | rural_b1m == 0
count
/* n = 1716*/

**********************************************************************
*** 2. EDR and health
**********************************************************************

tab edr1, m

// health factor analysis & reliability
factor FH1 FH2 FH3 FH4 FH5 FH6 FH7 FH8
alpha FH1 FH2 FH3 FH4 FH5 FH6 FH7 FH8

**********************************************************************
*** 3. Descriptives
**********************************************************************

tab1 cohort male pre famine post FH1-FH lninc edu edr1 fmparty age, m
tabstat FH FH1-FH8 edu lninc edr1 fmparty age, by(male) stat(mean sd count) long nototal format(%9.2f)

tabstat FH FH1-FH8 edu lninc fmparty age if cohort == 1, by(male) stat(mean sd count) long nototal format(%9.2f)
tabstat FH FH1-FH8 edu lninc fmparty age if cohort == 2, by(male) stat(mean sd count) long nototal format(%9.2f)
tabstat FH FH1-FH8 edu lninc fmparty age if cohort == 3, by(male) stat(mean sd count) long nototal format(%9.2f)

foreach var of varlist FH FH1-FH8 edu lninc fmparty age {
  display "variable == `var'"
  ttest `var', by(male)
}

foreach c of varlist pre famine post {
  foreach v of varlist FH FH1-FH8 edu lninc fmparty age {
    display "cohort = `c'; variable = `v'"
    ttest `v' if `c' == 1, by(male)
  }
}

foreach v of varlist FH FH1-FH8 edu lninc fmparty age {
  display "==> `v'"
  anova `v' cohort
}