library(postscriptfonts)
ps.options(setfont=ps.setfont.latin1)
ps.options( reset=T, colors = ps.colors.rgb[c("black", "green", "blue"),] )
options(mixed.text.screen=T)

postscript("C:/mikefolder/mikeFolder/survivalmethods/programsforweb/Figures_revisedpaper/GGfig1.ps", horizontal=T, color.p=T)
par(mfrow=c(1,1), tck=-0.03, mgp=c(1.5, 0.5, 0), mar=c(4,3,1,5), cex=1.5 )

ttx<-seq(0,3,.01)
plot(ttx, ttx, xlab="",ylab="",xlim=c(0,3),ylim=c(-1,3),axes=F,type="n",lwd=3,cex=1.0)
lines(c(0,0),c(-1.05,3),lwd=1, col=1)
#gamma
lines(ttx,ttx,lwd=4,col=1)
#ammag
lines(ttx,1/ttx,lwd=4,col=1)
#lognormal
lines(c(0,3),c(0,0),lty=2,lwd=1,col=16)
#weibull
lines(c(0,3),c(1,1),lty=2,lwd=1,col=16)
#inverse weibull
lines(c(0,3),c(-1,-1),lty=2,lwd=1,col=16)
#inverse gamma
lines(c(0,1.05),c(0,-1.05),lty=2,lwd=1,col=16)
#inverse ammag
lines(ttx[ttx>0.95],-1/ttx[ttx>0.95],lty=2,lwd=1,col=16)
axis(1,at=seq(0,3,1),label=c("0","1","2","3"),cex=1.0,tck=-0.01,lwd=1, col=1)


par(col=16, cex=0.8)
mixed.text(x=3, y=0.07, "log normal", adj=1 )
mixed.text(x=3, y=1.07, "Weibull", adj=1 )
mixed.text(x=3, y=-0.93, "Inverse Weibull", adj=1 )
mixed.text.vector(0.35, -0.27, c(1, 13, 1, 13, 1),c("Inverse Gamma (", "\l", " = - ", "\s", ")"), slope= -1)
mixed.text.vector(2.32, -0.33, c(1, 13, 1, 13, 1),c("Inverse Ammag (", "\l", " = - 1/ ", "\s", ")"), slope= 0.10)

par(col=1, cex=2)
mixed.text.vector(1.2, -1.5, c(1, 13, 1), c("Scale (", "\s", ")"))
mixed.text.vector(3.35, 0.5, c(1, 13, 1), c("Shape (", "\l", ")"), srt=90)
par(col=1, cex=2.4)
mixed.text(0.8, 0.2, "(0, 0)", adj=0 )
mixed.text.vector(0.8, 1.9, c(1, 13, 1, 13, 1), c("(", "\245", ", ", "\245", ")")  )
mixed.text.vector(0.2, 1.2, c(1, 13, 1), c("(0, ", "\245", ")")  )
mixed.text.vector(1.5, 1.2, c(1, 13, 1), c("(", "\245", ", 0)")  )

arrows(1.025, 0.6, 1.125, 0.4, size=0.2, open=T, lwd=6)
arrows(0.875, 0.4, 0.975, 0.6, size=0.2, open=T, lwd=6)
arrows(0.9, 1.7, 1.0, 1.5, size=0.2, open=T, lwd=6)
arrows(1.05, 1.5, 1.15, 1.7, size=0.2, open=T, lwd=6)
arrows(1.7, 1.55, 1.8, 1.35, size=0.2, open=T, lwd=6)
arrows(0.3, 1.35, 0.4, 1.55, size=0.2, open=T, lwd=6)

par(col=1, cex=0.9)
mixed.text(0.04, 0.27, "(0, d)")
mixed.text.vector(0.17, 2.7, c(1, 13, 1), c("(c, ", "\245", ")"))
arrows(0.09, 0.34, 0.15, 0.42, size=0.08, open=T, lwd=2)
arrows(0.22, 2.76, 0.28, 2.84, size=0.08, open=T, lwd=2)
mixed.text(2.85, 0.45, "(c, 0)")
mixed.text.vector(2.85, 2.7, c(1, 13, 1), c("(", "\245", ", d)"))
arrows(2.92, 2.84, 2.98, 2.76, size=0.08, open=T, lwd=2)
arrows(2.90, 0.62, 2.96, 0.54, size=0.08, open=T, lwd=2)

par(col=1, cex=1.2)
mixed.text.vector(2, 2.15, c(1, 13, 1, 13, 1), c("Gamma (", "\l", " = ", "\s", ")"), slope=1)
mixed.text.vector(2, 0.65, c(1, 13, 1, 13, 1), c("Ammag (", "\l", " = 1/", "\s", ")"), slope=-0.2)

par(mgp=c(1.5, 0.9, 0), col=1)
axis(4,at=seq(-1,3,1),label=c("-1","0","1","2","3"),cex=1.0,tck=-0.01, lwd=1)

dev.off()