search for: datx2

...of predict() in order to generate the same predicted probabilities as above. The same code as above does not work: p2a=predict(mod2, type="response", newdata=as.data.frame(cbind(x1, x2=0))) unique(p2a) Nor does creating a data frame that has the names "datx1" and "datx2," which is how the variables appear if you run a summary() on mod2. Looking at the model matrix of mod2 shows that the fitted model only shows two variables, the dependent variable y and one independent variable called "dat." It is as if my two variables x1 and x2 have be...