# Load libraries
library(tidyverse)
## ── Attaching core tidyverse packages ──────────────────────── tidyverse 2.0.0 ──
## ✔ dplyr     1.1.4     ✔ readr     2.1.5
## ✔ forcats   1.0.0     ✔ stringr   1.5.1
## ✔ ggplot2   3.5.1     ✔ tibble    3.2.1
## ✔ lubridate 1.9.4     ✔ tidyr     1.3.1
## ✔ purrr     1.0.2     
## ── Conflicts ────────────────────────────────────────── tidyverse_conflicts() ──
## ✖ dplyr::filter() masks stats::filter()
## ✖ dplyr::lag()    masks stats::lag()
## ℹ Use the conflicted package (<http://conflicted.r-lib.org/>) to force all conflicts to become errors
library(sjPlot)
library(broom)
library(modelsummary)
## `modelsummary` 2.0.0 now uses `tinytable` as its default table-drawing
##   backend. Learn more at: https://vincentarelbundock.github.io/tinytable/
## 
## Revert to `kableExtra` for one session:
## 
##   options(modelsummary_factory_default = 'kableExtra')
##   options(modelsummary_factory_latex = 'kableExtra')
##   options(modelsummary_factory_html = 'kableExtra')
## 
## Silence this message forever:
## 
##   config_modelsummary(startup_message = FALSE)
library(gt)
library(knitr)
library(scales)
## 
## Attaching package: 'scales'
## 
## The following object is masked from 'package:purrr':
## 
##     discard
## 
## The following object is masked from 'package:readr':
## 
##     col_factor
library(kableExtra)
## 
## Attaching package: 'kableExtra'
## 
## The following object is masked from 'package:dplyr':
## 
##     group_rows
library(flextable)
## 
## Attaching package: 'flextable'
## 
## The following objects are masked from 'package:kableExtra':
## 
##     as_image, footnote
## 
## The following object is masked from 'package:purrr':
## 
##     compose
library(car)
## Loading required package: carData
## 
## Attaching package: 'car'
## 
## The following object is masked from 'package:dplyr':
## 
##     recode
## 
## The following object is masked from 'package:purrr':
## 
##     some
library(fst)
getwd()
## [1] "/Users/emmanuelchigozie/Documents"
chs_data <- read.csv("CHS2021ECL_PUMF_cc13f1e41bae9850ad651344a7e9be09.csv")
# Define Atlantic Canada region codes (adjust if needed)
atlantic_codes <- c(10, 11, 12, 13)

# Subset for Atlantic Canada
chs_atlantic <- chs_data %>% filter(PGEOGR %in% atlantic_codes)

# Clean life satisfaction variable (remove invalid values)
chs_atlantic <- chs_atlantic %>% 
  mutate(PLIS_05 = ifelse(PLIS_05 < 0, NA, PLIS_05)) 

# Convert region codes to descriptive names
chs_atlantic <- chs_atlantic %>% 
  mutate(Region = factor(PGEOGR, 
                         levels = atlantic_codes, 
                         labels = c("Newfoundland and Labrador", 
                                    "Prince Edward Island", 
                                    "Nova Scotia", 
                                    "New Brunswick")))
ggplot(chs_atlantic, aes(x = Region, y = PLIS_05, fill = Region)) +
  geom_boxplot() +
  labs(
    title = "Life Satisfaction Across Atlantic Canada",
    x = "Region",
    y = "Life Satisfaction Score",
    fill = "Region"
  ) +
  theme_minimal() +
  theme(axis.text.x = element_text(angle = 45, hjust = 1))

# Run linear regression model
lm_model <- lm(PLIS_05 ~ Region, data = chs_atlantic)

# Plot coefficient visualization
plot_model(lm_model, type = "est", vline.color = "red")

# Create regression table
tab_model(lm_model, show.ci = FALSE, title = "Regression of Life Satisfaction on Region")
Regression of Life Satisfaction on Region
  PLIS_05
Predictors Estimates p
(Intercept) 7.61 <0.001
Region [Prince Edward
Island]
-0.83 0.011
Region [Nova Scotia] -0.86 0.007
Region [New Brunswick] -0.58 0.067
Observations 5171
R2 / R2 adjusted 0.002 / 0.002

Analysis of Regression for Life Satisfaction In Atlantic Canada

Our regression sought to understand whether life satisfaction varied by region for those respondents from Atlantic Canada. Thus, by using Newfoundland and Labrador as the reference group, we assess life satisfaction with respect to Prince Edward Island, Nova Scotia, and New Brunswick to assess whether differences exist.

Life Satisfaction Differs by Region

The results indicate that life satisfaction is lower in Nova Scotia and New Brunswick relative to Newfoundland and Labrador. Therefore, these provinces are on average less satisfied with their lives; the difference is significant, therefore, by p<0.05.

Yet Prince Edward Island does not show a significant difference when assessing life satisfaction relative to Newfoundland and Labrador. Therefore, it appears that those in these two regions feel similarly about their lives; yet those in Nova Scotia and New Brunswick are significantly lower in life satisfaction relative to those in Newfoundland and Labrador.

Magnitude of the Difference and Importance The regression coefficients assess the level of change. A negative coefficient indicates a lower life satisfaction compared to Newfoundland and Labrador. For instance, Nova Scotia has a coefficient of -0.8, meaning that the average life satisfaction score in Nova Scotia is 0.8 points lower than that in Newfoundland and Labrador.

Since Prince Edward Island’s coefficient is statistically not different from zero, we can conclude that life satisfaction in PEI is nearly the same as that in Newfoundland and Labrador.

R² Value – Variance Explained

The R² value is an indicator of variance in life satisfaction explained by region alone. Therefore, the R² value is consistent with regions accounting for about 0.05 or 5% of variance in life satisfaction. This means that regions only account for a small percentage of why some people have more life satisfaction than others. This indicates that there are other highly relevant factors such as income, health, chronic disease, employment status, quality and quantity of relationships. Since region does not strongly predict much of the variance in life satisfaction, our results suggest that where one lives is not as pertinent as the other more intimate and social factors that allow for happiness.

Is Region a Strong Predictor of Life Satisfaction?

In summary, region is a weak to moderate predictor of life satisfaction. For instance, Nova Scotia and New Brunswick have significantly lower life satisfaction, but region overall is not a strong predictor and the aggregate effect is small. Thus, where one lives does help determine happiness but not as much as other factors that contribute to more happiness over time.

This is a key takeaway for municipal leadership and governing bodies. For example, if Nova Scotia and New Brunswick have lower life satisfaction, perhaps the municipal governments could take some time to see why. Is there a regional economic downturn? Are there less social services? Are there quality of life issues that differentiate from this situation? This could help facilitate change for regional sustainability. Ultimately, we discovered that Nova Scotia and New Brunswick are less satisfied with life than Newfoundland and Labrador, yet they have comparable satisfaction to PEI. However, the region by itself was not a significant predictor of life satisfaction. Therefore, it seems that people and their situations are more important than any projected average value of life satisfaction by provincial boundaries.

