## ----echo=FALSE,eval=FALSE----------------------------------------------------
# rm(list=ls(all=TRUE))
# options(width=60)

## ----results='hide',eval=FALSE------------------------------------------------
# library(catdata)
# data(addiction)
# attach(addiction)

## ----eval=FALSE---------------------------------------------------------------
# ill01 <- ill
# ill01[ill==0] <- 1
# ill01[ill==2] <- 0
# 
# age2 <- age^2

## ----eval=FALSE---------------------------------------------------------------
# m01vs2 <- glm(ill01 ~ as.factor(gender) + as.factor(university) + age + age2,
# family=binomial())
# summary(m01vs2)

## ----results='hide',eval=FALSE------------------------------------------------
# detach(addiction)
# addiction2 <- addiction[addiction$ill!=2,]
# attach(addiction2)
# age2 <- age^2

## ----eval=FALSE---------------------------------------------------------------
# m0vs1 <- glm(ill ~ as.factor(gender) + as.factor(university) + age + age2,
# family=binomial())
# summary(m0vs1)