Displaying 2 results from an estimated 2 matches for "leafmin".
Did you mean:
leadmon
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)
plot(LeafDig,ylim=c(0,20),main="Leafminers",ylab="Observed
numbers",xlab="",adj=0,bty="l")
points(Diglyphus,pch=2)
lines(Diglyphus,lty=2)
legend(1,18,c("...
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 ~...