SAS两组产品

Question

我正在尝试按两组创建产品：id1和id2。

数据集如下所示：

data test;
input id1 id2 value;
datalines;
1 199001 1.762681948
1 199001 1.775245162
1 199001 1.428673376
1 199002 1.175974146
1 199002 1.236166022
1 199002 1.608842974
1 199003 1.673956674
1 199003 1.05879051
1 199003 1.565500916
2 199002 1.838999925
2 199002 .
2 199002 1.984415322
2 199003 1.096820927
2 199003 1.734215557
2 199003 1.157566337
;
run;

我想在第一列和第二列的基础上创建最后一列value的乘积。即组id1=1和id2=199001的输出应为：1.7626x1.7775x1.429=4.47。I.e输出应该看起来像：

1   199001  4.47059416
1   199002  2.338773875
1   199003  2.774645982
2   199002  .
2   199003  2.201834613

Answer 1

这是SAS中最佳三元组的完美示例：first./last.、retain和按组处理。

data want;
set test;
by id1 id2;
retain product;
if(first.id2) then product = value;
else product = product * value;
if(last.id2) then output;
run;

Answer 2

如果使用值的LOG((，您可以使用乘法转换为加法的事实。然后可以使用SQL的SUM((聚合函数。

proc sql ;
create table want as
select id1,id2,exp(sum(log(value))) as product
from have
group by 1,2
;
quit;

如果我们需要一个缺失值的存在来返回一个缺失的值，那么添加一个CASE语句。您可能还需要添加逻辑来处理零。

create table want as
select id1,id2
, case
when (sum(value=0)) then 0
when (sum(missing(value))) then .
else exp(sum(log(case when (value) then value else 1 end)))
end as product
from have
group by 1,2
;

内部CASE将消除有关尝试获取零或缺失值的LOG((的注释。

该代码使用SAS布尔逻辑为true生成1，为false生成0的事实，并将0或缺失的值评估为false，而任何其他值都被视为true。