rm(list=ls())

set.seed(364)

n = 10
k = 5
s = .2

X1 = runif(n,1,2)
X1 = sort(X1)
Y1 = 0 - X1 + rnorm(n,0,s)

X2 = runif(n,1.5,2.5)
X2 = sort(X2)
Y2 = 1 - X2 + rnorm(n,0,s)

X3 = runif(n,2,3)
X3 = sort(X3)
Y3 = 2 - X3 + rnorm(n,0,s)

X4 = runif(n,2.5,3.5)
X4 = sort(X4)
Y4 = 3 - X4 + rnorm(n,0,s)

X5 = runif(n,3,4)
X5 = sort(X5)
Y5 = 4 - X5 + rnorm(n,0,s)


X = c(X1,X2,X3,X4,X5)
Y = c(Y1,Y2,Y3,Y4,Y5)

pdf("Simpson0.pdf")
plot(X,Y,pch=21,cex=2, bg="lightpink",xlab="",ylab="",ylim=c(-2.2,1),xlim=c(.5,4.5))
dev.off()


pdf("Simpson1.pdf")

plot(X,Y,pch=21,cex=2,bg="lightpink",xlab="",ylab="",ylim=c(-2.2,1),xlim=c(.5,4.5))
out = lm(Y ~ X)
abline(out,col="purple",lwd=4)
dev.off()


pdf("Simpson2.pdf")
plot(X,Y,pch=21,cex=2,bg="lightpink",xlab="",ylab="",ylim=c(-2.2,1),xlim=c(.5,4.5))
                   # , xaxt="n",yaxt="n")
lines(X1,0-X1, col="purple",lwd=4)
lines(X2,1-X2, col="purple",lwd=4)
lines(X3,2-X3, col="purple",lwd=4)
lines(X4,3-X4, col="purple",lwd=4)
lines(X5,4-X5, col="purple",lwd=4)



pdf("Simpson3.pdf")
grid = seq(1,4,length=100)
plot(grid,2-grid,type="l",lwd=3,xlab="",ylab="",ylim=c(-2.2,1),xlim=c(.5,4.5),col="purple")
points(X,Y,pch=21,cex=2,bg="lightpink")



dev.off()