Displaying 1 result from an estimated 1 matches for "herbivoryedit".
2012 Mar 03
1
interpreting the output of a glm with an ordered categorical predictor.
...d a binomial response (presence/absence of herbivore damage).
My predictor has four categories: high, medium, low, and none. I used the
"ordered" function to sort my categories for a glm.
ah <-
read.csv("http://depot.northwestern.edu/class/2012WI_PBC_435-0_AND_BIOL_SCI_313/muller/herbivoryEdit.csv")
ah1<- ah[ah$date=="110810",]
ah2<-ah[ah$date=="110904",]
aphidOrder <- ordered(ah2$aphidLevelMax,levels=c("none", "low", "med",
"high"))
ordAph <- glm(chewholebinom~aphidOrder,family=binomial,data=ah2)
When I ra...