Mos<-read.csv(".//MS-Data//MOSQUITO.CSV")

xtabs(Ratio~Insecticide+Location,data=Mos)


summary(aov(Ratio~Insecticide,data=Mos))
summary(aov(Ratio~Location,data=Mos))

summary(aov(Ratio~Insecticide+Location,data=Mos))


stripchart(Ratio~Insecticide,data=Mos,vertical=T,pch=19,xlim=range(0,6))

stripchart(Ratio~Location,data=Mos,vertical=T,pch=19,xlim=range(0,8))


stripchart(Ratio~Insecticide,data=Mos,vertical=T,pch=19,xlim=range(0.5,5.5),
subset=(Location=="Anguilla"),col="red",ylim=range(0,10))

stripchart(Ratio~Insecticide,data=Mos,vertical=T,pch=15,
subset=(Location=="Antigua"),col="blue",add=T)

stripchart(Ratio~Insecticide,data=Mos,vertical=T,pch=17,
subset=(Location=="Dominica"),col="green",add=T)

stripchart(Ratio~Insecticide,data=Mos,vertical=T,pch=18,
subset=(Location=="Guyana"),col="cyan",add=T)

stripchart(Ratio~Insecticide,data=Mos,vertical=T,pch=8,
subset=(Location=="Jamaica"),col="purple",add=T)


#####################################################################
stripchart(log(Ratio)~Insecticide,data=Mos,vertical=T,pch=19,xlim=range(0,6))

stripchart(log(Ratio)~Location,data=Mos,vertical=T,pch=19,xlim=range(0,8))


summary(aov(log(Ratio)~Insecticide,data=Mos))
summary(aov(log(Ratio)~Location,data=Mos))

summary(aov(log(Ratio)~Insecticide+Location,data=Mos))