### SCRIPT PRODUCTION CURVES 

PHPB=read.csv("WinterProductionCurves.csv", header=TRUE, sep=";", quote="\"")
View(PHPB)

### Pmax ###
kruskal.test(PmaxP[Espece=="LC"]~Tratamiento[Espece=="LC"], data=PHPB)
kruskal.test(PmaxC[Espece=="LC"]~Tratamiento[Espece=="LC"], data=PHPB)

kruskal.test(PmaxP[Espece=="PC"]~Tratamiento[Espece=="PC"], data=PHPB)
kruskal.test(PmaxC[Espece=="PC"]~Tratamiento[Espece=="PC"], data=PHPB)

kruskal.test(PmaxP[Espece=="LI"]~Tratamiento[Espece=="LI"], data=PHPB)
kruskal.test(PmaxC[Espece=="LI"]~Tratamiento[Espece=="LI"], data=PHPB)

### Ek ###
kruskal.test(EkP[Espece=="LC"]~Tratamiento[Espece=="LC"], data=PHPB)
kruskal.test(EkC[Espece=="LC"]~Tratamiento[Espece=="LC"], data=PHPB)

kruskal.test(EkP[Espece=="PC"]~Tratamiento[Espece=="PC"], data=PHPB)
kruskal.test(EkC[Espece=="PC"]~Tratamiento[Espece=="PC"], data=PHPB)

kruskal.test(EkP[Espece=="LI"]~Tratamiento[Espece=="LI"], data=PHPB)
kruskal.test(EkC[Espece=="LI"]~Tratamiento[Espece=="LI"], data=PHPB)

### TEST HOMOGENEITY OF VARIANCES ###

fligner.test(PmaxC[Espece=="PC"] ~ Tratamiento[Espece=="PC"], data=PHPB)
fligner.test(PmaxP[Espece=="PC"] ~ Tratamiento[Espece=="PC"], data=PHPB)
fligner.test(EkC[Espece=="PC"] ~ Tratamiento[Espece=="PC"], data=PHPB)
fligner.test(EkP[Espece=="PC"] ~ Tratamiento[Espece=="PC"], data=PHPB)

fligner.test(PmaxC[Espece=="LI"] ~ Tratamiento[Espece=="LI"], data=PHPB)
fligner.test(PmaxP[Espece=="LI"] ~ Tratamiento[Espece=="LI"], data=PHPB)
fligner.test(EkP[Espece=="LI"] ~ Tratamiento[Espece=="LI"], data=PHPB)
fligner.test(EkC[Espece=="LI"] ~ Tratamiento[Espece=="LI"], data=PHPB)

PCT2=read.csv("SummerProductionCurves.csv", header=TRUE, sep=";", quote="\"")
View(PCT2)

### Photosynthesis ###

kruskal.test(Pmax[X=="LC"]~X.1[X=="LC"], data=PCT2)
kruskal.test(Pmax[X=="PC"]~X.1[X=="PC"], data=PCT2)
kruskal.test(Pmax[X=="LI"]~X.1[X=="LI"], data=PCT2)

kruskal.test(Ek[X=="LC"]~X.1[X=="LC"], data=PCT2)
kruskal.test(Ek[X=="PC"]~X.1[X=="PC"], data=PCT2)
kruskal.test(Ek[X=="LI"]~X.1[X=="LI"], data=PCT2)

### Calcification ###

kruskal.test(PmaxC[X=="LC"]~X.1[X=="LC"], data=PCT2)
kruskal.test(PmaxC[X=="PC"]~X.1[X=="PC"], data=PCT2)
kruskal.test(PmaxC[X=="LI"]~X.1[X=="LI"], data=PCT2)

kruskal.test(EkC[X=="LC"]~X.1[X=="LC"], data=PCT2)
kruskal.test(EkC[X=="PC"]~X.1[X=="PC"], data=PCT2)
kruskal.test(EkC[X=="LI"]~X.1[X=="LI"], data=PCT2)