diff --git a/R/lfs_clean_health.R b/R/lfs_clean_health.R
index ef567f6..8250acc 100644
--- a/R/lfs_clean_health.R
+++ b/R/lfs_clean_health.R
@@ -8,53 +8,188 @@ lfs_clean_health <- function(
     data
 ) {
 
+  ################################################
   # Disability (under the Equality Act 2010)
   if ("disea" %in% colnames(data)) {
   data[disea == 1 , disability := "disabled"]
   data[disea == 2 , disability := "not_disabled"]
   } else {
 
-  data[, disability := NULL]
+  data[, disability := NA]
   }
 
-  # Main health problem/disability
+  #############################################################################################
+  # Main health problem/disability (that limits the kind of paid work the respondent can do)
+
+  if ("health" %in% colnames(data)){
+  data[ , main_health_condition := "none"]
+  data[health == 1 , main_health_condition := "msk_arms_hands"]
+  data[health == 2 , main_health_condition := "msk_legs_feet"]
+  data[health == 3 , main_health_condition := "msk_back_neck"]
+  data[health == 4 , main_health_condition := "difficulty_seeing"]
+  data[health == 5 , main_health_condition := "difficulty_hearing"]
+  data[health == 6 , main_health_condition := "speech_impediment"]
+  data[health == 7 , main_health_condition := "skin_conditions_disfigurement"]
+  data[health == 8 , main_health_condition := "respiratory_problems"]
+  data[health == 9 , main_health_condition := "heart_blood_circulatory_problems"]
+  data[health == 10 , main_health_condition := "digestive_problems"]
+  data[health == 11 , main_health_condition := "diabetes"]
+  data[health == 12 , main_health_condition := "depression_anxiety"]
+  data[health == 13 , main_health_condition := "epilepsy"]
+  data[health == 18 , main_health_condition := "autism"]
+  data[health == 14 , main_health_condition := "learning_difficulties"]
+  data[health == 15 , main_health_condition := "mental_illness"]
+  data[health == 16 , main_health_condition := "progressive_illness_nec"]
+  data[health == 17 , main_health_condition := "other_problems"]
+
+
+  ## only applies to those who are working age so set NA if not otherwise
+  ## the "none" category will be incorrect
+  data[(sex == 1 & !(age %in% 16:64)) | (sex == 2 & !(age %in% 16:59)), main_health_condition := NA]
+  } else {
+
+  data[, main_health_condition := NA]
+  }
 
   if ("health" %in% colnames(data)){
-  data[ , health_condition := "none"]
-  data[health == 1 , health_condition := "msk_arms_hands"]
-  data[health == 2 , health_condition := "msk_legs_feet"]
-  data[health == 3 , health_condition := "msk_back_neck"]
-  data[health == 4 , health_condition := "difficulty_seeing"]
-  data[health == 5 , health_condition := "difficulty_hearing"]
-  data[health == 6 , health_condition := "speech_impediment"]
-  data[health == 7 , health_condition := "skin_conditions_disfigurement"]
-  data[health == 8 , health_condition := "respiratory_problems"]
-  data[health == 9 , health_condition := "heart_blood_circulatory_problems"]
-  data[health == 10 , health_condition := "digestive_problems"]
-  data[health == 11 , health_condition := "diabetes"]
-  data[health == 12 , health_condition := "depression_anxiety"]
-  data[health == 13 , health_condition := "epilepsy"]
-  data[health == 18 , health_condition := "autism"]
-  data[health == 14 , health_condition := "learning_difficulties"]
-  data[health == 15 , health_condition := "mental_illness"]
-  data[health == 16 , health_condition := "progressive_illness_nec"]
-  data[health == 17 , health_condition := "other_problems"]
+    data[ , main_health_condition_8cat := "none"]
+    data[health %in% 1:3 , main_health_condition_8cat := "musculoskeletal"]
+    data[health %in% 4:5 , main_health_condition_8cat := "sight_hearing"]
+    data[health %in% 8 , main_health_condition_8cat := "respiratory"]
+    data[health %in% 9:11 , main_health_condition_8cat := "blood_digestive_diabetes"]
+    data[health %in% c(12,15) , main_health_condition_8cat := "mental_health"]
+    data[health %in% c(14,18) , main_health_condition_8cat := "lda"]
+    data[health %in% c(16,17,13,6,7) , main_health_condition_8cat := "other_health_cond"]
+
+    ## only applies to those who are working age so set NA if not otherwise
+    ## the "none" category will be incorrect
+    data[(sex == 1 & !(age %in% 16:64)) | (sex == 2 & !(age %in% 16:59)), main_health_condition_8cat := NA]
   } else {
 
-  data[, health_condition := NULL]
+    data[, main_health_condition_8cat := NA]
   }
 
+  ############################################################################
+  # Dummies for each health condition (in the 8 category specification above)
+
+  if ("health" %in% colnames(data)){
+    data[, musculoskeletal := ifelse(heal01 %in% 1:3 |
+                                     heal02 %in% 1:3 |
+                                     heal03 %in% 1:3 |
+                                     heal04 %in% 1:3 |
+                                     heal05 %in% 1:3 |
+                                     heal06 %in% 1:3 |
+                                     heal07 %in% 1:3 |
+                                     heal08 %in% 1:3 |
+                                     heal09 %in% 1:3 |
+                                     heal10 %in% 1:3, 1, 0)]
+
+    data[, sight_hearing := ifelse(heal01 %in% 4:5 |
+                                   heal02 %in% 4:5 |
+                                   heal03 %in% 4:5 |
+                                   heal04 %in% 4:5 |
+                                   heal05 %in% 4:5 |
+                                   heal06 %in% 4:5 |
+                                   heal07 %in% 4:5 |
+                                   heal08 %in% 4:5 |
+                                   heal09 %in% 4:5 |
+                                   heal10 %in% 4:5, 1, 0)]
+
+    data[, respiratory := ifelse(heal01 %in% 8 |
+                                 heal02 %in% 8 |
+                                 heal03 %in% 8 |
+                                 heal04 %in% 8 |
+                                 heal05 %in% 8 |
+                                 heal06 %in% 8 |
+                                 heal07 %in% 8 |
+                                 heal08 %in% 8 |
+                                 heal09 %in% 8 |
+                                 heal10 %in% 8, 1, 0)]
+
+    data[, blood_digestive_diabetes := ifelse(heal01 %in% 9:11 |
+                                              heal02 %in% 9:11 |
+                                              heal03 %in% 9:11 |
+                                              heal04 %in% 9:11 |
+                                              heal05 %in% 9:11 |
+                                              heal06 %in% 9:11 |
+                                              heal07 %in% 9:11 |
+                                              heal08 %in% 9:11 |
+                                              heal09 %in% 9:11 |
+                                              heal10 %in% 9:11, 1, 0)]
+
+    data[, mental_health := ifelse(heal01 %in% c(12,15) |
+                                   heal02 %in% c(12,15) |
+                                   heal03 %in% c(12,15) |
+                                   heal04 %in% c(12,15) |
+                                   heal05 %in% c(12,15) |
+                                   heal06 %in% c(12,15) |
+                                   heal07 %in% c(12,15) |
+                                   heal08 %in% c(12,15) |
+                                   heal09 %in% c(12,15) |
+                                   heal10 %in% c(12,15), 1, 0)]
+
+    data[, lda := ifelse(heal01 %in% c(14,18) |
+                         heal02 %in% c(14,18) |
+                         heal03 %in% c(14,18) |
+                         heal04 %in% c(14,18) |
+                         heal05 %in% c(14,18) |
+                         heal06 %in% c(14,18) |
+                         heal07 %in% c(14,18) |
+                         heal08 %in% c(14,18) |
+                         heal09 %in% c(14,18) |
+                         heal10 %in% c(14,18), 1, 0)]
+
+
+    data[, other_health_cond := ifelse(heal01 %in% c(16,17,13,6,7) |
+                                       heal02 %in% c(16,17,13,6,7) |
+                                       heal03 %in% c(16,17,13,6,7) |
+                                       heal04 %in% c(16,17,13,6,7) |
+                                       heal05 %in% c(16,17,13,6,7) |
+                                       heal06 %in% c(16,17,13,6,7) |
+                                       heal07 %in% c(16,17,13,6,7) |
+                                       heal08 %in% c(16,17,13,6,7) |
+                                       heal09 %in% c(16,17,13,6,7) |
+                                       heal10 %in% c(16,17,13,6,7), 1, 0)]
+
+    data[(sex == 1 & !(age %in% 16:64)) | (sex == 2 & !(age %in% 16:59)), musculoskeletal := NA]
+    data[(sex == 1 & !(age %in% 16:64)) | (sex == 2 & !(age %in% 16:59)), sight_hearing := NA]
+    data[(sex == 1 & !(age %in% 16:64)) | (sex == 2 & !(age %in% 16:59)), respiratory := NA]
+    data[(sex == 1 & !(age %in% 16:64)) | (sex == 2 & !(age %in% 16:59)), blood_digestive_diabetes := NA]
+    data[(sex == 1 & !(age %in% 16:64)) | (sex == 2 & !(age %in% 16:59)), mental_health := NA]
+    data[(sex == 1 & !(age %in% 16:64)) | (sex == 2 & !(age %in% 16:59)), lda := NA]
+    data[(sex == 1 & !(age %in% 16:64)) | (sex == 2 & !(age %in% 16:59)), other_health_cond := NA]
+
+  } else {
+
+    data[, musculoskeletal := NA]
+    data[, sight_hearing := NA]
+    data[, respiratory := NA]
+    data[, blood_digestive_diabetes := NA]
+    data[, mental_health := NA]
+    data[, lda := NA]
+    data[, other_health_cond := NA]
+
+  }
+
+  ####################################
+  # Main health problem/disability
+
   ### convert the variables to factors
-  data$disability          <- as.factor(data$disability)
-  data$health_condition    <- as.factor(data$health_condition)
+  data$disability                 <- as.factor(data$disability)
+  data$main_health_condition      <- as.factor(data$main_health_condition)
+  data$main_health_condition_8cat <- as.factor(data$main_health_condition_8cat)
 
   ###############################
   ### RETAIN CLEANED VARIABLES
 
   final_data <- data[, c("obs_id",
-                         "disability", "health_condition")]
+                         "disability", "main_health_condition", "main_health_condition_8cat",
+                         "musculoskeletal","sight_hearing","respiratory","blood_digestive_diabetes",
+                         "mental_health","lda","other_health_cond")]
 
