Let there be $n$ assets and $k$ factors in the market. We assume multivariate normally distributed factor returns
with $k\times k$ factor covariance ...
This is not exactly an answer, but I hope to discuss a few things.
LinearFactorModel is definitely not Fama-Macbeth regression. The first stage seems the same, but the second stage only performs one ...