Tilapia<-5
Mystus_singaringan<-c(1,2,2,2)
Osteochilus_vittatus<-c(26,30,11,2,28,6,19,12,4,4,18,13,29,21,9,8,21,19,28)
Henicorhynchus_siamensis<-c(46,22,33,28,50,42,31,4,13,8,11,8,18)
Barbodes_gonionotus<-c(8,11,66)
Anabas_testudineus<-c(12,26)
Trichogaster_microlepis<-c(9,5)
Trichogaster_pectoralis<-c(21,12,6,3)
counts<-c(Tilapia,Mystus_singaringan,Osteochilus_vittatus,Henicorhynchus_siamensis, Barbodes_gonionotus,Anabas_testudineus,Trichogaster_microlepis,Trichogaster_pectoralis)
sample_sizes<-c(130.0,28.95,27.2,17.92,26.12,81.65,78.31,59.46,62.84,58.28,45.03,35.02,46.46,43.54,24.64,31.52,57.05,58.87,59.7,39.02,28.64,65.03,43.28,36.2,64.57,67.43,70.42,89.3,75.44,51.38,78.92,30.69,21.69,21.84,23.45,18.21,49.42,110,105.63,56.25,51.55,91.42, 32.78,39.83,45.76,62.03,34.32,14.77) 
species<-c("Tilapia",rep("Mystus singaringan",4),rep("Osteochilus vittatus",19),rep("Henicorhynchus siamensis",13),rep("Barbodes gonionotus",3),rep("Anabas testudineus",2),rep("Trichogaster microlepis",2), rep("Trichogaster pectoralis",4))
counts_per_gram<-counts/sample_sizes
dataf<-cbind(as.factor(species),counts_per_gram)
boxplot(counts_per_gram ~ factor(species),ylab="Metacercaria per gram")
###
boxplot(counts_per_gram ~ factor(species),ylab="Metacercaria per gram",las=2)
###
par(mar=c(13,4,4,2))
boxplot(counts_per_gram ~ factor(species),ylab="Metacercaria per gram",las=2,cex=0.8)
###


# sort(tapply(counts/sample_sizes,species,mean))
#                 Tilapia       Mystus singaringan  Trichogaster microlepis  Trichogaster pectoralis 
#              0.03846154               0.07406214               0.20004559               0.25757761 
#      Anabas testudineus     Osteochilus vittatus      Barbodes gonionotus Henicorhynchus siamensis 
#              0.25859268               0.33030379               0.45006590               0.46634414 
#
Fish<-rownames(as.matrix(sort(tapply(counts/sample_sizes,species,mean))))
boxplot(counts_per_gram ~ factor(species,Fish),ylab="Metacercaria per gram",las=2,cex=0.8)
###
boxplot(counts_per_gram ~ factor(species,Fish),notch=TRUE, ylab="Metacercaria per gram", las=2, cex=0.8)


plot(TukeyHSD(aov(counts_per_gram ~ factor(species))))
levels(as.factor(species))
LETTERS[1:length(levels(as.factor(species)))]
x<-factor(species,labels=LETTERS[1:length(levels(as.factor(species)))])
plot(TukeyHSD(aov(counts_per_gram ~ factor(species))))
levels(as.factor(species))

LETTERS[1:length(levels(as.factor(species)))]
species<-factor(species,labels=LETTERS[1:length(levels(as.factor(species)))])
x<-factor(species,labels=LETTERS[1:length(levels(as.factor(species)))])
plot(TukeyHSD(aov(counts_per_gram ~ factor(x))))
TukeyHSD(aov(counts_per_gram ~ factor(x)))[1]