# ANOVA script by Elie Maza and Christian Chervin, AgroToulouse
# Loading packages

library(agricolae)

# Loading data

HedoScore=read.table("DataForANOVA.csv",sep=";",header=TRUE)

# Descriptive analysis

summary(HedoScore)

HedoScore$Product=as.factor(HedoScore$Product)
HedoScore$NumPanel=as.factor(HedoScore$NumPanel)
summary(HedoScore,maxsum=Inf)

# ANOVA

res=aov(HedoScore~Product+NumPanel,data=HedoScore)
res

anova(res)

# Analysis of residuals

par(mfrow=c(2,2))
plot(res)
par(mfrow=c(1,1))

# Multiple comparison with Tuckey HSD test

groups_product=HSD.test(res,"Product")
groups_product

pdf(paste("fig_ANOVA.pdf",sep=""),width=10,height=4)
par(mar=c(3.1,4.1,1.1,0.1))
plot(groups_product,variation="SE",ylab="Hedonic score",main="",cex=1.25,cex.names=1.25,cex.axis=1.25,cex.lab=1.25,ylim=c(3.5,7))
dev.off()

#---------
# The End
#---------