Skip to content
Prev 30269 / 398506 Next

log-linear

The presence/absence nature of the outcome variable strongly supports
using logistic regression and nothing else.  I strongly encourage you
to stick with logistic regression.  The model formula and interaction
term capabilities in R are just the same for logistic regression as for
log-linear models.  (In some textbooks, log-linear models are used as
the motivation and example for introducing the ideas of interaction
terms, but once introduced, the ideas apply very generally.)

I would set up the data as you have, as a data frame or a matrix with
columns representing the number of landslide presence cells, the number
of landslide absence cells, and then one column for each predictor.

Then use  glm() with a call something like:

result <- glm(cbind(present, absent) ~ (a+b+c+d)^3,  family=binomial,
                 data = name.of.data.frame)

In  help("glm"), there's a sentence under "Details" which describes
the cbind() syntax I've used above, and  help("formula")  explains
the (.)^3 syntax.

-  tom blackwell  -  u michigan medical school  -  ann arbor  -
On Mon, 7 Apr 2003, orkun wrote: