0

我很想知道我是否可以在不使用 SAS IML 的情况下在 SAS Base 中处理这种情况;

假设我有向量

            a    b     c     d    e        f
           1001 JPN 10,000  50%  JPN    2,000
           1001 EUR 12,648  100% EUR    3,000
           1001 USD 15,997  50%  USD    5,000
           1001 JPN 20,233  20%  JPN    8,000
           1001 EUR 25,591  20%  EUR    9,000
           1001 USD 32,368  50%  USD    4,000
           1002 JPN 28,393  50%  JPN    6,000
           1002 EUR 24,906  100% EUR    4,000
           1002 USD 21,847  50%  USD    8,000
           1002 TRY 19,164  20%  JPN    6,000
           1002 EUR 16,811  50%  EUR    15,000
           1002 USD 14,746  100% USD    52,000
           1003 USD 10,000  50%  XVN    8,000



 %macro;
 % let i = 1;
 data want;
 set have;
 %do %while a[&i]=a[eval(&i+1)] ;
 b = &i;
 &i=eval(&i+1);
 %end
 %mend

我想要做的是 a 与 b=e 取 max(c) 和 max(f) 的差并将这个差与 d 相乘,然后对每个不同的 a 求和这些结果。这将是迭代的。我在这里组成的表格只是案例的一个小代表。

谢谢

4

3 回答 3

1

听起来您所说的“for a with b=e”的意思是您想通过 A 处理 WHERE B=E 的记录。
因此,在 SAS 中,如果您想过滤要使用的记录,您将使用 WHERE 语句或选项或 IF 语句。要处理具有一组变量的相同值的记录组,您将使用 BY 语句(或 PROC SQL 代码中的 GROUP BY)。所以这是您的示例数据集。

data have ;
  length a 8 b $3 c d 8 e $3 f 8;
  informat c f comma32. d percent. ;
  input a b c d e f ;
cards;
1001 JPN 10,000  50%  JPN    2,000
1001 EUR 12,648  100% EUR    3,000
1001 USD 15,997  50%  USD    5,000
1001 JPN 20,233  20%  JPN    8,000
1001 EUR 25,591  20%  EUR    9,000
1001 USD 32,368  50%  USD    4,000
1002 JPN 28,393  50%  JPN    6,000
1002 EUR 24,906  100% EUR    4,000
1002 USD 21,847  50%  USD    8,000
1002 TRY 19,164  20%  JPN    6,000
1002 EUR 16,811  50%  EUR    15,000
1002 USD 14,746  100% USD    52,000
1003 USD 10,000  50%  XVN    8,000
run;

首先,您可以使用 PROC Summary 找到 C 和 F 变量的最大值。

proc summary data=have nway ;
  by a ;
  where b=e ;
  var c f ;
  output out=summary max= / autoname ;
run;

然后您可以与详细记录合并,以应用计算差异并乘以百分比。

data middle ;
  merge have summary ;
  by a ;
  if b=e ;
  diff = c_max - f_max ;
  product = diff * d ;
run;

然后您可以再次使用 PROC Summary 来查找每个组内的总和。

proc summary data=middle nway ;
  by a ;
  var product ;
  output out=want sum=new_var ;
run;

总计

Obs      a     _TYPE_    _FREQ_     new_var
1     1001       0         6       67767.2
2     1002       0         5      -82624.5
于 2015-09-18T13:00:59.227 回答
1

感谢您的发表!我假设您在b = e时想要最大值,但为了以防万一,我会给您两种可能的解决方案。请注意,您只需要proc sql在任何 SQL 计算组中使用一次,但为了清楚起见,我们只需要使用两次。

第 1 步:在b = e时获取cf的最大值

proc sql noprint;
    create table maxes as
        select b, e, 
               max(c) as max_c, max(f) as max_f
        from have
        where upcase(b) = upcase(e)
        group by b, e
    ;
quit;

步骤2:将b = e时的最大值放入表中,将所有情况总结为a where b = e

proc sql noprint;
    create table want as
        select a, sum(result) as result

        /* Get max values into the table. 
           Only interested in cases where a = b */
        from(select st1.a, st1.b, st1.e, 
                    max_c, max_f, 
                    (max_c - max_f)*st1.d as result
             from have as st1
             INNER JOIN
                  maxes as st2
             ON upcase(st1.b) = upcase(st2.b)
                AND upcase(st1.e) = upcase(st2.e)
             where upcase(st1.b) = upcase(st1.e) )
        group by a
    ;
quit;

现在,如果要对cf的表最大值执行此操作,可以使用宏变量、按a排序的表、按组处理Sum 语句来执行此操作:

第 1 步:将cf的最大值读入宏变量

proc sql noprint;
    select max(c), max(f)
    into :max_c, :max_f
    from have;
quit;

第 2 步:将每个a的结果相加

data want;
    set have;
    by a;
    where upcase(b) = upcase(e);

    result+( (&max_c - &max_f)*d );

    if(last.a) then do;
        output;
        call missing(result); *Reset the sum for the next a group;
    end;

    keep a result;
run;

我希望这有帮助!这只是解决问题的一种方法,但还有许多其他好方法。这完全取决于您的目标、环境和编程风格。

于 2015-09-18T14:17:03.357 回答
0

因此,基本数据步骤或 proc 调用已经是对输入数据的循环。如果您想使用数据集作为宏调用的参数值的来源,那么您需要一个简单的驱动程序。

假设您已经定义了一个宏,它执行一系列复杂的数据步骤和过程步骤,并采用一个或多个输入参数。然后,您可以创建一个参数值数据集,并在一个简单的数据步骤中使用 CALL EXECUTE 生成对传递参数值的宏的调用。因此,如果宏 %MYMACRO() 定义了 5 个位置参数,并且您的数据集 HAVE 具有五个变量 PARM1 到 PARM5 和 10 个观察值,那么这个数据步骤将生成 10 个宏调用。

 data _null_;
    set have;
    call execute(cats('%nrstr(%mymacro)(',catx(',',of parm1-parm5),')'));
 run;
于 2015-09-18T11:35:37.437 回答