portugal_data <- read_fst("portugal_data.fst")
colnames(portugal_data)
##    [1] "name"      "essround"  "edition"   "proddate"  "cntry"     "idno"     
##    [7] "dweight"   "pspwght"   "pweight"   "anweight"  "prob"      "stratum"  
##   [13] "psu"       "netuse"    "netusoft"  "netustm"   "nwspol"    "nwsppol"  
##   [19] "nwsptot"   "pplfair"   "pplhlp"    "ppltrst"   "rdpol"     "rdtot"    
##   [25] "tvpol"     "tvtot"     "actrolg"   "actrolga"  "badge"     "bctprd"   
##   [31] "bghtprd"   "clsprty"   "contplt"   "cptppol"   "cptppola"  "dclagr"   
##   [37] "dclaid"    "dclcrm"    "dcldef"    "dclenv"    "dclintr"   "dclmig"   
##   [43] "dclwlfr"   "dmcntov"   "dntmny"    "ecohenv"   "etapapl"   "euftf"    
##   [49] "freehms"   "gincdif"   "ginveco"   "hmsacld"   "hmsfmlsh"  "ilglpst"  
##   [55] "lawobey"   "lrscale"   "mmbprty"   "needtru"   "pbldmn"    "pbldmna"  
##   [61] "pltcare"   "pltinvt"   "polactiv"  "polcmpl"   "poldcs"    "polintr"  
##   [67] "prtdgcl"   "prtyban"   "psppipl"   "psppipla"  "psppsgv"   "psppsgva" 
##   [73] "pstplonl"  "ptcpplt"   "scnsenv"   "sgnptit"   "stfdem"    "stfeco"   
##   [79] "stfedu"    "stfgov"    "stfhlth"   "stflife"   "trstep"    "trstlgl"  
##   [85] "trstplc"   "trstplt"   "trstprl"   "trstprt"   "trstun"    "vote"     
##   [91] "wrkorg"    "wrkprty"   "imsmetn"   "imdfetn"   "impcntr"   "imbgeco"  
##   [97] "imueclt"   "imwbcnt"   "volunfp"   "lrnobed"   "loylead"   "donprty"  
##  [103] "prtvde1"   "prtvade1"  "prtvbde1"  "prtvcde1"  "prtvdde1"  "prtvede1" 
##  [109] "prtvfde1"  "prtvde2"   "prtvade2"  "prtvbde2"  "prtvcde2"  "prtvdde2" 
##  [115] "prtvede2"  "prtvfde2"  "prtvlt1"   "prtvalt1"  "prtvblt1"  "prtvclt1" 
##  [121] "prtvlt2"   "prtvalt2"  "prtvblt2"  "prtvclt2"  "prtvlt3"   "prtvalt3" 
##  [127] "prtvblt3"  "prtvclt3"  "prtvtal"   "prtvtat"   "prtvtaat"  "prtvtbat" 
##  [133] "prtvtcat"  "prtvtbe"   "prtvtabe"  "prtvtbbe"  "prtvtcbe"  "prtvtdbe" 
##  [139] "prtvtebe"  "prtvtbg"   "prtvtabg"  "prtvtbbg"  "prtvtcbg"  "prtvtdbg" 
##  [145] "prtvtebg"  "prtvtch"   "prtvtach"  "prtvtbch"  "prtvtcch"  "prtvtdch" 
##  [151] "prtvtech"  "prtvtfch"  "prtvtgch"  "prtvthch"  "prtvtcy"   "prtvtacy" 
##  [157] "prtvtbcy"  "prtvtccy"  "prtvtcz"   "prtvtacz"  "prtvtbcz"  "prtvtccz" 
##  [163] "prtvtdcz"  "prtvtecz"  "prtvtdk"   "prtvtadk"  "prtvtbdk"  "prtvtcdk" 
##  [169] "prtvtddk"  "prtvtee"   "prtvtaee"  "prtvtbee"  "prtvtcee"  "prtvtdee" 
##  [175] "prtvteee"  "prtvtfee"  "prtvtgee"  "prtvthee"  "prtvtes"   "prtvtaes" 
##  [181] "prtvtbes"  "prtvtces"  "prtvtdes"  "prtvtees"  "prtvtfes"  "prtvtfi"  
##  [187] "prtvtafi"  "prtvtbfi"  "prtvtcfi"  "prtvtdfi"  "prtvtefi"  "prtvtfr"  
##  [193] "prtvtafr"  "prtvtbfr"  "prtvtcfr"  "prtvtdfr"  "prtvtefr"  "prtvtgb"  
##  [199] "prtvtagb"  "prtvtbgb"  "prtvtcgb"  "prtvtdgb"  "prtvtgr"   "prtvtagr" 
##  [205] "prtvtbgr"  "prtvtcgr"  "prtvtdgr"  "prtvthr"   "prtvtahr"  "prtvtbhr" 
##  [211] "prtvthu"   "prtvtahu"  "prtvtbhu"  "prtvtchu"  "prtvtdhu"  "prtvtehu" 
##  [217] "prtvtfhu"  "prtvtghu"  "prtvtie"   "prtvtaie"  "prtvtbie"  "prtvtcie" 
##  [223] "prtvtdie"  "prtvtil"   "prtvtail"  "prtvtbil"  "prtvtcil"  "prtvtdil" 
##  [229] "prtvtis"   "prtvtais"  "prtvtbis"  "prtvtcis"  "prtvtdis"  "prtvtit"  
##  [235] "prtvtait"  "prtvtbit"  "prtvtcit"  "prtvtdit"  "prtvtlu"   "prtvtlv"  
##  [241] "prtvtalv"  "prtvtnl"   "prtvtanl"  "prtvtbnl"  "prtvtcnl"  "prtvtdnl" 
##  [247] "prtvtenl"  "prtvtfnl"  "prtvtgnl"  "prtvthnl"  "prtvtno"   "prtvtano" 
##  [253] "prtvtbno"  "prtvtpl"   "prtvtapl"  "prtvtbpl"  "prtvtcpl"  "prtvtdpl" 
##  [259] "prtvtepl"  "prtvtpt"   "prtvtapt"  "prtvtbpt"  "prtvtcpt"  "prtvtdpt" 
##  [265] "prtvtro"   "prtvtaro"  "prtvtru"   "prtvtaru"  "prtvtbru"  "prtvtcru" 
##  [271] "prtvtdru"  "prtvtse"   "prtvtase"  "prtvtbse"  "prtvtcse"  "prtvtdse" 
##  [277] "prtvtsi"   "prtvtasi"  "prtvtbsi"  "prtvtcsi"  "prtvtdsi"  "prtvtesi" 
##  [283] "prtvtfsi"  "prtvtsk"   "prtvtask"  "prtvtbsk"  "prtvtcsk"  "prtvtdsk" 
##  [289] "prtvtesk"  "prtvttr"   "prtvtatr"  "prtvtua"   "prtvtaua"  "prtvtbua" 
##  [295] "prtvtcua"  "prtvtxk"   "prtvtrs"   "prtvtars"  "prtvtme"   "prtvtame" 
##  [301] "prtclal"   "prtclat"   "prtclaat"  "prtclcat"  "prtcldat"  "prtclbe"  
##  [307] "prtclabe"  "prtclbbe"  "prtclcbe"  "prtcldbe"  "prtclebe"  "prtclbg"  
##  [313] "prtclabg"  "prtclbbg"  "prtclcbg"  "prtcldbg"  "prtclebg"  "prtclch"  
##  [319] "prtclach"  "prtclbch"  "prtclcch"  "prtcldch"  "prtclech"  "prtclfch" 
##  [325] "prtclgch"  "prtclhch"  "prtclsk"   "prtclask"  "prtclbsk"  "prtclcsk" 
##  [331] "prtcldsk"  "prtclesk"  "prtclcy"   "prtclacy"  "prtclbcy"  "prtclccy" 
##  [337] "prtclcz"   "prtclacz"  "prtclbcz"  "prtclccz"  "prtcldcz"  "prtclecz" 
##  [343] "prtclde"   "prtclade"  "prtclbde"  "prtclcde"  "prtcldde"  "prtclede" 
##  [349] "prtclfde"  "prtcldk"   "prtcladk"  "prtclbdk"  "prtclcdk"  "prtclddk" 
##  [355] "prtclee"   "prtclaee"  "prtclbee"  "prtclcee"  "prtcldee"  "prtcleee" 
##  [361] "prtclfee"  "prtclgee"  "prtclhee"  "prtcles"   "prtclaes"  "prtclbes" 
##  [367] "prtclces"  "prtcldes"  "prtclees"  "prtclfes"  "prtclges"  "prtclfi"  
##  [373] "prtclafi"  "prtclbfi"  "prtclcfi"  "prtcldfi"  "prtclefi"  "prtclffi" 
##  [379] "prtclfr"   "prtclafr"  "prtclbfr"  "prtclcfr"  "prtcldfr"  "prtclefr" 
##  [385] "prtclffr"  "prtclgb"   "prtclagb"  "prtclbgb"  "prtclcgb"  "prtcldgb" 
##  [391] "prtclgr"   "prtclagr"  "prtclbgr"  "prtclcgr"  "prtcldgr"  "prtclhr"  
##  [397] "prtclahr"  "prtclbhr"  "prtclhu"   "prtclahu"  "prtclbhu"  "prtclchu" 
##  [403] "prtcldhu"  "prtclehu"  "prtclfhu"  "prtclghu"  "prtclhhu"  "prtclie"  
##  [409] "prtclaie"  "prtclbie"  "prtclcie"  "prtcldie"  "prtcleie"  "prtclfie" 
##  [415] "prtclil"   "prtclail"  "prtclbil"  "prtclcil"  "prtcldil"  "prtcleil" 
##  [421] "prtclis"   "prtclais"  "prtclbis"  "prtclcis"  "prtcldis"  "prtclit"  
##  [427] "prtclait"  "prtclbit"  "prtclcit"  "prtcldit"  "prtcleit"  "prtcllt"  
##  [433] "prtclalt"  "prtclblt"  "prtclclt"  "prtcllu"   "prtcllv"   "prtclalv" 
##  [439] "prtclnl"   "prtclanl"  "prtclbnl"  "prtclcnl"  "prtcldnl"  "prtclenl" 
##  [445] "prtclfnl"  "prtclgnl"  "prtclno"   "prtclano"  "prtclbno"  "prtclpl"  
##  [451] "prtclapl"  "prtclbpl"  "prtclcpl"  "prtcldpl"  "prtclepl"  "prtclfpl" 
##  [457] "prtclgpl"  "prtclhpl"  "prtclipl"  "prtclpt"   "prtclapt"  "prtclbpt" 
##  [463] "prtclcpt"  "prtcldpt"  "prtclept"  "prtclfpt"  "prtclro"   "prtclaro" 
##  [469] "prtclru"   "prtclaru"  "prtclbru"  "prtclcru"  "prtcldru"  "prtclse"  
##  [475] "prtclase"  "prtclbse"  "prtclcse"  "prtcldse"  "prtclsi"   "prtclasi" 
##  [481] "prtclbsi"  "prtclcsi"  "prtcldsi"  "prtclesi"  "prtclfsi"  "prtcltr"  
##  [487] "prtclatr"  "prtclua"   "prtclaua"  "prtclbua"  "prtclcua"  "prtcldua" 
##  [493] "prtclxk"   "prtclrs"   "prtclars"  "prtclme"   "prtclame"  "prtmbat"  
##  [499] "prtmbaat"  "prtmbbe"   "prtmbabe"  "prtmbbbe"  "prtmbcbe"  "prtmbbg"  
##  [505] "prtmbabg"  "prtmbbbg"  "prtmbch"   "prtmbach"  "prtmbbch"  "prtmbcch" 
##  [511] "prtmbcy"   "prtmbcz"   "prtmbacz"  "prtmbbcz"  "prtmbde"   "prtmbade" 
##  [517] "prtmbbde"  "prtmbcde"  "prtmbdk"   "prtmbadk"  "prtmbbdk"  "prtmbee"  
##  [523] "prtmbaee"  "prtmbbee"  "prtmbcee"  "prtmbes"   "prtmbaes"  "prtmbbes" 
##  [529] "prtmbfi"   "prtmbafi"  "prtmbbfi"  "prtmbfr"   "prtmbafr"  "prtmbbfr" 
##  [535] "prtmbcfr"  "prtmbgb"   "prtmbagb"  "prtmbgr"   "prtmbagr"  "prtmbbgr" 
##  [541] "prtmbcgr"  "prtmbhr"   "prtmbhu"   "prtmbahu"  "prtmbbhu"  "prtmbchu" 
##  [547] "prtmbie"   "prtmbaie"  "prtmbbie"  "prtmbil"   "prtmbail"  "prtmbbil" 
##  [553] "prtmbis"   "prtmbit"   "prtmbait"  "prtmblt"   "prtmblu"   "prtmblv"  
##  [559] "prtmbnl"   "prtmbanl"  "prtmbbnl"  "prtmbcnl"  "prtmbno"   "prtmbano" 
##  [565] "prtmbpl"   "prtmbapl"  "prtmbbpl"  "prtmbcpl"  "prtmbdpl"  "prtmbpt"  
##  [571] "prtmbapt"  "prtmbbpt"  "prtmbro"   "prtmbaro"  "prtmbru"   "prtmbaru" 
##  [577] "prtmbbru"  "prtmbse"   "prtmbase"  "prtmbsi"   "prtmbasi"  "prtmbbsi" 
##  [583] "prtmbcsi"  "prtmbsk"   "prtmbask"  "prtmbbsk"  "prtmbtr"   "prtmbatr" 
##  [589] "prtmbua"   "prtmbaua"  "prtmbbua"  "prtmbcua"  "aesfdrk"   "atchctr"  
##  [595] "atcherp"   "blgetmg"   "brghmef"   "brghmwr"   "brncntr"   "cntbrth"  
##  [601] "cntbrtha"  "cntbrthb"  "cntbrthc"  "cntbrthd"  "crmvct"    "crvctef"  
##  [607] "crvctwr"   "ctzcntr"   "ctzship"   "ctzshipa"  "ctzshipb"  "ctzshipc" 
##  [613] "ctzshipd"  "dscrage"   "dscrdk"    "dscrdsb"   "dscretn"   "dscrgnd"  
##  [619] "dscrgrp"   "dscrlng"   "dscrna"    "dscrnap"   "dscrntn"   "dscroth"  
##  [625] "dscrrce"   "dscrref"   "dscrrlg"   "dscrsex"   "facntn"    "facntr"   
##  [631] "fbrncnt"   "fbrncnta"  "fbrncntb"  "fbrncntc"  "happy"     "health"   
##  [637] "hlthhmp"   "inmdisc"   "inprdsc"   "livecntr"  "livecnta"  "lnghoma"  
##  [643] "lnghom1"   "lnghomb"   "lnghom2"   "mbrncnt"   "mbrncnta"  "mbrncntb" 
##  [649] "mbrncntc"  "mmbrn"     "mocntn"    "mocntr"    "pray"      "rlgatnd"  
##  [655] "rlgblg"    "rlgblge"   "rlgdgr"    "rlgdnm"    "rlgdnme"   "sclact"   
##  [661] "sclmeet"   "trrcnyr"   "trrenyr"   "trrprsn"   "trrtort"   "feethngr" 
##  [667] "vteurmmb"  "vteubcmb"  "vteumbgb"  "rlgdnal"   "rlgdnat"   "rlgdnbat" 
##  [673] "rlgdnbe"   "rlgdnch"   "rlgdnach"  "rlgdncy"   "rlgdnacy"  "rlgdnde"  
##  [679] "rlgdnade"  "rlgdndk"   "rlgdnfi"   "rlgdnafi"  "rlgdngb"   "rlgdngr"  
##  [685] "rlgdnagr"  "rlgdnhu"   "rlgdnie"   "rlgdnil"   "rlgdnis"   "rlgdnais" 
##  [691] "rlgdnlt"   "rlgdnlv"   "rlgdnnl"   "rlgdnanl"  "rlgdnno"   "rlgdnpl"  
##  [697] "rlgdnapl"  "rlgdnbpl"  "rlgdnpt"   "rlgdnro"   "rlgdnru"   "rlgdnaru" 
##  [703] "rlgdnse"   "rlgdnase"  "rlgdnsi"   "rlgdnsk"   "rlgdnask"  "rlgdnbsk" 
##  [709] "rlgdnua"   "rlgdnrs"   "rlgdme"    "rlgdeal"   "rlgdeat"   "rlgdebat" 
##  [715] "rlgdebe"   "rlgdech"   "rlgdeach"  "rlgdecy"   "rlgdeacy"  "rlgdede"  
##  [721] "rlgdeade"  "rlgdedk"   "rlgdefi"   "rlgdeafi"  "rlgdegb"   "rlgdegr"  
##  [727] "rlgdeagr"  "rlgdehu"   "rlgdeie"   "rlgdeil"   "rlgdeis"   "rlgdeais" 
##  [733] "rlgdelt"   "rlgdelv"   "rlgdenl"   "rlgdeanl"  "rlgdeno"   "rlgdepl"  
##  [739] "rlgdeapl"  "rlgdebpl"  "rlgdept"   "rlgdero"   "rlgderu"   "rlgdearu" 
##  [745] "rlgdese"   "rlgdease"  "rlgdesi"   "rlgdesk"   "rlgdeask"  "rlgdebsk" 
##  [751] "rlgdeua"   "rlgders"   "rlgdeme"   "hhmmb"     "gndr"      "gndr2"    
##  [757] "gndr3"     "gndr4"     "gndr5"     "gndr6"     "gndr7"     "gndr8"    
##  [763] "gndr9"     "gndr10"    "gndr11"    "gndr12"    "gndr13"    "gndr14"   
##  [769] "gndr15"    "gndr16"    "gndr17"    "gndr18"    "gndr19"    "gndr20"   
##  [775] "gndr21"    "gndr22"    "gndr23"    "gndr24"    "rship2"    "rshipa2"  
##  [781] "rship3"    "rshipa3"   "rship4"    "rshipa4"   "rship5"    "rshipa5"  
##  [787] "rship6"    "rshipa6"   "rship7"    "rshipa7"   "rship8"    "rshipa8"  
##  [793] "rship9"    "rshipa9"   "rship10"   "rshipa10"  "rship11"   "rshipa11" 
##  [799] "rship12"   "rshipa12"  "rship13"   "rshipa13"  "rship14"   "rshipa14" 
##  [805] "rship15"   "rshipa15"  "rshipa16"  "rshipa17"  "rshipa18"  "rshipa19" 
##  [811] "rshipa20"  "rshipa21"  "rshipa22"  "rshipa23"  "rshipa24"  "yrbrn2"   
##  [817] "yrbrn3"    "yrbrn4"    "yrbrn5"    "yrbrn6"    "yrbrn7"    "yrbrn8"   
##  [823] "yrbrn9"    "yrbrn10"   "yrbrn11"   "yrbrn12"   "yrbrn13"   "yrbrn14"  
##  [829] "yrbrn15"   "yrbrn16"   "yrbrn17"   "yrbrn18"   "yrbrn19"   "yrbrn20"  
##  [835] "yrbrn21"   "yrbrn22"   "yrbrn23"   "yrbrn24"   "age"       "agea"     
##  [841] "ablrtr"    "anctry1"   "anctry2"   "atncrse"   "brwmny"    "chldhhe"  
##  [847] "chldhm"    "cmsrv"     "cmsrvp"    "crpdwk"    "crpdwkp"   "dngdk"    
##  [853] "dngdkp"    "dngna"     "dngnap"    "dngnapp"   "dngoth"    "dngothp"  
##  [859] "dngref"    "dngrefp"   "domicil"   "dsbld"     "dsbldp"    "dvrcdev"  
##  [865] "dvrcdeva"  "edctn"     "edctnp"    "edufld"    "edulvla"   "edulvlb"  
##  [871] "edulvlfa"  "edulvlfb"  "edulvlma"  "edulvlmb"  "edulvlpa"  "edulvlpb" 
##  [877] "eduyrs"    "eisced"    "eiscedf"   "eiscedm"   "eiscedp"   "emplno"   
##  [883] "emplnof"   "emplnom"   "emplnop"   "emplrel"   "emprelp"   "emprf14"  
##  [889] "emprm14"   "estsz"     "fxltph"    "faalv"     "hhmodwl"   "hincfel"  
##  [895] "hincsrc"   "hincsrca"  "hinctnt"   "hinctnta"  "hswrk"     "hswrkp"   
##  [901] "iccohbt"   "icmnart"   "icmsw"     "icomdng"   "icomdnp"   "icpdwk2"  
##  [907] "icpdwrk"   "icppdwk"   "icwhct"    "inttph"    "ioactp"    "iorgact"  
##  [913] "isco08"    "isco08p"   "iscoco"    "iscocop"   "jbspv"     "jbspvf"   
##  [919] "jbspvm"    "jbspvp"    "lvghw"     "lvghwa"    "lvgoptn"   "lvgptn"   
##  [925] "lvgptna"   "lvgptne"   "lvgptnea"  "mainact"   "marital"   "maritala" 
##  [931] "maritalb"  "marsts"    "mbltph"    "mbtru"     "mnactic"   "mnactp"   
##  [937] "moalv"     "njbspv"    "njbspvp"   "occf14"    "occf14a"   "occf14b"  
##  [943] "occf14ie"  "occm14"    "occm14a"   "occm14b"   "occm14ie"  "orgwrk"   
##  [949] "partner"   "pdjobev"   "pdjobyr"   "pdwrk"     "pdwrkp"    "rmhhus"   
##  [955] "rshpsts"   "rtdsbyr"   "rtrd"      "rtrdp"     "tporgwk"   "uemp12m"  
##  [961] "uemp3m"    "uemp5yr"   "uempla"    "uemplap"   "uempli"    "uemplip"  
##  [967] "wkdcorp"   "wkhtot"    "wkhtotp"   "wrkac6m"   "wrkctr"    "wrkctra"  
##  [973] "wrkctrhu"  "wkdcorg"   "wkdcorga"  "nacer1"    "nacer11"   "nacer2"   
##  [979] "edagegb"   "edlvat"    "edlveat"   "eduat1"    "eduat2"    "edlvbe"   
##  [985] "edlvabe"   "edlvdbe"   "edlvebe"   "edlvbg"    "edlvdbg"   "edlvebg"  
##  [991] "edlvch"    "edlvach"   "edlvbch"   "edlvcch"   "edlvdch"   "edlvcy"   
##  [997] "edlvacy"   "edlvdcy"   "edlvecy"   "edlvgcy"   "edlvcz"    "edlvdcz"  
## [1003] "edlvdal"   "edlvde"    "edlvade"   "edlvdfi"   "edlvdis"   "edlvdk"   
## [1009] "edlvadk"   "edlvddk"   "edlvdxk"   "edlvee"    "edlvaee"   "edlvbee"  
## [1015] "edlvdee"   "edlves"    "edlvaes"   "edlvdes"   "edlvees"   "edlvges"  
## [1021] "edlvhes"   "edlvfr"    "edlvafr"   "edlvbfr"   "edlvdfr"   "edlvgb"   
## [1027] "edlvagb"   "edlvgr"    "edlvagr"   "edlvdgr"   "edlvegr"   "edlvhr"   
## [1033] "edlvdhr"   "edlvehr"   "edlvhu"    "edlvahu"   "edlvbhu"   "edlvdhu"  
## [1039] "edlvdahu"  "edlvie"    "edlvaie"   "edlvbie"   "edlvdie"   "edlvil"   
## [1045] "edlvail"   "eduil1"    "eduail1"   "edubil1"   "edlvit"    "edlvait"  
## [1051] "edlvdit"   "edlveit"   "edlvlt"    "edlvdlt"   "edlvlu"    "edlvlv"   
## [1057] "edlvdlv"   "edlvelv"   "edlvnl"    "edlvdnl"   "edlvenl"   "edlvno"   
## [1063] "edlvano"   "edlvdno"   "edlveno"   "edlvpl"    "edlvapl"   "edlvbpl"  
## [1069] "edlvdpl"   "edlvepl"   "edlvgpl"   "edlvhpl"   "edlvpt"    "edlvapt"  
## [1075] "edlvbpt"   "edlvdpt"   "edlvro"    "edlvru"    "edlvdru"   "edlvse"   
## [1081] "edlvase"   "edlvdse"   "edlvsi"    "edlvasi"   "edlvdsi"   "edlvesi"  
## [1087] "edlvsk"    "edlvask"   "edlvdsk"   "edlvtr"    "edlvua"    "edlvaua"  
## [1093] "edlvdua"   "edude1"    "eduade1"   "edubde1"   "educde1"   "edude2"   
## [1099] "eduade2"   "edubde2"   "edude3"    "eduade3"   "edugb2"    "eduagb2"  
## [1105] "edubgb2"   "eduil2"    "eduail2"   "edupl2"    "eduyrpl"   "edugb1"   
## [1111] "eduagb1"   "edubgb1"   "educgb1"   "edlvdrs"   "edlvdme"   "edagepgb" 
## [1117] "edlvpat"   "edlvpeat"  "edlvpfat"  "edlvpbe"   "edlvpdbe"  "edlvpebe" 
## [1123] "edlvpch"   "edlvpdch"  "edlvpcy"   "edlvpdcy"  "edlvpecy"  "edlvpgcy" 
## [1129] "edlvpcz"   "edlvpdcz"  "edlvpdbg"  "edlvpebg"  "edlvpdfi"  "edlvpdis" 
## [1135] "edlvpdit"  "edlvpeit"  "edlvpdk"   "edlvpddk"  "edlvpdxk"  "edlvpee"  
## [1141] "edlvpdee"  "edlvpes"   "edlvpdes"  "edlvpees"  "edlvpfes"  "edlvpges" 
## [1147] "edlvpfr"   "edlvpdfr"  "edlvpgb"   "edlvpgr"   "edlvpdgr"  "edlvpegr" 
## [1153] "edlvphr"   "edlvpdhr"  "edlvpehr"  "edlvphu"   "edlvpdhu"  "edlvpdahu"
## [1159] "edlvpie"   "edlvpdie"  "edlvpil"   "edupil1"   "edupail1"  "edupbil1" 
## [1165] "edlvplt"   "edlvpdlt"  "edlvplv"   "edlvpdlv"  "edlvpelv"  "edlvpnl"  
## [1171] "edlvpdnl"  "edlvpenl"  "edlvpno"   "edlvpdno"  "edlvpeno"  "edlvppl"  
## [1177] "edlvpdpl"  "edlvpepl"  "edlvpfpl"  "edlvpgpl"  "edlvppt"   "edlvpdpt" 
## [1183] "edlvpro"   "edlvpru"   "edlvpdru"  "edlvpse"   "edlvpdse"  "edlvpsi"  
## [1189] "edlvpdsi"  "edlvpesi"  "edlvpsk"   "edlvpdsk"  "edlvptr"   "edlvpua"  
## [1195] "edlvpdua"  "edupde1"   "edupade1"  "edupbde1"  "edupcde1"  "edupde2"  
## [1201] "edupade2"  "edupbde2"  "edupde3"   "edupade3"  "edupgb2"   "edupagb2" 
## [1207] "edupbgb2"  "edupil2"   "edupail2"  "eduppl2"   "eduyrppl"  "edlvpdal" 
## [1213] "edupgb1"   "edupagb1"  "edupbgb1"  "edupcgb1"  "edlvpdrs"  "edlvpdme" 
## [1219] "edlvfdal"  "edlvfat"   "edlvfeat"  "edlvfbe"   "edlvfdbe"  "edlvfebe" 
## [1225] "edlvfdbg"  "edlvfebg"  "edlvfhr"   "edlvfdhr"  "edlvfehr"  "edlvfcy"  
## [1231] "edlvfdcy"  "edlvfecy"  "edlvfgcy"  "edlvfcz"   "edlvfdcz"  "edlvfdk"  
## [1237] "edlvfddk"  "edlvfee"   "edlvfdee"  "edlvfdfi"  "edlvffr"   "edlvfdfr" 
## [1243] "edufde1"   "edufade1"  "edufbde1"  "edufde2"   "edufade2"  "edufde3"  
## [1249] "edufade3"  "edlvfgr"   "edlvfdgr"  "edlvfegr"  "edlvfhu"   "edlvfdhu" 
## [1255] "edlvfdahu" "edlvfdis"  "edlvfie"   "edlvfdie"  "edlvfil"   "edufil1"  
## [1261] "edufail1"  "edufbil1"  "edufil2"   "edufail2"  "edlvfdit"  "edlvfeit" 
## [1267] "edlvfdxk"  "edlvflv"   "edlvfdlv"  "edlvflt"   "edlvfdlt"  "edlvfnl"  
## [1273] "edlvfdnl"  "edlvfenl"  "edlvfno"   "edlvfdno"  "edlvfeno"  "edlvfpl"  
## [1279] "edlvfdpl"  "edlvfepl"  "edlvffpl"  "edlvfpt"   "edlvfdpt"  "edlvfro"  
## [1285] "edlvfru"   "edlvfdru"  "edlvfsk"   "edlvfdsk"  "edlvfsi"   "edlvfdsi" 
## [1291] "edlvfesi"  "edlvfes"   "edlvfdes"  "edlvfees"  "edlvffes"  "edlvfse"  
## [1297] "edlvfdse"  "edlvfch"   "edlvfdch"  "edlvftr"   "edlvfua"   "edlvfdua" 
## [1303] "edlvfgb"   "edufgb1"   "edufagb1"  "edufbgb1"  "edufcgb1"  "edufgb2"  
## [1309] "edufagb2"  "edufbgb2"  "edagefgb"  "edlvfdrs"  "edlvfdme"  "edagemgb" 
## [1315] "edlvmat"   "edlvmeat"  "edlvmbe"   "edlvmdbe"  "edlvmebe"  "edlvmch"  
## [1321] "edlvmdch"  "edlvmcy"   "edlvmdcy"  "edlvmecy"  "edlvmgcy"  "edlvmcz"  
## [1327] "edlvmdcz"  "edlvmdal"  "edlvmdfi"  "edlvmdis"  "edlvmdit"  "edlvmeit" 
## [1333] "edlvmdk"   "edlvmddk"  "edlvmdxk"  "edlvmee"   "edlvmdee"  "edlvmes"  
## [1339] "edlvmdes"  "edlvmees"  "edlvmfes"  "edlvmfr"   "edlvmdfr"  "edlvmgb"  
## [1345] "edlvmgr"   "edlvmdgr"  "edlvmegr"  "edlvmhr"   "edlvmdhr"  "edlvmehr" 
## [1351] "edlvmhu"   "edlvmdhu"  "edlvmdahu" "edlvmie"   "edlvmdie"  "edlvmil"  
## [1357] "edumil1"   "edumail1"  "edumbil1"  "edlvmlt"   "edlvmdlt"  "edlvmlv"  
## [1363] "edlvmdlv"  "edlvmnl"   "edlvmdnl"  "edlvmenl"  "edlvmno"   "edlvmdno" 
## [1369] "edlvmeno"  "edlvmpl"   "edlvmdpl"  "edlvmepl"  "edlvmfpl"  "edlvmpt"  
## [1375] "edlvmdpt"  "edlvmro"   "edlvmru"   "edlvmdru"  "edlvmse"   "edlvmdse" 
## [1381] "edlvmsi"   "edlvmdsi"  "edlvmesi"  "edlvmsk"   "edlvmdsk"  "edlvmtr"  
## [1387] "edlvmua"   "edlvmdua"  "edumde1"   "edumade1"  "edumbde1"  "edumde2"  
## [1393] "edumade2"  "edumde3"   "edumade3"  "edumgb1"   "edumagb1"  "edumbgb1" 
## [1399] "edumcgb1"  "edumgb2"   "edumagb2"  "edumbgb2"  "edumil2"   "edumail2" 
## [1405] "edlvmdbg"  "edlvmebg"  "edlvmdrs"  "edlvmdme"  "marstfi"   "marstgb"  
## [1411] "marstie"   "martlfr"   "rshpsfi"   "rshpsgb"   "impdiff"   "impenv"   
## [1417] "impfree"   "impfun"    "imprich"   "impsafe"   "imptrad"   "ipadvnt"  
## [1423] "ipbhprp"   "ipcrtiv"   "ipeqopt"   "ipfrule"   "ipgdtim"   "iphlppl"  
## [1429] "iplylfr"   "ipmodst"   "iprspot"   "ipshabt"   "ipstrgv"   "ipsuces"  
## [1435] "ipudrst"   "regionat"  "regionbe"  "regionbg"  "regionch"  "regioach" 
## [1441] "regioncy"  "regioncz"  "regioacz"  "regionde"  "regiondk"  "regioadk" 
## [1447] "regionee"  "regiones"  "regioaes"  "regionfi"  "regioafi"  "regionfr" 
## [1453] "regiongb"  "regiongr"  "regioagr"  "regionhr"  "regionhu"  "regionie" 
## [1459] "regioaie"  "regiobie"  "regionil"  "regionis"  "regionit"  "regionlt" 
## [1465] "regionlu"  "regionlv"  "regionnl"  "regionno"  "regionpl"  "regionpt" 
## [1471] "regioapt"  "regionro"  "regionru"  "regionse"  "regionsi"  "regionsk" 
## [1477] "regiontr"  "regionua"  "regunit"   "inwtm"     "inwshh"    "inwsmm"   
## [1483] "inwyys"    "inwdd"     "inwdds"    "inwmm"     "inwmms"    "inwyr"    
## [1489] "inwyye"    "inwds"     "inwehh"    "inwemm"    "inwdde"    "inwmme"   
## [1495] "inwde"     "spltadm"   "spltadma"  "spltadmb"  "spltadmc"  "spltadmd" 
## [1501] "spltadme"  "spltadmf"  "supqad1"   "supqadm"   "supqad2"   "supqdd"   
## [1507] "supqmm"    "supqyr"    "intewde"   "mode"      "actvgrs"   "agertr"   
## [1513] "bseftwk"   "bsmw"      "chdohwk"   "cld12hh"   "cldcrmr"   "cldcrot"  
## [1519] "cldnhh"    "cldnhhd"   "cldnhhg"   "cldnhhn"   "clfncrc"   "clfncsp"  
## [1525] "clhwkrc"   "clhwksp"   "clmrlx"    "clnhhbo"   "clnhhby"   "clnhhyb"  
## [1531] "crspce"    "cuthheq"   "dfcnswk"   "dfcnswka"  "dsdclve"   "dsgrhwk"  
## [1537] "dsgrmny"   "dsgrmnya"  "dsgrpwk"   "edul12m"   "edupdem"   "emplrjb"  
## [1543] "flthmcc"   "fmpvtst"   "frshrst"   "frsptjb"   "fstjbyr"   "fthcncr"  
## [1549] "gdsprt"    "grspay"    "grspaya"   "gvprppv"   "gwdvhwk"   "gwhhprc"  
## [1555] "hlpcowk"   "hlthrwk"   "hmeqphw"   "hwkmono"   "hwkpwd1"   "hwkpwd2"  
## [1561] "hwkpwdp"   "hwkpwe1"   "hwkpwe2"   "hwkpwep"   "hwkstrs"   "hwktwd1"  
## [1567] "hwktwd2"   "hwktwe1"   "hwktwe2"   "hwwkhs"    "icag45y"   "icago45"  
## [1573] "icagu70"   "icagu70a"  "icchld"    "iccldnh"   "icempl"    "icemplr"  
## [1579] "icgndr"    "ichwk1"    "ichwk2"    "icmnac"    "icmnact"   "icnopfm"  
## [1585] "icnopfma"  "icpart1"   "icpart2"   "icpart3"   "icptn"     "icptnwk"  
## [1591] "icptnwka"  "indisdc"   "ipjbhin"   "ipjbini"   "ipjbprm"   "ipjbscr"  
## [1597] "ipjbtro"   "ipjbwfm"   "jbcoedu"   "jbedyrs"   "jblrn"     "jbrqlrn"  
## [1603] "jbscr"     "jbtmppm"   "jbtsktm"   "lfintr"    "linwk3y"   "lkafohh"  
## [1609] "lscjb3y"   "lsintjb"   "mlohinc"   "mnrsefw"   "mnrsphm"   "nbsrsp"   
## [1615] "netpay"    "nevdnjb"   "npemp3y"   "oprtad"    "orgfd3y"   "payprd"   
## [1621] "payprda"   "pdjbndm"   "phwwkhs"   "plnchld"   "pphincr"   "ppwwkp"   
## [1627] "prmpls"    "prntghr"   "prrfmly"   "pthcncr"   "ptmhmcc"   "ptnwkwe"  
## [1633] "puemp3y"   "pwkhsch"   "quclss"    "rdpay3y"   "rpljbde"   "rsnlvem"  
## [1639] "rtryr"     "scrsefw"   "smblvjb"   "stdask"    "stdhrsw"   "stdlvl"   
## [1645] "stdmcdo"   "stdtruf"   "stpvtfm"   "tchints"   "tchlcrt"   "tchlp"    
## [1651] "tchtruf"   "tmtowrk"   "tngdohm"   "trndnjb"   "truinwk"   "uemp3y"   
## [1657] "updhlrl"   "updhlrp"   "useknwl"   "vrtywrk"   "wgdpeft"   "wkdcpce"  
## [1663] "wkengtp"   "wkhsch"    "wkjbndm"   "wkovrtm"   "wkovtmp"   "wkshr3y"  
## [1669] "wmcpwrk"   "wntrtr"    "wpmwkcp"   "wrkengt"   "wrkhrd"    "wrkspv"   
## [1675] "wrkwe"     "yrcremp"   "yrlvptn"   "yrskdwk"   "yrspdwk"   "yrspdwka" 
## [1681] "absddng"   "accdng"    "atnoact"   "chlrnnw"   "dclvlf"    "deaimpp"  
## [1687] "dngval"    "enjstm"    "enrglot"   "flapppl"   "flclpla"   "flrms"    
## [1693] "fltanx"    "fltbrd"    "fltpcfl"   "fltrstm"   "flttrd"    "fmlenj"   
## [1699] "fmlstrs"   "fstwttv"   "hlpoth"    "hlprtrn"   "inccmp"    "inccmpw"  
## [1705] "jbintr"    "jbstrs"    "lchshcp"   "lfcllk"    "lfwrs"     "lotsgot"  
## [1711] "lrnnew"    "lrnntlf"   "nhpftr"    "optftr"    "pactlot"   "pdwrkcr"  
## [1717] "physact"   "plinsoc"   "plprftr"   "pplahlp"   "ppllfcr"   "prhlppl"  
## [1723] "pstvms"    "rcndsrv"   "rehlppl"   "sedirlf"   "stfjb"     "stflfsf"  
## [1729] "stfsdlv"   "tmabdng"   "tmdotwa"   "tmendng"   "tmimdng"   "tnapsur"  
## [1735] "trtrsp"    "trtunf"    "uempnyr"   "wkvlorg"   "wrbknrm"   "acetalv"  
## [1741] "admaimg"   "algyplv"   "aljewlv"   "allbpe"    "allbpne"   "almuslv"  
## [1747] "alpfpe"    "alpfpne"   "alwspsc"   "blncmig"   "comnlng"   "cpimpop"  
## [1753] "ctbfsmv"   "dfegcf"    "dfegcon"   "dfeghbg"   "eimgrpc"   "eimpcnt"  
## [1759] "eimrcnt"   "fclcntr"   "gvrfgap"   "gvtrimg"   "icbrnct"   "idetalv"  
## [1765] "imacrlv"   "imbghct"   "imbleco"   "imdetbs"   "imdetmr"   "imfljob"  
## [1771] "imgclg"    "imgetn"    "imgfrnd"   "imgrpc"    "imhecop"   "imrcntr"  
## [1777] "imrsprc"   "imscrlv"   "imsetbs"   "imsetmr"   "imsmrgt"   "imtcjob"  
## [1783] "imunplv"   "imwbcrm"   "imwgdwn"   "lwdscwp"   "lwpeth"    "noimbro"  
## [1789] "pplstrd"   "qfimchr"   "qfimcmt"   "qfimedu"   "qfimfml"   "qfimlng"  
## [1795] "qfimwht"   "qfimwlt"   "qfimwsk"   "rfgawrk"   "rfgbfml"   "rfgdtcn"  
## [1801] "rfgfrpc"   "rfggvfn"   "rlgueim"   "shrrfg"    "smctmbe"   "smegbhw"  
## [1807] "smegbli"   "stimrdt"   "vrtrlg"    "cltodm"    "cltofrd"   "cltommb"  
## [1813] "cltona"    "cltonn"    "cltoptp"   "cltoref"   "cltovw"    "cnsodm"   
## [1819] "cnsofrd"   "cnsommb"   "cnsona"    "cnsonn"    "cnsoptp"   "cnsoref"  
## [1825] "cnsovw"    "discpol"   "empl"      "epaodm"    "epaofrd"   "epaommb"  
## [1831] "epaona"    "epaonn"    "epaoptp"   "epaoref"   "epaovw"    "hlpppl"   
## [1837] "hmnodm"    "hmnofrd"   "hmnommb"   "hmnona"    "hmnonn"    "hmnoptp"  
## [1843] "hmnoref"   "hmnovw"    "impapol"   "impavo"    "impfml"    "impfrds"  
## [1849] "implsrt"   "impoblw"   "impopin"   "imppol"    "imprlg"    "imprwcr"  
## [1855] "imprwct"   "imprwkc"   "impsppl"   "impvo"     "impvote"   "impwrk"   
## [1861] "othvdm"    "othvfrd"   "othvmmb"   "othvna"    "othvnn"    "othvptp"  
## [1867] "othvref"   "othvvw"    "prfodm"    "prfofrd"   "prfommb"   "prfona"   
## [1873] "prfonn"    "prfoptp"   "prforef"   "prfovw"    "prtydm"    "prtyfrd"  
## [1879] "prtymmb"   "prtyna"    "prtynn"    "prtyptp"   "prtyref"   "prtyvw"   
## [1885] "rlgodm"    "rlgofrd"   "rlgommb"   "rlgona"    "rlgonn"    "rlgoptp"  
## [1891] "rlgoref"   "rlgovw"    "sclcdm"    "sclcfrd"   "sclcmmb"   "sclcna"   
## [1897] "sclcnn"    "sclcptp"   "sclcref"   "sclcvw"    "setodm"    "setofrd"  
## [1903] "setommb"   "setona"    "setonn"    "setoptp"   "setoref"   "setovw"   
## [1909] "sptcdm"    "sptcfrd"   "sptcmmb"   "sptcna"    "sptcnn"    "sptcptp"  
## [1915] "sptcref"   "sptcvw"    "stfhwkp"   "strtbsn"   "trudm"     "trufrd"   
## [1921] "truiwkp"   "trummb"    "truna"     "trunn"     "truptp"    "truref"   
## [1927] "trusay"    "truvw"     "truwrkp"   "wkchtsk"   "wkdcsin"   "wkenvin"  
## [1933] "wrkflex"   "yrlvdae"   "advbach"   "advhach"   "advslep"   "advsthr"  
## [1939] "chsrgp"    "dcadmms"   "dcdfcwr"   "dcdisc"    "dckptrt"   "dctreql"  
## [1945] "ddprsmd"   "dsplvpr"   "follwdc"   "hltherb"   "illcure"   "mdhair"   
## [1951] "mdhappy"   "mdlswgt"   "mdmemo"    "mdsdeff"   "mdsexlf"   "pplcure"  
## [1957] "pprlydc"   "prcbach"   "prchach"   "prcslep"   "prcsthr"   "prfmddc"  
## [1963] "prfsmdc"   "psmdcpr"   "ptnrlcq"   "tmcnsdc"   "usmdprs"   "usmprse"  
## [1969] "bnkfldl"   "bsnprft"   "cmprcti"   "ctzchtx"   "ctzhlpo"   "flgvbnf"  
## [1975] "flinsr"    "flinsrw"   "fodcncl"   "frdbnft"   "frmwktg"   "icsbfm"   
## [1981] "ignrlaw"   "kptchng"   "mnyacth"   "musdocm"   "olwmsop"   "payavtx"  
## [1987] "pboafvr"   "pbofvr"    "pbofvrw"   "pyavtxw"   "rprochg"   "scbevts"  
## [1993] "scndhfl"   "slcnflw"   "slcnsfl"   "tstfnch"   "tstpboh"   "tstrprh"  
## [1999] "wrytrdh"   "acldnmr"   "adftjob"   "adllvhm"   "adlvptn"   "adpnt"    
## [2005] "advcyc"    "aftjbyc"   "ageadlt"   "agemage"   "ageoage"   "alvgptn"  
## [2011] "anvcld"    "brnmm"     "bthcld"    "evlvptn"   "evmar"     "evpdemp"  
## [2017] "fcldbrn"   "ggchld"    "iaglptn"   "iaglvmr"   "iagpnt"    "iagrtr"   
## [2023] "lvpntyr"   "lvptnyr"   "maryr"     "nbthcld"   "ngchld"    "oldfrl"   
## [2029] "oldgpnt"   "oldhlpo"   "pdempyr"   "plnftr"    "rcldnmr"   "rdvcyc"   
## [2035] "rftjbyc"   "rlvgptn"   "rnvcld"    "rpntyng"   "rspslvo"   "rwkold"   
## [2041] "svclvo"    "tochld"    "tolvpnt"   "towkht"    "tygledu"   "tyglvmr"  
## [2047] "tyglvp"    "tygpnt"    "tygrtr"    "tygsexi"   "wrinco"    "ycldbyr"  
## [2053] "ygcdbyr"   "admub"     "basinc"    "bennent"   "bnlwinc"   "brnocnt"  
## [2059] "cldcrsv"   "dcndleq"   "dfincac"   "ditxssp"   "earnpen"   "earnueb"  
## [2065] "eduunmp"   "eudcnbf"   "eusclbf"   "gvcldcr"   "gvhlthc"   "gvjbevn"  
## [2071] "gvpdlwk"   "gvslvol"   "gvslvue"   "hlthcef"   "imrccon"   "imsclbn"  
## [2077] "insfben"   "lbenent"   "lklpwcf"   "lknemny"   "lknhlcn"   "lkuemp"   
## [2083] "lvoapen"   "lvpbhlt"   "nmnybsc"   "prtsick"   "sbbsntx"   "sbcwkfm"  
## [2089] "sbenccm"   "sbeqsoc"   "sblazy"    "sblwcoa"   "sblwlka"   "sbprvpv"  
## [2095] "sbstrec"   "schtaut"   "sickdsb"   "slvpens"   "slvuemp"   "smdfslv"  
## [2101] "txadleq"   "txautef"   "txearn"    "ub20edu"   "ub20pay"   "ub20unp"  
## [2107] "ub50edu"   "ub50pay"   "ub50unp"   "ubedu"     "ubpay"     "ubspedu"  
## [2113] "ubsppay"   "ubspunp"   "ubunp"     "uemplwk"   "uentrjb"   "wrkprbf"  
## [2119] "yngfnjb"   "acq30bs"   "acq70bs"   "agdcold"   "agdcyng"   "agegrbl"  
## [2125] "ageyo30"   "agrpblg"   "buproag"   "chg1530"   "dprschg"   "dprsf70"  
## [2131] "dprso70"   "dprsy30"   "frndo70"   "frndy30"   "lkrspag"   "mbfmo70"  
## [2137] "oafl20"    "oafl70"    "p20cnec"   "p20cust"   "p70cnec"   "p70cust"  
## [2143] "p70hltb"   "predage"   "predetn"   "predsex"   "pvolwrk"   "sgr2070"  
## [2149] "srdscag"   "stsp20"    "stsp40"    "stsp70"    "suproag"   "trtbdag"  
## [2155] "twcol20"   "twcol70"   "v20adm"    "v20cntm"   "v20comp"   "v20envy"  
## [2161] "v20frnd"   "v20mrst"   "v20pity"   "v20resp"   "v70adm"    "v70cntm"  
## [2167] "v70comp"   "v70envy"   "v70frnd"   "v70mrst"   "v70pity"   "v70resp"  
## [2173] "wr20crm"   "wr20emp"   "bplcdc"    "bstln5y"   "bystlct"   "bystlwr"  
## [2179] "caplcst"   "ctfrdc"    "ctinplt"   "ctjbcnt"   "ctmtgfr"   "ctprpwr"  
## [2185] "dbctvrd"   "doplcsy"   "dpcstrb"   "flsin5y"   "gsupplc"   "insclct"  
## [2191] "insclwr"   "jdgcbrb"   "lwstrob"   "plcarcr"   "plccbrb"   "plccbrg"  
## [2197] "plccont"   "plcexdc"   "plcfrdc"   "plciplt"   "plcipvl"   "plcjbcn"  
## [2203] "plcpvcr"   "plcrgwr"   "plcrspc"   "plcstf"    "plcvcrc"   "plcvcrp"  
## [2209] "rgbrklw"   "stcbg2t"   "tmprs"     "trfoct"    "trfowr"    "troff5y"  
## [2215] "wevdct"    "widprsn"   "wraccrc"   "wraccrp"   "chpldm"    "chpldmc"  
## [2221] "chpldmi"   "ctstogv"   "cttresa"   "cttresac"  "dfprtal"   "dfprtalc" 
## [2227] "dspplvt"   "dspplvtc"  "fairelc"   "fairelcc"  "fplvdm"    "fplvdmc"  
## [2233] "fplvdmi"   "gptpelc"   "gptpelcc"  "grdfinc"   "grdfincc"  "gvcodmc"  
## [2239] "gvcodmi"   "gvctzpv"   "gvctzpvc"  "gvexpdc"   "gvexpdcc"  "gvspcdm"  
## [2245] "gvspdmc"   "gvspdmi"   "imvtctz"   "medcrgv"   "medcrgvc"  "meprinf"  
## [2251] "meprinfc"  "oppcrgv"   "oppcrgvc"  "pltavie"   "pltaviec"  "pplvdmc"  
## [2257] "pplvdmi"   "rghmgpr"   "rghmgprc"  "stpldmc"   "stpldmi"   "votedir"  
## [2263] "votedirc"  "viepol"    "wpestop"   "keydec"    "viepolc"   "wpestopc" 
## [2269] "keydecc"   "admit"     "showcv"    "accalaw"   "impdema"   "impdemb"  
## [2275] "impdemc"   "impdemd"   "impdeme"   "alcbnge"   "alcfreq"   "alcwkdy"  
## [2281] "alcwknd"   "anypacc"   "cgtsday"   "cgtsmke"   "cnfpplh"   "dosprt"   
## [2287] "dshltdk"   "dshltgp"   "dshltms"   "dshltna"   "dshltnt"   "dshltref" 
## [2293] "eatveg"    "etfruit"   "fnsdfml"   "height"    "hlpfmhr"   "hlpfmly"  
## [2299] "hltphal"   "hltphbn"   "hltphbp"   "hltphdi"   "hltphdk"   "hltphhb"  
## [2305] "hltphhc"   "hltphna"   "hltphnap"  "hltphnt"   "hltphpa"   "hltphpf"  
## [2311] "hltphref"  "hltphsc"   "hltphsd"   "hltphsh"   "hltpral"   "hltprbn"  
## [2317] "hltprbp"   "hltprca"   "hltprdi"   "hltprdk"   "hltprhb"   "hltprhc"  
## [2323] "hltprna"   "hltprnt"   "hltprpa"   "hltprpf"   "hltprref"  "hltprsc"  
## [2329] "hltprsd"   "hltprsh"   "jbexebs"   "jbexecp"   "jbexedk"   "jbexena"  
## [2335] "jbexenap"  "jbexent"   "jbexera"   "jbexeref"  "jbexevc"   "jbexevh"  
## [2341] "jbexevl"   "jbexpdk"   "jbexpmc"   "jbexpml"   "jbexpna"   "jbexpnap" 
## [2347] "jbexpnt"   "jbexpref"  "jbexpti"   "jbexpvi"   "medtrdk"   "medtrna"  
## [2353] "medtrnap"  "medtrnl"   "medtrnp"   "medtrnt"   "medtrnu"   "medtrnaa" 
## [2359] "medtroc"   "medtroth"  "medtrref"  "medtrun"   "medtrwl"   "trhltacp" 
## [2365] "trhltacu"  "trhltch"   "trhltcm"   "trhltdk"   "trhltho"   "trhltht"  
## [2371] "trhlthy"   "trhltmt"   "trhltna"   "trhltnt"   "trhltos"   "trhltpt"  
## [2377] "trhltre"   "trhltref"  "trhltsh"   "weight"    "icgndra"   "banhhap"  
## [2383] "ccgdbd"    "cflsenr"   "clmchng"   "clmthgt1"  "clmthgt2"  "elgbio"   
## [2389] "elgcoal"   "elghydr"   "elgngas"   "elgnuc"    "elgsun"    "elgwind"  
## [2395] "eneffap"   "gvsrdcc"   "inctxff"   "lklmten"   "lkredcc"   "ownrdcc"  
## [2401] "rdcenr"    "sbsrnen"   "wrdpfos"   "wrdpimp"   "wrenexp"   "wrinspw"  
## [2407] "wrntdis"   "wrpwrct"   "wrtcfl"    "wrtratc"   "btminfr"   "evfredu"  
## [2413] "evfrjob"   "frlgrsp"   "frlneti"   "frprtpl"   "fvgabc"    "grspfr"   
## [2419] "grsplet"   "grspnum"   "gvintcz"   "ifredu"    "ifrjob"    "infqbst"  
## [2425] "jstprev"   "netifr"    "netilet"   "netinum"   "occinfr"   "pcmpinj"  
## [2431] "poltran"   "ppldsrv"   "recexp"    "recgndr"   "recimg"    "recknow"  
## [2437] "recskil"   "sofrdst"   "sofrpr"    "sofrprv"   "sofrwrk"   "topinfr"  
## [2443] "wltdffr"   "iincsrc"   "acchome"   "accwrk"    "accmove"   "accoth"   
## [2449] "accnone"   "accref"    "accdk"     "accna"     "fampref"   "famadvs"  
## [2455] "fampdf"    "mcclose"   "mcinter"   "mccoord"   "mcpriv"    "mcmsinf"  
## [2461] "chldo12"   "gndro12a"  "gndro12b"  "ageo12"    "hhlio12"   "closeo12" 
## [2467] "ttmino12"  "speako12"  "como12"    "scrno12"   "phoneo12"  "livpnt"   
## [2473] "pntmofa"   "agepnt"    "hhlipnt"   "closepnt"  "ttminpnt"  "speakpnt" 
## [2479] "compnt"    "scrnpnt"   "phonepnt"  "wrkhome"   "wrklong"   "wrkresp"  
## [2485] "scrhlm"    "sclmsup"   "mansupp"   "manhlp"    "manwrkpl"  "manspeak" 
## [2491] "mancom"    "manscrn"   "manphone"  "teamfeel"  "wrkextra"  "colprop"  
## [2497] "colhlp"    "colspeak"  "colcom"    "colscrn"   "colphone"  "mcwrkhom" 
## [2503] "c19spo12"  "c19mco12"  "c19sppnt"  "c19mcpnt"  "c19whome"  "c19wplch" 
## [2509] "c19whacc"  "c19spwrk"  "c19mcwrk"  "secgrdec"  "scidecpb"  "admc19"   
## [2515] "panpriph"  "panmonpb"  "govpriph"  "govmonpb"  "panfolru"  "panclobo" 
## [2521] "panresmo"  "gvhanc19"  "gvjobc19"  "gveldc19"  "gvfamc19"  "hscopc19" 
## [2527] "gvbalc19"  "gvimpc19"  "gvconc19"  "respc19"   "reshhc19"  "hapljc19" 
## [2533] "hapirc19"  "hapwrc19"  "hapfuc19"  "hapfoc19"  "hapnoc19"  "hapnwc19" 
## [2539] "hapnpc19"  "haprec19"  "hapdkc19"  "hapnac19"  "icvacc19"  "getavc19" 
## [2545] "getnvc19"  "vdcond"    "vdovexre"  "vdtype"    "vdtpsvre"  "vdtpitre" 
## [2551] "vdtpscre"  "vdtpaure"  "vdtpvire"  "vdtpoire"  "vdtpntre"  "vdtpapre" 
## [2557] "vdtprere"  "vdtpdkre"  "vdtpnare"  "cldgng"    "enjlf"     "fltdpr"   
## [2563] "flteeff"   "fltlnl"    "fltlnla"   "fltsd"     "hrshsnt"   "hrshsnta" 
## [2569] "mnrgtjb"   "pdaprp"    "pdaprpa"   "slprl"     "smbtjob"   "smbtjoba" 
## [2575] "stfjbot"   "wrhpp"     "stfmjob"   "wrywprb"   "trdawrk"   "jbprtfp"  
## [2581] "pfmfdjb"   "pfmfdjba"  "dcsfwrk"   "dcsfwrka"  "implvdm"   "ccrdprs"  
## [2587] "ccnthum"   "wrclmch"   "yrbrn"     "wkhct"     "region"    "admdw"    
## [2593] "admge"     "iagmr"     "tygmr"     "domain"    "trstsci"   "prtvtmk"  
## [2599] "prtclmk"   "rlgdnmk"   "admrclc"   "testic34"  "testic35"  "testic36" 
## [2605] "testic37"  "testic38"  "testic39"  "testic40"  "testic41"  "testic42" 
## [2611] "rshpsmk"   "marstmk"   "edlvdmk"   "edlvpdmk"  "edlvfdmk"  "edlvmdmk" 
## [2617] "testii1"   "testii2"   "testii3"   "testii4"   "testii5"   "testii6"  
## [2623] "testii7"   "testii8"   "testii9"   "ainws"     "ainwe"     "binwe"    
## [2629] "cinwe"     "dinwe"     "finwe"     "ginwe"     "hinwe"     "iinwe"    
## [2635] "kinwe"     "vinwe"     "jinws"     "jinwe"     "scrlgblg"  "scchpldm" 
## [2641] "scwflt"    "scrwiat"   "sctmfeel"  "scrwwc"    "scclprp"   "scwsds"   
## [2647] "scwass"    "scwa5se"   "scwa11se"  "scwa59se"  "scwa89se"  "scwb30se" 
## [2653] "scwc55se"  "scwd62se"  "scwd78se"  "scwrecse"  "scwtype"   "uagent"   
## [2659] "smplresp"  "questcmp"  "intagea"   "intgndr"   "intnum"    "intlnga"  
## [2665] "rlgdemk"
colnames(portugal_data)
##    [1] "name"      "essround"  "edition"   "proddate"  "cntry"     "idno"     
##    [7] "dweight"   "pspwght"   "pweight"   "anweight"  "prob"      "stratum"  
##   [13] "psu"       "netuse"    "netusoft"  "netustm"   "nwspol"    "nwsppol"  
##   [19] "nwsptot"   "pplfair"   "pplhlp"    "ppltrst"   "rdpol"     "rdtot"    
##   [25] "tvpol"     "tvtot"     "actrolg"   "actrolga"  "badge"     "bctprd"   
##   [31] "bghtprd"   "clsprty"   "contplt"   "cptppol"   "cptppola"  "dclagr"   
##   [37] "dclaid"    "dclcrm"    "dcldef"    "dclenv"    "dclintr"   "dclmig"   
##   [43] "dclwlfr"   "dmcntov"   "dntmny"    "ecohenv"   "etapapl"   "euftf"    
##   [49] "freehms"   "gincdif"   "ginveco"   "hmsacld"   "hmsfmlsh"  "ilglpst"  
##   [55] "lawobey"   "lrscale"   "mmbprty"   "needtru"   "pbldmn"    "pbldmna"  
##   [61] "pltcare"   "pltinvt"   "polactiv"  "polcmpl"   "poldcs"    "polintr"  
##   [67] "prtdgcl"   "prtyban"   "psppipl"   "psppipla"  "psppsgv"   "psppsgva" 
##   [73] "pstplonl"  "ptcpplt"   "scnsenv"   "sgnptit"   "stfdem"    "stfeco"   
##   [79] "stfedu"    "stfgov"    "stfhlth"   "stflife"   "trstep"    "trstlgl"  
##   [85] "trstplc"   "trstplt"   "trstprl"   "trstprt"   "trstun"    "vote"     
##   [91] "wrkorg"    "wrkprty"   "imsmetn"   "imdfetn"   "impcntr"   "imbgeco"  
##   [97] "imueclt"   "imwbcnt"   "volunfp"   "lrnobed"   "loylead"   "donprty"  
##  [103] "prtvde1"   "prtvade1"  "prtvbde1"  "prtvcde1"  "prtvdde1"  "prtvede1" 
##  [109] "prtvfde1"  "prtvde2"   "prtvade2"  "prtvbde2"  "prtvcde2"  "prtvdde2" 
##  [115] "prtvede2"  "prtvfde2"  "prtvlt1"   "prtvalt1"  "prtvblt1"  "prtvclt1" 
##  [121] "prtvlt2"   "prtvalt2"  "prtvblt2"  "prtvclt2"  "prtvlt3"   "prtvalt3" 
##  [127] "prtvblt3"  "prtvclt3"  "prtvtal"   "prtvtat"   "prtvtaat"  "prtvtbat" 
##  [133] "prtvtcat"  "prtvtbe"   "prtvtabe"  "prtvtbbe"  "prtvtcbe"  "prtvtdbe" 
##  [139] "prtvtebe"  "prtvtbg"   "prtvtabg"  "prtvtbbg"  "prtvtcbg"  "prtvtdbg" 
##  [145] "prtvtebg"  "prtvtch"   "prtvtach"  "prtvtbch"  "prtvtcch"  "prtvtdch" 
##  [151] "prtvtech"  "prtvtfch"  "prtvtgch"  "prtvthch"  "prtvtcy"   "prtvtacy" 
##  [157] "prtvtbcy"  "prtvtccy"  "prtvtcz"   "prtvtacz"  "prtvtbcz"  "prtvtccz" 
##  [163] "prtvtdcz"  "prtvtecz"  "prtvtdk"   "prtvtadk"  "prtvtbdk"  "prtvtcdk" 
##  [169] "prtvtddk"  "prtvtee"   "prtvtaee"  "prtvtbee"  "prtvtcee"  "prtvtdee" 
##  [175] "prtvteee"  "prtvtfee"  "prtvtgee"  "prtvthee"  "prtvtes"   "prtvtaes" 
##  [181] "prtvtbes"  "prtvtces"  "prtvtdes"  "prtvtees"  "prtvtfes"  "prtvtfi"  
##  [187] "prtvtafi"  "prtvtbfi"  "prtvtcfi"  "prtvtdfi"  "prtvtefi"  "prtvtfr"  
##  [193] "prtvtafr"  "prtvtbfr"  "prtvtcfr"  "prtvtdfr"  "prtvtefr"  "prtvtgb"  
##  [199] "prtvtagb"  "prtvtbgb"  "prtvtcgb"  "prtvtdgb"  "prtvtgr"   "prtvtagr" 
##  [205] "prtvtbgr"  "prtvtcgr"  "prtvtdgr"  "prtvthr"   "prtvtahr"  "prtvtbhr" 
##  [211] "prtvthu"   "prtvtahu"  "prtvtbhu"  "prtvtchu"  "prtvtdhu"  "prtvtehu" 
##  [217] "prtvtfhu"  "prtvtghu"  "prtvtie"   "prtvtaie"  "prtvtbie"  "prtvtcie" 
##  [223] "prtvtdie"  "prtvtil"   "prtvtail"  "prtvtbil"  "prtvtcil"  "prtvtdil" 
##  [229] "prtvtis"   "prtvtais"  "prtvtbis"  "prtvtcis"  "prtvtdis"  "prtvtit"  
##  [235] "prtvtait"  "prtvtbit"  "prtvtcit"  "prtvtdit"  "prtvtlu"   "prtvtlv"  
##  [241] "prtvtalv"  "prtvtnl"   "prtvtanl"  "prtvtbnl"  "prtvtcnl"  "prtvtdnl" 
##  [247] "prtvtenl"  "prtvtfnl"  "prtvtgnl"  "prtvthnl"  "prtvtno"   "prtvtano" 
##  [253] "prtvtbno"  "prtvtpl"   "prtvtapl"  "prtvtbpl"  "prtvtcpl"  "prtvtdpl" 
##  [259] "prtvtepl"  "prtvtpt"   "prtvtapt"  "prtvtbpt"  "prtvtcpt"  "prtvtdpt" 
##  [265] "prtvtro"   "prtvtaro"  "prtvtru"   "prtvtaru"  "prtvtbru"  "prtvtcru" 
##  [271] "prtvtdru"  "prtvtse"   "prtvtase"  "prtvtbse"  "prtvtcse"  "prtvtdse" 
##  [277] "prtvtsi"   "prtvtasi"  "prtvtbsi"  "prtvtcsi"  "prtvtdsi"  "prtvtesi" 
##  [283] "prtvtfsi"  "prtvtsk"   "prtvtask"  "prtvtbsk"  "prtvtcsk"  "prtvtdsk" 
##  [289] "prtvtesk"  "prtvttr"   "prtvtatr"  "prtvtua"   "prtvtaua"  "prtvtbua" 
##  [295] "prtvtcua"  "prtvtxk"   "prtvtrs"   "prtvtars"  "prtvtme"   "prtvtame" 
##  [301] "prtclal"   "prtclat"   "prtclaat"  "prtclcat"  "prtcldat"  "prtclbe"  
##  [307] "prtclabe"  "prtclbbe"  "prtclcbe"  "prtcldbe"  "prtclebe"  "prtclbg"  
##  [313] "prtclabg"  "prtclbbg"  "prtclcbg"  "prtcldbg"  "prtclebg"  "prtclch"  
##  [319] "prtclach"  "prtclbch"  "prtclcch"  "prtcldch"  "prtclech"  "prtclfch" 
##  [325] "prtclgch"  "prtclhch"  "prtclsk"   "prtclask"  "prtclbsk"  "prtclcsk" 
##  [331] "prtcldsk"  "prtclesk"  "prtclcy"   "prtclacy"  "prtclbcy"  "prtclccy" 
##  [337] "prtclcz"   "prtclacz"  "prtclbcz"  "prtclccz"  "prtcldcz"  "prtclecz" 
##  [343] "prtclde"   "prtclade"  "prtclbde"  "prtclcde"  "prtcldde"  "prtclede" 
##  [349] "prtclfde"  "prtcldk"   "prtcladk"  "prtclbdk"  "prtclcdk"  "prtclddk" 
##  [355] "prtclee"   "prtclaee"  "prtclbee"  "prtclcee"  "prtcldee"  "prtcleee" 
##  [361] "prtclfee"  "prtclgee"  "prtclhee"  "prtcles"   "prtclaes"  "prtclbes" 
##  [367] "prtclces"  "prtcldes"  "prtclees"  "prtclfes"  "prtclges"  "prtclfi"  
##  [373] "prtclafi"  "prtclbfi"  "prtclcfi"  "prtcldfi"  "prtclefi"  "prtclffi" 
##  [379] "prtclfr"   "prtclafr"  "prtclbfr"  "prtclcfr"  "prtcldfr"  "prtclefr" 
##  [385] "prtclffr"  "prtclgb"   "prtclagb"  "prtclbgb"  "prtclcgb"  "prtcldgb" 
##  [391] "prtclgr"   "prtclagr"  "prtclbgr"  "prtclcgr"  "prtcldgr"  "prtclhr"  
##  [397] "prtclahr"  "prtclbhr"  "prtclhu"   "prtclahu"  "prtclbhu"  "prtclchu" 
##  [403] "prtcldhu"  "prtclehu"  "prtclfhu"  "prtclghu"  "prtclhhu"  "prtclie"  
##  [409] "prtclaie"  "prtclbie"  "prtclcie"  "prtcldie"  "prtcleie"  "prtclfie" 
##  [415] "prtclil"   "prtclail"  "prtclbil"  "prtclcil"  "prtcldil"  "prtcleil" 
##  [421] "prtclis"   "prtclais"  "prtclbis"  "prtclcis"  "prtcldis"  "prtclit"  
##  [427] "prtclait"  "prtclbit"  "prtclcit"  "prtcldit"  "prtcleit"  "prtcllt"  
##  [433] "prtclalt"  "prtclblt"  "prtclclt"  "prtcllu"   "prtcllv"   "prtclalv" 
##  [439] "prtclnl"   "prtclanl"  "prtclbnl"  "prtclcnl"  "prtcldnl"  "prtclenl" 
##  [445] "prtclfnl"  "prtclgnl"  "prtclno"   "prtclano"  "prtclbno"  "prtclpl"  
##  [451] "prtclapl"  "prtclbpl"  "prtclcpl"  "prtcldpl"  "prtclepl"  "prtclfpl" 
##  [457] "prtclgpl"  "prtclhpl"  "prtclipl"  "prtclpt"   "prtclapt"  "prtclbpt" 
##  [463] "prtclcpt"  "prtcldpt"  "prtclept"  "prtclfpt"  "prtclro"   "prtclaro" 
##  [469] "prtclru"   "prtclaru"  "prtclbru"  "prtclcru"  "prtcldru"  "prtclse"  
##  [475] "prtclase"  "prtclbse"  "prtclcse"  "prtcldse"  "prtclsi"   "prtclasi" 
##  [481] "prtclbsi"  "prtclcsi"  "prtcldsi"  "prtclesi"  "prtclfsi"  "prtcltr"  
##  [487] "prtclatr"  "prtclua"   "prtclaua"  "prtclbua"  "prtclcua"  "prtcldua" 
##  [493] "prtclxk"   "prtclrs"   "prtclars"  "prtclme"   "prtclame"  "prtmbat"  
##  [499] "prtmbaat"  "prtmbbe"   "prtmbabe"  "prtmbbbe"  "prtmbcbe"  "prtmbbg"  
##  [505] "prtmbabg"  "prtmbbbg"  "prtmbch"   "prtmbach"  "prtmbbch"  "prtmbcch" 
##  [511] "prtmbcy"   "prtmbcz"   "prtmbacz"  "prtmbbcz"  "prtmbde"   "prtmbade" 
##  [517] "prtmbbde"  "prtmbcde"  "prtmbdk"   "prtmbadk"  "prtmbbdk"  "prtmbee"  
##  [523] "prtmbaee"  "prtmbbee"  "prtmbcee"  "prtmbes"   "prtmbaes"  "prtmbbes" 
##  [529] "prtmbfi"   "prtmbafi"  "prtmbbfi"  "prtmbfr"   "prtmbafr"  "prtmbbfr" 
##  [535] "prtmbcfr"  "prtmbgb"   "prtmbagb"  "prtmbgr"   "prtmbagr"  "prtmbbgr" 
##  [541] "prtmbcgr"  "prtmbhr"   "prtmbhu"   "prtmbahu"  "prtmbbhu"  "prtmbchu" 
##  [547] "prtmbie"   "prtmbaie"  "prtmbbie"  "prtmbil"   "prtmbail"  "prtmbbil" 
##  [553] "prtmbis"   "prtmbit"   "prtmbait"  "prtmblt"   "prtmblu"   "prtmblv"  
##  [559] "prtmbnl"   "prtmbanl"  "prtmbbnl"  "prtmbcnl"  "prtmbno"   "prtmbano" 
##  [565] "prtmbpl"   "prtmbapl"  "prtmbbpl"  "prtmbcpl"  "prtmbdpl"  "prtmbpt"  
##  [571] "prtmbapt"  "prtmbbpt"  "prtmbro"   "prtmbaro"  "prtmbru"   "prtmbaru" 
##  [577] "prtmbbru"  "prtmbse"   "prtmbase"  "prtmbsi"   "prtmbasi"  "prtmbbsi" 
##  [583] "prtmbcsi"  "prtmbsk"   "prtmbask"  "prtmbbsk"  "prtmbtr"   "prtmbatr" 
##  [589] "prtmbua"   "prtmbaua"  "prtmbbua"  "prtmbcua"  "aesfdrk"   "atchctr"  
##  [595] "atcherp"   "blgetmg"   "brghmef"   "brghmwr"   "brncntr"   "cntbrth"  
##  [601] "cntbrtha"  "cntbrthb"  "cntbrthc"  "cntbrthd"  "crmvct"    "crvctef"  
##  [607] "crvctwr"   "ctzcntr"   "ctzship"   "ctzshipa"  "ctzshipb"  "ctzshipc" 
##  [613] "ctzshipd"  "dscrage"   "dscrdk"    "dscrdsb"   "dscretn"   "dscrgnd"  
##  [619] "dscrgrp"   "dscrlng"   "dscrna"    "dscrnap"   "dscrntn"   "dscroth"  
##  [625] "dscrrce"   "dscrref"   "dscrrlg"   "dscrsex"   "facntn"    "facntr"   
##  [631] "fbrncnt"   "fbrncnta"  "fbrncntb"  "fbrncntc"  "happy"     "health"   
##  [637] "hlthhmp"   "inmdisc"   "inprdsc"   "livecntr"  "livecnta"  "lnghoma"  
##  [643] "lnghom1"   "lnghomb"   "lnghom2"   "mbrncnt"   "mbrncnta"  "mbrncntb" 
##  [649] "mbrncntc"  "mmbrn"     "mocntn"    "mocntr"    "pray"      "rlgatnd"  
##  [655] "rlgblg"    "rlgblge"   "rlgdgr"    "rlgdnm"    "rlgdnme"   "sclact"   
##  [661] "sclmeet"   "trrcnyr"   "trrenyr"   "trrprsn"   "trrtort"   "feethngr" 
##  [667] "vteurmmb"  "vteubcmb"  "vteumbgb"  "rlgdnal"   "rlgdnat"   "rlgdnbat" 
##  [673] "rlgdnbe"   "rlgdnch"   "rlgdnach"  "rlgdncy"   "rlgdnacy"  "rlgdnde"  
##  [679] "rlgdnade"  "rlgdndk"   "rlgdnfi"   "rlgdnafi"  "rlgdngb"   "rlgdngr"  
##  [685] "rlgdnagr"  "rlgdnhu"   "rlgdnie"   "rlgdnil"   "rlgdnis"   "rlgdnais" 
##  [691] "rlgdnlt"   "rlgdnlv"   "rlgdnnl"   "rlgdnanl"  "rlgdnno"   "rlgdnpl"  
##  [697] "rlgdnapl"  "rlgdnbpl"  "rlgdnpt"   "rlgdnro"   "rlgdnru"   "rlgdnaru" 
##  [703] "rlgdnse"   "rlgdnase"  "rlgdnsi"   "rlgdnsk"   "rlgdnask"  "rlgdnbsk" 
##  [709] "rlgdnua"   "rlgdnrs"   "rlgdme"    "rlgdeal"   "rlgdeat"   "rlgdebat" 
##  [715] "rlgdebe"   "rlgdech"   "rlgdeach"  "rlgdecy"   "rlgdeacy"  "rlgdede"  
##  [721] "rlgdeade"  "rlgdedk"   "rlgdefi"   "rlgdeafi"  "rlgdegb"   "rlgdegr"  
##  [727] "rlgdeagr"  "rlgdehu"   "rlgdeie"   "rlgdeil"   "rlgdeis"   "rlgdeais" 
##  [733] "rlgdelt"   "rlgdelv"   "rlgdenl"   "rlgdeanl"  "rlgdeno"   "rlgdepl"  
##  [739] "rlgdeapl"  "rlgdebpl"  "rlgdept"   "rlgdero"   "rlgderu"   "rlgdearu" 
##  [745] "rlgdese"   "rlgdease"  "rlgdesi"   "rlgdesk"   "rlgdeask"  "rlgdebsk" 
##  [751] "rlgdeua"   "rlgders"   "rlgdeme"   "hhmmb"     "gndr"      "gndr2"    
##  [757] "gndr3"     "gndr4"     "gndr5"     "gndr6"     "gndr7"     "gndr8"    
##  [763] "gndr9"     "gndr10"    "gndr11"    "gndr12"    "gndr13"    "gndr14"   
##  [769] "gndr15"    "gndr16"    "gndr17"    "gndr18"    "gndr19"    "gndr20"   
##  [775] "gndr21"    "gndr22"    "gndr23"    "gndr24"    "rship2"    "rshipa2"  
##  [781] "rship3"    "rshipa3"   "rship4"    "rshipa4"   "rship5"    "rshipa5"  
##  [787] "rship6"    "rshipa6"   "rship7"    "rshipa7"   "rship8"    "rshipa8"  
##  [793] "rship9"    "rshipa9"   "rship10"   "rshipa10"  "rship11"   "rshipa11" 
##  [799] "rship12"   "rshipa12"  "rship13"   "rshipa13"  "rship14"   "rshipa14" 
##  [805] "rship15"   "rshipa15"  "rshipa16"  "rshipa17"  "rshipa18"  "rshipa19" 
##  [811] "rshipa20"  "rshipa21"  "rshipa22"  "rshipa23"  "rshipa24"  "yrbrn2"   
##  [817] "yrbrn3"    "yrbrn4"    "yrbrn5"    "yrbrn6"    "yrbrn7"    "yrbrn8"   
##  [823] "yrbrn9"    "yrbrn10"   "yrbrn11"   "yrbrn12"   "yrbrn13"   "yrbrn14"  
##  [829] "yrbrn15"   "yrbrn16"   "yrbrn17"   "yrbrn18"   "yrbrn19"   "yrbrn20"  
##  [835] "yrbrn21"   "yrbrn22"   "yrbrn23"   "yrbrn24"   "age"       "agea"     
##  [841] "ablrtr"    "anctry1"   "anctry2"   "atncrse"   "brwmny"    "chldhhe"  
##  [847] "chldhm"    "cmsrv"     "cmsrvp"    "crpdwk"    "crpdwkp"   "dngdk"    
##  [853] "dngdkp"    "dngna"     "dngnap"    "dngnapp"   "dngoth"    "dngothp"  
##  [859] "dngref"    "dngrefp"   "domicil"   "dsbld"     "dsbldp"    "dvrcdev"  
##  [865] "dvrcdeva"  "edctn"     "edctnp"    "edufld"    "edulvla"   "edulvlb"  
##  [871] "edulvlfa"  "edulvlfb"  "edulvlma"  "edulvlmb"  "edulvlpa"  "edulvlpb" 
##  [877] "eduyrs"    "eisced"    "eiscedf"   "eiscedm"   "eiscedp"   "emplno"   
##  [883] "emplnof"   "emplnom"   "emplnop"   "emplrel"   "emprelp"   "emprf14"  
##  [889] "emprm14"   "estsz"     "fxltph"    "faalv"     "hhmodwl"   "hincfel"  
##  [895] "hincsrc"   "hincsrca"  "hinctnt"   "hinctnta"  "hswrk"     "hswrkp"   
##  [901] "iccohbt"   "icmnart"   "icmsw"     "icomdng"   "icomdnp"   "icpdwk2"  
##  [907] "icpdwrk"   "icppdwk"   "icwhct"    "inttph"    "ioactp"    "iorgact"  
##  [913] "isco08"    "isco08p"   "iscoco"    "iscocop"   "jbspv"     "jbspvf"   
##  [919] "jbspvm"    "jbspvp"    "lvghw"     "lvghwa"    "lvgoptn"   "lvgptn"   
##  [925] "lvgptna"   "lvgptne"   "lvgptnea"  "mainact"   "marital"   "maritala" 
##  [931] "maritalb"  "marsts"    "mbltph"    "mbtru"     "mnactic"   "mnactp"   
##  [937] "moalv"     "njbspv"    "njbspvp"   "occf14"    "occf14a"   "occf14b"  
##  [943] "occf14ie"  "occm14"    "occm14a"   "occm14b"   "occm14ie"  "orgwrk"   
##  [949] "partner"   "pdjobev"   "pdjobyr"   "pdwrk"     "pdwrkp"    "rmhhus"   
##  [955] "rshpsts"   "rtdsbyr"   "rtrd"      "rtrdp"     "tporgwk"   "uemp12m"  
##  [961] "uemp3m"    "uemp5yr"   "uempla"    "uemplap"   "uempli"    "uemplip"  
##  [967] "wkdcorp"   "wkhtot"    "wkhtotp"   "wrkac6m"   "wrkctr"    "wrkctra"  
##  [973] "wrkctrhu"  "wkdcorg"   "wkdcorga"  "nacer1"    "nacer11"   "nacer2"   
##  [979] "edagegb"   "edlvat"    "edlveat"   "eduat1"    "eduat2"    "edlvbe"   
##  [985] "edlvabe"   "edlvdbe"   "edlvebe"   "edlvbg"    "edlvdbg"   "edlvebg"  
##  [991] "edlvch"    "edlvach"   "edlvbch"   "edlvcch"   "edlvdch"   "edlvcy"   
##  [997] "edlvacy"   "edlvdcy"   "edlvecy"   "edlvgcy"   "edlvcz"    "edlvdcz"  
## [1003] "edlvdal"   "edlvde"    "edlvade"   "edlvdfi"   "edlvdis"   "edlvdk"   
## [1009] "edlvadk"   "edlvddk"   "edlvdxk"   "edlvee"    "edlvaee"   "edlvbee"  
## [1015] "edlvdee"   "edlves"    "edlvaes"   "edlvdes"   "edlvees"   "edlvges"  
## [1021] "edlvhes"   "edlvfr"    "edlvafr"   "edlvbfr"   "edlvdfr"   "edlvgb"   
## [1027] "edlvagb"   "edlvgr"    "edlvagr"   "edlvdgr"   "edlvegr"   "edlvhr"   
## [1033] "edlvdhr"   "edlvehr"   "edlvhu"    "edlvahu"   "edlvbhu"   "edlvdhu"  
## [1039] "edlvdahu"  "edlvie"    "edlvaie"   "edlvbie"   "edlvdie"   "edlvil"   
## [1045] "edlvail"   "eduil1"    "eduail1"   "edubil1"   "edlvit"    "edlvait"  
## [1051] "edlvdit"   "edlveit"   "edlvlt"    "edlvdlt"   "edlvlu"    "edlvlv"   
## [1057] "edlvdlv"   "edlvelv"   "edlvnl"    "edlvdnl"   "edlvenl"   "edlvno"   
## [1063] "edlvano"   "edlvdno"   "edlveno"   "edlvpl"    "edlvapl"   "edlvbpl"  
## [1069] "edlvdpl"   "edlvepl"   "edlvgpl"   "edlvhpl"   "edlvpt"    "edlvapt"  
## [1075] "edlvbpt"   "edlvdpt"   "edlvro"    "edlvru"    "edlvdru"   "edlvse"   
## [1081] "edlvase"   "edlvdse"   "edlvsi"    "edlvasi"   "edlvdsi"   "edlvesi"  
## [1087] "edlvsk"    "edlvask"   "edlvdsk"   "edlvtr"    "edlvua"    "edlvaua"  
## [1093] "edlvdua"   "edude1"    "eduade1"   "edubde1"   "educde1"   "edude2"   
## [1099] "eduade2"   "edubde2"   "edude3"    "eduade3"   "edugb2"    "eduagb2"  
## [1105] "edubgb2"   "eduil2"    "eduail2"   "edupl2"    "eduyrpl"   "edugb1"   
## [1111] "eduagb1"   "edubgb1"   "educgb1"   "edlvdrs"   "edlvdme"   "edagepgb" 
## [1117] "edlvpat"   "edlvpeat"  "edlvpfat"  "edlvpbe"   "edlvpdbe"  "edlvpebe" 
## [1123] "edlvpch"   "edlvpdch"  "edlvpcy"   "edlvpdcy"  "edlvpecy"  "edlvpgcy" 
## [1129] "edlvpcz"   "edlvpdcz"  "edlvpdbg"  "edlvpebg"  "edlvpdfi"  "edlvpdis" 
## [1135] "edlvpdit"  "edlvpeit"  "edlvpdk"   "edlvpddk"  "edlvpdxk"  "edlvpee"  
## [1141] "edlvpdee"  "edlvpes"   "edlvpdes"  "edlvpees"  "edlvpfes"  "edlvpges" 
## [1147] "edlvpfr"   "edlvpdfr"  "edlvpgb"   "edlvpgr"   "edlvpdgr"  "edlvpegr" 
## [1153] "edlvphr"   "edlvpdhr"  "edlvpehr"  "edlvphu"   "edlvpdhu"  "edlvpdahu"
## [1159] "edlvpie"   "edlvpdie"  "edlvpil"   "edupil1"   "edupail1"  "edupbil1" 
## [1165] "edlvplt"   "edlvpdlt"  "edlvplv"   "edlvpdlv"  "edlvpelv"  "edlvpnl"  
## [1171] "edlvpdnl"  "edlvpenl"  "edlvpno"   "edlvpdno"  "edlvpeno"  "edlvppl"  
## [1177] "edlvpdpl"  "edlvpepl"  "edlvpfpl"  "edlvpgpl"  "edlvppt"   "edlvpdpt" 
## [1183] "edlvpro"   "edlvpru"   "edlvpdru"  "edlvpse"   "edlvpdse"  "edlvpsi"  
## [1189] "edlvpdsi"  "edlvpesi"  "edlvpsk"   "edlvpdsk"  "edlvptr"   "edlvpua"  
## [1195] "edlvpdua"  "edupde1"   "edupade1"  "edupbde1"  "edupcde1"  "edupde2"  
## [1201] "edupade2"  "edupbde2"  "edupde3"   "edupade3"  "edupgb2"   "edupagb2" 
## [1207] "edupbgb2"  "edupil2"   "edupail2"  "eduppl2"   "eduyrppl"  "edlvpdal" 
## [1213] "edupgb1"   "edupagb1"  "edupbgb1"  "edupcgb1"  "edlvpdrs"  "edlvpdme" 
## [1219] "edlvfdal"  "edlvfat"   "edlvfeat"  "edlvfbe"   "edlvfdbe"  "edlvfebe" 
## [1225] "edlvfdbg"  "edlvfebg"  "edlvfhr"   "edlvfdhr"  "edlvfehr"  "edlvfcy"  
## [1231] "edlvfdcy"  "edlvfecy"  "edlvfgcy"  "edlvfcz"   "edlvfdcz"  "edlvfdk"  
## [1237] "edlvfddk"  "edlvfee"   "edlvfdee"  "edlvfdfi"  "edlvffr"   "edlvfdfr" 
## [1243] "edufde1"   "edufade1"  "edufbde1"  "edufde2"   "edufade2"  "edufde3"  
## [1249] "edufade3"  "edlvfgr"   "edlvfdgr"  "edlvfegr"  "edlvfhu"   "edlvfdhu" 
## [1255] "edlvfdahu" "edlvfdis"  "edlvfie"   "edlvfdie"  "edlvfil"   "edufil1"  
## [1261] "edufail1"  "edufbil1"  "edufil2"   "edufail2"  "edlvfdit"  "edlvfeit" 
## [1267] "edlvfdxk"  "edlvflv"   "edlvfdlv"  "edlvflt"   "edlvfdlt"  "edlvfnl"  
## [1273] "edlvfdnl"  "edlvfenl"  "edlvfno"   "edlvfdno"  "edlvfeno"  "edlvfpl"  
## [1279] "edlvfdpl"  "edlvfepl"  "edlvffpl"  "edlvfpt"   "edlvfdpt"  "edlvfro"  
## [1285] "edlvfru"   "edlvfdru"  "edlvfsk"   "edlvfdsk"  "edlvfsi"   "edlvfdsi" 
## [1291] "edlvfesi"  "edlvfes"   "edlvfdes"  "edlvfees"  "edlvffes"  "edlvfse"  
## [1297] "edlvfdse"  "edlvfch"   "edlvfdch"  "edlvftr"   "edlvfua"   "edlvfdua" 
## [1303] "edlvfgb"   "edufgb1"   "edufagb1"  "edufbgb1"  "edufcgb1"  "edufgb2"  
## [1309] "edufagb2"  "edufbgb2"  "edagefgb"  "edlvfdrs"  "edlvfdme"  "edagemgb" 
## [1315] "edlvmat"   "edlvmeat"  "edlvmbe"   "edlvmdbe"  "edlvmebe"  "edlvmch"  
## [1321] "edlvmdch"  "edlvmcy"   "edlvmdcy"  "edlvmecy"  "edlvmgcy"  "edlvmcz"  
## [1327] "edlvmdcz"  "edlvmdal"  "edlvmdfi"  "edlvmdis"  "edlvmdit"  "edlvmeit" 
## [1333] "edlvmdk"   "edlvmddk"  "edlvmdxk"  "edlvmee"   "edlvmdee"  "edlvmes"  
## [1339] "edlvmdes"  "edlvmees"  "edlvmfes"  "edlvmfr"   "edlvmdfr"  "edlvmgb"  
## [1345] "edlvmgr"   "edlvmdgr"  "edlvmegr"  "edlvmhr"   "edlvmdhr"  "edlvmehr" 
## [1351] "edlvmhu"   "edlvmdhu"  "edlvmdahu" "edlvmie"   "edlvmdie"  "edlvmil"  
## [1357] "edumil1"   "edumail1"  "edumbil1"  "edlvmlt"   "edlvmdlt"  "edlvmlv"  
## [1363] "edlvmdlv"  "edlvmnl"   "edlvmdnl"  "edlvmenl"  "edlvmno"   "edlvmdno" 
## [1369] "edlvmeno"  "edlvmpl"   "edlvmdpl"  "edlvmepl"  "edlvmfpl"  "edlvmpt"  
## [1375] "edlvmdpt"  "edlvmro"   "edlvmru"   "edlvmdru"  "edlvmse"   "edlvmdse" 
## [1381] "edlvmsi"   "edlvmdsi"  "edlvmesi"  "edlvmsk"   "edlvmdsk"  "edlvmtr"  
## [1387] "edlvmua"   "edlvmdua"  "edumde1"   "edumade1"  "edumbde1"  "edumde2"  
## [1393] "edumade2"  "edumde3"   "edumade3"  "edumgb1"   "edumagb1"  "edumbgb1" 
## [1399] "edumcgb1"  "edumgb2"   "edumagb2"  "edumbgb2"  "edumil2"   "edumail2" 
## [1405] "edlvmdbg"  "edlvmebg"  "edlvmdrs"  "edlvmdme"  "marstfi"   "marstgb"  
## [1411] "marstie"   "martlfr"   "rshpsfi"   "rshpsgb"   "impdiff"   "impenv"   
## [1417] "impfree"   "impfun"    "imprich"   "impsafe"   "imptrad"   "ipadvnt"  
## [1423] "ipbhprp"   "ipcrtiv"   "ipeqopt"   "ipfrule"   "ipgdtim"   "iphlppl"  
## [1429] "iplylfr"   "ipmodst"   "iprspot"   "ipshabt"   "ipstrgv"   "ipsuces"  
## [1435] "ipudrst"   "regionat"  "regionbe"  "regionbg"  "regionch"  "regioach" 
## [1441] "regioncy"  "regioncz"  "regioacz"  "regionde"  "regiondk"  "regioadk" 
## [1447] "regionee"  "regiones"  "regioaes"  "regionfi"  "regioafi"  "regionfr" 
## [1453] "regiongb"  "regiongr"  "regioagr"  "regionhr"  "regionhu"  "regionie" 
## [1459] "regioaie"  "regiobie"  "regionil"  "regionis"  "regionit"  "regionlt" 
## [1465] "regionlu"  "regionlv"  "regionnl"  "regionno"  "regionpl"  "regionpt" 
## [1471] "regioapt"  "regionro"  "regionru"  "regionse"  "regionsi"  "regionsk" 
## [1477] "regiontr"  "regionua"  "regunit"   "inwtm"     "inwshh"    "inwsmm"   
## [1483] "inwyys"    "inwdd"     "inwdds"    "inwmm"     "inwmms"    "inwyr"    
## [1489] "inwyye"    "inwds"     "inwehh"    "inwemm"    "inwdde"    "inwmme"   
## [1495] "inwde"     "spltadm"   "spltadma"  "spltadmb"  "spltadmc"  "spltadmd" 
## [1501] "spltadme"  "spltadmf"  "supqad1"   "supqadm"   "supqad2"   "supqdd"   
## [1507] "supqmm"    "supqyr"    "intewde"   "mode"      "actvgrs"   "agertr"   
## [1513] "bseftwk"   "bsmw"      "chdohwk"   "cld12hh"   "cldcrmr"   "cldcrot"  
## [1519] "cldnhh"    "cldnhhd"   "cldnhhg"   "cldnhhn"   "clfncrc"   "clfncsp"  
## [1525] "clhwkrc"   "clhwksp"   "clmrlx"    "clnhhbo"   "clnhhby"   "clnhhyb"  
## [1531] "crspce"    "cuthheq"   "dfcnswk"   "dfcnswka"  "dsdclve"   "dsgrhwk"  
## [1537] "dsgrmny"   "dsgrmnya"  "dsgrpwk"   "edul12m"   "edupdem"   "emplrjb"  
## [1543] "flthmcc"   "fmpvtst"   "frshrst"   "frsptjb"   "fstjbyr"   "fthcncr"  
## [1549] "gdsprt"    "grspay"    "grspaya"   "gvprppv"   "gwdvhwk"   "gwhhprc"  
## [1555] "hlpcowk"   "hlthrwk"   "hmeqphw"   "hwkmono"   "hwkpwd1"   "hwkpwd2"  
## [1561] "hwkpwdp"   "hwkpwe1"   "hwkpwe2"   "hwkpwep"   "hwkstrs"   "hwktwd1"  
## [1567] "hwktwd2"   "hwktwe1"   "hwktwe2"   "hwwkhs"    "icag45y"   "icago45"  
## [1573] "icagu70"   "icagu70a"  "icchld"    "iccldnh"   "icempl"    "icemplr"  
## [1579] "icgndr"    "ichwk1"    "ichwk2"    "icmnac"    "icmnact"   "icnopfm"  
## [1585] "icnopfma"  "icpart1"   "icpart2"   "icpart3"   "icptn"     "icptnwk"  
## [1591] "icptnwka"  "indisdc"   "ipjbhin"   "ipjbini"   "ipjbprm"   "ipjbscr"  
## [1597] "ipjbtro"   "ipjbwfm"   "jbcoedu"   "jbedyrs"   "jblrn"     "jbrqlrn"  
## [1603] "jbscr"     "jbtmppm"   "jbtsktm"   "lfintr"    "linwk3y"   "lkafohh"  
## [1609] "lscjb3y"   "lsintjb"   "mlohinc"   "mnrsefw"   "mnrsphm"   "nbsrsp"   
## [1615] "netpay"    "nevdnjb"   "npemp3y"   "oprtad"    "orgfd3y"   "payprd"   
## [1621] "payprda"   "pdjbndm"   "phwwkhs"   "plnchld"   "pphincr"   "ppwwkp"   
## [1627] "prmpls"    "prntghr"   "prrfmly"   "pthcncr"   "ptmhmcc"   "ptnwkwe"  
## [1633] "puemp3y"   "pwkhsch"   "quclss"    "rdpay3y"   "rpljbde"   "rsnlvem"  
## [1639] "rtryr"     "scrsefw"   "smblvjb"   "stdask"    "stdhrsw"   "stdlvl"   
## [1645] "stdmcdo"   "stdtruf"   "stpvtfm"   "tchints"   "tchlcrt"   "tchlp"    
## [1651] "tchtruf"   "tmtowrk"   "tngdohm"   "trndnjb"   "truinwk"   "uemp3y"   
## [1657] "updhlrl"   "updhlrp"   "useknwl"   "vrtywrk"   "wgdpeft"   "wkdcpce"  
## [1663] "wkengtp"   "wkhsch"    "wkjbndm"   "wkovrtm"   "wkovtmp"   "wkshr3y"  
## [1669] "wmcpwrk"   "wntrtr"    "wpmwkcp"   "wrkengt"   "wrkhrd"    "wrkspv"   
## [1675] "wrkwe"     "yrcremp"   "yrlvptn"   "yrskdwk"   "yrspdwk"   "yrspdwka" 
## [1681] "absddng"   "accdng"    "atnoact"   "chlrnnw"   "dclvlf"    "deaimpp"  
## [1687] "dngval"    "enjstm"    "enrglot"   "flapppl"   "flclpla"   "flrms"    
## [1693] "fltanx"    "fltbrd"    "fltpcfl"   "fltrstm"   "flttrd"    "fmlenj"   
## [1699] "fmlstrs"   "fstwttv"   "hlpoth"    "hlprtrn"   "inccmp"    "inccmpw"  
## [1705] "jbintr"    "jbstrs"    "lchshcp"   "lfcllk"    "lfwrs"     "lotsgot"  
## [1711] "lrnnew"    "lrnntlf"   "nhpftr"    "optftr"    "pactlot"   "pdwrkcr"  
## [1717] "physact"   "plinsoc"   "plprftr"   "pplahlp"   "ppllfcr"   "prhlppl"  
## [1723] "pstvms"    "rcndsrv"   "rehlppl"   "sedirlf"   "stfjb"     "stflfsf"  
## [1729] "stfsdlv"   "tmabdng"   "tmdotwa"   "tmendng"   "tmimdng"   "tnapsur"  
## [1735] "trtrsp"    "trtunf"    "uempnyr"   "wkvlorg"   "wrbknrm"   "acetalv"  
## [1741] "admaimg"   "algyplv"   "aljewlv"   "allbpe"    "allbpne"   "almuslv"  
## [1747] "alpfpe"    "alpfpne"   "alwspsc"   "blncmig"   "comnlng"   "cpimpop"  
## [1753] "ctbfsmv"   "dfegcf"    "dfegcon"   "dfeghbg"   "eimgrpc"   "eimpcnt"  
## [1759] "eimrcnt"   "fclcntr"   "gvrfgap"   "gvtrimg"   "icbrnct"   "idetalv"  
## [1765] "imacrlv"   "imbghct"   "imbleco"   "imdetbs"   "imdetmr"   "imfljob"  
## [1771] "imgclg"    "imgetn"    "imgfrnd"   "imgrpc"    "imhecop"   "imrcntr"  
## [1777] "imrsprc"   "imscrlv"   "imsetbs"   "imsetmr"   "imsmrgt"   "imtcjob"  
## [1783] "imunplv"   "imwbcrm"   "imwgdwn"   "lwdscwp"   "lwpeth"    "noimbro"  
## [1789] "pplstrd"   "qfimchr"   "qfimcmt"   "qfimedu"   "qfimfml"   "qfimlng"  
## [1795] "qfimwht"   "qfimwlt"   "qfimwsk"   "rfgawrk"   "rfgbfml"   "rfgdtcn"  
## [1801] "rfgfrpc"   "rfggvfn"   "rlgueim"   "shrrfg"    "smctmbe"   "smegbhw"  
## [1807] "smegbli"   "stimrdt"   "vrtrlg"    "cltodm"    "cltofrd"   "cltommb"  
## [1813] "cltona"    "cltonn"    "cltoptp"   "cltoref"   "cltovw"    "cnsodm"   
## [1819] "cnsofrd"   "cnsommb"   "cnsona"    "cnsonn"    "cnsoptp"   "cnsoref"  
## [1825] "cnsovw"    "discpol"   "empl"      "epaodm"    "epaofrd"   "epaommb"  
## [1831] "epaona"    "epaonn"    "epaoptp"   "epaoref"   "epaovw"    "hlpppl"   
## [1837] "hmnodm"    "hmnofrd"   "hmnommb"   "hmnona"    "hmnonn"    "hmnoptp"  
## [1843] "hmnoref"   "hmnovw"    "impapol"   "impavo"    "impfml"    "impfrds"  
## [1849] "implsrt"   "impoblw"   "impopin"   "imppol"    "imprlg"    "imprwcr"  
## [1855] "imprwct"   "imprwkc"   "impsppl"   "impvo"     "impvote"   "impwrk"   
## [1861] "othvdm"    "othvfrd"   "othvmmb"   "othvna"    "othvnn"    "othvptp"  
## [1867] "othvref"   "othvvw"    "prfodm"    "prfofrd"   "prfommb"   "prfona"   
## [1873] "prfonn"    "prfoptp"   "prforef"   "prfovw"    "prtydm"    "prtyfrd"  
## [1879] "prtymmb"   "prtyna"    "prtynn"    "prtyptp"   "prtyref"   "prtyvw"   
## [1885] "rlgodm"    "rlgofrd"   "rlgommb"   "rlgona"    "rlgonn"    "rlgoptp"  
## [1891] "rlgoref"   "rlgovw"    "sclcdm"    "sclcfrd"   "sclcmmb"   "sclcna"   
## [1897] "sclcnn"    "sclcptp"   "sclcref"   "sclcvw"    "setodm"    "setofrd"  
## [1903] "setommb"   "setona"    "setonn"    "setoptp"   "setoref"   "setovw"   
## [1909] "sptcdm"    "sptcfrd"   "sptcmmb"   "sptcna"    "sptcnn"    "sptcptp"  
## [1915] "sptcref"   "sptcvw"    "stfhwkp"   "strtbsn"   "trudm"     "trufrd"   
## [1921] "truiwkp"   "trummb"    "truna"     "trunn"     "truptp"    "truref"   
## [1927] "trusay"    "truvw"     "truwrkp"   "wkchtsk"   "wkdcsin"   "wkenvin"  
## [1933] "wrkflex"   "yrlvdae"   "advbach"   "advhach"   "advslep"   "advsthr"  
## [1939] "chsrgp"    "dcadmms"   "dcdfcwr"   "dcdisc"    "dckptrt"   "dctreql"  
## [1945] "ddprsmd"   "dsplvpr"   "follwdc"   "hltherb"   "illcure"   "mdhair"   
## [1951] "mdhappy"   "mdlswgt"   "mdmemo"    "mdsdeff"   "mdsexlf"   "pplcure"  
## [1957] "pprlydc"   "prcbach"   "prchach"   "prcslep"   "prcsthr"   "prfmddc"  
## [1963] "prfsmdc"   "psmdcpr"   "ptnrlcq"   "tmcnsdc"   "usmdprs"   "usmprse"  
## [1969] "bnkfldl"   "bsnprft"   "cmprcti"   "ctzchtx"   "ctzhlpo"   "flgvbnf"  
## [1975] "flinsr"    "flinsrw"   "fodcncl"   "frdbnft"   "frmwktg"   "icsbfm"   
## [1981] "ignrlaw"   "kptchng"   "mnyacth"   "musdocm"   "olwmsop"   "payavtx"  
## [1987] "pboafvr"   "pbofvr"    "pbofvrw"   "pyavtxw"   "rprochg"   "scbevts"  
## [1993] "scndhfl"   "slcnflw"   "slcnsfl"   "tstfnch"   "tstpboh"   "tstrprh"  
## [1999] "wrytrdh"   "acldnmr"   "adftjob"   "adllvhm"   "adlvptn"   "adpnt"    
## [2005] "advcyc"    "aftjbyc"   "ageadlt"   "agemage"   "ageoage"   "alvgptn"  
## [2011] "anvcld"    "brnmm"     "bthcld"    "evlvptn"   "evmar"     "evpdemp"  
## [2017] "fcldbrn"   "ggchld"    "iaglptn"   "iaglvmr"   "iagpnt"    "iagrtr"   
## [2023] "lvpntyr"   "lvptnyr"   "maryr"     "nbthcld"   "ngchld"    "oldfrl"   
## [2029] "oldgpnt"   "oldhlpo"   "pdempyr"   "plnftr"    "rcldnmr"   "rdvcyc"   
## [2035] "rftjbyc"   "rlvgptn"   "rnvcld"    "rpntyng"   "rspslvo"   "rwkold"   
## [2041] "svclvo"    "tochld"    "tolvpnt"   "towkht"    "tygledu"   "tyglvmr"  
## [2047] "tyglvp"    "tygpnt"    "tygrtr"    "tygsexi"   "wrinco"    "ycldbyr"  
## [2053] "ygcdbyr"   "admub"     "basinc"    "bennent"   "bnlwinc"   "brnocnt"  
## [2059] "cldcrsv"   "dcndleq"   "dfincac"   "ditxssp"   "earnpen"   "earnueb"  
## [2065] "eduunmp"   "eudcnbf"   "eusclbf"   "gvcldcr"   "gvhlthc"   "gvjbevn"  
## [2071] "gvpdlwk"   "gvslvol"   "gvslvue"   "hlthcef"   "imrccon"   "imsclbn"  
## [2077] "insfben"   "lbenent"   "lklpwcf"   "lknemny"   "lknhlcn"   "lkuemp"   
## [2083] "lvoapen"   "lvpbhlt"   "nmnybsc"   "prtsick"   "sbbsntx"   "sbcwkfm"  
## [2089] "sbenccm"   "sbeqsoc"   "sblazy"    "sblwcoa"   "sblwlka"   "sbprvpv"  
## [2095] "sbstrec"   "schtaut"   "sickdsb"   "slvpens"   "slvuemp"   "smdfslv"  
## [2101] "txadleq"   "txautef"   "txearn"    "ub20edu"   "ub20pay"   "ub20unp"  
## [2107] "ub50edu"   "ub50pay"   "ub50unp"   "ubedu"     "ubpay"     "ubspedu"  
## [2113] "ubsppay"   "ubspunp"   "ubunp"     "uemplwk"   "uentrjb"   "wrkprbf"  
## [2119] "yngfnjb"   "acq30bs"   "acq70bs"   "agdcold"   "agdcyng"   "agegrbl"  
## [2125] "ageyo30"   "agrpblg"   "buproag"   "chg1530"   "dprschg"   "dprsf70"  
## [2131] "dprso70"   "dprsy30"   "frndo70"   "frndy30"   "lkrspag"   "mbfmo70"  
## [2137] "oafl20"    "oafl70"    "p20cnec"   "p20cust"   "p70cnec"   "p70cust"  
## [2143] "p70hltb"   "predage"   "predetn"   "predsex"   "pvolwrk"   "sgr2070"  
## [2149] "srdscag"   "stsp20"    "stsp40"    "stsp70"    "suproag"   "trtbdag"  
## [2155] "twcol20"   "twcol70"   "v20adm"    "v20cntm"   "v20comp"   "v20envy"  
## [2161] "v20frnd"   "v20mrst"   "v20pity"   "v20resp"   "v70adm"    "v70cntm"  
## [2167] "v70comp"   "v70envy"   "v70frnd"   "v70mrst"   "v70pity"   "v70resp"  
## [2173] "wr20crm"   "wr20emp"   "bplcdc"    "bstln5y"   "bystlct"   "bystlwr"  
## [2179] "caplcst"   "ctfrdc"    "ctinplt"   "ctjbcnt"   "ctmtgfr"   "ctprpwr"  
## [2185] "dbctvrd"   "doplcsy"   "dpcstrb"   "flsin5y"   "gsupplc"   "insclct"  
## [2191] "insclwr"   "jdgcbrb"   "lwstrob"   "plcarcr"   "plccbrb"   "plccbrg"  
## [2197] "plccont"   "plcexdc"   "plcfrdc"   "plciplt"   "plcipvl"   "plcjbcn"  
## [2203] "plcpvcr"   "plcrgwr"   "plcrspc"   "plcstf"    "plcvcrc"   "plcvcrp"  
## [2209] "rgbrklw"   "stcbg2t"   "tmprs"     "trfoct"    "trfowr"    "troff5y"  
## [2215] "wevdct"    "widprsn"   "wraccrc"   "wraccrp"   "chpldm"    "chpldmc"  
## [2221] "chpldmi"   "ctstogv"   "cttresa"   "cttresac"  "dfprtal"   "dfprtalc" 
## [2227] "dspplvt"   "dspplvtc"  "fairelc"   "fairelcc"  "fplvdm"    "fplvdmc"  
## [2233] "fplvdmi"   "gptpelc"   "gptpelcc"  "grdfinc"   "grdfincc"  "gvcodmc"  
## [2239] "gvcodmi"   "gvctzpv"   "gvctzpvc"  "gvexpdc"   "gvexpdcc"  "gvspcdm"  
## [2245] "gvspdmc"   "gvspdmi"   "imvtctz"   "medcrgv"   "medcrgvc"  "meprinf"  
## [2251] "meprinfc"  "oppcrgv"   "oppcrgvc"  "pltavie"   "pltaviec"  "pplvdmc"  
## [2257] "pplvdmi"   "rghmgpr"   "rghmgprc"  "stpldmc"   "stpldmi"   "votedir"  
## [2263] "votedirc"  "viepol"    "wpestop"   "keydec"    "viepolc"   "wpestopc" 
## [2269] "keydecc"   "admit"     "showcv"    "accalaw"   "impdema"   "impdemb"  
## [2275] "impdemc"   "impdemd"   "impdeme"   "alcbnge"   "alcfreq"   "alcwkdy"  
## [2281] "alcwknd"   "anypacc"   "cgtsday"   "cgtsmke"   "cnfpplh"   "dosprt"   
## [2287] "dshltdk"   "dshltgp"   "dshltms"   "dshltna"   "dshltnt"   "dshltref" 
## [2293] "eatveg"    "etfruit"   "fnsdfml"   "height"    "hlpfmhr"   "hlpfmly"  
## [2299] "hltphal"   "hltphbn"   "hltphbp"   "hltphdi"   "hltphdk"   "hltphhb"  
## [2305] "hltphhc"   "hltphna"   "hltphnap"  "hltphnt"   "hltphpa"   "hltphpf"  
## [2311] "hltphref"  "hltphsc"   "hltphsd"   "hltphsh"   "hltpral"   "hltprbn"  
## [2317] "hltprbp"   "hltprca"   "hltprdi"   "hltprdk"   "hltprhb"   "hltprhc"  
## [2323] "hltprna"   "hltprnt"   "hltprpa"   "hltprpf"   "hltprref"  "hltprsc"  
## [2329] "hltprsd"   "hltprsh"   "jbexebs"   "jbexecp"   "jbexedk"   "jbexena"  
## [2335] "jbexenap"  "jbexent"   "jbexera"   "jbexeref"  "jbexevc"   "jbexevh"  
## [2341] "jbexevl"   "jbexpdk"   "jbexpmc"   "jbexpml"   "jbexpna"   "jbexpnap" 
## [2347] "jbexpnt"   "jbexpref"  "jbexpti"   "jbexpvi"   "medtrdk"   "medtrna"  
## [2353] "medtrnap"  "medtrnl"   "medtrnp"   "medtrnt"   "medtrnu"   "medtrnaa" 
## [2359] "medtroc"   "medtroth"  "medtrref"  "medtrun"   "medtrwl"   "trhltacp" 
## [2365] "trhltacu"  "trhltch"   "trhltcm"   "trhltdk"   "trhltho"   "trhltht"  
## [2371] "trhlthy"   "trhltmt"   "trhltna"   "trhltnt"   "trhltos"   "trhltpt"  
## [2377] "trhltre"   "trhltref"  "trhltsh"   "weight"    "icgndra"   "banhhap"  
## [2383] "ccgdbd"    "cflsenr"   "clmchng"   "clmthgt1"  "clmthgt2"  "elgbio"   
## [2389] "elgcoal"   "elghydr"   "elgngas"   "elgnuc"    "elgsun"    "elgwind"  
## [2395] "eneffap"   "gvsrdcc"   "inctxff"   "lklmten"   "lkredcc"   "ownrdcc"  
## [2401] "rdcenr"    "sbsrnen"   "wrdpfos"   "wrdpimp"   "wrenexp"   "wrinspw"  
## [2407] "wrntdis"   "wrpwrct"   "wrtcfl"    "wrtratc"   "btminfr"   "evfredu"  
## [2413] "evfrjob"   "frlgrsp"   "frlneti"   "frprtpl"   "fvgabc"    "grspfr"   
## [2419] "grsplet"   "grspnum"   "gvintcz"   "ifredu"    "ifrjob"    "infqbst"  
## [2425] "jstprev"   "netifr"    "netilet"   "netinum"   "occinfr"   "pcmpinj"  
## [2431] "poltran"   "ppldsrv"   "recexp"    "recgndr"   "recimg"    "recknow"  
## [2437] "recskil"   "sofrdst"   "sofrpr"    "sofrprv"   "sofrwrk"   "topinfr"  
## [2443] "wltdffr"   "iincsrc"   "acchome"   "accwrk"    "accmove"   "accoth"   
## [2449] "accnone"   "accref"    "accdk"     "accna"     "fampref"   "famadvs"  
## [2455] "fampdf"    "mcclose"   "mcinter"   "mccoord"   "mcpriv"    "mcmsinf"  
## [2461] "chldo12"   "gndro12a"  "gndro12b"  "ageo12"    "hhlio12"   "closeo12" 
## [2467] "ttmino12"  "speako12"  "como12"    "scrno12"   "phoneo12"  "livpnt"   
## [2473] "pntmofa"   "agepnt"    "hhlipnt"   "closepnt"  "ttminpnt"  "speakpnt" 
## [2479] "compnt"    "scrnpnt"   "phonepnt"  "wrkhome"   "wrklong"   "wrkresp"  
## [2485] "scrhlm"    "sclmsup"   "mansupp"   "manhlp"    "manwrkpl"  "manspeak" 
## [2491] "mancom"    "manscrn"   "manphone"  "teamfeel"  "wrkextra"  "colprop"  
## [2497] "colhlp"    "colspeak"  "colcom"    "colscrn"   "colphone"  "mcwrkhom" 
## [2503] "c19spo12"  "c19mco12"  "c19sppnt"  "c19mcpnt"  "c19whome"  "c19wplch" 
## [2509] "c19whacc"  "c19spwrk"  "c19mcwrk"  "secgrdec"  "scidecpb"  "admc19"   
## [2515] "panpriph"  "panmonpb"  "govpriph"  "govmonpb"  "panfolru"  "panclobo" 
## [2521] "panresmo"  "gvhanc19"  "gvjobc19"  "gveldc19"  "gvfamc19"  "hscopc19" 
## [2527] "gvbalc19"  "gvimpc19"  "gvconc19"  "respc19"   "reshhc19"  "hapljc19" 
## [2533] "hapirc19"  "hapwrc19"  "hapfuc19"  "hapfoc19"  "hapnoc19"  "hapnwc19" 
## [2539] "hapnpc19"  "haprec19"  "hapdkc19"  "hapnac19"  "icvacc19"  "getavc19" 
## [2545] "getnvc19"  "vdcond"    "vdovexre"  "vdtype"    "vdtpsvre"  "vdtpitre" 
## [2551] "vdtpscre"  "vdtpaure"  "vdtpvire"  "vdtpoire"  "vdtpntre"  "vdtpapre" 
## [2557] "vdtprere"  "vdtpdkre"  "vdtpnare"  "cldgng"    "enjlf"     "fltdpr"   
## [2563] "flteeff"   "fltlnl"    "fltlnla"   "fltsd"     "hrshsnt"   "hrshsnta" 
## [2569] "mnrgtjb"   "pdaprp"    "pdaprpa"   "slprl"     "smbtjob"   "smbtjoba" 
## [2575] "stfjbot"   "wrhpp"     "stfmjob"   "wrywprb"   "trdawrk"   "jbprtfp"  
## [2581] "pfmfdjb"   "pfmfdjba"  "dcsfwrk"   "dcsfwrka"  "implvdm"   "ccrdprs"  
## [2587] "ccnthum"   "wrclmch"   "yrbrn"     "wkhct"     "region"    "admdw"    
## [2593] "admge"     "iagmr"     "tygmr"     "domain"    "trstsci"   "prtvtmk"  
## [2599] "prtclmk"   "rlgdnmk"   "admrclc"   "testic34"  "testic35"  "testic36" 
## [2605] "testic37"  "testic38"  "testic39"  "testic40"  "testic41"  "testic42" 
## [2611] "rshpsmk"   "marstmk"   "edlvdmk"   "edlvpdmk"  "edlvfdmk"  "edlvmdmk" 
## [2617] "testii1"   "testii2"   "testii3"   "testii4"   "testii5"   "testii6"  
## [2623] "testii7"   "testii8"   "testii9"   "ainws"     "ainwe"     "binwe"    
## [2629] "cinwe"     "dinwe"     "finwe"     "ginwe"     "hinwe"     "iinwe"    
## [2635] "kinwe"     "vinwe"     "jinws"     "jinwe"     "scrlgblg"  "scchpldm" 
## [2641] "scwflt"    "scrwiat"   "sctmfeel"  "scrwwc"    "scclprp"   "scwsds"   
## [2647] "scwass"    "scwa5se"   "scwa11se"  "scwa59se"  "scwa89se"  "scwb30se" 
## [2653] "scwc55se"  "scwd62se"  "scwd78se"  "scwrecse"  "scwtype"   "uagent"   
## [2659] "smplresp"  "questcmp"  "intagea"   "intgndr"   "intnum"    "intlnga"  
## [2665] "rlgdemk"
head(portugal_data[, c("proddate", "essround")])
##     proddate essround
## 1 01.12.2018        1
## 2 01.12.2018        1
## 3 01.12.2018        1
## 4 01.12.2018        1
## 5 01.12.2018        1
## 6 01.12.2018        1
portugal_data <- portugal_data %>%
  mutate(birth_year = 2018 - age)
