rm(list=ls())
c <- Miln_et_al_Dryad_data
str(c)
c$tank <- as.factor(c$tank)
c$day <- as.factor(c$day)
library(lmPerm)
d1 <- subset(c, day == "d1")
d4 <- subset(c, day == "d4")
d7 <- subset(c, day == "d7")

M1 <- aovp(maxrealised ~ cstart+day+tank, data = c, perm = "Exact")
summary(M1)

M12 <- aovp(asprint ~ cstart, data = c, perm = "Exact")
summary(M12)

M3 <- aovp(compsuccess ~ cstart+maxrealised+socialscore+day+tank, data = c, perm = "Exact")
summary(M3)


M4 <- aovp(socialscore ~ maxrealised+day+tank, data = c, perm = "Exact")
summary(M4)

M41 <- aovp(socialscore ~ maxrealised, data = d1, perm = "Exact")
summary(M41)

M42 <- aovp(socialscore ~ maxrealised, data = d4, perm = "Exact")
summary(M42)

M43 <- aovp(socialscore ~ maxrealised, data = d7, perm = "Exact")
summary(M43)