use "H:\529\WFA\c3 activities data no 014 in zemp.dta", clear

regress WAGES023 jwcw jsocw pecw jtcw cjcw jwrpre cjrpre werpre perpre jtrpre jsrpre totpre FRST9899 FRST9900 FRST0001 ednone edged edvotec edcoll lddum adddum felony healthpf wklimp mentalh CHILDLT3 childrn kidsickd eastr westr eastu AGELT20 AGE31_40 AGEOVR40 rlatino rafam rapi rnatam parwelf UNEMPL02 UE02MISS english teenpreg drugany alcohol JWC973 JWC974 JWC981 JWC982 JWC983 JWC984 JWC991 JWC992 JWC993 JWC994 JWC001 JWC002 JWC003 JWC004 JWC011 JWC012 JWC013 JSC973 JSC974 JSC981 JSC982 JSC983 JSC984 JSC991 JSC992 JSC993 JSC994 JSC001  JSC002 JSC003 JSC004 JSC011 JSC012 JSC013 PEC973 PEC974 PEC981 PEC982 PEC983 PEC984 PEC991 PEC992 PEC993 PEC994 PEC001 PEC002 PEC003 PEC004 PEC011 PEC012 PEC013 JTC973 JTC974 JTC981 JTC982 JTC983 JTC984 JTC991 JTC992 JTC993 JTC994 JTC001 JTC002 JTC003 JTC004 JTC011 JTC012 JTC013 WAGES991 WAGES992 WAGES993 WAGES994 WAGES001 WAGES002 WAGES003 WAGES004 WAGES011 WAGES012 WAGES013 EMP991 EMP992 EMP993 EMP994 EMP001 EMP002 EMP003 EMP004 EMP011 EMP012 EMP013 


predict yhatbase if e(sample), xb

generate jwcwreal=jwcw

replace jwcw=0

predict yhat0 if e(sample), xb

replace jwcw=1

predict yhat1 if e(sample), xb

summarize yhatbase yhat1 yhat0 if jwcwreal==1

twoway (scatter yhat1 yhat0 if jwcwreal==1)

*******************DON"T SAVE THE DATA IF YOU MADE CHANGES!!!!!!!!!!!!!!!********

*****now logit***********************************************

*******OPEN NEW DATA VERSION***********************************
use "H:\529\WFA\c3 activities data no 014 in zemp.dta", clear
logit EMP023 jwcw jsocw pecw jtcw cjcw jwrpre cjrpre werpre perpre jtrpre jsrpre totpre FRST9899 FRST9900 FRST0001 ednone edged edvotec edcoll lddum adddum felony healthpf wklimp mentalh CHILDLT3 childrn kidsickd eastr westr eastu AGELT20 AGE31_40 AGEOVR40 rlatino rafam rapi rnatam parwelf UNEMPL02 UE02MISS english teenpreg drugany alcohol JWC973 JWC974 JWC981 JWC982 JWC983 JWC984 JWC991 JWC992 JWC993 JWC994 JWC001 JWC002 JWC003 JWC004 JWC011 JWC012 JWC013 JSC973 JSC974 JSC981 JSC982 JSC983 JSC984 JSC991 JSC992 JSC993 JSC994 JSC001  JSC002 JSC003 JSC004 JSC011 JSC012 JSC013 PEC973 PEC974 PEC981 PEC982 PEC983 PEC984 PEC991 PEC992 PEC993 PEC994 PEC001 PEC002 PEC003 PEC004 PEC011 PEC012 PEC013 JTC973 JTC974 JTC981 JTC982 JTC983 JTC984 JTC991 JTC992 JTC993 JTC994 JTC001 JTC002 JTC003 JTC004 JTC011 JTC012 JTC013 WAGES991 WAGES992 WAGES993 WAGES994 WAGES001 WAGES002 WAGES003 WAGES004 WAGES011 WAGES012 WAGES013 EMP991 EMP992 EMP993 EMP994 EMP001 EMP002 EMP003 EMP004 EMP011 EMP012 EMP013 



predict phatbase if e(sample), pr

generate jwcwreal=jwcw

replace jwcw=0

predict phat0 if e(sample), pr

replace jwcw=1

predict phat1 if e(sample), pr

summarize phatbase phat1 phat0 if jwcwreal==1

twoway (scatter phat1 phat0 if jwcwreal==1)