CD = read.table('../../Data/Cracker.dat', header=TRUE, skip=0)
CD$Diet = as.factor(CD$Diet)
##print(head(CD))

T = table(CD$Diet, CD$Bloat)
print(T)

## Run a chi-squared test:
##
res = chisq.test(T)
print(res)

## Run an exact test:
##
res = fisher.test(T) #, alternative='greater')
print(res)


CD$CombDiet ='Diets_1_N_4'
CD[ (CD$Diet==2) | (CD$Diet==3), 'CombDiet'] = 'Diets_2_N_3'

CD$CombBloat ='Bloat_H_N_M'
CD[(CD$Bloat=='low') | (CD$Bloat=='none'), 'CombBloat'] = 'Bloat_L_N_N'

T = table(CD$CombDiet, CD$CombBloat)
print(T)

## Run a chi-squared test:
##
res = chisq.test(T)
print(res)

boxplot(Digested ~ Diet, data=CD, xlab='Diet', ylab='Digested')

model = aov(Digested ~ Diet, data=CD)
print(summary(model) )