求GLM->MIANALYZE过程指点

[请勿回复，暂时已解决]




例如我现在有如下数据表格，我想用MI（3次填补）->GLM（建模分析）->MIANALYZE（MI结果综合）三个步骤实现，请老师给我一点提示，我想先把数据模拟弄好，心里能踏实点。

我第一步得到的MI结果：
原表：

N	group	v1	v2
1	1	11	33
2	2	22	44

整理后：

_Imputation_	V(变量)	Group	Value	N(观测)
1	1	1	11	1
1	1	2	22	2
1	2	1	33	1
1	2	2	44	2
2	1	1	11	1
2	1	2	22	2
2	2	1	33	1
2	2	2	44	2
3	1	1	11	1
3	1	2	22	2
3	2	1	33	1
3	2	2	44	2

我的最终想法是这样：
得到v1和v2共两个P值，就是固定v看group组间的统计学意义，我不考虑v1与v2之间的交互作用。

其实我想用这个办法做：

1. proc glm data=mia;
   
2. class group visit;
   
3. model value=group visit group*visit;
   
4. lsmeans group*visit /pdiff;
   
5. run;

复制代码

直接去看group*visit的最小二乘均数，不过结果做出来很有问题。

后来看到MI过程是必须用MIANALYZE做结果整合的，MIANALYZE过程真的不知道是起啥作用的，我照猫画虎的写了如下，当然得到的不是想要的group组间比较的结果，请各位老师帮我查错解释一下，非常感谢，谢谢！

1. proc glm data=mitemp;
   
2. class visit;
   
3. model value=group visit/solution inverse ;
   
4. by _imputation_;
   
5. ods output ParameterEstimates=glmparms
   
6.            InvXPX=glmxpxi;
   
7. quit;
   
8. 
   
9. 
   
10. proc mianalyze parms=glmparms xpxi=glmxpxi;
    
11.    modeleffects Intercept group visit;
    
12. run;

复制代码

自己顶一下，拜托学过MI过程的朋友帮下忙，我现在最关键的问题是：这个MIANALYZE到底到底是干什么的？

我跳过MIANALYZE过程直接用GLM过程对MI填补过后的数据集进行统计分析，发现若是参照GLM过程的最小二乘法估计，多次模拟后的一类错误非常高，可是不知道怎么将GLM的结果和MIANALYZE过程进行衔接。

mi 过程是用于填补缺失值的；通常采用多重填补来填补缺失值；
缺失值填补完成后，采用SAS标准的过程步来分析每一次_impuation_的数据；
而MIANALYZE 作用则是：
The MIANALYZE procedure reads the parameter estimates and associated covariance matrix that are computed by the standard statistical procedure for each imputed
data set. The MIANALYZE procedure then derives valid univariate and multivariate
inferences for these parameters；
通过标准过程步，得到的参数估计和协方差矩阵等等，对这些变量派生出单一或者多元的
推断；简言之，就是综合一下上述每一次_imputation_的分析结果。