I have the following set up:
$y_i = \beta_0 + \beta_1 x_i + \beta_2 z_i + e_i$, where $e_i$ is extracted from a Normal (0,1) distribution independently of $x$ and $z$. The true values are $\beta_0 = \beta_1 = \beta_2 = 1$.
I want to compare the coefficients of regressing $y$ on $x$ and $z$ with the coefficients of regressing $y$ on the residuals of $x$ (obtained by regressing $x$ on $z$) and the residuals of $z$ (obtained by regressing $z$ on $x$).
I can see that the estimates are now biased but I can't find a way to understand the rational as to why they are biased. Any ideas?