x <- 0:1e6/1e6


png("link-sym.png", width=6, height=4, units="in", res=180)
plot(logit(x),x, type="l", las=1, ylab="Probability", xlab="x", xlim=c(-6,6), lwd=2, col=1)
lines(probit(x),x, lwd=1, col=6)
lines(cauchit(x),x, lwd=1, col=4)
legend("topleft", c("Cauchit","Probit","Logit"), lwd=c(1,1,2), col=c(4,6,1), bty="n" )
dev.off()


png("link-asym.png", width=6, height=4, units="in", res=180)
plot(logit(x),x, type="l", las=1, ylab="Probability", xlab="x", xlim=c(-6,6), lwd=2, col=1)
lines(cloglog(x),x, lwd=1, col=4)
lines(loglog(x),x, lwd=1, col=6)
legend("topleft", c("Complementary Log-Log","Log-Log","Logit"), lwd=c(1,1,2), col=c(4,6,1), bty="n" )
dev.off()