Merge pull request #307 from centreformicrosimulation/bugfix/302-assign-random-variation-to-mcs-pcs-life-satisfaction-scores

Bugfix/302 assign random variation to mcs pcs life satisfaction scores

Author: justin-ven

input/InitialPopulations/compile/01_prepare_UKHLS_pooled_data.do

@@ -198,10 +198,12 @@ preserve
 * Generate UC benefit marker
 gen benefits_uc=(ficode==40)
 label var benefits_uc "Universal Credit indicator"
+gen benefits_lb=(ficode == 15 | ficode == 16 | ficode == 19 | ficode == 20 | ficode == 22 | ficode == 33)
+label var benefits_lb "Legacy Benefits indicator"
 
 
-keep hidp pidp swv benefits_uc
-collapse (max) benefits_uc, by(hidp swv)
+keep hidp pidp swv benefits_uc benefits_lb
+collapse (max) benefits_uc benefits_lb, by(hidp swv)
 compress
 
 save "$dir_data/tmp_ucrcpt", replace

input/InitialPopulations/compile/02_create_UKHLS_variables.do

@@ -352,20 +352,43 @@ eststo predict_dhm: reg dhm c.dag i.dgn i.swv i.dhe, vce(robust) // Physical hea
 restore
 estimates restore predict_dhm
 predict dhm_prediction
-fre dhm_prediction
+// fre dhm_prediction
 
 gen dhm_flag = missing(dhm)
 replace dhm = round(dhm_prediction) if missing(dhm) 
 bys dhm_flag : sum dhm 
 
+/**************************Subjective well-being: GHQ 0-12 score *****************************/
+/*dhm_ghq   scghq2_dv    "DEMOGRAPHIC: Subjective wellbeing (GHQ): rescaled to 0-12, assumed to be continuous"
+This measure converts valid answers to 12 questions of the General Health Questionnaire (GHQ) to a 
+single scale by recoding the values of individual variables of 1 and 2 to 0, and values of 3 and 4 to 1 before summing them. 
+This produces a scale that ranges from 0 (indicating the least amount of distress) to 12 (indicating the greatest amount of distress).
+*/
+gen dhm_ghq	= scghq2_dv 
+replace dhm_ghq = . if scghq2_dv <0
+la var dhm_ghq "DEMOGRAPHIC: Subjective wellbeing (GHQ): 0-12 score"
+gen scghq2_dv_miss_flag = (scghq2_dv == .)
+
 
+preserve
+drop if dgn < 0 | dag<0 | dhe<0
+eststo predict_dhm_ghq: reg dhm_ghq c.dag i.dgn i.swv i.dhe, vce(robust) // Physical health has a big impact, so included as covariate.  
+restore
+estimates restore predict_dhm_ghq
+predict dhm_ghq_prediction
+
+gen dhm_ghq_flag = missing(dhm_ghq)
+replace dhm_ghq = round(dhm_ghq_prediction) if missing(dhm_ghq) 
+bys dhm_ghq_flag : sum dhm_ghq
+
+/* Alternative method of GHQ caseness - binary cutoff at scghq2_dv <4
 /**************************Subjective wellbeing (GHQ): Caseness ******************************
 0: not psychologically distressed, scghq2_dv < 4 
 1: psychologically distressed, scghq2_dv >= 4
 This measure converts valid answers to 12 questions of the General Health Questionnaire (GHQ) to a single scale by recoding 1 and 2 values 
 on individual variables to 0, and 3 and 4 values to 1, and then summing, giving a scale running from 0 (the least distressed) to 12 
 (the most distressed). A binary indicator is then created, equal to 1 for values >= 4.*/
-fre scghq2_dv
+// fre scghq2_dv
 recode scghq2_dv (-9/-1 . = .)
 gen scghq2_dv_miss_flag = (scghq2_dv == .)
 
@@ -451,12 +474,16 @@ replace dhe_pcssp=-9 if missing(dhe_pcssp) & idpartner>0
 //fre dhe_mcssp dhe_pcssp if idpartner>0
 
 
-/***************************** Life Satisfaction ***************************************************************************/
-/* Life satisfaction, self report. Continuous scale 0 to 7. */
+/***************************** Life Satisfaction 0-10 ****************************************************/
+/* Life satisfaction, self report. Continuous scale.
+  UKHLS variable sclfsato records scores 1-7.
+  Convert to 0-10 in line with ONS Life Satisfaction scale.
+ */
 
 
 gen dls = sclfsato
 replace dls = . if sclfsato < 0
+replace dls = (dls-1) * 10/6  // Change to 0-10 scale
 lab var dls "DEMOGRAPHIC: Life Satisfaction"
 // fre dls if dag>0 & dag<16
 
@@ -1751,6 +1778,11 @@ gen econ_benefits_uc=econ_benefits
 replace econ_benefits_uc=0 if benefits_uc==0
 label var econ_benefits_uc "Household income includes UC benefits"
 
+gen econ_benefits_lb=benefits_lb
+replace econ_benefits_lb=0 if benefits_lb==.
+replace econ_benefits_lb=0 if econ_benefits_uc==1
+label var econ_benefits_lb "Household income includes Legacy Benefits"
+
 
 /***************************** Financial Distress ***************************************************************************/
 // This is a measure of subjective financial distress, corresponding to answering 4 or 5 to the question below:

input/InitialPopulations/compile/09_finalise_input_data.do

@@ -205,7 +205,7 @@ rename dhe_mcs healthMentalMcs
 rename dhe_pcs healthPhysicalPcs
 rename dhe_mcssp healthMentalPartnerMcs
 rename dhe_pcssp healthPhysicalPartnerPcs
-rename dls demLifeSatScore1to7
+rename dls demLifeSatScore0to10
 rename dot demEthnC4
 rename dot01 demEthnC6
 rename financial_distress yFinDstrssFlag