众所周知，当为群体拟合回归模型时，群体之间的不等方差会导致麻烦。 R中是否有办法解释这个问题？
本演示文稿中提出了Stata的解决方案http://www.stata.com/meeting/fnasug08/gutierrez.pdf（幻灯片9-17）。幻灯片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中是否有类似的解决方案？