sort

n <- 100
m <- 1000
coef <- rep(0,m)
coef2 <- rep(0,m)
beta <- 0

for (j in seq(1,m) ){
  x <- rnorm(n, 0, 1)
  y <- rnorm(n, x* beta, 1)
  
  fit <- lm(y ~ x)
  coef[j] <- fit$coefficients[1]
  fit <- lm(sort(y) ~ sort(x))
  coef2[j] <- fit$coefficients[1]
}