The linear model is
$$y_{i}= x_{i}'\beta+u_{i}$$
When written in vector notation such that $y_{i}$ is a $1$ x $1$ matrix of outcomes, $x_{i}'$ is a $1$ x $k$ matrix of control variables, $\beta$ is a $k$ x $1$ matrix of population parameters for the coefficients of the control variables, and $u_{i}$ is a $1$ x $1$ matrix for the residual error term.
$F_{R}$ is the robust F statistic for the linear model.
$$F_{R}= (R\hat\beta_{ols}- c)'\Big[R[\sum{x_{i}x_{i}'}]^{-1}[\sum{\hat u^2_{i}x_{i}x_{i}'}][\sum{x_{i}x_{i}'}]^{-1}R'\Big]^{-1}(R\hat\beta_{ols}- c)/J$$
where $R$ is a $j$ x $k$ matrix of restrictions, $c$ is a $j$ x $1$ matrix of null hypothesis, $\hat\beta_{ols}$ is a $k$ x $1$ matrix of estimates of coefficients of $x_{i}$, $\hat u_{i}$ is a $1$ x $1$ matrix of the estimated residual error, $J$ is the degrees freedom, I think it is a scalar but I am not 100% certain.
How do you prove that $F_{R}$ is asymptotically approximately chi squared distributed with $J$ degrees of freedom.
$$ F_{R} \sim \chi^{2}(J) $$