var} ====================================================  Consider a measurement model $y = X\beta^\star + \epsilon$, where $X\in\mathbb{R}^{n\times p}$ is the [*design*]{} matrix and $\epsilon\in\mathbb{R}^n$ is 