portugal_data <- portugal_data %>%
  mutate(Generation = case_when(
    birth_year >= 1928 & birth_year <= 1945 ~ "Interwar",
    birth_year >= 1946 & birth_year <= 1964 ~ "Baby Boomers",
    birth_year >= 1965 & birth_year <= 1980 ~ "Generation X",
    birth_year >= 1981 & birth_year <= 1996 ~ "Millennials",
    birth_year >= 1997 & birth_year <= 2005 ~ "Gen Z",
    TRUE ~ NA_character_
  )) %>%
  mutate(Generation = factor(Generation, levels = c("Interwar", "Baby Boomers", "Generation X", "Millennials", "Gen Z")))
library(ggplot2)
ggplot(portugal_data, aes(x = Generation)) +
  geom_bar(fill = "steelblue") +
  labs(title = "Distribution of Generations", x = "Generation", y = "Count") +
  theme_minimal()

# Create populist_attitude from trust variables
portugal_data <- portugal_data %>%
  mutate(
    trust_scale = rowMeans(select(., trstplt, trstlgl, trstprl), na.rm = TRUE),
    trust_scale_100 = scales::rescale(trust_scale, to = c(0, 100), na.rm = TRUE),
    populist_attitude = 100 - trust_scale_100
  )
ggplot(portugal_data, aes(x = Generation, y = populist_attitude, fill = Generation)) +
  geom_boxplot() +
  labs(
    title = "Populist Attitudes Across Generations in Portugal",
    subtitle = "Based on Trust in Political Institutions",
    x = "Generation",
    y = "Populist Attitude (0-100 Scale)",
    fill = "Generation"
  ) +
  theme_minimal() +
  theme(axis.text.x = element_text(angle = 45, hjust = 1))

