attach(data)
plot(x, y, pch=21, col="blue", bg="orange")
var(x)
var(y)
var(x,y)
var(x,y)/sqrt(var(x)*var(y))
cor(x,y)