tj<-read.table(file="Tunjai_regeneration.txt", header=TRUE)
colnames(tj)[9]<- "PercentEstab"
zeros<-which(is.na(tj$PercentEstab))
tj$PercentEstab[zeros]<-0
min(tj$PercentEstab>0)
min_non_zero<-min(tj$PercentEstab[which(tj$PercentEstab>0)])
tj$PercentEstab[zeros]<-0.5*min_non_zero
tj$Size_categ<-ordered(tj$Size_categ, levels=c("S" ,"I", "L"))
tj$Coat<-ordered(tj$Coat, levels=c("Tn" ,"M", "Tk"))





par(mfrow=c(2,2))
tj$logit_establishment<-with(tj,log((PercentEstab/100)/(1-(PercentEstab/100))))
modelD<-aov(tj$logit_establishment ~ tj$Size_categ)
plot(modelD)


mof<-aov((tj$PercentEstab ~ tj$Size_categ))
plot(mof)

####