add married / non-married binary variable

STAPM · Mar 12, 2024 · 0de3611 · 0de3611
1 parent 88bc401
commit 0de3611
Show file tree

Hide file tree

Showing 18 changed files with 37 additions and 25 deletions.
diff --git a/R/lfs_clean_demographic.R b/R/lfs_clean_demographic.R
@@ -145,6 +145,18 @@ lfs_clean_demographic <- function(
 
   }
 
+  # married or not binary
+
+  if ("mardy6" %in% colnames(data)) {
+    data[mardy6 == 1, married_cohabiting := 1]
+    data[mardy6 == 2, married_cohabiting := 0]
+
+  } else {
+
+    data[, married_cohabiting := NA]
+
+  }
+
   # dependent children in the family
 
   data[is.na(fdpch16) | fdpch16 == 0, dep_children := "no"]
@@ -164,10 +176,10 @@ lfs_clean_demographic <- function(
 
   final_data <- data[, c("obs_id",
                          "age", "sex", "ethnicity_2cat", "ethnicity_4cat",
-                         "region", "country", "marstat", "dep_children")]
+                         "region", "country", "marstat", "married_cohabiting", "dep_children")]
 
   var_names <- c("age", "sex", "ethnicity_2cat", "ethnicity_4cat",
-                 "region", "country", "marstat", "dep_children")
+                 "region", "country", "marstat", "married_cohabiting", "dep_children")
 
   setnames(final_data, var_names, paste0("d_", var_names))
 

diff --git a/R/lfs_read_2007.R b/R/lfs_read_2007.R
@@ -57,7 +57,7 @@ lfs_read_2007 <- function(
     setnames(data, names(data), tolower(names(data)))
 
     weights_vars     <- Hmisc::Cs(pwt14,piwt14)
-    demographic_vars <- Hmisc::Cs(age,sex,govtof,eth01,marsta,fdpch16)
+    demographic_vars <- Hmisc::Cs(age,sex,govtof,eth01,marsta,mardy6,fdpch16)
     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)
