library(MASS)
data(snails)
spA<-snails[1:48,]; spB<-snails[48:96,]
par(mfrow=c(2,2))
boxplot(spA$Deaths~spA$Temp,main="Species A",xlab="Temperature")
boxplot(spA$Deaths~spA$Rel.Hum,main="Species A",xlab="Relative Humidity")
boxplot(spB$Deaths~spB$Temp,main="Species B",xlab="Temperature")
boxplot(spB$Deaths~spB$Rel.Hum,main="Species B",xlab="Relative Humidity")
snails$deadORalive<-cbind(snails$Deaths,20-snails$Deaths) 
# head(snails$deadORalive)
m1<-glm(deadORalive ~ Species/(Temp+Rel.Hum), family=binomial, data=snails)
m2<-glm(deadORalive ~ Species/(Temp+Rel.Hum), family=quasibinomial, data=snails)
summary(m1)
summary(m2)