fjbjor<-c(5,2,9,7,3,3,4,5,8,3,5,5,6,7,1,4)
amagn<-c(0.1,0.03,0.19,0.095,0.07,0.02,0.07,0.085,0.12,0.04,0.06,0.05,0.1,0.09,0.01,0.05)
n<-length(fjbjor)

plot(fjbjor,amagn,pch=16,col=1,cex=0.6,axes=F,xlim=c(0,10),ylim=c(0,0.2),
xaxs='i',yaxs='i',xlab='',ylab='',cex.lab=1.1)
abline(lm(amagn~fjbjor),col=colors()[209],lwd=1)
for(i in 1:4){
y1<-amagn+runif(n,-0.05,0.05)
x1<-fjbjor+runif(n,-0.7,0.7)
abline(lm(y1~x1),col=colors()[209],lwd=1)
i = i+1
}
axis(1,cex.axis=1.2)
axis(2,cex.axis=1.2)
mtext('Number of beers',1,cex=1.3,2.5)
mtext('Alcohol level',2,cex=1.3,2.5)