用proc mean过程

用by

分组汇总比较简单,关键是按收入的高低来分为三组的实现.

proc univariate data=temp;

var inc;

output out=temp1 pctlpts=33 67 pctlpre=pre;

run;

输出数据集将包含33%分位数和67%分位数.然后建立分组变量.分组别来计算平均值.

data a;

set temp;

if _n_=1 then set temp1;

if inc<=pre33 then gr="第一组";

else if inc<=pre67 then gr="第二组";

else gr="第三组";

run;

proc sql;

select avg(inc) as avginc,

avg(exp) as avgexp,

avg(tax) as avgtax

from a

group by gr;

quit;

楼上的贴 well-done! 我再加一些：

1. proc sql 部分,与以下code相同:

proc means data=a mean std median p25 p75 min max maxdec=2;

class gr;

var inc exp tax;

run;

Also you can use proc tabulate, proc summary, proc univariate to get 收入、支出、所得税的均值.

2. 如果你想算P 值， 比较三组的收入、支出、所得税的均值的不同，use proc glm:

proc glm data=a;

model inc=gr;

run;

proc glm data=a;

model exp=gr;

run;

proc glm data=a;

model tax=gr;

run;

inc 部分的P值没啥意思，因为你就是inc用分的三组。有点象自己跟自己比.

[此贴子已经被作者于2007-3-25 8:14:39编辑过]