#File syncrony2.xls
#load lme4, blme, multcomp


x1<-read.table(file="clipboard",header=T)
attach(Synchrony2R)


fit1<-lmer(ZCorrR~Treat*Time+(1|Pair))

fit5<-blmer(ZCorrRankR~Treat*Time+(1|Pair))
summary(fit2)
plot(fit2)
fit6<-blmer(ZCorrRankR~Treat+Time+(1|Pair))

anova(fit5,fit6)

fitnull<-blmer(ZCorr~1+(1|Pair))

anova(fitnull,fit1)

AB<-interaction(Treat,Time)
fit4<-blmer(ZCorr~AB+(1|Pair))

summary(glht(fit4, linfct=mcp(AB = "Tukey")))