proc import datafile="d:\spss\syntax\spss15-4.sav" out=work.bjcj;
/*特别提醒:spss的变量名用英文字符,不能用中文,否则,sas无法识别,本实例只要有score字段就可以了*/
run;
proc contents data=bjcj;
run;
proc sql;
create table m as
select score,count(score) as sum
from bjcj; /*把分数计数赋值给summ*/
quit;
proc sort;
by score;/*按照分数大小排序m*/
run;
data m;
set m;
length group1 $10;
pctn=_n_/sum;
if pctn<=0.27 then group1='1';
else if 0.27<pctn<=0.73 then group1='2';
else if 0.73<pctn then group1='3';
run;/*根据百分位分成高(group=3)、中、低(group=1)三组*/
proc sql;
create table qufendu as select avg(score), group1 from m
group by group1;
quit; /*按照分组情况,计算平均值,并把这些值储存至新创建的qufendu表*/
proc transpose data=qufendu out=trsqfd;
run;/*转置qufendu表,新生成为trsqfd表,这样原来的高、中、低三个组就变成了变量*/
data qfdout1;
set trsqfd;
if _n_=1 then qufendu=(col3-col1)/100;/*取第一行观测值,实质上只有一行,进行计算,用col3(高分组均值)减去col1(低分组均值),运用区分度公式计算得到区分度*/
run;
proc print data=qfdout1;/*输出区分度值*/
var qufendu;
run;