The empirical model was fitted based on previous literature using the STATA interaction expansion function (xi), which expands terms containing categorical variables into a dummy variable set. The independent variables explored in both the multinomial and OLS models include gender, ethnicity, marital status, employment status, occupation category, working sector, monthly household income, health status and household size.