diff --git a/analysis/employment analysis.do b/analysis/employment analysis.do
new file mode 100644
index 000000000..3e34105c3
--- /dev/null
+++ b/analysis/employment analysis.do	
@@ -0,0 +1,83 @@
+/**************************************************************************************
+*
+*	PROGRAM TO ANALYSE EMPLOYMENT TRANSATIONS
+*
+*	Last version:  Justin van de Ven, 06 May 2025
+*	First version: Justin van de Ven, 06 May 2025
+*
+**************************************************************************************/
+
+clear all
+global moddir = "C:\Justin\dev\SimPaths\output\20250506123915\csv"
+global outdir = "C:\Justin\dev\SimPaths\analysis\"
+cd "$outdir"
+
+
+/**************************************************************************************
+*	start
+**************************************************************************************/
+global year_start = 2019
+global year_end = 2030
+
+
+/**************************************************************************************
+*	load data
+**************************************************************************************/
+import delimited using "$moddir/BenefitUnit.csv", clear
+rename *, l
+rename id_benefitunit idbenefitunit
+gsort idbenefitunit time
+save "$outdir/temp0", replace
+import delimited using "$moddir/Person.csv", clear
+rename *, l
+rename id_person idperson
+rename socialcareprovision socialcareprovision_p
+gsort idbenefitunit time idperson
+merge m:1 idbenefitunit time using temp0
+gsort time idbenefitunit idperson
+gen refbenefitunit = 0
+replace refbenefitunit = 1 if (idbenefitunit != idbenefitunit[_n-1])
+
+destring hoursworkedweekly, replace force
+recode hoursworkedweekly (missing=0)
+gen idNotEmployedAdult = (hoursworkedweekly<0.1 & dag>17)
+
+gen led = (deh_c3=="Low")
+gen med = (deh_c3=="Medium")
+gen hed = (deh_c3=="High")
+
+gen male = (dgn=="Male")
+
+gen idna = (dag>17)
+gen idnk = (dag<18)
+bys time idbenefitunit: egen partnered = sum(idna)
+replace partnered = partnered - 1
+bys time idbenefitunit: egen nk = sum(idnk)
+
+save "$outdir/temp1", replace
+
+
+/**************************************************************************************
+*	analysis
+**************************************************************************************/
+use "$outdir/temp1", clear
+global year_ref = $year_start-1
+
+// block 1 statistics
+gen emp = (hoursworkedweekly>0.5)
+gen emp_to_nemp = 0
+gen nemp_to_emp = 0
+gsort idperson time
+replace emp_to_nemp = 1 if (idperson[_n-1]==idperson & time[_n-1]+1==time & emp[_n-1]==1 & emp==0)
+replace nemp_to_emp = 1 if (idperson[_n-1]==idperson & time[_n-1]+1==time & emp[_n-1]==0 & emp==1)
+order time idperson hoursworkedweekly emp emp_to_nemp
+matrix store1 = J($year_end-$year_ref,2,.)
+forvalues yy = $year_start/$year_end {
+	qui{
+		sum emp_to_nemp if (time==`yy' & dag>17 & dag<65), mean
+		mat store1[`yy'-${year_ref},1] = r(mean)
+		sum nemp_to_emp if (time==`yy' & dag>17 & dag<65), mean
+		mat store1[`yy'-${year_ref},2] = r(mean)
+	}
+}
+matlist store1
diff --git a/input/DatabaseCountryYear.xlsx b/input/DatabaseCountryYear.xlsx
index b44d82113..0c791de8a 100644
Binary files a/input/DatabaseCountryYear.xlsx and b/input/DatabaseCountryYear.xlsx differ
diff --git a/input/EUROMODpolicySchedule.xlsx b/input/EUROMODpolicySchedule.xlsx
index a5d81548d..ce8925d5b 100644
Binary files a/input/EUROMODpolicySchedule.xlsx and b/input/EUROMODpolicySchedule.xlsx differ
diff --git a/input/InitialPopulations/compile/00_master.do b/input/InitialPopulations/compile/00_master.do
index 674a624ce..369476d23 100644
--- a/input/InitialPopulations/compile/00_master.do
+++ b/input/InitialPopulations/compile/00_master.do
@@ -8,7 +8,7 @@
 * DATA:         	    UKHLS EUL version - UKDA-6614-stata [to wave n]
 *						WAS EUL version - UKDA-7215-stata [to wave 7]
 * AUTHORS: 				Daria Popova, Justin van de Ven
-* LAST UPDATE:          14 Jan 2025
+* LAST UPDATE:          30 Apr 2025
 ***************************************************************************************
 
 ***************************************************************************************
@@ -36,8 +36,7 @@ set matsize 1000
 **************************************************************************************/
 
 * Working directory
-//global dir_work "C:\MyFiles\99 DEV ENV\JAS-MINE\data work\initial_populations"
-global dir_work "D:\Dasha\ESSEX\ESPON 2024\UK\initial_populations"
+global dir_work "C:\MyFiles\99 DEV ENV\JAS-MINE\data work\initial_populations"
 
 * Directory which contains do files
 global dir_do "${dir_work}/do"
@@ -49,12 +48,12 @@ global dir_data "${dir_work}/data"
 global dir_log "${dir_work}/log"
 
 * Directory which contains UKHLS data
-//global dir_ukhls_data "J:\01 DATA\UK\ukhls\wave13\stata\stata13_se\ukhls"
-global dir_ukhls_data "D:\Dasha\UK-original-data\USoc\UKDA-6614-stata\stata\stata13_se\ukhls"
+global dir_ukhls_data "J:\01 DATA\UK\ukhls\wave14\stata\stata13_se\ukhls"
+//global dir_ukhls_data "D:\Dasha\UK-original-data\USoc\UKDA-6614-stata\stata\stata13_se\ukhls"
 
 * Directory which contains WAS data
-//global dir_was_data "J:\01 DATA\UK\was\wave7\stata\stata13_se"
-global dir_was_data "D:\Dasha\UK-original-data\WAS\UKDA-7215-stata\stata\stata13_se"
+global dir_was_data "J:\01 DATA\UK\was\wave7\stata\stata13_se"
+//global dir_was_data "D:\Dasha\UK-original-data\WAS\UKDA-7215-stata\stata\stata13_se"
 
 * Directory which contains original initial popultions 
 global dir_ipop_orig "${dir_work}/original_initial_populations"