[R] logit regression, test among groups
Martin Henry H. Stevens
HStevens at muohio.edu
Tue May 28 12:59:23 CEST 2002
Dear all:
My logistic regression model includes one qualitative and one quantitative
predictor variable,
aes <- glm(p.a ~ spp * log(light), family=binomial(link=logit)),
where spp is abundance of 3 species and light is subcanopy light
availability varying from 0 1.
I want to test differences among levels of the quantitative variable at a
value of x other than the current log(light)=0. I tried using
I(log(light+.5)) but this caused serious lack of fit (AIC went from 1574 to
1656, and diagnostic plots demonstrated poor fit). How would create a model
statement such that I test differences among groups at log(light)=approx. 2?
Many thanks.
Martin Henry H. Stevens, Assistant Professor
338 Pearson Hall
Botany Department
Miami University
Oxford, OH 45056
Office: (513) 529-4206
Lab: (513) 529-4262
FAX: (513) 529-4243
http://www.muohio.edu/~botcwis/bot/henry.html
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !) To: r-help-request at stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._
More information about the R-help
mailing list