@@ -68,7 +68,7 @@ lfs_read_2007 <- function(
                                   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_2008.R b/R/lfs_read_2008.R
@@ -57,7 +57,7 @@ lfs_read_2008 <- function(
     setnames(data, names(data), tolower(names(data)))
 
     weights_vars     <- Hmisc::Cs(pwt14,piwt14)
-    demographic_vars <- Hmisc::Cs(age,sex,govtof,eth01,marsta,fdpch16)
+    demographic_vars <- Hmisc::Cs(age,sex,govtof,eth01,marsta,mardy6,fdpch16)
     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)
@@ -68,7 +68,7 @@ lfs_read_2008 <- function(
                                   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_2009.R b/R/lfs_read_2009.R
@@ -56,7 +56,7 @@ lfs_read_2009 <- function(
     setnames(data, names(data), tolower(names(data)))
 
     weights_vars     <- Hmisc::Cs(pwt14,piwt14)
-    demographic_vars <- Hmisc::Cs(age,sex,govtof,eth01,marsta)
+    demographic_vars <- Hmisc::Cs(age,sex,govtof,eth01,marsta,mardy6,fdpch16)
     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)
@@ -107,7 +107,7 @@ lfs_read_2009 <- function(
     setnames(data, names(data), tolower(names(data)))
 
     weights_vars     <- Hmisc::Cs(pwt14,piwt14)
-    demographic_vars <- Hmisc::Cs(age,sex,govtof,eth01,marsta,fdpch16)
+    demographic_vars <- Hmisc::Cs(age,sex,govtof,eth01,marsta,mardy6,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,

diff --git a/R/lfs_read_2010.R b/R/lfs_read_2010.R
@@ -57,7 +57,7 @@ lfs_read_2010 <- function(
     setnames(data, names(data), tolower(names(data)))
 
     weights_vars     <- Hmisc::Cs(pwt14,piwt14)
-    demographic_vars <- Hmisc::Cs(age,sex,govtof,eth01,marsta,fdpch16)
+    demographic_vars <- Hmisc::Cs(age,sex,govtof,eth01,marsta,mardy6,fdpch16)
     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,
                                   quls1001,quls1002,quls1003,quls1004,quls1005,quls1006,quls1007,quls1008,quls1009,quls1010,

diff --git a/R/lfs_read_2011.R b/R/lfs_read_2011.R
@@ -172,7 +172,7 @@ lfs_read_2011 <- function(
     setnames(data, names(data), tolower(names(data)))
 
     weights_vars     <- Hmisc::Cs(pwt18,piwt18)
-    demographic_vars <- Hmisc::Cs(age,sex,govtof,ethukeul,marsta,fdpch16)
+    demographic_vars <- Hmisc::Cs(age,sex,govtof,ethukeul,marsta,mardy6,fdpch16)
     education_vars   <- Hmisc::Cs(hiqul11d,hiqual11,bte11,sctvc11,gnvq11,nvq11,rsa11,cag11,numol5,numal,numas,hst,advhst,typhst1,typhst2,typhst3,typhst4,typhst5,
                                   gcseful1,gcseful2,gcseful3,gcseful4,qgcse41,qgcse42,qgcse43,qgcse44,
                                   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,

diff --git a/R/lfs_read_2012.R b/R/lfs_read_2012.R
@@ -56,7 +56,7 @@ lfs_read_2012 <- function(
     setnames(data, names(data), tolower(names(data)))
 
     weights_vars     <- Hmisc::Cs(pwt18,piwt18)
-    demographic_vars <- Hmisc::Cs(age,sex,govtof,ethukeul,marsta,fdpch16)
+    demographic_vars <- Hmisc::Cs(age,sex,govtof,ethukeul,marsta,mardy6,fdpch16)
     education_vars   <- Hmisc::Cs(edage,hiqul11d,hiqual11,bte11,sctvc11,gnvq11,nvq11,rsa11,cag11,numol5,numal,numas,hst,advhst,typhst1,typhst2,typhst3,typhst4,typhst5,
                                   gcseful1,gcseful2,gcseful3,gcseful4,qgcse41,qgcse42,qgcse43,qgcse44,
                                   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,

diff --git a/R/lfs_read_2013.R b/R/lfs_read_2013.R
@@ -57,7 +57,7 @@ lfs_read_2013 <- function(
     setnames(data, names(data), tolower(names(data)))
 
     weights_vars     <- Hmisc::Cs(pwt18,piwt18)
-    demographic_vars <- Hmisc::Cs(age,sex,govtof,ethukeul,marsta,fdpch16)
+    demographic_vars <- Hmisc::Cs(age,sex,govtof,ethukeul,marsta,mardy6,fdpch16)
     education_vars   <- Hmisc::Cs(edage,hiqul11d,hiqual11,bte11,sctvc11,gnvq11,nvq11,rsa11,cag11,numol5,numal,numas,hst,advhst,typhst1,typhst2,typhst3,typhst4,typhst5,
                                   gcseful1,gcseful2,gcseful3,gcseful4,qgcse41,qgcse42,qgcse43,qgcse44,
                                   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,
@@ -97,7 +97,7 @@ lfs_read_2013 <- function(
     setnames(data, names(data), tolower(names(data)))
 
     weights_vars     <- Hmisc::Cs(pwt18,piwt18)
-    demographic_vars <- Hmisc::Cs(age,sex,govtof,ethukeul,marsta,fdpch16)
+    demographic_vars <- Hmisc::Cs(age,sex,govtof,ethukeul,marsta,mardy6,fdpch16)
     education_vars   <- Hmisc::Cs(edage,hiqul11d,hiqual11,bte11,sctvc11,gnvq11,nvq11,rsa11,cag11,numol5,numal,numas,hst,advhst,typhst1,typhst2,typhst3,typhst4,typhst5,
                                   gcseful1,gcseful2,gcseful3,gcseful4,qgcse41,qgcse42,qgcse43,qgcse44,
                                   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,

diff --git a/R/lfs_read_2014.R b/R/lfs_read_2014.R
@@ -56,7 +56,7 @@ lfs_read_2014 <- function(
     setnames(data, names(data), tolower(names(data)))
 
     weights_vars     <- Hmisc::Cs(pwt18,piwt18)
-    demographic_vars <- Hmisc::Cs(age,sex,gor9d,ethukeul,marsta,fdpch16)
+    demographic_vars <- Hmisc::Cs(age,sex,gor9d,ethukeul,marsta,mardy6,fdpch16)
     education_vars   <- Hmisc::Cs(edage,hiqul11d,hiqual11,bte11,sctvc11,gnvq11,nvq11,rsa11,cag11,numol5,numal,numas,hst,advhst,typhst1,typhst2,typhst3,typhst4,typhst5,
                                   gcseful1,gcseful2,gcseful3,gcseful4,qgcse41,qgcse42,qgcse43,qgcse44,
                                   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,

diff --git a/R/lfs_read_2015.R b/R/lfs_read_2015.R
@@ -56,7 +56,7 @@ lfs_read_2015 <- function(
     setnames(data, names(data), tolower(names(data)))
 
     weights_vars     <- Hmisc::Cs(pwt18,piwt18)
-    demographic_vars <- Hmisc::Cs(age,sex,gor9d,ethukeul,marsta,fdpch16)
+    demographic_vars <- Hmisc::Cs(age,sex,gor9d,ethukeul,marsta,mardy6,fdpch16)
     education_vars   <- Hmisc::Cs(edage,hiqul15d,hiqual15,bte11,sctvc11,gnvq11,nvq11,rsa11,cag11,numol5,numal,numas,hst,advhst,typhst1,typhst2,typhst3,typhst4,typhst5,
                                   gcseful1,gcseful2,gcseful3,gcseful4,qgcse41,qgcse42,qgcse43,qgcse44,
                                   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,

diff --git a/R/lfs_read_2016.R b/R/lfs_read_2016.R
@@ -56,7 +56,7 @@ lfs_read_2016 <- function(
     setnames(data, names(data), tolower(names(data)))
 
     weights_vars     <- Hmisc::Cs(pwt18,piwt18)
-    demographic_vars <- Hmisc::Cs(age,sex,gor9d,ethukeul,marsta,fdpch16)
+    demographic_vars <- Hmisc::Cs(age,sex,gor9d,ethukeul,marsta,mardy6,fdpch16)
     education_vars   <- Hmisc::Cs(edage,hiqul15d,hiqual15,bte11,sctvc11,gnvq11,nvq11,rsa11,cag11,numol5,numal,numas,hst,advhst,typhst1,typhst2,typhst3,typhst4,typhst5,
                                   gcseful1,gcseful2,gcseful3,gcseful4,qgcse41,qgcse42,qgcse43,qgcse44,
                                   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,

diff --git a/R/lfs_read_2017.R b/R/lfs_read_2017.R
@@ -56,7 +56,7 @@ lfs_read_2017 <- function(
     setnames(data, names(data), tolower(names(data)))
 
     weights_vars     <- Hmisc::Cs(pwt18,piwt18)
-    demographic_vars <- Hmisc::Cs(age,sex,gor9d,ethukeul,marsta,fdpch16)
+    demographic_vars <- Hmisc::Cs(age,sex,gor9d,ethukeul,marsta,mardy6,fdpch16)
     education_vars   <- Hmisc::Cs(edage,hiqul15d,hiqual15,bte11,sctvc11,gnvq11,nvq11,rsa11,cag11,numol5,numal,numas,hst,advhst,typhst1,typhst2,typhst3,typhst4,typhst5,
                                   gcseful1,gcseful2,gcseful3,gcseful4,qgcse41,qgcse42,qgcse43,qgcse44,
                                   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,

diff --git a/R/lfs_read_2018.R b/R/lfs_read_2018.R
@@ -56,7 +56,7 @@ lfs_read_2018 <- function(
     setnames(data, names(data), tolower(names(data)))
 
     weights_vars     <- Hmisc::Cs(pwt18,piwt18)
-    demographic_vars <- Hmisc::Cs(age,sex,gor9d,ethukeul,marsta,fdpch16)
+    demographic_vars <- Hmisc::Cs(age,sex,gor9d,ethukeul,marsta,mardy6,fdpch16)
     education_vars   <- Hmisc::Cs(edage,hiqul15d,hiqual15,bte11,sctvc11,gnvq11,nvq11,rsa11,cag11,numol5,numal,numas,hst,advhst,typhst1,typhst2,typhst3,typhst4,typhst5,
                                   gcseful1,gcseful2,gcseful3,gcseful4,qgcse41,qgcse42,qgcse43,qgcse44,
                                   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,

diff --git a/R/lfs_read_2019.R b/R/lfs_read_2019.R
@@ -56,7 +56,7 @@ lfs_read_2019 <- function(
     setnames(data, names(data), tolower(names(data)))
 
     weights_vars     <- Hmisc::Cs(pwt18,piwt18)
-    demographic_vars <- Hmisc::Cs(age,sex,gor9d,ethukeul,marsta,fdpch16)
+    demographic_vars <- Hmisc::Cs(age,sex,gor9d,ethukeul,marsta,mardy6,fdpch16)
     education_vars   <- Hmisc::Cs(edage,hiqul15d,hiqual15,bte11,sctvc11,gnvq11,nvq11,rsa11,cag11,numol5,numal,numas,hst,advhst,typhst1,typhst2,typhst3,typhst4,typhst5,
                                   gcseful1,gcseful2,gcseful3,gcseful4,qgcse41,qgcse42,qgcse43,qgcse44,
                                   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,

diff --git a/R/lfs_read_2020.R b/R/lfs_read_2020.R
@@ -56,7 +56,7 @@ lfs_read_2020 <- function(
     setnames(data, names(data), tolower(names(data)))
 
     weights_vars     <- Hmisc::Cs(pwt22,piwt22)
-    demographic_vars <- Hmisc::Cs(age,sex,gor9d,ethukeul,marsta,fdpch16)
+    demographic_vars <- Hmisc::Cs(age,sex,gor9d,ethukeul,marsta,mardy6,fdpch16)
     education_vars   <- Hmisc::Cs(edage,hiqul15d,hiqual15,bte11,sctvc11,gnvq11,nvq11,rsa11,cag11,numol5,numal,numas,hst,advhst,typhst1,typhst2,typhst3,typhst4,typhst5,
                                   gcseful1,gcseful2,gcseful3,gcseful4,qgcse41,qgcse42,qgcse43,qgcse44,
                                   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,

diff --git a/R/lfs_read_2021.R b/R/lfs_read_2021.R
@@ -56,7 +56,7 @@ lfs_read_2021 <- function(
     setnames(data, names(data), tolower(names(data)))
 
     weights_vars     <- Hmisc::Cs(pwt22,piwt22)
-    demographic_vars <- Hmisc::Cs(age,sex,gor9d,ethukeul,marsta,fdpch16)
+    demographic_vars <- Hmisc::Cs(age,sex,gor9d,ethukeul,marsta,mardy6,fdpch16)
     education_vars   <- Hmisc::Cs(edage,hiqul15d,hiqual15,bte11,sctvc11,gnvq11,nvq11,rsa11,cag11,numol5,numal,numas,hst,advhst,typhst1,typhst2,typhst3,typhst4,typhst5,
                                   gcseful1,gcseful2,gcseful3,gcseful4,qgcse41,qgcse42,qgcse43,qgcse44,
                                   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,

diff --git a/R/lfs_read_2022.R b/R/lfs_read_2022.R
@@ -56,7 +56,7 @@ lfs_read_2022 <- function(
     setnames(data, names(data), tolower(names(data)))
 
     weights_vars     <- Hmisc::Cs(pwt22,piwt22)
-    demographic_vars <- Hmisc::Cs(age,sex,gor9d,ethukeul,marsta,fdpch16)
+    demographic_vars <- Hmisc::Cs(age,sex,gor9d,ethukeul,marsta,mardy6,fdpch16)
     education_vars   <- Hmisc::Cs(edage,hiqul22d,hiqual22)
     health_vars      <- Hmisc::Cs(health20,discurr20,disea,
                                   heal2001,heal2002,heal2003,heal2004,heal2005,

diff --git a/R/lfs_read_2023.R b/R/lfs_read_2023.R
@@ -56,7 +56,7 @@ lfs_read_2023 <- function(
     setnames(data, names(data), tolower(names(data)))
 
     weights_vars     <- Hmisc::Cs(pwt22,piwt22)
-    demographic_vars <- Hmisc::Cs(age,sex,gor9d,ethukeul,marsta,fdpch16)
+    demographic_vars <- Hmisc::Cs(age,sex,gor9d,ethukeul,marsta,mardy6,fdpch16)
     education_vars   <- Hmisc::Cs(edage,hiqul22d,hiqual22)
     health_vars      <- Hmisc::Cs(health20,discurr20,disea,
                                   heal2001,heal2002,heal2003,heal2004,heal2005,
@@ -101,7 +101,7 @@ lfs_read_2023 <- function(
     setnames(data, names(data), tolower(names(data)))
 
     weights_vars     <- Hmisc::Cs(pwt23,piwt22)
-    demographic_vars <- Hmisc::Cs(age,sex,gor9d,ethukeul,marsta,fdpch16)
+    demographic_vars <- Hmisc::Cs(age,sex,gor9d,ethukeul,marsta,mardy6,fdpch16)
     education_vars   <- Hmisc::Cs(edage,hiqul22d,hiqual22)
     health_vars      <- Hmisc::Cs(health20,discurr20,disea,
                                   heal2001,heal2002,heal2003,heal2004,heal2005,
@@ -149,7 +149,7 @@ lfs_read_2023 <- function(
     setnames(data, names(data), tolower(names(data)))
 
     weights_vars     <- Hmisc::Cs(pwt23,piwt23)
-    demographic_vars <- Hmisc::Cs(age,sex,gor9d,ethukeul,marsta,fdpch16)
+    demographic_vars <- Hmisc::Cs(age,sex,gor9d,ethukeul,marsta,mardy6,fdpch16)
     education_vars   <- Hmisc::Cs(edage,hiqul22d,hiqual22)
     health_vars      <- Hmisc::Cs(health20,discurr20,disea,
                                   heal2001,heal2002,heal2003,heal2004,heal2005,