Displaying 2 results from an estimated 2 matches for "leafmines".
Did you mean:
leafminers
2003 Sep 02
1
Plotting dates
I'm trying to plot observations against observation dates and getting julian
dates along the x-axis:
library(date)
Week<-as.date(c("05/02/03","05/09/03","05/16/03","05/23/03","05/30/03","06/0
7/03","06/14/03"))
Leafminers<-c(0,2,5,10,4,6,5)
Diglyphus<-c(0,0,4,5,7,3,1)
LeafDig<-cbind(Week,Leafminers,Diglyphus)
2006 Apr 23
1
Comparing GLMMs and GLMs with quasi-binomial errors?
...from the MASS library(after using the 'cbind()' command on
the two columns with total number of leaves per seedling and number of
leaves attacked by that herbivory category) - and then compared these models
to GLMs without the random effect.
## model example1: leaf mines GLMM
proportion.leafmines <- cbind(leaves.affected, total.leaves -
leaves.affected)
leafminesGLMM <- glmmPQL(proportion.leafmines ~ PLATEAU * DISTANCE,
random=~1| TreeID, family=binomial(link=logit))
##AIC(leafminesGLMM) = 474.773
## model example2: leaf mines GLM
leafminesGLM <- glm(proportion.leafmines ~ PL...