# Remission Example -- source code file remission.r

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

# Compute logistic regression model.
fullmodel = glm(remiss ~ cell + smear + infil + li + blast + temp,
	 data=remission, family=binomial(link="logit"))
print(summary(fullmodel))

# Perform backward selection:
cat("\nBackward Selection:\n")
model = step(fullmodel, direction="backward")
print(summary(model))

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