各位大神,我的数据中存在缺失,我使用MI过程填补数据后,使用MIXED model 计算了5次填补的结果,最后使用mianalyze过程进行整合,但是我分析的变量里存在分类变量,mianalyze跑完后,赋值最小那一类被当做了参照,结果为0,但是我需要让每一类都有结果,不知道要怎么处理。
ods select none;
proc mixed data=XX;
class no1 period(ref="1");
by _Imputation_;
model gwg=trimester trimester*pollution/noint solution CL;
random int/subject=no1 type=un;
estimate 'period 1' period*pollution 0 0 1/cl;
estimate 'period 2' period*pollution 1 0 0/cl;
estimate 'period 3' period*pollution 0 1 0/cl;
ods output SolutionF=mixparms;
run;
ods select all;
proc mianalyze parms=mixparms;
class period;
modeleffects period period*pollution ;
run;
其中period为分类变量,赋值分别为1/2/3,pollution为连续型变量
最后结果为:
Parameter Estimates
ParameterperiodEstimateStd Error95% Confidence LimitsDFMinimumMaximumTheta0t for H0:
Parameter=Theta0Pr > |t|
pollution*period1.0000000.0908320.043826...0.0908320.0908320..
pollution*period2.0000000.0074310.053884-0.103820.11868623.822-0.0441380.03529700.140.8915
pollution*period3.0000000.0583150.055686-0.054340.17096738.8350.0382290.10857801.050.3015有没有什么办法让period为1的那组分类也能算出数据么?