exemplo onde y é explicado por x + z mas não por z

exemplo_regressao

x <- rnorm(1000)
z <- runif(1000)
d <- data.frame(x = x, z = z, y= x + z/10 + rnorm(1000, sd=.001))
cor(d)
summary(lm(y~x+z, data=d))
summary(lm(y~x, data=d))
summary(lm(y~z, data=d))

Outro exemplo: aqui eu gero o y como uma função determinística de x, z1, z2, z3, z4, z5.
Parece que a função glm dá bug, mas a função step acaba tirando as variáveis z...

x <- rnorm(1000)
z1 <- runif(1000)
z2 <- runif(1000)
z3 <- runif(1000)
z4 <- runif(1000)
z5 <- runif(1000)
z <-runif(1000)/1000
y=ifelse(x + (z1+z2+z3+z4+z5)/1000 > .2, 1, 0)
d <- data.frame(x=x, z=z, z1=z1, z2=z2, z3=z3, z4=z4, z5=z5, y=y)
cor(d)
summary(glm(y~., data=d, family=binomial))
step(glm(y~., data=d, family=binomial))