nests<- c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0,1,1,2,2,2,2,2,4,5,5,7,7,9,11,15,20,24,30) 
total_damage_percent<-c(3.72, 3.11, 2.15, 1.40, 1.12, 0.94, 0.55, 0.59, 0.73, 0.83, 0.19, 0.44, 0.19, 0.64, 0.66, 0.68, 1.03, 0.19, 0.20, 0.13, 0.10, 0.59, 0.09, 0.05, 0.06, 0.06, 0.04, 0.14)
leafbeetle_damage_percent<-c(0.44, 0.49, 0.59, 0.70, 0.73, 0.94, 1.03, 1.09, 1.42, 2.10, 0.19, 0.39, 0.03, 0.04, 0.19, 0.36, 0.42, 0.20, 0.21, 0.09, 0.05, 0.16, 0.07, 0.05, 0.06, 0.06, 0.01, 0.14)
plot(nests, total_damage_percent, xlab= "Number of weaver ant nests", ylab= "Percentage leaf damage", pch=1)
points(nests, leafbeetle_damage_percent, pch=4)

modelA<-lm(leafbeetle_damage_percent ~ nests)
plot(modelA)  # hit return until you get the Q-Q plot

modelB<-lm(leafbeetle_damage_percent ~ log(nests + 1)) 
plot(modelB) 

plot(log(nests+1), total_damage_percent^0.333, xlab= "ln(weaver ant nests+1)", ylab= "Percentage leaf damage ^0.33", pch=1, col="red")
points(log(nests+1), leafbeetle_damage_percent^0.333, pch=4,col="green")

modelC<-lm(total_damage_percent^(1/3) ~ log(nests + 1))
summary(modelC)

modelT<-lm(leafbeetle_damage_percent^(1/3) ~ log(nests + 1))
summary(modelT)

abline(lm(total_damage_percent^(1/3) ~ log(nests + 1)),lwd=2, col= "red")
abline(modelT,lwd=2, col= "green")