portugal_data <- portugal_data %>%
  mutate(
    # Clean invalid trust values
    trstplt = ifelse(trstplt %in% c(77, 88, 99), NA, trstplt),
    trstprl = ifelse(trstprl %in% c(77, 88, 99), NA, trstprl),
    trstprt = ifelse(trstprt %in% c(77, 88, 99), NA, trstprt),
    
    # Calculate trust average and scale it
    trust_scale = rowMeans(select(., trstplt, trstprl, trstprt), na.rm = TRUE),
    trust_scale_100 = scales::rescale(trust_scale, to = c(0, 100), na.rm = TRUE),
    
    # Create populist scale (reverse of trust)
    populism = 100 - trust_scale_100
  )
summary(portugal_data$populism)
##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
##    0.00   95.29   96.97   94.70   98.65  100.00
portugal_data$Gender <- ifelse(
  portugal_data$gndr == 1, "Male",
  ifelse(portugal_data$gndr == 2, "Female", NA)
)

portugal_data$Gender <- as.factor(portugal_data$Gender)
table(portugal_data$Gender, useNA = "ifany")
## 
## Female   Male 
##  10613   7268
model1 <- lm(populism ~ Gender, data = portugal_data)
table(portugal_data$Generation, useNA = "ifany")
## 
##     Interwar Baby Boomers Generation X  Millennials        Gen Z         <NA> 
##          351          656          484          479          120        15791
portugal_data <- portugal_data %>%
  mutate(
    Generation = case_when(
      yrbrn >= 1928 & yrbrn <= 1945 ~ "Interwar",
      yrbrn >= 1946 & yrbrn <= 1964 ~ "Baby Boomers",
      yrbrn >= 1965 & yrbrn <= 1980 ~ "Generation X",
      yrbrn >= 1981 & yrbrn <= 1996 ~ "Millennials",
      yrbrn >= 1997 & yrbrn <= 2005 ~ "Gen Z",
      TRUE ~ NA_character_
    ),
    Generation = factor(Generation, levels = c("Interwar", "Baby Boomers", "Generation X", "Millennials", "Gen Z"))
  )
