R中的g2反向

Question

在解决Proc混合中的正规方程时，SAS显然使用了一种叫做g2 inverse的东西（Pringle和Rayner，1971）。 R中的ginv给出了不同的moore-penrose逆。 R中是否有任何包可以作为SAS给我g2逆？对我来说非常重要的是，我得到了与SAS相同的广义逆，g2逆而不是moore-penrose。基本上，当（X＆＃39; X）矩阵不是满秩时，我需要在R中获得与SAS中相同的固定效果估计。 谢谢。 任何帮助或建议将不胜感激。