Setup: Consider a random sample of size n with binary outcome $Y_i\in\{0,1\}$. Assume $Y_i\sim Bern(\pi_i)$. Use a linear probability model so that $\pi_i=X_i^\intercal\beta$, where $X_i$ is a predictor vector of length P. Here consider $\beta$ by maximum likelihood.
Question: Show that the error term from the linear regression of $Y_i$ on $X_i$ is always heteroskedastic. If you were to use OLS for this model, would this problem be corrected?
Comment: The error term is heteroskedastic if it depends on $X_i$, but I'm not sure what approach I should use to get a formula for it.