在混合 Proc 中求解正规方程时，SAS 显然使用了一种称为 g2 逆的东西（Pringle 和 Rayner，1971）。 R 中的 ginv 给出了不同的 moore-penrose 逆。 R 中是否有任何包可以为我提供 SAS 的 g2 逆？对我来说非常重要的是，我得到了与 SAS 相同的广义逆，即 g2 逆，而不是 moore-penrose。基本上，当 (X'X) 矩阵不是满秩时，我需要在 R 中获得与 SAS 中相同的固定效应估计。 谢谢。 任何帮助或建议将不胜感激。
您可以使用 MATLAB 或 PROC IML 轻松获得 SAS 逆矩阵，遵循 Searle 在他的书《线性模型：我将尝试做一个总结》。 您可以通过获得系数矩阵的 rref 来实现线性相关的列，然后： a) 删除 SAS 将解向量中的效应置零的行和列。 b) 取a 中剩下的矩阵的正则逆。 c) 将b得到的逆元素放在原来的位置 d) 填写在执行步骤 a 时清零的行和列。
如果你写下增广矩阵，你可以使用PROC IML的SWEEP函数更容易地得到G2逆矩阵。
我希望这可以帮助别人