function(PDF=F){
if(PDF==T) pdf(file="shiftmodel.pdf",width=11)
par(mfrow=c(2,1),mar=c(5,4,0.2,1)+.1)
x=seq(-3,4,.01)
y1=pnorm(x)
y2=pnorm(x,1,1)
plot(x,y1,xlab="x",ylab="CDF",type="l",axes=F)
axis(1)
axis(2)
abline(h=c(0,1),col="grey")
lines(x,y2,lty=2)
legend(-3,.9,expression(F(x)==Phi(x)~"standard normal",
G(x)==Phi(x-1)),lty=c(1,2),bty="n",y.intersp=1.2)

x=seq(0,5,.01)
y1=punif(x,1,3)
y2=punif(x,2,4)
plot(x,y1,xlab="x",ylab="CDF",type="l",axes=F)
legend(0,.9,expression(list(F(x)==(x-1)/2, x %in%~"(1,3)")~"  Uniform(1,3)",
list(G(x)==F(x-1))~"  Uniform(2,4)"),bty="n",lty=c(1,2),y.intersp=1.2)
axis(1)
axis(2)
abline(h=c(0,1),col="grey")
lines(x,y1)
lines(x,y2,lty=2)
par(mfrow=c(1,1))
if(PDF==T) dev.off()
}