x <- rnorm(20)
y <- 2*x + rnorm(20)
plot(x,y)
rsquare(x,y)
r2 <- sapply(1:10000, function(i){
x <- rnorm(20); y <- 2*x + rnorm(20); rsquare(x,y) })
hist(r2, breaks=70, col=5,
main= "10'000 times x <- rnorm(20); y <- 2*x + rnorm(20); rsquare(x,y)")
Run the code above in your browser using DataLab