我可以使用proc sql
中的子查询来计算模式,但这是最简单的方法吗?此代码处理在依赖于计算模式时可能发生的关系,但max
中的proc sql
函数会断开关系。
ods html file = "sas_output.html";
data raw_data;
input cust_id $
category $
amount;
datalines;
A red 72.83
A red 80.22
A blue 0.2
A blue 33.62
A blue 30.63
A green 89.04
B blue 10.29
B red 97.07
B red 68.71
B red 98.2
B red 1.79
C green 92.94
C green 0.96
C red 15.23
D red 49.94
D blue 75.82
E blue 20.97
E blue 78.49
F green 87.92
F green 32.29
;
run;
proc sql;
create table modes as
select cust_id, mean(amount) as mean_amount, category as mode_category
from (
select *, count(1) as count from raw_data group by cust_id, category
)
group by cust_id
having count=max(count)
order by cust_id;
quit;
data modes;
set modes;
by cust_id;
if first.cust_id then output;
run;
data final_data;
merge raw_data modes;
by cust_id;
run;
proc print data=final_data noobs;
title "final data";
run;
ods html close;
我试图像这样使用proc means
proc means data=raw_data;
class cust_id;
var category;
output out=modes mode=mode_category;
run;
但是我收到错误“列表中的变量类别与此列表规定的类型不匹配”,因为proc means
不支持字符变量。
答案 0 :(得分:2)
SQL当然是一种很好的方法。这是一个双DOW loop的数据步骤解决方案。如果需要,您也可以在同一步骤中使用此方法计算平均值。
对数据进行排序以便按组使用。
proc sort data=raw_data;
by cust_id category;
run;
读取数据集一次,按cust_id计算每个类别的出现次数。变量maxfreq
存储最高计数,变量mode
保留具有最高计数的类别。因为数据是按类别变量排序的,所以如果是破坏者,这将返回最高的字母值。
第二个循环从第一个循环输出值和模式。
data want(drop=freq maxfreq);
do until (last.cust_id);
set raw_data;
by cust_id category;
if first.category then freq=0;
freq+1;
maxfreq=max(freq,maxfreq);
if freq=maxfreq then mode=category;
end;
do until (last.cust_id);
set raw_data;
by cust_id;
output;
end;
run;