rm(list = ls())
Trauma <- read.table(
url("http://stat.unm.edu/~fletcher/LLM/DATA/TRAUMAa.DAT"),
#"C:\\E-drive\\Books\\LOGLIN3\\DATA\\TRAUMAa.DAT",
		  sep="",col.names=c("ID","Death","ISS","TI","RTS","AGE"))
# Alternative file read
#"url("http://stat.unm.edu/~fletcher/LLM/DATA/TRAUMA.DAT"),
#		  sep="",header=T)
attach(Trauma)
summary(Trauma)
# From the summary, 0.93 percent "die" in the ER
# It is labeled backwards.
Death=1-Death

par(mfrow=c(2,2))
boxplot(ISS~Death,ylab="ISS")
boxplot(RTS~Death)
boxplot(AGE~Death)
par(mfrow=c(1,1))