table(portugal_data$Generation, useNA = "ifany")
## 
##     Interwar Baby Boomers Generation X  Millennials        Gen Z         <NA> 
##         4382         5371         4416         2645          307          760
model2 <- lm(populism ~ Gender + Generation, data = portugal_data)
library(modelsummary)

modelsummary(
  list("Gender Only" = model1, "Gender + Generation" = model2),
  stars = TRUE,
  coef_map = c("GenderMale" = "Male", 
               "GenerationBaby Boomers" = "Baby Boomers",
               "GenerationGeneration X" = "Generation X",
               "GenerationMillennials" = "Millennials",
               "GenerationGen Z" = "Gen Z"),
  gof_map = c("nobs", "r.squared", "adj.r.squared", "AIC", "BIC"),
  title = "Regression Results: Predicting Populist Attitudes in Portugal"
)
Regression Results: Predicting Populist Attitudes in Portugal
Gender Only Gender + Generation
+ p < 0.1, * p < 0.05, ** p < 0.01, *** p < 0.001
Male 1.603*** 1.348***
(0.181) (0.174)
Baby Boomers 2.300***
(0.227)
Generation X 2.552***
(0.238)
Millennials 1.785***
(0.275)
Gen Z 0.625
(0.659)
Num.Obs. 17881 17121
R2 0.004 0.012
R2 Adj. 0.004 0.012
library(modelsummary)

