rm(list=ls())

set.seed(4251)

pdf("linear_class.pdf")
plot(c(0,1),c(0,1),type="n",xlab="",ylab="",axes=FALSE)   ## xaxt="n",yaxt="n")   ##
abline(.9,-.8,lwd=4,col="purple")
x1 = runif(12,.2,.9)
y1 = .9 - .8*x1
y1 = y1 +.6*rnorm(12,.9-.8*x1,.8)
## y1 = y1 +rnorm(12,.9-.8*x1,.8)
points(x1,y1,pch="") 
text(x1,y1,label=0,col="blue",cex=2)

x2 = runif(15,.2,.9)
y2 = .9 - .8*x2
y2 = y2 -.3*rnorm(15,.9-.8*x1,.8)

points(x2,y2,pch="")
text(x2,y2,label=1,col='red',cex=2)

dev.off()