ch = glm(y ~ Ag+S+D+Ch+H+W,family=binomial)
chstep <- step(ch, direction="backward")
chstep