全部版块 我的主页
论坛 休闲区 十二区 灌水吧
991 0
2018-10-09
各位大神,我的数据中存在缺失,我使用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的那组分类也能算出数据么?
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群