因为混合截面数据的每一年所含样本不太一样那么每一年对应的空间矩阵应该也不一样，这种情况下分析的代码如何实现？