islands<-c("Yang","Butong","Adang","Rawi","Tarutao")
areas<- c(1.0,3.98,29.8,28.2,151)
rodents<-c(1,1,2,3,4)
par(mfrow=c(2,2)) # to plot two square graphs side by side
plot(areas, rodents, xlab="Island area (square km)", ylab="Number of murid species", pch=15)
plot(log(areas,10),log(rodents,10), xlab="Log (island area (square km))", ylab="Log(Number of murid species)", pch=15)
modelt<-lm(log(rodents,10)~log(areas,10))
abline(modelt,col="blue")
summary(modelt)