Message-ID: <Pine.SOL.4.44.0304070910570.10983-100000@zektor.gpcc.itd.umich.edu>
Date: 2003-04-07T13:36:26Z
From: Thomas W Blackwell
Subject: log-linear
In-Reply-To: <3E913835.10308@deprem.gov.tr>
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:
> hello
>
> I have spatial data which contain number of landslide presence cells
> with respect to landslide predictors and number of landslide absence
> cells with respect to same predictors.
>
> predictors are essentially categorical data.
>
> I tried logistic regression. But because of providing interaction
> capability of predictors, I want to use log-linear method.
> I hesitate the way I should use landslide count as response variable.
> only landslide presence data should be regarded ? or both landslide
> presence and absent data should be regarded as response variable ?
>
> I will appreciate if anyone can supply information
>
> thanks in advance
>
> Ahmet Temiz
> Gen Dir of Disaster of Affairs
> TURKEY