************
***** Invoked by: bhps_educ_master*.do 
************






*******************
** Descriptive analysis: 

use $path9\bh_temp.dta, clear


table wave , c(mean educ1 mean educ2 n educ1)
table wave if !missing(midage) , c(mean educ1 mean educ2 n educ1)
* Year-to-year, proportions with higher education increase; 
*   age matters to these proportions => restrict to ages 30-50 in the below. 

table wave [aweight=zxrwght] if !missing(midage) ///
     , c(mean educ1 mean educ2 n educ1)
* Weighting makes a bit of a difference (not surprising given differences in samples
*    but actually the magnitude of difference isn't huge) 



** Does education impact upon desire to move?

pwcorr loc1 educ1 educ2, star(05) sig
bysort wave: pwcorr loc1 educ1 educ2, star(05) 
table wave loc1 [aweight=zxrwght] if !missing(midage) ///
    , c(mean educ1 mean educ2 n educ1)
* By and large, some weak correlations, people with higher levels of education
*   are slightly more likely to want to move home 


summarize loc1 midage fem educ1 educ2 
xtlogit loc1 if !missing(educ1), i(pid)
est store nulla
xtlogit loc1 if !missing(educ1) & !missing(midage) & !missing(fem), i(pid)
est store nullb

xtlogit loc1 educ1 , i(pid)
est store mod1a
xtlogit loc1 educ2 , i(pid)
est store mod2a
xtlogit loc1 midage fem educ1 , i(pid)
est store mod1b
xtlogit loc1 midage fem educ2 , i(pid)
est store mod2b

est table nulla mod1a mod2a nullb  mod1b mod2b, stats(N ll bic rho) star b(%9.3g)



*******************