# Discrim Example -- source code file discrim.r

# Read and print discrim data frame.
discrim = read.table("discrim.txt", header=T)
cat("discrim data frame:\n")
print(discrim)

# Compute logistic regression model.
model = glm(hire ~ educ + exp + gender,
            data=discrim, family=binomial(link="logit"))
print(summary(model))

# Print predicted values and residuals.
cat("Predicted values:\n")
print(fitted(model, type="response"))
cat("Residuals:\n")
print(residuals(model, type="deviance"))