[R] nontabular logistic regression

Jeffrey Stratford stratja at auburn.edu
Fri Oct 13 16:28:23 CEST 2006


Hi.  I'm attempting to fit a logistic/binomial model so I can determine
the influence of landscape on the probability that a box gets used by a
bird.  I've looked at a few sources (MASS text, Dalgaard, Fox and
google) and the examples are almost always based on tabular predictor
variables.  My data, however are not.  I'm not sure if that is the
source of the problems or not because the one example that includes a
continuous predictor looks to be coded exactly the same way.  Looking at
the output, I get estimates for each case when I should get a single
estimate for purbank.  Any suggestions?

Many thanks,

Jeff


THE DATA: (200 boxes total, used [0 if unoccupied, 1 occupied], the rest
are landscape variables).  

box	use	purbank	purban2	purban1	pgrassk	pgrass2	pgrass1	grassdist	grasspatchk
1	1	0.003813435	0.02684564	0.06896552	0.3282487	0.6845638	0.7586207	0	3.73
2	1	0.04429451	0.1610738	0.1724138	0.1534174	0.3825503	0.6551724	0	1.023261
3	1	0.04458785	0.06040268	0	0.1628043	0.557047	0.7586207	0	0.9605769
4	1	0.06072162	0.2080537	0.06896552	0.01936052	0	0	323.1099	0.2284615
5	0	0.6080962	0.6979866	0.6896552	0.03168084	0.1275168	0.2413793	30	0.2627027
6	1	0.6060428	0.6107383	0.3448276	0.04077442	0.2885906	0.4482759	30	0.2978571
7	1	0.3807568	0.4362416	0.6896552	0.06864183	0.03355705	0	94.86833	0.468
8	0	0.3649164	0.3154362	0.4137931	0.06277501	0.1275168	0	120	0.4585714

THE CODE:

box.use<- read.csv("c:\\eabl\\2004\\use_logistic2.csv", header=TRUE)
attach(box.use)
box.use <- na.omit(box.use)
use <- factor(use, levels=0:1)
levels(use) <- c("unused", "used")
glm1 <- glm(use ~ purbank, binomial)

THE OUTPUT:

Coefficients:
                     Estimate Std. Error   z value Pr(>|z|)
(Intercept)        -4.544e-16  1.414e+00 -3.21e-16    1.000
purbank0            2.157e+01  2.923e+04     0.001    0.999
purbank0.001173365  2.157e+01  2.067e+04     0.001    0.999
purbank0.001466706  2.157e+01  2.923e+04     0.001    0.999
purbank0.001760047  6.429e-16  2.000e+00  3.21e-16    1.000
purbank0.002346729  2.157e+01  2.923e+04     0.001    0.999
purbank0.003813435  2.157e+01  2.923e+04     0.001    0.999
purbank0.004106776  2.157e+01  2.067e+04     0.001    0.999
purbank0.004693458  2.157e+01  2.067e+04     0.001    0.999


****************************************
Jeffrey A. Stratford, Ph.D.
Postdoctoral Associate
331 Funchess Hall
Department of Biological Sciences
Auburn University
Auburn, AL 36849
334-329-9198
FAX 334-844-9234
http://www.auburn.edu/~stratja



More information about the R-help mailing list