x <- c(2,3,7,10) y <- c(7,1,2,30) R <- cor(x,y) R R^2