modelsummary(
  list("Gender Only" = model1, "Gender + Generation" = model2),
  stars = TRUE,
  coef_map = c(
    "GenderMale" = "Male", 
    "GenerationBaby Boomers" = "Baby Boomers",
    "GenerationGeneration X" = "Generation X",
    "GenerationMillennials" = "Millennials",
    "GenerationGen Z" = "Gen Z"
  ),
  gof_map = c("nobs", "r.squared", "adj.r.squared", "AIC", "BIC"),
  title = "Regression Results: Predicting Populist Attitudes in Portugal"
)
Regression Results: Predicting Populist Attitudes in Portugal
Gender Only Gender + Generation
+ p < 0.1, * p < 0.05, ** p < 0.01, *** p < 0.001
Male 1.603*** 1.348***
(0.181) (0.174)
Baby Boomers 2.300***
(0.227)
Generation X 2.552***
(0.238)
Millennials 1.785***
(0.275)
Gen Z 0.625
(0.659)
Num.Obs. 17881 17121
R2 0.004 0.012
R2 Adj. 0.004 0.012
library(sjPlot)

plot_model(model2, 
           type = "est", 
           vline.color = "red", 
           title = "Effect of Gender and Generation on Populist Attitudes")

The first regression model examined the effect of gender on populist attitudes. In this model, the intercept represents the average populist score for females (the reference category). The coefficient for “Male” indicates the difference in average populist attitudes between men and women. The result suggests that men exhibit slightly higher populist attitudes than women, but the size of this effect is small and may not be statistically significant depending on the p-value. The model’s R² is also low, indicating that gender alone does not explain much of the variation in populist attitudes.

The second model expanded on this by adding generation as an additional predictor. In this model, the intercept represents the predicted populist score for females in the “Interwar” generation (the reference category). Each generational coefficient shows how much more or less populist the other generations are in comparison to Interwar, controlling for gender.

From the results, we can see notable generational differences. For example, “Gen Z” has significantly higher populist attitudes compared to the Interwar generation, suggesting a generational shift in political trust. Other generations (like Millennials and Generation X) also show differences, though not all may be statistically significant. These coefficients reflect real generational patterns in political trust — younger people, especially Gen Z, may be more skeptical of institutions, translating into higher populist scores.

The gender coefficient in Model 2 remains similar in direction but may slightly change in magnitude. This shift reflects how adding generation as a predictor helps to explain part of the variance previously attributed to gender. In other words, some differences between men and women in Model 1 were actually due to generational differences.

When comparing model fit, Model 2 clearly outperforms Model 1. It has a higher R², indicating it explains more of the variance in populist attitudes. Additionally, the AIC and BIC values are lower in Model 2, which supports its superiority in model selection. This shows that including both gender and generation provides a more complete understanding of populist attitudes in Portugal.

Overall, while gender shows a mild influence, generational identity plays a stronger role in shaping populist attitudes. Younger generations are more likely to express skepticism or distrust in political institutions, and these attitudes appear to be systematically patterned by birth cohort. These findings highlight the importance of considering both demographic and generational context when analyzing political attitudes.

BONUS QUESTION Why the Intercept and Gender Coefficient Changed

The intercept changed between Models 1 and 2 because it always reflects the expected value of the outcome for the reference categories of the predictors. In Model 1, the intercept was the average populist score for females. In Model 2, it is now the average for females in the Interwar generation, since we added Generation to the model.

The coefficient for “Male” also changed slightly between models because in Model 1, it captured all gender-related differences in populism. But in Model 2, part of what was originally attributed to gender may now be better explained by generational differences. This is a common shift when you add new predictors — earlier coefficients often change as overlapping sources of variation are separated out.

This highlights how gender and generation are interrelated in shaping attitudes, and why including all relevant predictors is crucial for an accurate model.

Bonus Question (Task 1)

Q: Looking only at the coefficient visualization from sjPlot, how can you determine which coefficients are statistically significant? Why is this helpful?

Answer:

In the coefficient plot from sjPlot, each coefficient is shown with a horizontal line representing its confidence interval. If this line does not cross the red vertical line at zero, the coefficient is statistically significant. If the line does cross zero, the effect is not statistically significant.

This visual approach is helpful because it allows you to quickly identify which predictors have meaningful effects without digging through tables of p-values. It also shows you the size and direction of effects (positive or negative), making interpretation easier at a glance.

# Load required packages if not already loaded
library(tidyverse)
library(fst)

