众所周知的事实是，在为组拟合回归模型时，组之间的不等方差会导致麻烦。R中有没有办法解决这个问题？
本演示文稿http://www.stata.com/meeting/fnasug08/gutierrez.pdf（幻灯片 9-17）建议了 Stata 的解决方案。幻灯片 14 中显示的实际代码如下：
# where I might normally specify
. xtmixed . . . || id: boy ageXboy girl ageXgirl, nocons cov(un)
# instead I want
. xtmixed . . . || id: boy ageXboy, nocons cov(un) || id: girl ageXgirl, nocons cov(un)
R中有类似的解决方案吗？