# Depress Example -- source code file depress.r

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

# Compute logistic regression model.
model = glm(isdep ~ obs + age + educ + income + score,
	 data=depress, family=binomial(link="logit"))
print(summary(model))

cat("Predicted Values:\n")
print(fitted(model, type="response"))
	
cat("Residuals:\n")
print(residuals(model, type="deviance"))