-  var_names <- c("disability", "health_condition")
+  var_names <- c("disability", "main_health_condition", "main_health_condition_8cat",
+                 "musculoskeletal","sight_hearing","respiratory","blood_digestive_diabetes",
+                 "mental_health","lda","other_health_cond")
 
   setnames(final_data, var_names, paste0("h_", var_names))
 
diff --git a/R/lfs_read_2000.R b/R/lfs_read_2000.R
index ad2daed..29800bf 100644
--- a/R/lfs_read_2000.R
+++ b/R/lfs_read_2000.R
@@ -60,11 +60,13 @@ lfs_read_2000 <- function(
     demographic_vars <- Hmisc::Cs(age,sex,govtor,ethcen,marstt,fdpch16)
     education_vars   <- Hmisc::Cs(edage,hiqual,hiquald,btec,sctvec,gnvq,nvqlev,rsa,candg,numal,numas,numol,qgcse,
                                   quals01,quals02,quals03,quals04,quals05,quals06,quals07,quals08,quals09,quals10,quals11,qualch1)
-    health_vars      <- Hmisc::Cs(health,discurr)
+    health_vars      <- Hmisc::Cs(health,discurr,
+                                  heal01,heal02,heal03,heal04,heal05,
+                                  heal06,heal07,heal08,heal09,heal10)
     work_vars        <- Hmisc::Cs(inecacr,grsswk,ftptwk,ttachr,ttushr,mpnor,publicr,indm92m,indd92m,inds92m,socmain,socminm)
     other_vars       <- Hmisc::Cs(refwkm,thiswv)
 
-    names <- c(demographic_vars,education_vars,work_vars, weights_vars,other_vars)
+    names <- c(demographic_vars,education_vars,work_vars,health_vars,weights_vars,other_vars)
     names <- tolower(names)
 
     data <- data[ ,names, with=F]
diff --git a/R/lfs_read_2001.R b/R/lfs_read_2001.R
index 84061c9..b7461a0 100644
--- a/R/lfs_read_2001.R
+++ b/R/lfs_read_2001.R
@@ -60,7 +60,9 @@ lfs_read_2001 <- function(
     demographic_vars <- Hmisc::Cs(age,sex,govtof,fdpch16)
     education_vars   <- Hmisc::Cs(edage,btec,sctvec,gnvq,nvqlev,rsa,candg,numal,numas,numol,gcse,qgcse,
                                   quals01,quals02,quals03,quals04,quals05,quals06,quals07,quals08,quals09,quals10,quals11,qualch1)
-    health_vars      <- Hmisc::Cs(health,discurr)
+    health_vars      <- Hmisc::Cs(health,discurr,
+                                  heal01,heal02,heal03,heal04,heal05,
+                                  heal06,heal07,heal08,heal09,heal10)
     work_vars        <- Hmisc::Cs(inecacr,     ftptwk,ttachr,ttushr,      publicr,indm92m,indd92m,inds92m)
     other_vars       <- Hmisc::Cs(refwkm,thiswv)
 
@@ -103,7 +105,9 @@ lfs_read_2001 <- function(
     demographic_vars <- Hmisc::Cs(age,sex,govtof,eth01,marstt,fdpch16)
     education_vars   <- Hmisc::Cs(edage,btec,sctvec,gnvq,nvqlev,rsa,candg,numal,numas,numol,
                                   quals01,quals02,quals03,quals04,quals05,quals06,quals07,quals08,quals09,quals10,quals11,qualch1)
-    health_vars      <- Hmisc::Cs(health,discurr)
+    health_vars      <- Hmisc::Cs(health,discurr,
+                                  heal01,heal02,heal03,heal04,heal05,
+                                  heal06,heal07,heal08,heal09,heal10)
     work_vars        <- Hmisc::Cs(inecacr,grsswk,ftptwk,ttachr,ttushr,mpnr01,publicr,indm92m,indd92m,inds92m,soc2km,sc2kmmn)
     other_vars       <- Hmisc::Cs(refwkm)
 
@@ -147,10 +151,13 @@ lfs_read_2001 <- function(
     demographic_vars <- Hmisc::Cs(age,sex,govtof,eth01,marstt,fdpch16)
     education_vars   <- Hmisc::Cs(hiqual,hiquald,btec,sctvec,gnvq,nvqlev,rsa,candg,numal,numas,numol,
                                   quals01,quals02,quals03,quals04,quals05,quals06,quals07,quals08,quals09,quals10,quals11,qualch1)
+    health_vars      <- Hmisc::Cs(health,discurr,
+                                  heal01,heal02,heal03,heal04,heal05,
+                                  heal06,heal07,heal08,heal09,heal10)
     work_vars        <- Hmisc::Cs(inecacr,grsswk,ftptwk,ttachr,ttushr,mpnr01,publicr,indm92m,indd92m,inds92m,soc2km,sc2kmmn)
     other_vars       <- Hmisc::Cs(refwkm)
 
-    names <- c(demographic_vars,education_vars,work_vars, weights_vars,other_vars)
+    names <- c(demographic_vars,education_vars,work_vars,health_vars,weights_vars,other_vars)
     names <- tolower(names)
 
     data <- data[ ,names, with=F]
diff --git a/R/lfs_read_2002.R b/R/lfs_read_2002.R
index d1276c3..7c7876e 100644
--- a/R/lfs_read_2002.R
+++ b/R/lfs_read_2002.R
@@ -60,11 +60,13 @@ lfs_read_2002 <- function(
     demographic_vars <- Hmisc::Cs(age,sex,govtof,eth01,marstt,fdpch16)
     education_vars   <- Hmisc::Cs(edage,hiqual,hiquald,btec,sctvec,gnvq,nvqlev,rsa,candg,numal,numas,numol,gcse,qgcse,
                                   quals01,quals02,quals03,quals04,quals05,quals06,quals07,quals08,quals09,quals10,quals11,qualch1)
-    health_vars      <- Hmisc::Cs(health,discurr)
+    health_vars      <- Hmisc::Cs(health,discurr,
+                                  heal01,heal02,heal03,heal04,heal05,
+                                  heal06,heal07,heal08,heal09,heal10)
     work_vars        <- Hmisc::Cs(inecacr,grsswk,ftptwk,ttachr,ttushr,mpnr01,publicr,indm92m,indd92m,inds92m,soc2km,sc2kmmn)
     other_vars       <- Hmisc::Cs(refwkm,thiswv)
 
-    names <- c(demographic_vars,education_vars,work_vars, weights_vars,other_vars)
+    names <- c(demographic_vars,education_vars,work_vars,health_vars,weights_vars,other_vars)
     names <- tolower(names)
 
     data <- data[ ,names, with=F]
@@ -104,12 +106,14 @@ lfs_read_2002 <- function(
     demographic_vars <- Hmisc::Cs(age,sex,govtof,eth01,marstt,fdpch16)
     education_vars   <- Hmisc::Cs(edage,btec,sctvec,gnvq,nvqlev,rsa,candg,numal,numas,numol,
                                   quals01,quals02,quals03,quals04,quals05,quals06,quals07,quals08,quals09,quals10,quals11,qualch1)
-    health_vars      <- Hmisc::Cs(health,discurr)
+    health_vars      <- Hmisc::Cs(health,discurr,
+                                  heal01,heal02,heal03,heal04,heal05,
+                                  heal06,heal07,heal08,heal09,heal10)
     work_vars        <- Hmisc::Cs(inecacr,grsswk,ftptwk,ttachr,ttushr,mpnr02,publicr,indm92m,indd92m,inds92m,soc2km,sc2kmmn,
                                   undemp,undhrs,ovhrs,lespay2)
     other_vars       <- Hmisc::Cs(refwkm)
 
-    names <- c(demographic_vars,education_vars,work_vars, weights_vars,other_vars)
+    names <- c(demographic_vars,education_vars,work_vars,health_vars,weights_vars,other_vars)
     names <- tolower(names)
 
     data <- data[ ,names, with=F]
diff --git a/R/lfs_read_2003.R b/R/lfs_read_2003.R
index fc2c8de..e039363 100644
--- a/R/lfs_read_2003.R
+++ b/R/lfs_read_2003.R
@@ -60,12 +60,14 @@ lfs_read_2003 <- function(
     demographic_vars <- Hmisc::Cs(age,sex,govtof,eth01,marstt,fdpch16)
     education_vars   <- Hmisc::Cs(edage,hiqual,hiquald,btec,sctvec,gnvq,nvqlev,rsa,candg,numal,numas,numol,gcse,qgcse,
                                   quals01,quals02,quals03,quals04,quals05,quals06,quals07,quals08,quals09,quals10,quals11,qualch1)
-    health_vars      <- Hmisc::Cs(health,discurr)
+    health_vars      <- Hmisc::Cs(health,discurr,
+                                  heal01,heal02,heal03,heal04,heal05,
+                                  heal06,heal07,heal08,heal09,heal10)
     work_vars        <- Hmisc::Cs(inecacr,grsswk,ftptwk,ttachr,ttushr,mpnr02,publicr,indm92m,indd92m,inds92m,soc2km,sc2kmmn,
                                   undemp,undhrs,ovhrs,lespay2)
     other_vars       <- Hmisc::Cs(refwkm,thiswv)
 
-    names <- c(demographic_vars,education_vars,work_vars, weights_vars,other_vars)
+    names <- c(demographic_vars,education_vars,work_vars,health_vars,weights_vars,other_vars)
     names <- tolower(names)
 
     data <- data[ ,names, with=F]
diff --git a/R/lfs_read_2004.R b/R/lfs_read_2004.R
index f83fbab..7887e95 100644
--- a/R/lfs_read_2004.R
+++ b/R/lfs_read_2004.R
@@ -96,12 +96,14 @@ lfs_read_2004 <- function(
                                   gcseful1,gcseful2,gcseful3,gcseful4,qgcse41,qgcse42,qgcse43,qgcse44,
                                   quals401,quals402,quals403,quals404,quals405,quals406,quals407,quals408,quals409,quals410,quals411,
                                   subcod1,subcod2,subcod3,subcod4,subcod5,subcod6,subcod7,subcod8)
-    health_vars      <- Hmisc::Cs(health,discurr)
+    health_vars      <- Hmisc::Cs(health,discurr,
+                                  heal01,heal02,heal03,heal04,heal05,
+                                  heal06,heal07,heal08,heal09,heal10)
     work_vars        <- Hmisc::Cs(inecacr,grsswk,ftptwk,ttachr,ttushr,mpnr02,publicr,indm92m,indd92m,inds92m,soc2km,sc2kmmn,
                                   undemp,undhrs,ovhrs,lespay2)
     other_vars       <- Hmisc::Cs(refwkm)
 
-    names <- c(demographic_vars,education_vars,work_vars, weights_vars,other_vars)
+    names <- c(demographic_vars,education_vars,work_vars,health_vars,weights_vars,other_vars)
     names <- tolower(names)
 
     data <- data[ ,names, with=F]
diff --git a/R/lfs_read_2005.R b/R/lfs_read_2005.R
index cfd2ef4..9b7cf88 100644
--- a/R/lfs_read_2005.R
+++ b/R/lfs_read_2005.R
@@ -60,12 +60,14 @@ lfs_read_2005 <- function(
     demographic_vars <- Hmisc::Cs(age,sex,govtof,eth01,marstt,fdpch16)
     education_vars   <- Hmisc::Cs(edage,hiqual4d,hiqual4,btec,sctvec,gnvq4,nvqlev,rsa,candg,numal,numas,hst,advhst,
                                   quals401,quals402,quals403,quals404,quals405,quals406,quals407,quals408,quals409,quals410,quals411)
-    health_vars      <- Hmisc::Cs(health,discurr)
+    health_vars      <- Hmisc::Cs(health,discurr,
+                                  heal01,heal02,heal03,heal04,heal05,
+                                  heal06,heal07,heal08,heal09,heal10)
     work_vars        <- Hmisc::Cs(inecacr,grsswk,ftptwk,ttachr,ttushr,mpnr02,publicr,indm92m,indd92m,inds92m,soc2km,sc2kmmn,
                                   undemp,undhrs,ovhrs,lespay2)
     other_vars       <- Hmisc::Cs(refwkm,thiswv)
 
-    names <- c(demographic_vars,education_vars,work_vars, weights_vars,other_vars)
+    names <- c(demographic_vars,education_vars,work_vars,health_vars,weights_vars,other_vars)
     names <- tolower(names)
 
     data <- data[ ,names, with=F]
@@ -106,11 +108,13 @@ lfs_read_2005 <- function(
     education_vars   <- Hmisc::Cs(hiqual5d,hiqual5,btec,sctvec,gnvq4,nvqlev,rsa,candg,numol5,numal,numas,hst,advhst,
                                   gcseful1,gcseful2,gcseful3,gcseful4,qgcse41,qgcse42,qgcse43,qgcse44,
                                   quals401,quals402,quals403,quals404,quals405,quals406,quals407,quals408,quals409,quals410,quals411,qualch51)
-    health_vars      <- Hmisc::Cs(health,discurr)
+    health_vars      <- Hmisc::Cs(health,discurr,
+                                  heal01,heal02,heal03,heal04,heal05,
+                                  heal06,heal07,heal08,heal09,heal10)
     work_vars        <- Hmisc::Cs(inecac05,grsswk,ftptwk,ttachr,ttushr,mpnr02,publicr,indm92m,indd92m,inds92m,soc2km,sc2kmmn)
     other_vars       <- Hmisc::Cs(refwkm)
 
-    names <- c(demographic_vars,education_vars,work_vars, weights_vars,other_vars)
+    names <- c(demographic_vars,education_vars,work_vars,health_vars,weights_vars,other_vars)
     names <- tolower(names)
 
     data <- data[ ,names, with=F]
diff --git a/R/lfs_read_2006.R b/R/lfs_read_2006.R
index a2f5891..1bd0874 100644
--- a/R/lfs_read_2006.R
+++ b/R/lfs_read_2006.R
@@ -61,12 +61,14 @@ lfs_read_2006 <- function(
     education_vars   <- Hmisc::Cs(edage,hiqual5d,hiqual5,btec,sctvec,gnvq4,nvqlev,rsa,candg,numol5,numal,numas,hst,advhst,typhst1,typhst2,typhst3,typhst4,typhst5,
                                   gcseful1,gcseful2,gcseful3,gcseful4,qgcse41,qgcse42,qgcse43,qgcse44,
                                   quals601,quals602,quals603,quals604,quals605,quals606,quals607,quals608,quals609,quals610,quals611,qualch51)
-    health_vars      <- Hmisc::Cs(health,discurr)
+    health_vars      <- Hmisc::Cs(health,discurr,
+                                  heal01,heal02,heal03,heal04,heal05,
+                                  heal06,heal07,heal08,heal09,heal10)
     work_vars        <- Hmisc::Cs(inecac05,grsswk,ftptwk,ttachr,ttushr,mpnr02,publicr,indm92m,indd92m,inds92m,soc2km,sc2kmmn,
                                   undemp,undhrs,ovhrs,lespay2)
     other_vars       <- Hmisc::Cs(refwkm,thiswv)
 
-    names <- c(demographic_vars,education_vars,work_vars, weights_vars,other_vars)
+    names <- c(demographic_vars,education_vars,work_vars,health_vars,weights_vars,other_vars)
     names <- tolower(names)
 
     data <- data[ ,names, with=F]
diff --git a/R/lfs_read_2007.R b/R/lfs_read_2007.R
index 5e5409d..cd2db61 100644
--- a/R/lfs_read_2007.R
+++ b/R/lfs_read_2007.R
@@ -61,7 +61,9 @@ lfs_read_2007 <- function(
     education_vars   <- Hmisc::Cs(edage,hiqual5d,hiqual5,btec7,sctvec,gnvq4,nvqlev,rsa,candg,numol5,numal,numas,hst,advhst,typhst1,typhst2,typhst3,typhst4,typhst5,
                                   gcseful1,gcseful2,gcseful3,gcseful4,qgcse41,qgcse42,qgcse43,qgcse44,
                                   quals601,quals602,quals603,quals604,quals605,quals606,quals607,quals608,quals609,quals610,quals611,qualch51)
-    health_vars      <- Hmisc::Cs(health,discurr)
+    health_vars      <- Hmisc::Cs(health,discurr,
+                                  heal01,heal02,heal03,heal04,heal05,
+                                  heal06,heal07,heal08,heal09,heal10)
     work_vars        <- Hmisc::Cs(inecac05,grsswk,ftptwk,ttachr,ttushr,mpnr02,publicr,indm92m,indd92m,inds92m,soc2km,sc2kmmn,
                                   undemp,undhrs,ovhrs,lespay2)
     other_vars       <- Hmisc::Cs(refwkm,thiswv)
diff --git a/R/lfs_read_2008.R b/R/lfs_read_2008.R
index 285b671..dd6a678 100644
--- a/R/lfs_read_2008.R
+++ b/R/lfs_read_2008.R
@@ -61,7 +61,9 @@ lfs_read_2008 <- function(
     education_vars   <- Hmisc::Cs(edage,hiqual8d,hiqual8,btec8,sctvec,gnvq4,nvqlev,rsa,candg,numol5,numal,numas,hst,advhst,typhst1,typhst2,typhst3,typhst4,typhst5,
                                   gcseful1,gcseful2,gcseful3,gcseful4,qgcse41,qgcse42,qgcse43,qgcse44,
                                   quals601,quals602,quals603,quals604,quals605,quals606,quals607,quals608,quals609,quals610,quals611,qualch51)
-    health_vars      <- Hmisc::Cs(health,discurr)
+    health_vars      <- Hmisc::Cs(health,discurr,
+                                  heal01,heal02,heal03,heal04,heal05,
+                                  heal06,heal07,heal08,heal09,heal10)
     work_vars        <- Hmisc::Cs(inecac05,grsswk,ftptwk,ttachr,ttushr,mpnr02,publicr,indm92m,indd92m,inds92m,soc2km,sc2kmmn,
                                   undemp,undhrs,ovhrs,lespay2)
     other_vars       <- Hmisc::Cs(refwkm,thiswv)
diff --git a/R/lfs_read_2009.R b/R/lfs_read_2009.R
index 3702ddd..c4313a5 100644
--- a/R/lfs_read_2009.R
+++ b/R/lfs_read_2009.R
@@ -60,12 +60,14 @@ lfs_read_2009 <- function(
     education_vars   <- Hmisc::Cs(edage,hiqual8d,hiqual8,btec8,sctvec,gnvq4,nvqlev,rsa,candg,numol5,numal,numas,hst,advhst,typhst1,typhst2,typhst3,typhst4,typhst5,
                                   gcseful1,gcseful2,gcseful3,gcseful4,qgcse41,qgcse42,qgcse43,qgcse44,
                                   quals601,quals602,quals603,quals604,quals605,quals606,quals607,quals608,quals609,quals610,quals611,qualch51)
-    health_vars      <- Hmisc::Cs(health,discurr)
+    health_vars      <- Hmisc::Cs(health,discurr,
+                                  heal01,heal02,heal03,heal04,heal05,
+                                  heal06,heal07,heal08,heal09,heal10)
     work_vars        <- Hmisc::Cs(inecac05,grsswk,ftptwk,ttachr,ttushr,mpnr02,publicr,indc07m,indd07m,inds07m,soc2km,sc2kmmn,
                                   undemp,undhrs,ovhrs,lespay2)
     other_vars       <- Hmisc::Cs(refwkm,thiswv)
 
-    names <- c(demographic_vars,education_vars,work_vars, weights_vars,other_vars)
+    names <- c(demographic_vars,education_vars,work_vars,health_vars,weights_vars,other_vars)
     names <- tolower(names)
 
     data <- data[ ,names, with=F]
@@ -108,11 +110,13 @@ lfs_read_2009 <- function(
     demographic_vars <- Hmisc::Cs(age,sex,govtof,eth01,marsta,fdpch16)
     education_vars   <- Hmisc::Cs(hiqual8d,hiqual8,btec8,sctvec,gnvq4,nvqlev,rsa,candg,numol5,numal,numas,hst,advhst,
                                   quals601,quals602,quals603,quals604,quals605,quals606,quals607,quals608,quals609,quals610,quals611,qualch91)
-    health_vars      <- Hmisc::Cs(health,discurr)
+    health_vars      <- Hmisc::Cs(health,discurr,
+                                  heal01,heal02,heal03,heal04,heal05,
+                                  heal06,heal07,heal08,heal09,heal10)
     work_vars        <- Hmisc::Cs(inecac05,grsswk,ftptwk,ttachr,ttushr,mpnr02,publicr,indc07m,indd07m,inds07m,soc2km,sc2kmmn)
     other_vars       <- Hmisc::Cs(refwkm,thiswv)
 
-    names <- c(demographic_vars,education_vars,work_vars, weights_vars,other_vars)
+    names <- c(demographic_vars,education_vars,work_vars,health_vars,weights_vars,other_vars)
     names <- tolower(names)
 
     data <- data[ ,names, with=F]
diff --git a/R/lfs_read_2010.R b/R/lfs_read_2010.R
index 11a835a..99b96eb 100644
--- a/R/lfs_read_2010.R
+++ b/R/lfs_read_2010.R
@@ -62,7 +62,9 @@ lfs_read_2010 <- function(
                                   gcseful1,gcseful2,gcseful3,gcseful4,qgcse41,qgcse42,qgcse43,qgcse44,
                                   quls1001,quls1002,quls1003,quls1004,quls1005,quls1006,quls1007,quls1008,quls1009,quls1010,
                                   quls1011,quls1012,quls1013,quls1014,quls1015,quls1016,quls1017,quls1018,quls1019,quls1020,qualch91)
-    health_vars      <- Hmisc::Cs(health,discurr)
+    health_vars      <- Hmisc::Cs(health,discurr,
+                                  heal01,heal02,heal03,heal04,heal05,
+                                  heal06,heal07,heal08,heal09,heal10)
     work_vars        <- Hmisc::Cs(inecac05,grsswk,ftptwk,ttachr,ttushr,mpnr02,publicr,indc07m,indd07m,inds07m,soc2km,sc2kmmn,
                                   undemp,undhrs,ovhrs,lespay2)
     other_vars       <- Hmisc::Cs(refwkm,thiswv)
diff --git a/R/lfs_read_2011.R b/R/lfs_read_2011.R
index f0ef4d8..f1be454 100644
--- a/R/lfs_read_2011.R
+++ b/R/lfs_read_2011.R
@@ -62,7 +62,9 @@ lfs_read_2011 <- function(
                                   qual_1,qual_2,qual_3,qual_4,qual_5,qual_6,qual_7,qual_8,qual_9,qual_10,qual_11,qual_12,qual_13,qual_14,
                                   qual_15,qual_16,qual_17,qual_18,qual_19,qual_20,qual_21,qual_22,qual_23,qual_24,qual_25,qual_26,qual_27,
                                   qual_28,qual_29,qual_30,qual_31)
-    health_vars      <- Hmisc::Cs(health,discurr)
+    health_vars      <- Hmisc::Cs(health,discurr,
+                                  heal01,heal02,heal03,heal04,heal05,
+                                  heal06,heal07,heal08,heal09,heal10)
     work_vars        <- Hmisc::Cs(inecac05,grsswk,ftptwk,ttachr,ttushr,mpnr02,publicr,indc07m,indd07m,inds07m,soc10m,sc10mmn,
                                   undemp,undhrs,ovhrs,lespay2)
     other_vars       <- Hmisc::Cs(refwkm,thiswv)
@@ -117,7 +119,9 @@ lfs_read_2011 <- function(
                                   qual_1,qual_2,qual_3,qual_4,qual_5,qual_6,qual_7,qual_8,qual_9,qual_10,qual_11,qual_12,qual_13,qual_14,
                                   qual_15,qual_16,qual_17,qual_18,qual_19,qual_20,qual_21,qual_22,qual_23,qual_24,qual_25,qual_26,qual_27,
                                   qual_28,qual_29,qual_30,qual_31)
-    health_vars      <- Hmisc::Cs(health,discurr)
+    health_vars      <- Hmisc::Cs(health,discurr,
+                                  heal01,heal02,heal03,heal04,heal05,
+                                  heal06,heal07,heal08,heal09,heal10)
     work_vars        <- Hmisc::Cs(inecac05,grsswk,ftptwk,ttachr,ttushr,mpnr02,publicr,indc07m,indd07m,inds07m,soc10m,sc10mmn,
                                   undemp,undhrs,ovhrs,lespay2)
     other_vars       <- Hmisc::Cs(refwkm)
diff --git a/R/lfs_read_2012.R b/R/lfs_read_2012.R
index 2f65b84..ce78b94 100644
--- a/R/lfs_read_2012.R
+++ b/R/lfs_read_2012.R
@@ -62,7 +62,9 @@ lfs_read_2012 <- function(
                                   qual_1,qual_2,qual_3,qual_4,qual_5,qual_6,qual_7,qual_8,qual_9,qual_10,qual_11,qual_12,qual_13,qual_14,
                                   qual_15,qual_16,qual_17,qual_18,qual_19,qual_20,qual_21,qual_22,qual_23,qual_24,qual_25,qual_26,qual_27,
                                   qual_28,qual_29,qual_30,qual_31)
-    health_vars      <- Hmisc::Cs(health,discurr)
+    health_vars      <- Hmisc::Cs(health,discurr,
+                                  heal01,heal02,heal03,heal04,heal05,
+                                  heal06,heal07,heal08,heal09,heal10)
     work_vars        <- Hmisc::Cs(inecac05,grsswk,ftptwk,ttachr,ttushr,mpnr02,publicr,indc07m,indd07m,inds07m,soc10m,sc10mmn,
                                   undemp,undhrs,ovhrs,lespay2)
     other_vars       <- Hmisc::Cs(refwkm,thiswv)
diff --git a/R/lfs_read_2013.R b/R/lfs_read_2013.R
index bba8b7b..17677ed 100644
--- a/R/lfs_read_2013.R
+++ b/R/lfs_read_2013.R
@@ -63,7 +63,9 @@ lfs_read_2013 <- function(
                                   qual_1,qual_2,qual_3,qual_4,qual_5,qual_6,qual_7,qual_8,qual_9,qual_10,qual_11,qual_12,qual_13,qual_14,
                                   qual_15,qual_16,qual_17,qual_18,qual_19,qual_20,qual_21,qual_22,qual_23,qual_24,qual_25,qual_26,qual_27,
                                   qual_28,qual_29,qual_30,qual_31)
-    health_vars      <- Hmisc::Cs(health,discurr)
+    health_vars      <- Hmisc::Cs(health,discurr,
+                                  heal01,heal02,heal03,heal04,heal05,
+                                  heal06,heal07,heal08,heal09,heal10)
     work_vars        <- Hmisc::Cs(inecac05,grsswk,ftptwk,ttachr,ttushr,mpnr02,publicr,indc07m,indd07m,inds07m,soc10m,sc10mmn,
                                   undemp,undhrs,ovhrs,lespay2)
     other_vars       <- Hmisc::Cs(refwkm,thiswv)
@@ -101,7 +103,9 @@ lfs_read_2013 <- function(
                                   qual_1,qual_2,qual_3,qual_4,qual_5,qual_6,qual_7,qual_8,qual_9,qual_10,qual_11,qual_12,qual_13,qual_14,
                                   qual_15,qual_16,qual_17,qual_18,qual_19,qual_20,qual_21,qual_22,qual_23,qual_24,qual_25,qual_26,qual_27,
                                   qual_28,qual_29,qual_30,qual_31)
-    health_vars      <- Hmisc::Cs(health,discurr13,disea)
+    health_vars      <- Hmisc::Cs(health,discurr13,disea,
+                                  heal01,heal02,heal03,heal04,heal05,
+                                  heal06,heal07,heal08,heal09,heal10)
     work_vars        <- Hmisc::Cs(inecac05,grsswk,ftptwk,ttachr,ttushr,mpnr02,publicr,indc07m,indd07m,inds07m,soc10m,sc10mmn,
                                   undemp,undhrs,ovhrs,lespay2)
     other_vars       <- Hmisc::Cs(refwkm,thiswv)
diff --git a/R/lfs_read_2014.R b/R/lfs_read_2014.R
index 84aa1bc..72c560a 100644
--- a/R/lfs_read_2014.R
+++ b/R/lfs_read_2014.R
@@ -62,7 +62,9 @@ lfs_read_2014 <- function(
                                   qual_1,qual_2,qual_3,qual_4,qual_5,qual_6,qual_7,qual_8,qual_9,qual_10,qual_11,qual_12,qual_13,qual_14,
                                   qual_15,qual_16,qual_17,qual_18,qual_19,qual_20,qual_21,qual_22,qual_23,qual_24,qual_25,qual_26,qual_27,
                                   qual_28,qual_29,qual_30,qual_31)
-    health_vars      <- Hmisc::Cs(health,discurr13,disea)
+    health_vars      <- Hmisc::Cs(health,discurr13,disea,
+                                  heal01,heal02,heal03,heal04,heal05,
+                                  heal06,heal07,heal08,heal09,heal10)
     work_vars        <- Hmisc::Cs(inecac05,grsswk,ftptwk,ttachr,ttushr,mpnr02,publicr,indc07m,indd07m,inds07m,soc10m,sc10mmn,
                                   undemp,undhrs,ovhrs,lespay2)
     other_vars       <- Hmisc::Cs(refwkm,thiswv)
diff --git a/R/lfs_read_2015.R b/R/lfs_read_2015.R
index abcbab2..fa27b26 100644
--- a/R/lfs_read_2015.R
+++ b/R/lfs_read_2015.R
@@ -62,7 +62,9 @@ lfs_read_2015 <- function(
                                   qual_1,qual_2,qual_3,qual_4,qual_5,qual_6,qual_7,qual_8,qual_9,qual_10,qual_11,qual_12,qual_13,qual_14,
                                   qual_15,qual_16,qual_17,qual_18,qual_19,qual_20,qual_21,qual_22,qual_23,qual_24,qual_25,qual_26,qual_27,
                                   qual_28,qual_29,qual_30,qual_31)
-    health_vars      <- Hmisc::Cs(health,discurr13,disea)
+    health_vars      <- Hmisc::Cs(health,discurr13,disea,
+                                  heal01,heal02,heal03,heal04,heal05,
+                                  heal06,heal07,heal08,heal09,heal10)
     work_vars        <- Hmisc::Cs(inecac05,grsswk,ftptwk,ttachr,ttushr,mpnr02,publicr,indc07m,indd07m,inds07m,soc10m,sc10mmn,
                                   undemp,undhrs,ovhrs,lespay2)
     benefit_vars     <- Hmisc::Cs(ooben)
diff --git a/R/lfs_read_2016.R b/R/lfs_read_2016.R
index af46828..e040ea1 100644
--- a/R/lfs_read_2016.R
+++ b/R/lfs_read_2016.R
@@ -62,7 +62,9 @@ lfs_read_2016 <- function(
                                   qual_1,qual_2,qual_3,qual_4,qual_5,qual_6,qual_7,qual_8,qual_9,qual_10,qual_11,qual_12,qual_13,qual_14,
                                   qual_15,qual_16,qual_17,qual_18,qual_19,qual_20,qual_21,qual_22,qual_23,qual_24,qual_25,qual_26,qual_27,
                                   qual_28,qual_29,qual_30,qual_31)
-    health_vars      <- Hmisc::Cs(health,discurr13,disea)
+    health_vars      <- Hmisc::Cs(health,discurr13,disea,
+                                  heal01,heal02,heal03,heal04,heal05,
+                                  heal06,heal07,heal08,heal09,heal10)
     work_vars        <- Hmisc::Cs(inecac05,grsswk,ftptwk,ttachr,ttushr,mpnr02,publicr,indc07m,indd07m,inds07m,soc10m,sc10mmn,
                                   undemp,undhrs,ovhrs,lespay2)
     benefit_vars     <- Hmisc::Cs(ooben)
diff --git a/R/lfs_read_2017.R b/R/lfs_read_2017.R
index bfeba88..c93aac5 100644
--- a/R/lfs_read_2017.R
+++ b/R/lfs_read_2017.R
@@ -62,7 +62,9 @@ lfs_read_2017 <- function(
                                   qual_1,qual_2,qual_3,qual_4,qual_5,qual_6,qual_7,qual_8,qual_9,qual_10,qual_11,qual_12,qual_13,qual_14,
                                   qual_15,qual_16,qual_17,qual_18,qual_19,qual_20,qual_21,qual_22,qual_23,qual_24,qual_25,qual_26,qual_27,
                                   qual_28,qual_29,qual_30,qual_31)
-    health_vars      <- Hmisc::Cs(health,discurr13,disea)
+    health_vars      <- Hmisc::Cs(health,discurr13,disea,
+                                  heal01,heal02,heal03,heal04,heal05,
+                                  heal06,heal07,heal08,heal09,heal10)
     work_vars        <- Hmisc::Cs(inecac05,grsswk,ftptwk,ttachr,ttushr,mpnr02,publicr,indc07m,indd07m,inds07m,soc10m,sc10mmn,
                                   undemp,undhrs,ovhrs,lespay2)
     benefit_vars     <- Hmisc::Cs(ooben)
diff --git a/R/lfs_read_2018.R b/R/lfs_read_2018.R
index 4a2fb74..d0acf3e 100644
--- a/R/lfs_read_2018.R
+++ b/R/lfs_read_2018.R
@@ -62,7 +62,9 @@ lfs_read_2018 <- function(
                                   qual_1,qual_2,qual_3,qual_4,qual_5,qual_6,qual_7,qual_8,qual_9,qual_10,qual_11,qual_12,qual_13,qual_14,
                                   qual_15,qual_16,qual_17,qual_18,qual_19,qual_20,qual_21,qual_22,qual_23,qual_24,qual_25,qual_26,qual_27,
                                   qual_28,qual_29,qual_30,qual_31)
-    health_vars      <- Hmisc::Cs(health,discurr13,disea)
+    health_vars      <- Hmisc::Cs(health,discurr13,disea,
+                                  heal01,heal02,heal03,heal04,heal05,
+                                  heal06,heal07,heal08,heal09,heal10)
     work_vars        <- Hmisc::Cs(inecac05,grsswk,ftptwk,ttachr,ttushr,mpnr02,publicr,indc07m,indd07m,inds07m,soc10m,sc10mmn,
                                   undemp,undhrs,ovhrs,lespay2)
     benefit_vars     <- Hmisc::Cs(ooben)
diff --git a/R/lfs_read_2019.R b/R/lfs_read_2019.R
index 5bfba76..05759c5 100644
--- a/R/lfs_read_2019.R
+++ b/R/lfs_read_2019.R
@@ -62,7 +62,9 @@ lfs_read_2019 <- function(
                                   qual_1,qual_2,qual_3,qual_4,qual_5,qual_6,qual_7,qual_8,qual_9,qual_10,qual_11,qual_12,qual_13,qual_14,
                                   qual_15,qual_16,qual_17,qual_18,qual_19,qual_20,qual_21,qual_22,qual_23,qual_24,qual_25,qual_26,qual_27,
                                   qual_28,qual_29,qual_30,qual_31)
-    health_vars      <- Hmisc::Cs(health,discurr13,disea)
+    health_vars      <- Hmisc::Cs(health,discurr13,disea,
+                                  heal01,heal02,heal03,heal04,heal05,
+                                  heal06,heal07,heal08,heal09,heal10)
     work_vars        <- Hmisc::Cs(inecac05,grsswk,ftptwk,ttachr,ttushr,mpnr02,publicr,indc07m,indd07m,inds07m,soc10m,sc10mmn,
                                   undemp,undhrs,ovhrs,lespay2)
     benefit_vars     <- Hmisc::Cs(ooben)
diff --git a/R/lfs_read_2020.R b/R/lfs_read_2020.R
index 665a8f2..b31fb04 100644
--- a/R/lfs_read_2020.R
+++ b/R/lfs_read_2020.R
@@ -62,7 +62,9 @@ lfs_read_2020 <- function(
                                   qual_1,qual_2,qual_3,qual_4,qual_5,qual_6,qual_7,qual_8,qual_9,qual_10,qual_11,qual_12,qual_13,qual_14,
                                   qual_15,qual_16,qual_17,qual_18,qual_19,qual_20,qual_21,qual_22,qual_23,qual_24,qual_25,qual_26,qual_27,
                                   qual_28,qual_29,qual_30,qual_31)
-    health_vars      <- Hmisc::Cs(health20,discurr20,disea)
+    health_vars      <- Hmisc::Cs(health20,discurr20,disea,
+                                  heal2001,heal2002,heal2003,heal2004,heal2005,
+                                  heal2006,heal2007,heal2008,heal2009,heal2010)
     work_vars        <- Hmisc::Cs(inecac05,grsswk,ftptwk,ttachr,ttushr,mpnr02,publicr,indc07m,indd07m,inds07m,soc10m,sc10mmn)
     benefit_vars     <- Hmisc::Cs(ooben)
     other_vars       <- Hmisc::Cs(refwkm,thiswv)
@@ -89,7 +91,13 @@ lfs_read_2020 <- function(
                            "qf11","qf12","qf13","qf14", "qf15","qf16","qf17","qf18","qf19","qf20","qf21","qf22",
                            "qf23","qf24","qf25","qf26","qf27","qf28","qf29","qf30","qf31") )
 
-    data.table::setnames(data, c("health20","discurr20"), c("health","discurr"))
+    data.table::setnames(data,
+                         c("health20","discurr20",
+                           "heal2001","heal2002","heal2003","heal2004","heal2005",
+                           "heal2006","heal2007","heal2008","heal2009","heal2010"),
+                         c("health","discurr",
+                           "heal01","heal02","heal03","heal04","heal05",
+                           "heal06","heal07","heal08","heal09","heal10"))
 
     #preliminary cleaning of the vocational qualification variables
     data[btec   %in% c(5,6), btec := NA]
diff --git a/R/lfs_read_2021.R b/R/lfs_read_2021.R
index 34fa563..75b0dd0 100644
--- a/R/lfs_read_2021.R
+++ b/R/lfs_read_2021.R
@@ -62,7 +62,9 @@ lfs_read_2021 <- function(
                                   qual_1,qual_2,qual_3,qual_4,qual_5,qual_6,qual_7,qual_8,qual_9,qual_10,qual_11,qual_12,qual_13,qual_14,
                                   qual_15,qual_16,qual_17,qual_18,qual_19,qual_20,qual_21,qual_22,qual_23,qual_24,qual_25,qual_26,qual_27,
                                   qual_28,qual_29,qual_30,qual_31)
-    health_vars      <- Hmisc::Cs(health20,discurr20,disea)
+    health_vars      <- Hmisc::Cs(health20,discurr20,disea,
+                                  heal2001,heal2002,heal2003,heal2004,heal2005,
+                                  heal2006,heal2007,heal2008,heal2009,heal2010)
     work_vars        <- Hmisc::Cs(inecac05,grsswk,ftptwk,ttachr,ttushr,mpnr02,publicr,indc07m,indd07m,inds07m,soc20m,sc20mmn,
                                   undemp,undhrs,ovhrs,lespay2)
     benefit_vars     <- Hmisc::Cs(ooben)
@@ -90,7 +92,13 @@ lfs_read_2021 <- function(
                            "qf11","qf12","qf13","qf14", "qf15","qf16","qf17","qf18","qf19","qf20","qf21","qf22",
                            "qf23","qf24","qf25","qf26","qf27","qf28","qf29","qf30","qf31") )
 
-    data.table::setnames(data, c("health20","discurr20"), c("health","discurr"))
+    data.table::setnames(data,
+                         c("health20","discurr20",
+                           "heal2001","heal2002","heal2003","heal2004","heal2005",
+                           "heal2006","heal2007","heal2008","heal2009","heal2010"),
+                         c("health","discurr",
+                           "heal01","heal02","heal03","heal04","heal05",
+                           "heal06","heal07","heal08","heal09","heal10"))
 
     #preliminary cleaning of the vocational qualification variables
     data[btec   %in% c(5,6), btec := NA]
diff --git a/R/lfs_read_2022.R b/R/lfs_read_2022.R
index e845dc1..bc01808 100644
--- a/R/lfs_read_2022.R
+++ b/R/lfs_read_2022.R
@@ -58,7 +58,9 @@ lfs_read_2022 <- function(
     weights_vars     <- Hmisc::Cs(pwt22,piwt22)
     demographic_vars <- Hmisc::Cs(age,sex,gor9d,ethukeul,marsta,fdpch16)
     education_vars   <- Hmisc::Cs(edage,hiqul22d,hiqual22)
-    health_vars      <- Hmisc::Cs(health20,discurr20,disea)
+    health_vars      <- Hmisc::Cs(health20,discurr20,disea,
+                                  heal2001,heal2002,heal2003,heal2004,heal2005,
+                                  heal2006,heal2007,heal2008,heal2009,heal2010)
     work_vars        <- Hmisc::Cs(inecac05,grsswk,ftptwk,ttachr,ttushr,mpnr02,publicr,indc07m,indd07m,inds07m,soc20m,sc20mmn,
                                   undemp,undhrs,ovhrs,lespay2)
     benefit_vars     <- Hmisc::Cs(ooben)
@@ -81,8 +83,13 @@ lfs_read_2022 <- function(
                          c("refwkm", "pwt22","piwt22"),
                          c("month", "pwt", "piwt") )
 
-    data.table::setnames(data, c("health20","discurr20"), c("health","discurr"))
-
+    data.table::setnames(data,
+                         c("health20","discurr20",
+                           "heal2001","heal2002","heal2003","heal2004","heal2005",
+                           "heal2006","heal2007","heal2008","heal2009","heal2010"),
+                         c("health","discurr",
+                           "heal01","heal02","heal03","heal04","heal05",
+                           "heal06","heal07","heal08","heal09","heal10"))
 
     clean.data.list[[l]] <- data
   }
diff --git a/R/lfs_read_2023.R b/R/lfs_read_2023.R
index c7fce24..3fbbd80 100644
--- a/R/lfs_read_2023.R
+++ b/R/lfs_read_2023.R
@@ -58,7 +58,9 @@ lfs_read_2023 <- function(
     weights_vars     <- Hmisc::Cs(pwt22,piwt22)
     demographic_vars <- Hmisc::Cs(age,sex,gor9d,ethukeul,marsta,fdpch16)
     education_vars   <- Hmisc::Cs(edage,hiqul22d,hiqual22)
-    health_vars      <- Hmisc::Cs(health20,discurr20,disea)
+    health_vars      <- Hmisc::Cs(health20,discurr20,disea,
+                                  heal2001,heal2002,heal2003,heal2004,heal2005,
+                                  heal2006,heal2007,heal2008,heal2009,heal2010)
     work_vars        <- Hmisc::Cs(inecac05,grsswk,ftptwk,ttachr,ttushr,mpnr02,publicr,indc07m,indd07m,inds07m,soc20m,sc20mmn,
                                   undemp,undhrs,ovhrs,lespay2)
     benefit_vars     <- Hmisc::Cs(ooben)
@@ -81,7 +83,13 @@ lfs_read_2023 <- function(
                          c("refwkm", "pwt22","piwt22"),
                          c("month", "pwt", "piwt") )
 
-    data.table::setnames(data, c("health20","discurr20"), c("health","discurr"))
+    data.table::setnames(data,
+                         c("health20","discurr20",
+                           "heal2001","heal2002","heal2003","heal2004","heal2005",
+                           "heal2006","heal2007","heal2008","heal2009","heal2010"),
+                         c("health","discurr",
+                           "heal01","heal02","heal03","heal04","heal05",
+                           "heal06","heal07","heal08","heal09","heal10"))
 
     clean.data.list[[l]] <- data
   }
@@ -95,7 +103,9 @@ lfs_read_2023 <- function(
     weights_vars     <- Hmisc::Cs(pwt23,piwt22)
     demographic_vars <- Hmisc::Cs(age,sex,gor9d,ethukeul,marsta,fdpch16)
     education_vars   <- Hmisc::Cs(edage,hiqul22d,hiqual22)
-    health_vars      <- Hmisc::Cs(health20,discurr20,disea)
+    health_vars      <- Hmisc::Cs(health20,discurr20,disea,
+                                  heal2001,heal2002,heal2003,heal2004,heal2005,
+                                  heal2006,heal2007,heal2008,heal2009,heal2010)
     work_vars        <- Hmisc::Cs(inecac05,grsswk,ftptwk,ttachr,ttushr,mpnr02,publicr,indc07m,indd07m,inds07m,soc20m,sc20mmn,
                                   undemp,undhrs,ovhrs,lespay2)
     benefit_vars     <- Hmisc::Cs(ooben)
@@ -118,6 +128,13 @@ lfs_read_2023 <- function(
                          c("refwkm", "pwt23","piwt22"),
                          c("month", "pwt", "piwt") )
 
+    data.table::setnames(data,
+                         c("health20","discurr20",
+                           "heal2001","heal2002","heal2003","heal2004","heal2005",
+                           "heal2006","heal2007","heal2008","heal2009","heal2010"),
+                         c("health","discurr",
+                           "heal01","heal02","heal03","heal04","heal05",
+                           "heal06","heal07","heal08","heal09","heal10"))
 
     clean.data.list[[l]] <- data
   }
@@ -134,7 +151,9 @@ lfs_read_2023 <- function(
     weights_vars     <- Hmisc::Cs(pwt23,piwt23)
     demographic_vars <- Hmisc::Cs(age,sex,gor9d,ethukeul,marsta,fdpch16)
     education_vars   <- Hmisc::Cs(edage,hiqul22d,hiqual22)
-    health_vars      <- Hmisc::Cs(health20,discurr20,disea)
+    health_vars      <- Hmisc::Cs(health20,discurr20,disea,
+                                  heal2001,heal2002,heal2003,heal2004,heal2005,
+                                  heal2006,heal2007,heal2008,heal2009,heal2010)
     work_vars        <- Hmisc::Cs(inecac05,grsswk,ftptwk,ttachr,ttushr,mpnr02,publicr,indc07m,indd07m,inds07m,soc20m,sc20mmn,
                                   undemp,undhrs,ovhrs,lespay2)
     benefit_vars     <- Hmisc::Cs(ooben)
@@ -157,6 +176,14 @@ lfs_read_2023 <- function(
                          c("refwkm", "pwt23","piwt23"),
                          c("month", "pwt", "piwt") )
 
+    data.table::setnames(data,
+                         c("health20","discurr20",
+                           "heal2001","heal2002","heal2003","heal2004","heal2005",
+                           "heal2006","heal2007","heal2008","heal2009","heal2010"),
+                         c("health","discurr",
+                           "heal01","heal02","heal03","heal04","heal05",
+                           "heal06","heal07","heal08","heal09","heal10"))
+
 
     clean.data.list[[l]] <- data
   }