library(foreign); library(dplyr)
hten <- read.dta("https://teaching.sociology.ul.ie/so5032/hten.dta")
hten$tenure <- droplevels(hten$tenure); hten$jbstat <- droplevels(hten$jbstat)
hten$ten3 <- hten$tenure
hten$ten3 <- recode(hten$ten3,
                    "owned outright" = "Owner",
                    "owned with mortgage" = "Owner",
                    "local authority rented" = "Social renter",
                    "housing assoc. rented" = "Social renter",
                    "rented from employer" = "Private renter",
                    "rented private unfurnished" = "Private renter",
                    "rented private furnished" = "Private renter",
                    "other rented" = "Private renter")
hten$eun <- hten$jbstat
hten$eun <- recode(hten$eun,
                    "self-employed" = "Employed",
                   "employed" = "Employed",
                   "unemployed" = "Unemployed",
                   "retired" = "Retired",
                   "maternity leave" = "Employed",
                   "family care" = "Not in LabMkt",
                   "ft studt, school" = "Not in LabMkt",
                   "lt sick, disabld" = "Not in LabMkt",
                   "gvt trng scheme" = "Not in LabMkt",
                   "other" = "Not in LabMkt")

library(nnet)
summary(modm1 <- multinom(ten3 ~ age + eun, hten))

bhps <- read.dta("https://teaching.sociology.ul.ie/so5032/bhpsologit.dta")
bhps$rsex <- droplevels(bhps$rsex)
bhps$ropfamr <- droplevels(bhps$ropfamr)

summary(modb1 <- multinom(ropfamr ~ rsex + rage, bhps))

library(MASS)
summary(modb1 <- polr(ropfamr ~ rsex + rage, bhps))

summary(modb1)

summary(modlm <- lm(as.numeric(ropfamr) ~ rsex + rage, bhps))