# Wright-y-Muller-Landau-2006 future of tropical forests

d<-read.csv("Wright-y-Muller-Landau-2006.csv",header=TRUE)
head(d)
par(mfrow=c(2,2))
plot(d[,1],d[,2],pch=as.numeric(as.factor(as.numeric(d[,3]))),col=as.numeric(as.factor(as.numeric(d[,3]))), xlab=colnames(d)[1], ylab=colnames(d)[2])

plot(log10(d[,1]),d[,2],pch=as.numeric(as.factor(as.numeric(d[,3])))+14,col=as.numeric(as.factor(as.numeric(d[,3])))+3,type="p",xlab="Log(rural popn density)", ylab=colnames(d)[2])


plot(d[,1],d[,2],log="x",pch=as.numeric(as.factor(as.numeric(d[,3])))+14,col=as.numeric(as.factor(as.numeric(d[,3])))+3,axes=FALSE,xlim=c(0.1,1000),ylim=c(0,100),xlab="Rural popn density", ylab="% potential forest remaining")
axis(1,c("0.1","1.0","10","100","1000"),at=c(10^-1,10^0,10^1,10^2,10^3),cex.axis=0.8)
axis(2,at=c(0:10)*10,cex.axis=0.8)
legend("bottomleft",pch=as.numeric(levels(as.factor(as.numeric(d[,3]))))+14,col=as.numeric(levels(as.factor(as.numeric(d[,3]))))+3,levels((d[,3])),bty="n")


install.packages("sfsmisc")
library(sfsmisc)
plot(d[,1],d[,2],log="x",pch=as.numeric(as.factor(as.numeric(d[,3])))+14,col=as.numeric(as.factor(as.numeric(d[,3])))+3,axes=FALSE,xlim=c(0.1,1000),ylim=c(0,100),xlab="Rural popn density", ylab="% potential forest remaining")
eaxis(1,at=c(10^-1,10^0,10^1,10^2,10^3),cex.axis=0.8)
axis(2,at=c(0:10)*10,cex.axis=0.8)
legend("bottomleft",pch=as.numeric(levels(as.factor(as.numeric(d[,3]))))+14,col=as.numeric(levels(as.factor(as.numeric(d[,3]))))+3,levels((d[,3])),bty="n")