# R-Code zu GRM-Tutorium 5 

# Datensatz befindet sich im Paket AER 
install.packages("AER")
library("AER")
data("BankWages")
data<- BankWages

# Alternativ von der Homepage 
data <- read.table("BankWages.txt",header=TRUE)
# stelle ursprüngliche Kodierungen des Datensatzes wieder her! 
data$job <- relevel(data$job,ref="custodial")
#Referenzkategorie der Variablen Gender
data$gender <- relevel(data$gender,ref="male")

# Modell zu Aufgabe 1
library("nnet")
multinomial <- multinom(job ~ education + minority + gender, data = data)
summary(multinomial)



# Modell zu Aufgabe 2 
library("MASS")
proportional <- polr(formula = job ~ education + minority, data = data,
     subset = gender == "male", Hess = "TRUE")
summary(proportional)