Note that omitting the intercept the estimate for the coefficient $\beta$ is given by <source>
$\hat{\beta} = \frac{\sum^n_{i=1} x_iy_i}{\sum^n_{i=1}x_i^2}$
in other words dividing the sum of $x$ times $y$ by the sum of squares of $x$ $-$ can also take the mean.
Due to the random ordering matching $x$ and $y$ values in the linear model you end up with a smaller numerator and as a result your coefficient will be smaller than 1.
I must admit that this is probably not an entirely satisfactory answer as it is based mainly on observed patterns; I don't know if there is a better mathematical explanation for it. You can check the former for yourself though.
Just take the sum of $xy$ as they were generated and compare it with the sum when you order both variables (dividing them by the sum of squares of $x$).
obs
. Indeed, withobs
set to1
, the expected value of either regression coefficient is infinite! $\endgroup$