我有一个proc glm SAS代码,可以删除大约260个观察值,这很好。但是,我想知道对于某个变量的每个级别使用了多少个观察值。
我如何获得一些观测表,但按变量而不是整体列出?
由于您要排除具有一个缺失值的观测值,因此只需从模型之外获得结果即可(proc glm),进入一个简单的数据集。
data want;
length misfl $1;
set have;
array mvars[*] var1 var2 var3;
do i=1 to dim(mvars);
if mvars[i]=. and misfl='' then misfl='Y';
end;
drop i;
run;
proc freq data=want noprint;
by byvar;
table misfl / out=frwant;
run;
如果要在模型规范内获得此结果,则以下是一个有效的替代方法。由于对于排除的观测值,SAS不会计算残差(观测值在第一步被删除,因此不包括在估计中),因此您可以输出残差数据集,并将其添加到模型规范中
output out=want00 r=resid;
然后您将像我在上一个代码块中一样,使用简单的proc freq搜索缺失值(排除观察值->估计未完成->残差未计算)。
本文收集自互联网,转载请注明来源。
如有侵权,请联系[email protected] 删除。
我来说两句