# Load France data
france_data <- read_fst("france_data.fst")
colnames(france_data)
##    [1] "name"      "essround"  "edition"   "proddate"  "cntry"     "idno"     
##    [7] "dweight"   "pspwght"   "pweight"   "anweight"  "prob"      "stratum"  
##   [13] "psu"       "netuse"    "netusoft"  "netustm"   "nwspol"    "nwsppol"  
##   [19] "nwsptot"   "pplfair"   "pplhlp"    "ppltrst"   "rdpol"     "rdtot"    
##   [25] "tvpol"     "tvtot"     "actrolg"   "actrolga"  "badge"     "bctprd"   
##   [31] "bghtprd"   "clsprty"   "contplt"   "cptppol"   "cptppola"  "dclagr"   
##   [37] "dclaid"    "dclcrm"    "dcldef"    "dclenv"    "dclintr"   "dclmig"   
##   [43] "dclwlfr"   "dmcntov"   "dntmny"    "ecohenv"   "etapapl"   "euftf"    
##   [49] "freehms"   "gincdif"   "ginveco"   "hmsacld"   "hmsfmlsh"  "ilglpst"  
##   [55] "lawobey"   "lrscale"   "mmbprty"   "needtru"   "pbldmn"    "pbldmna"  
##   [61] "pltcare"   "pltinvt"   "polactiv"  "polcmpl"   "poldcs"    "polintr"  
##   [67] "prtdgcl"   "prtyban"   "psppipl"   "psppipla"  "psppsgv"   "psppsgva" 
##   [73] "pstplonl"  "ptcpplt"   "scnsenv"   "sgnptit"   "stfdem"    "stfeco"   
##   [79] "stfedu"    "stfgov"    "stfhlth"   "stflife"   "trstep"    "trstlgl"  
##   [85] "trstplc"   "trstplt"   "trstprl"   "trstprt"   "trstun"    "vote"     
##   [91] "wrkorg"    "wrkprty"   "imsmetn"   "imdfetn"   "impcntr"   "imbgeco"  
##   [97] "imueclt"   "imwbcnt"   "volunfp"   "lrnobed"   "loylead"   "donprty"  
##  [103] "prtvde1"   "prtvade1"  "prtvbde1"  "prtvcde1"  "prtvdde1"  "prtvede1" 
##  [109] "prtvfde1"  "prtvde2"   "prtvade2"  "prtvbde2"  "prtvcde2"  "prtvdde2" 
##  [115] "prtvede2"  "prtvfde2"  "prtvlt1"   "prtvalt1"  "prtvblt1"  "prtvclt1" 
##  [121] "prtvlt2"   "prtvalt2"  "prtvblt2"  "prtvclt2"  "prtvlt3"   "prtvalt3" 
##  [127] "prtvblt3"  "prtvclt3"  "prtvtal"   "prtvtat"   "prtvtaat"  "prtvtbat" 
##  [133] "prtvtcat"  "prtvtbe"   "prtvtabe"  "prtvtbbe"  "prtvtcbe"  "prtvtdbe" 
##  [139] "prtvtebe"  "prtvtbg"   "prtvtabg"  "prtvtbbg"  "prtvtcbg"  "prtvtdbg" 
##  [145] "prtvtebg"  "prtvtch"   "prtvtach"  "prtvtbch"  "prtvtcch"  "prtvtdch" 
##  [151] "prtvtech"  "prtvtfch"  "prtvtgch"  "prtvthch"  "prtvtcy"   "prtvtacy" 
##  [157] "prtvtbcy"  "prtvtccy"  "prtvtcz"   "prtvtacz"  "prtvtbcz"  "prtvtccz" 
##  [163] "prtvtdcz"  "prtvtecz"  "prtvtdk"   "prtvtadk"  "prtvtbdk"  "prtvtcdk" 
##  [169] "prtvtddk"  "prtvtee"   "prtvtaee"  "prtvtbee"  "prtvtcee"  "prtvtdee" 
##  [175] "prtvteee"  "prtvtfee"  "prtvtgee"  "prtvthee"  "prtvtes"   "prtvtaes" 
##  [181] "prtvtbes"  "prtvtces"  "prtvtdes"  "prtvtees"  "prtvtfes"  "prtvtfi"  
##  [187] "prtvtafi"  "prtvtbfi"  "prtvtcfi"  "prtvtdfi"  "prtvtefi"  "prtvtfr"  
##  [193] "prtvtafr"  "prtvtbfr"  "prtvtcfr"  "prtvtdfr"  "prtvtefr"  "prtvtgb"  
##  [199] "prtvtagb"  "prtvtbgb"  "prtvtcgb"  "prtvtdgb"  "prtvtgr"   "prtvtagr" 
##  [205] "prtvtbgr"  "prtvtcgr"  "prtvtdgr"  "prtvthr"   "prtvtahr"  "prtvtbhr" 
##  [211] "prtvthu"   "prtvtahu"  "prtvtbhu"  "prtvtchu"  "prtvtdhu"  "prtvtehu" 
##  [217] "prtvtfhu"  "prtvtghu"  "prtvtie"   "prtvtaie"  "prtvtbie"  "prtvtcie" 
##  [223] "prtvtdie"  "prtvtil"   "prtvtail"  "prtvtbil"  "prtvtcil"  "prtvtdil" 
##  [229] "prtvtis"   "prtvtais"  "prtvtbis"  "prtvtcis"  "prtvtdis"  "prtvtit"  
##  [235] "prtvtait"  "prtvtbit"  "prtvtcit"  "prtvtdit"  "prtvtlu"   "prtvtlv"  
##  [241] "prtvtalv"  "prtvtnl"   "prtvtanl"  "prtvtbnl"  "prtvtcnl"  "prtvtdnl" 
##  [247] "prtvtenl"  "prtvtfnl"  "prtvtgnl"  "prtvthnl"  "prtvtno"   "prtvtano" 
##  [253] "prtvtbno"  "prtvtpl"   "prtvtapl"  "prtvtbpl"  "prtvtcpl"  "prtvtdpl" 
##  [259] "prtvtepl"  "prtvtpt"   "prtvtapt"  "prtvtbpt"  "prtvtcpt"  "prtvtdpt" 
##  [265] "prtvtro"   "prtvtaro"  "prtvtru"   "prtvtaru"  "prtvtbru"  "prtvtcru" 
##  [271] "prtvtdru"  "prtvtse"   "prtvtase"  "prtvtbse"  "prtvtcse"  "prtvtdse" 
##  [277] "prtvtsi"   "prtvtasi"  "prtvtbsi"  "prtvtcsi"  "prtvtdsi"  "prtvtesi" 
##  [283] "prtvtfsi"  "prtvtsk"   "prtvtask"  "prtvtbsk"  "prtvtcsk"  "prtvtdsk" 
##  [289] "prtvtesk"  "prtvttr"   "prtvtatr"  "prtvtua"   "prtvtaua"  "prtvtbua" 
##  [295] "prtvtcua"  "prtvtxk"   "prtvtrs"   "prtvtars"  "prtvtme"   "prtvtame" 
##  [301] "prtclal"   "prtclat"   "prtclaat"  "prtclcat"  "prtcldat"  "prtclbe"  
##  [307] "prtclabe"  "prtclbbe"  "prtclcbe"  "prtcldbe"  "prtclebe"  "prtclbg"  
##  [313] "prtclabg"  "prtclbbg"  "prtclcbg"  "prtcldbg"  "prtclebg"  "prtclch"  
##  [319] "prtclach"  "prtclbch"  "prtclcch"  "prtcldch"  "prtclech"  "prtclfch" 
##  [325] "prtclgch"  "prtclhch"  "prtclsk"   "prtclask"  "prtclbsk"  "prtclcsk" 
##  [331] "prtcldsk"  "prtclesk"  "prtclcy"   "prtclacy"  "prtclbcy"  "prtclccy" 
##  [337] "prtclcz"   "prtclacz"  "prtclbcz"  "prtclccz"  "prtcldcz"  "prtclecz" 
##  [343] "prtclde"   "prtclade"  "prtclbde"  "prtclcde"  "prtcldde"  "prtclede" 
##  [349] "prtclfde"  "prtcldk"   "prtcladk"  "prtclbdk"  "prtclcdk"  "prtclddk" 
##  [355] "prtclee"   "prtclaee"  "prtclbee"  "prtclcee"  "prtcldee"  "prtcleee" 
##  [361] "prtclfee"  "prtclgee"  "prtclhee"  "prtcles"   "prtclaes"  "prtclbes" 
##  [367] "prtclces"  "prtcldes"  "prtclees"  "prtclfes"  "prtclges"  "prtclfi"  
##  [373] "prtclafi"  "prtclbfi"  "prtclcfi"  "prtcldfi"  "prtclefi"  "prtclffi" 
##  [379] "prtclfr"   "prtclafr"  "prtclbfr"  "prtclcfr"  "prtcldfr"  "prtclefr" 
##  [385] "prtclffr"  "prtclgb"   "prtclagb"  "prtclbgb"  "prtclcgb"  "prtcldgb" 
##  [391] "prtclgr"   "prtclagr"  "prtclbgr"  "prtclcgr"  "prtcldgr"  "prtclhr"  
##  [397] "prtclahr"  "prtclbhr"  "prtclhu"   "prtclahu"  "prtclbhu"  "prtclchu" 
##  [403] "prtcldhu"  "prtclehu"  "prtclfhu"  "prtclghu"  "prtclhhu"  "prtclie"  
##  [409] "prtclaie"  "prtclbie"  "prtclcie"  "prtcldie"  "prtcleie"  "prtclfie" 
##  [415] "prtclil"   "prtclail"  "prtclbil"  "prtclcil"  "prtcldil"  "prtcleil" 
##  [421] "prtclis"   "prtclais"  "prtclbis"  "prtclcis"  "prtcldis"  "prtclit"  
##  [427] "prtclait"  "prtclbit"  "prtclcit"  "prtcldit"  "prtcleit"  "prtcllt"  
##  [433] "prtclalt"  "prtclblt"  "prtclclt"  "prtcllu"   "prtcllv"   "prtclalv" 
##  [439] "prtclnl"   "prtclanl"  "prtclbnl"  "prtclcnl"  "prtcldnl"  "prtclenl" 
##  [445] "prtclfnl"  "prtclgnl"  "prtclno"   "prtclano"  "prtclbno"  "prtclpl"  
##  [451] "prtclapl"  "prtclbpl"  "prtclcpl"  "prtcldpl"  "prtclepl"  "prtclfpl" 
##  [457] "prtclgpl"  "prtclhpl"  "prtclipl"  "prtclpt"   "prtclapt"  "prtclbpt" 
##  [463] "prtclcpt"  "prtcldpt"  "prtclept"  "prtclfpt"  "prtclro"   "prtclaro" 
##  [469] "prtclru"   "prtclaru"  "prtclbru"  "prtclcru"  "prtcldru"  "prtclse"  
##  [475] "prtclase"  "prtclbse"  "prtclcse"  "prtcldse"  "prtclsi"   "prtclasi" 
##  [481] "prtclbsi"  "prtclcsi"  "prtcldsi"  "prtclesi"  "prtclfsi"  "prtcltr"  
##  [487] "prtclatr"  "prtclua"   "prtclaua"  "prtclbua"  "prtclcua"  "prtcldua" 
##  [493] "prtclxk"   "prtclrs"   "prtclars"  "prtclme"   "prtclame"  "prtmbat"  
##  [499] "prtmbaat"  "prtmbbe"   "prtmbabe"  "prtmbbbe"  "prtmbcbe"  "prtmbbg"  
##  [505] "prtmbabg"  "prtmbbbg"  "prtmbch"   "prtmbach"  "prtmbbch"  "prtmbcch" 
##  [511] "prtmbcy"   "prtmbcz"   "prtmbacz"  "prtmbbcz"  "prtmbde"   "prtmbade" 
##  [517] "prtmbbde"  "prtmbcde"  "prtmbdk"   "prtmbadk"  "prtmbbdk"  "prtmbee"  
##  [523] "prtmbaee"  "prtmbbee"  "prtmbcee"  "prtmbes"   "prtmbaes"  "prtmbbes" 
##  [529] "prtmbfi"   "prtmbafi"  "prtmbbfi"  "prtmbfr"   "prtmbafr"  "prtmbbfr" 
##  [535] "prtmbcfr"  "prtmbgb"   "prtmbagb"  "prtmbgr"   "prtmbagr"  "prtmbbgr" 
##  [541] "prtmbcgr"  "prtmbhr"   "prtmbhu"   "prtmbahu"  "prtmbbhu"  "prtmbchu" 
##  [547] "prtmbie"   "prtmbaie"  "prtmbbie"  "prtmbil"   "prtmbail"  "prtmbbil" 
##  [553] "prtmbis"   "prtmbit"   "prtmbait"  "prtmblt"   "prtmblu"   "prtmblv"  
##  [559] "prtmbnl"   "prtmbanl"  "prtmbbnl"  "prtmbcnl"  "prtmbno"   "prtmbano" 
##  [565] "prtmbpl"   "prtmbapl"  "prtmbbpl"  "prtmbcpl"  "prtmbdpl"  "prtmbpt"  
##  [571] "prtmbapt"  "prtmbbpt"  "prtmbro"   "prtmbaro"  "prtmbru"   "prtmbaru" 
##  [577] "prtmbbru"  "prtmbse"   "prtmbase"  "prtmbsi"   "prtmbasi"  "prtmbbsi" 
##  [583] "prtmbcsi"  "prtmbsk"   "prtmbask"  "prtmbbsk"  "prtmbtr"   "prtmbatr" 
##  [589] "prtmbua"   "prtmbaua"  "prtmbbua"  "prtmbcua"  "aesfdrk"   "atchctr"  
##  [595] "atcherp"   "blgetmg"   "brghmef"   "brghmwr"   "brncntr"   "cntbrth"  
##  [601] "cntbrtha"  "cntbrthb"  "cntbrthc"  "cntbrthd"  "crmvct"    "crvctef"  
##  [607] "crvctwr"   "ctzcntr"   "ctzship"   "ctzshipa"  "ctzshipb"  "ctzshipc" 
##  [613] "ctzshipd"  "dscrage"   "dscrdk"    "dscrdsb"   "dscretn"   "dscrgnd"  
##  [619] "dscrgrp"   "dscrlng"   "dscrna"    "dscrnap"   "dscrntn"   "dscroth"  
##  [625] "dscrrce"   "dscrref"   "dscrrlg"   "dscrsex"   "facntn"    "facntr"   
##  [631] "fbrncnt"   "fbrncnta"  "fbrncntb"  "fbrncntc"  "happy"     "health"   
##  [637] "hlthhmp"   "inmdisc"   "inprdsc"   "livecntr"  "livecnta"  "lnghoma"  
##  [643] "lnghom1"   "lnghomb"   "lnghom2"   "mbrncnt"   "mbrncnta"  "mbrncntb" 
##  [649] "mbrncntc"  "mmbrn"     "mocntn"    "mocntr"    "pray"      "rlgatnd"  
##  [655] "rlgblg"    "rlgblge"   "rlgdgr"    "rlgdnm"    "rlgdnme"   "sclact"   
##  [661] "sclmeet"   "trrcnyr"   "trrenyr"   "trrprsn"   "trrtort"   "feethngr" 
##  [667] "vteurmmb"  "vteubcmb"  "vteumbgb"  "rlgdnal"   "rlgdnat"   "rlgdnbat" 
##  [673] "rlgdnbe"   "rlgdnch"   "rlgdnach"  "rlgdncy"   "rlgdnacy"  "rlgdnde"  
##  [679] "rlgdnade"  "rlgdndk"   "rlgdnfi"   "rlgdnafi"  "rlgdngb"   "rlgdngr"  
##  [685] "rlgdnagr"  "rlgdnhu"   "rlgdnie"   "rlgdnil"   "rlgdnis"   "rlgdnais" 
##  [691] "rlgdnlt"   "rlgdnlv"   "rlgdnnl"   "rlgdnanl"  "rlgdnno"   "rlgdnpl"  
##  [697] "rlgdnapl"  "rlgdnbpl"  "rlgdnpt"   "rlgdnro"   "rlgdnru"   "rlgdnaru" 
##  [703] "rlgdnse"   "rlgdnase"  "rlgdnsi"   "rlgdnsk"   "rlgdnask"  "rlgdnbsk" 
##  [709] "rlgdnua"   "rlgdnrs"   "rlgdme"    "rlgdeal"   "rlgdeat"   "rlgdebat" 
##  [715] "rlgdebe"   "rlgdech"   "rlgdeach"  "rlgdecy"   "rlgdeacy"  "rlgdede"  
##  [721] "rlgdeade"  "rlgdedk"   "rlgdefi"   "rlgdeafi"  "rlgdegb"   "rlgdegr"  
##  [727] "rlgdeagr"  "rlgdehu"   "rlgdeie"   "rlgdeil"   "rlgdeis"   "rlgdeais" 
##  [733] "rlgdelt"   "rlgdelv"   "rlgdenl"   "rlgdeanl"  "rlgdeno"   "rlgdepl"  
##  [739] "rlgdeapl"  "rlgdebpl"  "rlgdept"   "rlgdero"   "rlgderu"   "rlgdearu" 
##  [745] "rlgdese"   "rlgdease"  "rlgdesi"   "rlgdesk"   "rlgdeask"  "rlgdebsk" 
##  [751] "rlgdeua"   "rlgders"   "rlgdeme"   "hhmmb"     "gndr"      "gndr2"    
##  [757] "gndr3"     "gndr4"     "gndr5"     "gndr6"     "gndr7"     "gndr8"    
##  [763] "gndr9"     "gndr10"    "gndr11"    "gndr12"    "gndr13"    "gndr14"   
##  [769] "gndr15"    "gndr16"    "gndr17"    "gndr18"    "gndr19"    "gndr20"   
##  [775] "gndr21"    "gndr22"    "gndr23"    "gndr24"    "rship2"    "rshipa2"  
##  [781] "rship3"    "rshipa3"   "rship4"    "rshipa4"   "rship5"    "rshipa5"  
##  [787] "rship6"    "rshipa6"   "rship7"    "rshipa7"   "rship8"    "rshipa8"  
##  [793] "rship9"    "rshipa9"   "rship10"   "rshipa10"  "rship11"   "rshipa11" 
##  [799] "rship12"   "rshipa12"  "rship13"   "rshipa13"  "rship14"   "rshipa14" 
##  [805] "rship15"   "rshipa15"  "rshipa16"  "rshipa17"  "rshipa18"  "rshipa19" 
##  [811] "rshipa20"  "rshipa21"  "rshipa22"  "rshipa23"  "rshipa24"  "yrbrn2"   
##  [817] "yrbrn3"    "yrbrn4"    "yrbrn5"    "yrbrn6"    "yrbrn7"    "yrbrn8"   
##  [823] "yrbrn9"    "yrbrn10"   "yrbrn11"   "yrbrn12"   "yrbrn13"   "yrbrn14"  
##  [829] "yrbrn15"   "yrbrn16"   "yrbrn17"   "yrbrn18"   "yrbrn19"   "yrbrn20"  
##  [835] "yrbrn21"   "yrbrn22"   "yrbrn23"   "yrbrn24"   "age"       "agea"     
##  [841] "ablrtr"    "anctry1"   "anctry2"   "atncrse"   "brwmny"    "chldhhe"  
##  [847] "chldhm"    "cmsrv"     "cmsrvp"    "crpdwk"    "crpdwkp"   "dngdk"    
##  [853] "dngdkp"    "dngna"     "dngnap"    "dngnapp"   "dngoth"    "dngothp"  
##  [859] "dngref"    "dngrefp"   "domicil"   "dsbld"     "dsbldp"    "dvrcdev"  
##  [865] "dvrcdeva"  "edctn"     "edctnp"    "edufld"    "edulvla"   "edulvlb"  
##  [871] "edulvlfa"  "edulvlfb"  "edulvlma"  "edulvlmb"  "edulvlpa"  "edulvlpb" 
##  [877] "eduyrs"    "eisced"    "eiscedf"   "eiscedm"   "eiscedp"   "emplno"   
##  [883] "emplnof"   "emplnom"   "emplnop"   "emplrel"   "emprelp"   "emprf14"  
##  [889] "emprm14"   "estsz"     "fxltph"    "faalv"     "hhmodwl"   "hincfel"  
##  [895] "hincsrc"   "hincsrca"  "hinctnt"   "hinctnta"  "hswrk"     "hswrkp"   
##  [901] "iccohbt"   "icmnart"   "icmsw"     "icomdng"   "icomdnp"   "icpdwk2"  
##  [907] "icpdwrk"   "icppdwk"   "icwhct"    "inttph"    "ioactp"    "iorgact"  
##  [913] "isco08"    "isco08p"   "iscoco"    "iscocop"   "jbspv"     "jbspvf"   
##  [919] "jbspvm"    "jbspvp"    "lvghw"     "lvghwa"    "lvgoptn"   "lvgptn"   
##  [925] "lvgptna"   "lvgptne"   "lvgptnea"  "mainact"   "marital"   "maritala" 
##  [931] "maritalb"  "marsts"    "mbltph"    "mbtru"     "mnactic"   "mnactp"   
##  [937] "moalv"     "njbspv"    "njbspvp"   "occf14"    "occf14a"   "occf14b"  
##  [943] "occf14ie"  "occm14"    "occm14a"   "occm14b"   "occm14ie"  "orgwrk"   
##  [949] "partner"   "pdjobev"   "pdjobyr"   "pdwrk"     "pdwrkp"    "rmhhus"   
##  [955] "rshpsts"   "rtdsbyr"   "rtrd"      "rtrdp"     "tporgwk"   "uemp12m"  
##  [961] "uemp3m"    "uemp5yr"   "uempla"    "uemplap"   "uempli"    "uemplip"  
##  [967] "wkdcorp"   "wkhtot"    "wkhtotp"   "wrkac6m"   "wrkctr"    "wrkctra"  
##  [973] "wrkctrhu"  "wkdcorg"   "wkdcorga"  "nacer1"    "nacer11"   "nacer2"   
##  [979] "edagegb"   "edlvat"    "edlveat"   "eduat1"    "eduat2"    "edlvbe"   
##  [985] "edlvabe"   "edlvdbe"   "edlvebe"   "edlvbg"    "edlvdbg"   "edlvebg"  
##  [991] "edlvch"    "edlvach"   "edlvbch"   "edlvcch"   "edlvdch"   "edlvcy"   
##  [997] "edlvacy"   "edlvdcy"   "edlvecy"   "edlvgcy"   "edlvcz"    "edlvdcz"  
## [1003] "edlvdal"   "edlvde"    "edlvade"   "edlvdfi"   "edlvdis"   "edlvdk"   
## [1009] "edlvadk"   "edlvddk"   "edlvdxk"   "edlvee"    "edlvaee"   "edlvbee"  
## [1015] "edlvdee"   "edlves"    "edlvaes"   "edlvdes"   "edlvees"   "edlvges"  
## [1021] "edlvhes"   "edlvfr"    "edlvafr"   "edlvbfr"   "edlvdfr"   "edlvgb"   
## [1027] "edlvagb"   "edlvgr"    "edlvagr"   "edlvdgr"   "edlvegr"   "edlvhr"   
## [1033] "edlvdhr"   "edlvehr"   "edlvhu"    "edlvahu"   "edlvbhu"   "edlvdhu"  
## [1039] "edlvdahu"  "edlvie"    "edlvaie"   "edlvbie"   "edlvdie"   "edlvil"   
## [1045] "edlvail"   "eduil1"    "eduail1"   "edubil1"   "edlvit"    "edlvait"  
## [1051] "edlvdit"   "edlveit"   "edlvlt"    "edlvdlt"   "edlvlu"    "edlvlv"   
## [1057] "edlvdlv"   "edlvelv"   "edlvnl"    "edlvdnl"   "edlvenl"   "edlvno"   
## [1063] "edlvano"   "edlvdno"   "edlveno"   "edlvpl"    "edlvapl"   "edlvbpl"  
## [1069] "edlvdpl"   "edlvepl"   "edlvgpl"   "edlvhpl"   "edlvpt"    "edlvapt"  
## [1075] "edlvbpt"   "edlvdpt"   "edlvro"    "edlvru"    "edlvdru"   "edlvse"   
## [1081] "edlvase"   "edlvdse"   "edlvsi"    "edlvasi"   "edlvdsi"   "edlvesi"  
## [1087] "edlvsk"    "edlvask"   "edlvdsk"   "edlvtr"    "edlvua"    "edlvaua"  
## [1093] "edlvdua"   "edude1"    "eduade1"   "edubde1"   "educde1"   "edude2"   
## [1099] "eduade2"   "edubde2"   "edude3"    "eduade3"   "edugb2"    "eduagb2"  
## [1105] "edubgb2"   "eduil2"    "eduail2"   "edupl2"    "eduyrpl"   "edugb1"   
## [1111] "eduagb1"   "edubgb1"   "educgb1"   "edlvdrs"   "edlvdme"   "edagepgb" 
## [1117] "edlvpat"   "edlvpeat"  "edlvpfat"  "edlvpbe"   "edlvpdbe"  "edlvpebe" 
## [1123] "edlvpch"   "edlvpdch"  "edlvpcy"   "edlvpdcy"  "edlvpecy"  "edlvpgcy" 
## [1129] "edlvpcz"   "edlvpdcz"  "edlvpdbg"  "edlvpebg"  "edlvpdfi"  "edlvpdis" 
## [1135] "edlvpdit"  "edlvpeit"  "edlvpdk"   "edlvpddk"  "edlvpdxk"  "edlvpee"  
## [1141] "edlvpdee"  "edlvpes"   "edlvpdes"  "edlvpees"  "edlvpfes"  "edlvpges" 
## [1147] "edlvpfr"   "edlvpdfr"  "edlvpgb"   "edlvpgr"   "edlvpdgr"  "edlvpegr" 
## [1153] "edlvphr"   "edlvpdhr"  "edlvpehr"  "edlvphu"   "edlvpdhu"  "edlvpdahu"
## [1159] "edlvpie"   "edlvpdie"  "edlvpil"   "edupil1"   "edupail1"  "edupbil1" 
## [1165] "edlvplt"   "edlvpdlt"  "edlvplv"   "edlvpdlv"  "edlvpelv"  "edlvpnl"  
## [1171] "edlvpdnl"  "edlvpenl"  "edlvpno"   "edlvpdno"  "edlvpeno"  "edlvppl"  
## [1177] "edlvpdpl"  "edlvpepl"  "edlvpfpl"  "edlvpgpl"  "edlvppt"   "edlvpdpt" 
## [1183] "edlvpro"   "edlvpru"   "edlvpdru"  "edlvpse"   "edlvpdse"  "edlvpsi"  
## [1189] "edlvpdsi"  "edlvpesi"  "edlvpsk"   "edlvpdsk"  "edlvptr"   "edlvpua"  
## [1195] "edlvpdua"  "edupde1"   "edupade1"  "edupbde1"  "edupcde1"  "edupde2"  
## [1201] "edupade2"  "edupbde2"  "edupde3"   "edupade3"  "edupgb2"   "edupagb2" 
## [1207] "edupbgb2"  "edupil2"   "edupail2"  "eduppl2"   "eduyrppl"  "edlvpdal" 
## [1213] "edupgb1"   "edupagb1"  "edupbgb1"  "edupcgb1"  "edlvpdrs"  "edlvpdme" 
## [1219] "edlvfdal"  "edlvfat"   "edlvfeat"  "edlvfbe"   "edlvfdbe"  "edlvfebe" 
## [1225] "edlvfdbg"  "edlvfebg"  "edlvfhr"   "edlvfdhr"  "edlvfehr"  "edlvfcy"  
## [1231] "edlvfdcy"  "edlvfecy"  "edlvfgcy"  "edlvfcz"   "edlvfdcz"  "edlvfdk"  
## [1237] "edlvfddk"  "edlvfee"   "edlvfdee"  "edlvfdfi"  "edlvffr"   "edlvfdfr" 
## [1243] "edufde1"   "edufade1"  "edufbde1"  "edufde2"   "edufade2"  "edufde3"  
## [1249] "edufade3"  "edlvfgr"   "edlvfdgr"  "edlvfegr"  "edlvfhu"   "edlvfdhu" 
## [1255] "edlvfdahu" "edlvfdis"  "edlvfie"   "edlvfdie"  "edlvfil"   "edufil1"  
## [1261] "edufail1"  "edufbil1"  "edufil2"   "edufail2"  "edlvfdit"  "edlvfeit" 
## [1267] "edlvfdxk"  "edlvflv"   "edlvfdlv"  "edlvflt"   "edlvfdlt"  "edlvfnl"  
## [1273] "edlvfdnl"  "edlvfenl"  "edlvfno"   "edlvfdno"  "edlvfeno"  "edlvfpl"  
## [1279] "edlvfdpl"  "edlvfepl"  "edlvffpl"  "edlvfpt"   "edlvfdpt"  "edlvfro"  
## [1285] "edlvfru"   "edlvfdru"  "edlvfsk"   "edlvfdsk"  "edlvfsi"   "edlvfdsi" 
## [1291] "edlvfesi"  "edlvfes"   "edlvfdes"  "edlvfees"  "edlvffes"  "edlvfse"  
## [1297] "edlvfdse"  "edlvfch"   "edlvfdch"  "edlvftr"   "edlvfua"   "edlvfdua" 
## [1303] "edlvfgb"   "edufgb1"   "edufagb1"  "edufbgb1"  "edufcgb1"  "edufgb2"  
## [1309] "edufagb2"  "edufbgb2"  "edagefgb"  "edlvfdrs"  "edlvfdme"  "edagemgb" 
## [1315] "edlvmat"   "edlvmeat"  "edlvmbe"   "edlvmdbe"  "edlvmebe"  "edlvmch"  
## [1321] "edlvmdch"  "edlvmcy"   "edlvmdcy"  "edlvmecy"  "edlvmgcy"  "edlvmcz"  
## [1327] "edlvmdcz"  "edlvmdal"  "edlvmdfi"  "edlvmdis"  "edlvmdit"  "edlvmeit" 
## [1333] "edlvmdk"   "edlvmddk"  "edlvmdxk"  "edlvmee"   "edlvmdee"  "edlvmes"  
## [1339] "edlvmdes"  "edlvmees"  "edlvmfes"  "edlvmfr"   "edlvmdfr"  "edlvmgb"  
## [1345] "edlvmgr"   "edlvmdgr"  "edlvmegr"  "edlvmhr"   "edlvmdhr"  "edlvmehr" 
## [1351] "edlvmhu"   "edlvmdhu"  "edlvmdahu" "edlvmie"   "edlvmdie"  "edlvmil"  
## [1357] "edumil1"   "edumail1"  "edumbil1"  "edlvmlt"   "edlvmdlt"  "edlvmlv"  
## [1363] "edlvmdlv"  "edlvmnl"   "edlvmdnl"  "edlvmenl"  "edlvmno"   "edlvmdno" 
## [1369] "edlvmeno"  "edlvmpl"   "edlvmdpl"  "edlvmepl"  "edlvmfpl"  "edlvmpt"  
## [1375] "edlvmdpt"  "edlvmro"   "edlvmru"   "edlvmdru"  "edlvmse"   "edlvmdse" 
## [1381] "edlvmsi"   "edlvmdsi"  "edlvmesi"  "edlvmsk"   "edlvmdsk"  "edlvmtr"  
## [1387] "edlvmua"   "edlvmdua"  "edumde1"   "edumade1"  "edumbde1"  "edumde2"  
## [1393] "edumade2"  "edumde3"   "edumade3"  "edumgb1"   "edumagb1"  "edumbgb1" 
## [1399] "edumcgb1"  "edumgb2"   "edumagb2"  "edumbgb2"  "edumil2"   "edumail2" 
## [1405] "edlvmdbg"  "edlvmebg"  "edlvmdrs"  "edlvmdme"  "marstfi"   "marstgb"  
## [1411] "marstie"   "martlfr"   "rshpsfi"   "rshpsgb"   "impdiff"   "impenv"   
## [1417] "impfree"   "impfun"    "imprich"   "impsafe"   "imptrad"   "ipadvnt"  
## [1423] "ipbhprp"   "ipcrtiv"   "ipeqopt"   "ipfrule"   "ipgdtim"   "iphlppl"  
## [1429] "iplylfr"   "ipmodst"   "iprspot"   "ipshabt"   "ipstrgv"   "ipsuces"  
## [1435] "ipudrst"   "regionat"  "regionbe"  "regionbg"  "regionch"  "regioach" 
## [1441] "regioncy"  "regioncz"  "regioacz"  "regionde"  "regiondk"  "regioadk" 
## [1447] "regionee"  "regiones"  "regioaes"  "regionfi"  "regioafi"  "regionfr" 
## [1453] "regiongb"  "regiongr"  "regioagr"  "regionhr"  "regionhu"  "regionie" 
## [1459] "regioaie"  "regiobie"  "regionil"  "regionis"  "regionit"  "regionlt" 
## [1465] "regionlu"  "regionlv"  "regionnl"  "regionno"  "regionpl"  "regionpt" 
## [1471] "regioapt"  "regionro"  "regionru"  "regionse"  "regionsi"  "regionsk" 
## [1477] "regiontr"  "regionua"  "regunit"   "inwtm"     "inwshh"    "inwsmm"   
## [1483] "inwyys"    "inwdd"     "inwdds"    "inwmm"     "inwmms"    "inwyr"    
## [1489] "inwyye"    "inwds"     "inwehh"    "inwemm"    "inwdde"    "inwmme"   
## [1495] "inwde"     "spltadm"   "spltadma"  "spltadmb"  "spltadmc"  "spltadmd" 
## [1501] "spltadme"  "spltadmf"  "supqad1"   "supqadm"   "supqad2"   "supqdd"   
## [1507] "supqmm"    "supqyr"    "intewde"   "mode"      "actvgrs"   "agertr"   
## [1513] "bseftwk"   "bsmw"      "chdohwk"   "cld12hh"   "cldcrmr"   "cldcrot"  
## [1519] "cldnhh"    "cldnhhd"   "cldnhhg"   "cldnhhn"   "clfncrc"   "clfncsp"  
## [1525] "clhwkrc"   "clhwksp"   "clmrlx"    "clnhhbo"   "clnhhby"   "clnhhyb"  
## [1531] "crspce"    "cuthheq"   "dfcnswk"   "dfcnswka"  "dsdclve"   "dsgrhwk"  
## [1537] "dsgrmny"   "dsgrmnya"  "dsgrpwk"   "edul12m"   "edupdem"   "emplrjb"  
## [1543] "flthmcc"   "fmpvtst"   "frshrst"   "frsptjb"   "fstjbyr"   "fthcncr"  
## [1549] "gdsprt"    "grspay"    "grspaya"   "gvprppv"   "gwdvhwk"   "gwhhprc"  
## [1555] "hlpcowk"   "hlthrwk"   "hmeqphw"   "hwkmono"   "hwkpwd1"   "hwkpwd2"  
## [1561] "hwkpwdp"   "hwkpwe1"   "hwkpwe2"   "hwkpwep"   "hwkstrs"   "hwktwd1"  
## [1567] "hwktwd2"   "hwktwe1"   "hwktwe2"   "hwwkhs"    "icag45y"   "icago45"  
## [1573] "icagu70"   "icagu70a"  "icchld"    "iccldnh"   "icempl"    "icemplr"  
## [1579] "icgndr"    "ichwk1"    "ichwk2"    "icmnac"    "icmnact"   "icnopfm"  
## [1585] "icnopfma"  "icpart1"   "icpart2"   "icpart3"   "icptn"     "icptnwk"  
## [1591] "icptnwka"  "indisdc"   "ipjbhin"   "ipjbini"   "ipjbprm"   "ipjbscr"  
## [1597] "ipjbtro"   "ipjbwfm"   "jbcoedu"   "jbedyrs"   "jblrn"     "jbrqlrn"  
## [1603] "jbscr"     "jbtmppm"   "jbtsktm"   "lfintr"    "linwk3y"   "lkafohh"  
## [1609] "lscjb3y"   "lsintjb"   "mlohinc"   "mnrsefw"   "mnrsphm"   "nbsrsp"   
## [1615] "netpay"    "nevdnjb"   "npemp3y"   "oprtad"    "orgfd3y"   "payprd"   
## [1621] "payprda"   "pdjbndm"   "phwwkhs"   "plnchld"   "pphincr"   "ppwwkp"   
## [1627] "prmpls"    "prntghr"   "prrfmly"   "pthcncr"   "ptmhmcc"   "ptnwkwe"  
## [1633] "puemp3y"   "pwkhsch"   "quclss"    "rdpay3y"   "rpljbde"   "rsnlvem"  
## [1639] "rtryr"     "scrsefw"   "smblvjb"   "stdask"    "stdhrsw"   "stdlvl"   
## [1645] "stdmcdo"   "stdtruf"   "stpvtfm"   "tchints"   "tchlcrt"   "tchlp"    
## [1651] "tchtruf"   "tmtowrk"   "tngdohm"   "trndnjb"   "truinwk"   "uemp3y"   
## [1657] "updhlrl"   "updhlrp"   "useknwl"   "vrtywrk"   "wgdpeft"   "wkdcpce"  
## [1663] "wkengtp"   "wkhsch"    "wkjbndm"   "wkovrtm"   "wkovtmp"   "wkshr3y"  
## [1669] "wmcpwrk"   "wntrtr"    "wpmwkcp"   "wrkengt"   "wrkhrd"    "wrkspv"   
## [1675] "wrkwe"     "yrcremp"   "yrlvptn"   "yrskdwk"   "yrspdwk"   "yrspdwka" 
## [1681] "absddng"   "accdng"    "atnoact"   "chlrnnw"   "dclvlf"    "deaimpp"  
## [1687] "dngval"    "enjstm"    "enrglot"   "flapppl"   "flclpla"   "flrms"    
## [1693] "fltanx"    "fltbrd"    "fltpcfl"   "fltrstm"   "flttrd"    "fmlenj"   
## [1699] "fmlstrs"   "fstwttv"   "hlpoth"    "hlprtrn"   "inccmp"    "inccmpw"  
## [1705] "jbintr"    "jbstrs"    "lchshcp"   "lfcllk"    "lfwrs"     "lotsgot"  
## [1711] "lrnnew"    "lrnntlf"   "nhpftr"    "optftr"    "pactlot"   "pdwrkcr"  
## [1717] "physact"   "plinsoc"   "plprftr"   "pplahlp"   "ppllfcr"   "prhlppl"  
## [1723] "pstvms"    "rcndsrv"   "rehlppl"   "sedirlf"   "stfjb"     "stflfsf"  
## [1729] "stfsdlv"   "tmabdng"   "tmdotwa"   "tmendng"   "tmimdng"   "tnapsur"  
## [1735] "trtrsp"    "trtunf"    "uempnyr"   "wkvlorg"   "wrbknrm"   "acetalv"  
## [1741] "admaimg"   "algyplv"   "aljewlv"   "allbpe"    "allbpne"   "almuslv"  
## [1747] "alpfpe"    "alpfpne"   "alwspsc"   "blncmig"   "comnlng"   "cpimpop"  
## [1753] "ctbfsmv"   "dfegcf"    "dfegcon"   "dfeghbg"   "eimgrpc"   "eimpcnt"  
## [1759] "eimrcnt"   "fclcntr"   "gvrfgap"   "gvtrimg"   "icbrnct"   "idetalv"  
## [1765] "imacrlv"   "imbghct"   "imbleco"   "imdetbs"   "imdetmr"   "imfljob"  
## [1771] "imgclg"    "imgetn"    "imgfrnd"   "imgrpc"    "imhecop"   "imrcntr"  
## [1777] "imrsprc"   "imscrlv"   "imsetbs"   "imsetmr"   "imsmrgt"   "imtcjob"  
## [1783] "imunplv"   "imwbcrm"   "imwgdwn"   "lwdscwp"   "lwpeth"    "noimbro"  
## [1789] "pplstrd"   "qfimchr"   "qfimcmt"   "qfimedu"   "qfimfml"   "qfimlng"  
## [1795] "qfimwht"   "qfimwlt"   "qfimwsk"   "rfgawrk"   "rfgbfml"   "rfgdtcn"  
## [1801] "rfgfrpc"   "rfggvfn"   "rlgueim"   "shrrfg"    "smctmbe"   "smegbhw"  
## [1807] "smegbli"   "stimrdt"   "vrtrlg"    "cltodm"    "cltofrd"   "cltommb"  
## [1813] "cltona"    "cltonn"    "cltoptp"   "cltoref"   "cltovw"    "cnsodm"   
## [1819] "cnsofrd"   "cnsommb"   "cnsona"    "cnsonn"    "cnsoptp"   "cnsoref"  
## [1825] "cnsovw"    "discpol"   "empl"      "epaodm"    "epaofrd"   "epaommb"  
## [1831] "epaona"    "epaonn"    "epaoptp"   "epaoref"   "epaovw"    "hlpppl"   
## [1837] "hmnodm"    "hmnofrd"   "hmnommb"   "hmnona"    "hmnonn"    "hmnoptp"  
## [1843] "hmnoref"   "hmnovw"    "impapol"   "impavo"    "impfml"    "impfrds"  
## [1849] "implsrt"   "impoblw"   "impopin"   "imppol"    "imprlg"    "imprwcr"  
## [1855] "imprwct"   "imprwkc"   "impsppl"   "impvo"     "impvote"   "impwrk"   
## [1861] "othvdm"    "othvfrd"   "othvmmb"   "othvna"    "othvnn"    "othvptp"  
## [1867] "othvref"   "othvvw"    "prfodm"    "prfofrd"   "prfommb"   "prfona"   
## [1873] "prfonn"    "prfoptp"   "prforef"   "prfovw"    "prtydm"    "prtyfrd"  
## [1879] "prtymmb"   "prtyna"    "prtynn"    "prtyptp"   "prtyref"   "prtyvw"   
## [1885] "rlgodm"    "rlgofrd"   "rlgommb"   "rlgona"    "rlgonn"    "rlgoptp"  
## [1891] "rlgoref"   "rlgovw"    "sclcdm"    "sclcfrd"   "sclcmmb"   "sclcna"   
## [1897] "sclcnn"    "sclcptp"   "sclcref"   "sclcvw"    "setodm"    "setofrd"  
## [1903] "setommb"   "setona"    "setonn"    "setoptp"   "setoref"   "setovw"   
## [1909] "sptcdm"    "sptcfrd"   "sptcmmb"   "sptcna"    "sptcnn"    "sptcptp"  
## [1915] "sptcref"   "sptcvw"    "stfhwkp"   "strtbsn"   "trudm"     "trufrd"   
## [1921] "truiwkp"   "trummb"    "truna"     "trunn"     "truptp"    "truref"   
## [1927] "trusay"    "truvw"     "truwrkp"   "wkchtsk"   "wkdcsin"   "wkenvin"  
## [1933] "wrkflex"   "yrlvdae"   "advbach"   "advhach"   "advslep"   "advsthr"  
## [1939] "chsrgp"    "dcadmms"   "dcdfcwr"   "dcdisc"    "dckptrt"   "dctreql"  
## [1945] "ddprsmd"   "dsplvpr"   "follwdc"   "hltherb"   "illcure"   "mdhair"   
## [1951] "mdhappy"   "mdlswgt"   "mdmemo"    "mdsdeff"   "mdsexlf"   "pplcure"  
## [1957] "pprlydc"   "prcbach"   "prchach"   "prcslep"   "prcsthr"   "prfmddc"  
## [1963] "prfsmdc"   "psmdcpr"   "ptnrlcq"   "tmcnsdc"   "usmdprs"   "usmprse"  
## [1969] "bnkfldl"   "bsnprft"   "cmprcti"   "ctzchtx"   "ctzhlpo"   "flgvbnf"  
## [1975] "flinsr"    "flinsrw"   "fodcncl"   "frdbnft"   "frmwktg"   "icsbfm"   
## [1981] "ignrlaw"   "kptchng"   "mnyacth"   "musdocm"   "olwmsop"   "payavtx"  
## [1987] "pboafvr"   "pbofvr"    "pbofvrw"   "pyavtxw"   "rprochg"   "scbevts"  
## [1993] "scndhfl"   "slcnflw"   "slcnsfl"   "tstfnch"   "tstpboh"   "tstrprh"  
## [1999] "wrytrdh"   "acldnmr"   "adftjob"   "adllvhm"   "adlvptn"   "adpnt"    
## [2005] "advcyc"    "aftjbyc"   "ageadlt"   "agemage"   "ageoage"   "alvgptn"  
## [2011] "anvcld"    "brnmm"     "bthcld"    "evlvptn"   "evmar"     "evpdemp"  
## [2017] "fcldbrn"   "ggchld"    "iaglptn"   "iaglvmr"   "iagpnt"    "iagrtr"   
## [2023] "lvpntyr"   "lvptnyr"   "maryr"     "nbthcld"   "ngchld"    "oldfrl"   
## [2029] "oldgpnt"   "oldhlpo"   "pdempyr"   "plnftr"    "rcldnmr"   "rdvcyc"   
## [2035] "rftjbyc"   "rlvgptn"   "rnvcld"    "rpntyng"   "rspslvo"   "rwkold"   
## [2041] "svclvo"    "tochld"    "tolvpnt"   "towkht"    "tygledu"   "tyglvmr"  
## [2047] "tyglvp"    "tygpnt"    "tygrtr"    "tygsexi"   "wrinco"    "ycldbyr"  
## [2053] "ygcdbyr"   "admub"     "basinc"    "bennent"   "bnlwinc"   "brnocnt"  
## [2059] "cldcrsv"   "dcndleq"   "dfincac"   "ditxssp"   "earnpen"   "earnueb"  
## [2065] "eduunmp"   "eudcnbf"   "eusclbf"   "gvcldcr"   "gvhlthc"   "gvjbevn"  
## [2071] "gvpdlwk"   "gvslvol"   "gvslvue"   "hlthcef"   "imrccon"   "imsclbn"  
## [2077] "insfben"   "lbenent"   "lklpwcf"   "lknemny"   "lknhlcn"   "lkuemp"   
## [2083] "lvoapen"   "lvpbhlt"   "nmnybsc"   "prtsick"   "sbbsntx"   "sbcwkfm"  
## [2089] "sbenccm"   "sbeqsoc"   "sblazy"    "sblwcoa"   "sblwlka"   "sbprvpv"  
## [2095] "sbstrec"   "schtaut"   "sickdsb"   "slvpens"   "slvuemp"   "smdfslv"  
## [2101] "txadleq"   "txautef"   "txearn"    "ub20edu"   "ub20pay"   "ub20unp"  
## [2107] "ub50edu"   "ub50pay"   "ub50unp"   "ubedu"     "ubpay"     "ubspedu"  
## [2113] "ubsppay"   "ubspunp"   "ubunp"     "uemplwk"   "uentrjb"   "wrkprbf"  
## [2119] "yngfnjb"   "acq30bs"   "acq70bs"   "agdcold"   "agdcyng"   "agegrbl"  
## [2125] "ageyo30"   "agrpblg"   "buproag"   "chg1530"   "dprschg"   "dprsf70"  
## [2131] "dprso70"   "dprsy30"   "frndo70"   "frndy30"   "lkrspag"   "mbfmo70"  
## [2137] "oafl20"    "oafl70"    "p20cnec"   "p20cust"   "p70cnec"   "p70cust"  
## [2143] "p70hltb"   "predage"   "predetn"   "predsex"   "pvolwrk"   "sgr2070"  
## [2149] "srdscag"   "stsp20"    "stsp40"    "stsp70"    "suproag"   "trtbdag"  
## [2155] "twcol20"   "twcol70"   "v20adm"    "v20cntm"   "v20comp"   "v20envy"  
## [2161] "v20frnd"   "v20mrst"   "v20pity"   "v20resp"   "v70adm"    "v70cntm"  
## [2167] "v70comp"   "v70envy"   "v70frnd"   "v70mrst"   "v70pity"   "v70resp"  
## [2173] "wr20crm"   "wr20emp"   "bplcdc"    "bstln5y"   "bystlct"   "bystlwr"  
## [2179] "caplcst"   "ctfrdc"    "ctinplt"   "ctjbcnt"   "ctmtgfr"   "ctprpwr"  
## [2185] "dbctvrd"   "doplcsy"   "dpcstrb"   "flsin5y"   "gsupplc"   "insclct"  
## [2191] "insclwr"   "jdgcbrb"   "lwstrob"   "plcarcr"   "plccbrb"   "plccbrg"  
## [2197] "plccont"   "plcexdc"   "plcfrdc"   "plciplt"   "plcipvl"   "plcjbcn"  
## [2203] "plcpvcr"   "plcrgwr"   "plcrspc"   "plcstf"    "plcvcrc"   "plcvcrp"  
## [2209] "rgbrklw"   "stcbg2t"   "tmprs"     "trfoct"    "trfowr"    "troff5y"  
## [2215] "wevdct"    "widprsn"   "wraccrc"   "wraccrp"   "chpldm"    "chpldmc"  
## [2221] "chpldmi"   "ctstogv"   "cttresa"   "cttresac"  "dfprtal"   "dfprtalc" 
## [2227] "dspplvt"   "dspplvtc"  "fairelc"   "fairelcc"  "fplvdm"    "fplvdmc"  
## [2233] "fplvdmi"   "gptpelc"   "gptpelcc"  "grdfinc"   "grdfincc"  "gvcodmc"  
## [2239] "gvcodmi"   "gvctzpv"   "gvctzpvc"  "gvexpdc"   "gvexpdcc"  "gvspcdm"  
## [2245] "gvspdmc"   "gvspdmi"   "imvtctz"   "medcrgv"   "medcrgvc"  "meprinf"  
## [2251] "meprinfc"  "oppcrgv"   "oppcrgvc"  "pltavie"   "pltaviec"  "pplvdmc"  
## [2257] "pplvdmi"   "rghmgpr"   "rghmgprc"  "stpldmc"   "stpldmi"   "votedir"  
## [2263] "votedirc"  "viepol"    "wpestop"   "keydec"    "viepolc"   "wpestopc" 
## [2269] "keydecc"   "admit"     "showcv"    "accalaw"   "impdema"   "impdemb"  
## [2275] "impdemc"   "impdemd"   "impdeme"   "alcbnge"   "alcfreq"   "alcwkdy"  
## [2281] "alcwknd"   "anypacc"   "cgtsday"   "cgtsmke"   "cnfpplh"   "dosprt"   
## [2287] "dshltdk"   "dshltgp"   "dshltms"   "dshltna"   "dshltnt"   "dshltref" 
## [2293] "eatveg"    "etfruit"   "fnsdfml"   "height"    "hlpfmhr"   "hlpfmly"  
## [2299] "hltphal"   "hltphbn"   "hltphbp"   "hltphdi"   "hltphdk"   "hltphhb"  
## [2305] "hltphhc"   "hltphna"   "hltphnap"  "hltphnt"   "hltphpa"   "hltphpf"  
## [2311] "hltphref"  "hltphsc"   "hltphsd"   "hltphsh"   "hltpral"   "hltprbn"  
## [2317] "hltprbp"   "hltprca"   "hltprdi"   "hltprdk"   "hltprhb"   "hltprhc"  
## [2323] "hltprna"   "hltprnt"   "hltprpa"   "hltprpf"   "hltprref"  "hltprsc"  
## [2329] "hltprsd"   "hltprsh"   "jbexebs"   "jbexecp"   "jbexedk"   "jbexena"  
## [2335] "jbexenap"  "jbexent"   "jbexera"   "jbexeref"  "jbexevc"   "jbexevh"  
## [2341] "jbexevl"   "jbexpdk"   "jbexpmc"   "jbexpml"   "jbexpna"   "jbexpnap" 
## [2347] "jbexpnt"   "jbexpref"  "jbexpti"   "jbexpvi"   "medtrdk"   "medtrna"  
## [2353] "medtrnap"  "medtrnl"   "medtrnp"   "medtrnt"   "medtrnu"   "medtrnaa" 
## [2359] "medtroc"   "medtroth"  "medtrref"  "medtrun"   "medtrwl"   "trhltacp" 
## [2365] "trhltacu"  "trhltch"   "trhltcm"   "trhltdk"   "trhltho"   "trhltht"  
## [2371] "trhlthy"   "trhltmt"   "trhltna"   "trhltnt"   "trhltos"   "trhltpt"  
## [2377] "trhltre"   "trhltref"  "trhltsh"   "weight"    "icgndra"   "banhhap"  
## [2383] "ccgdbd"    "cflsenr"   "clmchng"   "clmthgt1"  "clmthgt2"  "elgbio"   
## [2389] "elgcoal"   "elghydr"   "elgngas"   "elgnuc"    "elgsun"    "elgwind"  
## [2395] "eneffap"   "gvsrdcc"   "inctxff"   "lklmten"   "lkredcc"   "ownrdcc"  
## [2401] "rdcenr"    "sbsrnen"   "wrdpfos"   "wrdpimp"   "wrenexp"   "wrinspw"  
## [2407] "wrntdis"   "wrpwrct"   "wrtcfl"    "wrtratc"   "btminfr"   "evfredu"  
## [2413] "evfrjob"   "frlgrsp"   "frlneti"   "frprtpl"   "fvgabc"    "grspfr"   
## [2419] "grsplet"   "grspnum"   "gvintcz"   "ifredu"    "ifrjob"    "infqbst"  
## [2425] "jstprev"   "netifr"    "netilet"   "netinum"   "occinfr"   "pcmpinj"  
## [2431] "poltran"   "ppldsrv"   "recexp"    "recgndr"   "recimg"    "recknow"  
## [2437] "recskil"   "sofrdst"   "sofrpr"    "sofrprv"   "sofrwrk"   "topinfr"  
## [2443] "wltdffr"   "iincsrc"   "acchome"   "accwrk"    "accmove"   "accoth"   
## [2449] "accnone"   "accref"    "accdk"     "accna"     "fampref"   "famadvs"  
## [2455] "fampdf"    "mcclose"   "mcinter"   "mccoord"   "mcpriv"    "mcmsinf"  
## [2461] "chldo12"   "gndro12a"  "gndro12b"  "ageo12"    "hhlio12"   "closeo12" 
## [2467] "ttmino12"  "speako12"  "como12"    "scrno12"   "phoneo12"  "livpnt"   
## [2473] "pntmofa"   "agepnt"    "hhlipnt"   "closepnt"  "ttminpnt"  "speakpnt" 
## [2479] "compnt"    "scrnpnt"   "phonepnt"  "wrkhome"   "wrklong"   "wrkresp"  
## [2485] "scrhlm"    "sclmsup"   "mansupp"   "manhlp"    "manwrkpl"  "manspeak" 
## [2491] "mancom"    "manscrn"   "manphone"  "teamfeel"  "wrkextra"  "colprop"  
## [2497] "colhlp"    "colspeak"  "colcom"    "colscrn"   "colphone"  "mcwrkhom" 
## [2503] "c19spo12"  "c19mco12"  "c19sppnt"  "c19mcpnt"  "c19whome"  "c19wplch" 
## [2509] "c19whacc"  "c19spwrk"  "c19mcwrk"  "secgrdec"  "scidecpb"  "admc19"   
## [2515] "panpriph"  "panmonpb"  "govpriph"  "govmonpb"  "panfolru"  "panclobo" 
## [2521] "panresmo"  "gvhanc19"  "gvjobc19"  "gveldc19"  "gvfamc19"  "hscopc19" 
## [2527] "gvbalc19"  "gvimpc19"  "gvconc19"  "respc19"   "reshhc19"  "hapljc19" 
## [2533] "hapirc19"  "hapwrc19"  "hapfuc19"  "hapfoc19"  "hapnoc19"  "hapnwc19" 
## [2539] "hapnpc19"  "haprec19"  "hapdkc19"  "hapnac19"  "icvacc19"  "getavc19" 
## [2545] "getnvc19"  "vdcond"    "vdovexre"  "vdtype"    "vdtpsvre"  "vdtpitre" 
## [2551] "vdtpscre"  "vdtpaure"  "vdtpvire"  "vdtpoire"  "vdtpntre"  "vdtpapre" 
## [2557] "vdtprere"  "vdtpdkre"  "vdtpnare"  "cldgng"    "enjlf"     "fltdpr"   
## [2563] "flteeff"   "fltlnl"    "fltlnla"   "fltsd"     "hrshsnt"   "hrshsnta" 
## [2569] "mnrgtjb"   "pdaprp"    "pdaprpa"   "slprl"     "smbtjob"   "smbtjoba" 
## [2575] "stfjbot"   "wrhpp"     "stfmjob"   "wrywprb"   "trdawrk"   "jbprtfp"  
## [2581] "pfmfdjb"   "pfmfdjba"  "dcsfwrk"   "dcsfwrka"  "implvdm"   "ccrdprs"  
## [2587] "ccnthum"   "wrclmch"   "yrbrn"     "wkhct"     "region"    "admdw"    
## [2593] "admge"     "iagmr"     "tygmr"     "domain"    "trstsci"   "prtvtmk"  
## [2599] "prtclmk"   "rlgdnmk"   "admrclc"   "testic34"  "testic35"  "testic36" 
## [2605] "testic37"  "testic38"  "testic39"  "testic40"  "testic41"  "testic42" 
## [2611] "rshpsmk"   "marstmk"   "edlvdmk"   "edlvpdmk"  "edlvfdmk"  "edlvmdmk" 
## [2617] "testii1"   "testii2"   "testii3"   "testii4"   "testii5"   "testii6"  
## [2623] "testii7"   "testii8"   "testii9"   "ainws"     "ainwe"     "binwe"    
## [2629] "cinwe"     "dinwe"     "finwe"     "ginwe"     "hinwe"     "iinwe"    
## [2635] "kinwe"     "vinwe"     "jinws"     "jinwe"     "scrlgblg"  "scchpldm" 
## [2641] "scwflt"    "scrwiat"   "sctmfeel"  "scrwwc"    "scclprp"   "scwsds"   
## [2647] "scwass"    "scwa5se"   "scwa11se"  "scwa59se"  "scwa89se"  "scwb30se" 
## [2653] "scwc55se"  "scwd62se"  "scwd78se"  "scwrecse"  "scwtype"   "uagent"   
## [2659] "smplresp"  "questcmp"  "intagea"   "intgndr"   "intnum"    "intlnga"  
## [2665] "rlgdemk"
# View relevant variable names
grep("trad|rule|safe|gov|obey|behav", names(france_data), value = TRUE)
## [1] "lawobey"  "stfgov"   "impsafe"  "imptrad"  "ipfrule"  "rlgovw"   "govpriph"
## [8] "govmonpb"
france_data <- france_data %>%
  mutate(across(
    c(lawobey, stfgov, impsafe, imptrad, ipfrule),
    ~ ifelse(.x %in% c(77, 88, 99), NA, .x)
  ))
