datafilename="nitrate.csv"
data.ex2=read.table(datafilename,header=T,sep=";")   #read the data into a table
data.ex2                                      #show the data
aov.ex2 = aov(Nitrate~pH+temperature+Turbulence,data=data.ex2)         #do the analysis of variance
summary(aov.ex2)                                    #show the summary table
print(model.tables(aov.ex2,"means"),digits=3)       #report the means and the number of subjects/cell
boxplot(Nitrate~pH+temperature+Turbulence,data=data.ex2)

boxplot(Nitrate~pH,data=data.ex2)
boxplot(Nitrate~temperature,data=data.ex2)
boxplot(Nitrate~pH,data=data.ex2)


==============================================================================

datafilename="nitrate.csv"
data.ex2=read.table(datafilename,header=T,sep=";")   #read the data into a table
data.ex2                                      #show the data
aov.ex2 = aov(Nitrate~pH+temperature+Turbulence,data=data.ex2)         #do the analysis of variance
summary(aov.ex2)                                    #show the summary table
print(model.tables(aov.ex2,"means"),digits=3)       #report the means and the number of subjects/cell
boxplot(Nitrate~temperature*Turbulence,data=data.ex2)

boxplot(Nitrate~temperature,data=data.ex2,ylab='acc. mass(mg)',xlab='Temperature (C)')

boxplot(Nitrate~Turbulence,data=data.ex2,ylab='acc. mass(mg)',xlab='Turbulence (rpm)')
boxplot(Nitrate~pH,data=data.ex2,ylab='acc. mass(mg)',xlab='pH')


===============================================================================

datafilename="phosphate.csv"
data.ex2=read.table(datafilename,header=T,sep=";")   #read the data into a table
data.ex2                                      #show the data
aov.ex2 = aov(Phosphate~pH+temperature+Turbulence,data=data.ex2)         #do the analysis of variance
summary(aov.ex2)                                    #show the summary table
print(model.tables(aov.ex2,"means"),digits=3)       #report the means and the number of subjects/cell

boxplot(Phosphate~temperature,data=data.ex2,ylab='acc. mass(mg)',xlab='Temperature (C)')

boxplot(Phosphate~Turbulence,data=data.ex2,ylab='acc. mass(mg)',xlab='Turbulence (rpm)')
boxplot(Phosphate~pH,data=data.ex2,ylab='acc. mass(mg)',xlab='pH')
boxplot(Phosphate~pH+Turbulence,data=data.ex2,ylab='acc. mass(mg)',xlab='pH')