# mimicry Banci #PAT=pattern(brown/coral/diad/intermediate), LOC=locality (Juquitiba/Sta Virginia), SEA=season (dry/rainy), SIZ=size(small/big) mim<-read.table("mim.csv",header=T,sep=";") head(mim) summary(mim) mim$PAT<-factor(mim$PAT) # analysis a,c = factors without interaction # analysis b,d = interaction between PAT paired to other factors #BIR=bird a<-glm(BIR~PAT+LOC+SEA+SIZ, data = mim, family = "binomial") b<-glm(BIR~(PAT+LOC+SEA+SIZ)*PAT-(PAT+LOC+SEA+SIZ),data = mim, family = "binomial") summary(a) summary(b) #MAM=mammal c<-glm(MAM~PAT+LOC+SEA+SIZ, data = mim, family = "binomial") d<-glm(MAM~(PAT+LOC+SEA+SIZ)*PAT-(PAT+LOC+SEA+SIZ),data = mim, family = "binomial") summary(c) summary(d) with(mim, table(BIR,PAT)) with(mim, table(BIR,LOC)) with(mim, table(BIR,SEA)) with(mim, table(BIR,SIZ)) with(mim, table(BIR,PAT,LOC)) with(mim, table(BIR,PAT,SEA)) with(mim, table(BIR,PAT,SIZ)) with(mim, table(MAM,PAT)) with(mim, table(MAM,LOC)) with(mim, table(MAM,SEA)) with(mim, table(MAM,SIZ)) with(mim, table(MAM,PAT,LOC)) with(mim, table(MAM,PAT,SEA)) with(mim, table(MAM,PAT,SIZ))