# Example reverse-coding (if necessary)
france_data <- france_data %>%
  mutate(
    stfgov = max(stfgov, na.rm = TRUE) + 1 - stfgov
  )
france_data <- france_data %>%
  mutate(
    authoritarian_raw = rowMeans(select(., lawobey, stfgov, impsafe, imptrad, ipfrule), na.rm = TRUE),
    authoritarian = scales::rescale(authoritarian_raw, to = c(0, 100), na.rm = TRUE)
  )
summary(france_data$authoritarian)
##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
##    0.00   29.41   38.24   38.78   47.06  100.00
hist(france_data$authoritarian, breaks = 30, main = "Distribution of Authoritarian Values", xlab = "Score (0–100)")

france_data <- france_data %>%
  mutate(
    Education = case_when(
      eisced %in% 0:2 ~ "Lower Secondary or Less",
      eisced %in% 3:4 ~ "Upper Secondary",
      eisced %in% 5:6 ~ "Tertiary Education",
      TRUE ~ NA_character_
    ),
    Education = factor(Education, levels = c("Lower Secondary or Less", "Upper Secondary", "Tertiary Education"))
  )
library(ggplot2)

france_data %>%
  filter(!is.na(Education), !is.na(authoritarian)) %>%
  ggplot(aes(x = Education, y = authoritarian, fill = Education)) +
  geom_violin(trim = FALSE) +
  geom_boxplot(width = 0.1, fill = "white", outlier.shape = NA) +
  labs(
    title = "Distribution of Authoritarian Values by Education Level",
    subtitle = "Based on Schwartz Value Indicators (ESS, France 2018)",
    x = "Education Level",
    y = "Authoritarian Values (0–100)",
    caption = "Source: European Social Survey, France"
  ) +
  theme_minimal() +
  theme(legend.position = "none", axis.text.x = element_text(angle = 30, hjust = 1))

Interpretation of Violin Plot The violin plot displays the distribution of authoritarian values across three education levels in France based on data from the 2018 European Social Survey. The most noticeable trend is that individuals with lower levels of education tend to score higher on authoritarian values, while those with tertiary education show the lowest scores.

The median score is highest for those with Lower Secondary or Less, slightly lower for Upper Secondary, and notably lower for Tertiary Education. This suggests a negative relationship between education level and authoritarian attitudes.

The spread (interquartile range) is fairly similar across all groups, but the Lower Secondary or Less group displays a slightly wider distribution, indicating greater variability in authoritarian views. Additionally, this group shows the most pronounced tail in the upper range, suggesting a subset of individuals with especially strong authoritarian tendencies.

Meanwhile, the Tertiary Education group has the narrowest spread and lowest overall values, suggesting more consistent and less authoritarian outlooks.

Overall, the plot suggests that higher education may be associated with lower authoritarian values and possibly more openness or liberal attitudes. This aligns with theories that education enhances critical thinking and tolerance, reducing tendencies toward authoritarianism.

model1 <- lm(authoritarian ~ Education, data = france_data)
france_data <- france_data %>%
  mutate(
    Residence = case_when(
      domicil %in% c(1, 2) ~ "Urban",   # big city or suburbs
      domicil %in% c(3, 4, 5) ~ "Rural", # town or village/farm
      TRUE ~ NA_character_
    ),
    Residence = factor(Residence, levels = c("Urban", "Rural"))
  )
model2 <- lm(authoritarian ~ Education + Residence, data = france_data)
# Create Economic Status variable
france_data <- france_data %>%
  mutate(
    EconomicStatus = case_when(
      hincfel == 1 ~ "Economic Comfort",  # Living comfortably
      hincfel %in% 2:4 ~ "Economic Strain",  # Coping to very difficult
      TRUE ~ NA_character_
    ),
    EconomicStatus = factor(EconomicStatus, levels = c("Economic Comfort", "Economic Strain"))
  )
# Model 3: Education + Residence + Economic Status
model3 <- lm(authoritarian ~ Education + Residence + EconomicStatus, data = france_data)
library(sjPlot)

tab_model(
  model1, model2, model3,
  show.ci = FALSE,
  show.se = TRUE,
  show.aic = TRUE,
  show.obs = TRUE,
  dv.labels = c("Model 1: Education", "Model 2: + Residence", "Model 3: + Economic Status"),
  title = "Regression Models Predicting Authoritarian Values (France 2018)"
)
Regression Models Predicting Authoritarian Values (France 2018)
  Model 1: Education Model 2: + Residence Model 3: + Economic Status
Predictors Estimates std. Error p Estimates std. Error p Estimates std. Error p
(Intercept) 36.54 0.15 <0.001 35.83 0.21 <0.001 35.16 0.31 <0.001
Education [Upper
Secondary]
3.29 0.22 <0.001 3.26 0.22 <0.001 2.58 0.24 <0.001
Education [Tertiary
Education]
4.37 0.27 <0.001 4.43 0.27 <0.001 4.02 0.29 <0.001
Residence [Rural] 1.00 0.22 <0.001 1.23 0.23 <0.001
EconomicStatus [Economic
Strain]
1.56 0.23 <0.001
Observations 17090 17085 13737
R2 / R2 adjusted 0.020 / 0.020 0.022 / 0.022 0.018 / 0.018
AIC 135521.103 135454.015 107445.932
library(knitr)

aic_bic_table <- data.frame(
  Model = c("Model 1", "Model 2", "Model 3"),
  AIC = c(AIC(model1), AIC(model2), AIC(model3)),
  BIC = c(BIC(model1), BIC(model2), BIC(model3))
)

kable(aic_bic_table, caption = "AIC and BIC Values for All Models")
AIC and BIC Values for All Models
Model AIC BIC
Model 1 135521.1 135552.1
Model 2 135454.0 135492.7
Model 3 107445.9 107491.1
plot_model(model3, 
           show.values = TRUE, 
           value.offset = 0.3, 
           vline.color = "red",
           title = "Predictors of Authoritarian Values in France (Final Model)")

This analysis explored how education, place of residence, and economic status predict authoritarian values in France using data from the 2018 European Social Survey. We built three linear regression models: Model 1 included education only, Model 2 added residence, and Model 3 added economic status.

Model 1 shows that education is a strong predictor of authoritarian values. Individuals with higher education levels (especially tertiary education) are significantly less likely to hold authoritarian views compared to those with lower secondary education or less (the reference group). This supports previous findings that education promotes critical thinking and tolerance, reducing authoritarian tendencies. The model’s R² is relatively low, suggesting education alone explains a modest portion of the variance.

In Model 2, we introduced residential area. Urban residents (big cities or suburbs) had significantly lower authoritarian values than rural residents, even after accounting for education. This aligns with literature suggesting rural communities may lean more conservative or traditional. The R² improved slightly compared to Model 1, and both AIC and BIC decreased, indicating a better model fit. Notably, the education coefficients remained stable, suggesting that residence and education are independently associated with authoritarian values.

Model 3 added economic status. Here, individuals experiencing economic strain (coping, difficult, or very difficult) showed significantly higher authoritarian values than those living comfortably. This variable had a strong and statistically significant effect. Its inclusion improved the model further, with the highest R² and lowest AIC/BIC among all three models. This means Model 3 explains the most variation and is the best-fitting model overall. Importantly, education and residence remained significant predictors even after accounting for economic stress, showing these are distinct drivers of authoritarian attitudes.

When comparing across models, we see that as additional predictors are introduced, the coefficients for education remain relatively consistent. This suggests the effect of education on authoritarianism is robust and not confounded by place or economic status. The model fit improvements—especially the drop in AIC and BIC—show that residence and financial hardship meaningfully enhance our understanding of what drives authoritarian values.

In conclusion, Model 3 offers the most complete picture. It shows that authoritarian attitudes are shaped by a combination of educational attainment, living environment, and financial stress. Those with less education, living in rural areas, and under financial strain are more likely to hold authoritarian beliefs. These insights highlight the importance of structural conditions in shaping political attitudes and may help inform social policy and civic education initiatives.

model4 <- lm(authoritarian ~ Education * Residence + EconomicStatus, data = france_data)
library(modelsummary)

model_list <- list(
  "Model 3: Add Economic Status" = model3,
  "Model 4: Interaction (Education * Residence)" = model4
)

modelsummary(
  model_list,
  statistic = "({p.value})",
  stars = TRUE,
  coef_map = c(
    "(Intercept)" = "Intercept",
    "EducationUpper Secondary" = "Upper Secondary",
    "EducationTertiary Education" = "Tertiary Education",
    "ResidenceRural" = "Rural",
    "EconomicStatusEconomic Strain" = "Economic Strain",
    "EducationUpper Secondary:ResidenceRural" = "Upper Sec. x Rural",
    "EducationTertiary Education:ResidenceRural" = "Tertiary x Rural"
  ),
  gof_map = c("nobs", "r.squared", "adj.r.squared", "AIC", "BIC"),
  title = "Comparing Models With and Without Education-Residence Interaction",
  notes = "Source: European Social Survey, France 2018"
)
Comparing Models With and Without Education-Residence Interaction
Model 3: Add Economic Status Model 4: Interaction (Education * Residence)
+ p < 0.1, * p < 0.05, ** p < 0.01, *** p < 0.001
Source: European Social Survey, France 2018
Intercept 35.163*** 34.527***
(<0.001) (<0.001)
Upper Secondary 2.582*** 3.177***
(<0.001) (<0.001)
Tertiary Education 4.018*** 5.410***
(<0.001) (<0.001)
Rural 1.225*** 2.090***
(<0.001) (<0.001)
Economic Strain 1.557*** 1.553***
(<0.001) (<0.001)
Upper Sec. x Rural -0.805
(0.138)
Tertiary x Rural -2.006**
(0.001)
Num.Obs. 13737 13737
R2 0.018 0.019
R2 Adj. 0.018 0.018

Bonus Question Interpretation (Model 4 – Interaction Effects):

The interaction model explores whether the effect of education on authoritarian values differs between urban and rural areas. In Model 4, interaction terms were added between education and residence, allowing us to assess whether the impact of education depends on whether someone lives in a rural or urban area.

The interaction term Tertiary x Rural is statistically significant and negative, suggesting that the effect of having a tertiary education on lowering authoritarian values is even stronger in rural areas than in urban ones. This implies that rural individuals with higher education levels show notably less authoritarianism than their peers with lower education levels—possibly because higher education in rural areas contrasts more sharply with the local norms or environment.

In contrast, the Upper Sec. x Rural interaction was not statistically significant, meaning the difference between upper secondary and lower education in rural areas does not significantly differ from that in urban areas.

Adding the interaction terms slightly improves the model fit (check the adjusted R², AIC, and BIC), but not dramatically. This tells us that while there is some variation in how education operates across locations, it’s not a dominant source of difference in authoritarian attitudes overall.