data(Cica1)
boxplot(Au_Bonitur~Treatment, data=Cica1)
ZikaNP<-pairwiseCI(Au_Bonitur~Treatment, data=Cica1,
method="HL.ratio", control="Iso", alternative="two.sided", conf.level=0.9)
plotCI(ZikaNP, lines=c(0.5,1,2), lineslty=c(3,1,3), lineslwd=c(1,2,1),
ylim=c(0.5,2.5), main="Auchenorhyncha visual assessment" )
####################Run the code above in your browser using DataLab