popbio (version 1.0.1)

logi.hist.plot: Plot logistic regression

Description

Plot combined graphs for logistic regressions

Usage

logi.hist.plot(independ, depend, logi.mod = 1, type = "dit", 
boxp = TRUE, rug = FALSE, ylabel = "Probability", ylabel2 = "Frequency", 
xlabel = "", mainlabel = "", las.h = 1, counts = FALSE, ...)

Arguments

independ

explanatory variable

depend

dependent variable, typically a logical vector

logi.mod

type of fitting, 1 = logistic; 2 = "gaussian" logistic

type

type of representation, "dit" = dit plot; "hist" = histogram

boxp

TRUE = with box plots, FALSE = without

rug

TRUE = with rug plots, FALSE = without

ylabel

y-axis label

ylabel2

2nd y-axis label

xlabel

x-axix label

mainlabel

overall title for plot

las.h

orientation of axes labels (0 = vertical, 1 = horizontal

counts

add counts above histogram bars

...

additional options passed to logi.hist

Value

• A logistic regression plot

References

de la Cruz Rot, M. 2005. Improving the Presentation of Results of Logistic Regression with R. ESA Bulletin 86:41-48.

http://esapubs.org/bulletin/backissues/086-1/bulletinjan2005.htm

Examples

data(aq.trans)

aq.trans$survived<-aq.trans$fate!="dead"

a<-subset(aq.trans, leaf<50 & stage!="recruit", c(leaf,survived))

logi.hist.plot(a$leaf,  a$survived, 
type="hist", boxp=FALSE, counts=TRUE, las.h=0, int=10, 
ylabel="Survival probability", ylabel2="Number of plants", 
 xlab="Number of leaves" )



b<-glm(survived ~ leaf, binomial, data=a)
 summary(b)