求助：怎么将一行数据和多行数据进行匹配，然后生成新变量？

shiyiming

data a;
input ID$ Date yymmdd8. Num Act;
format Date yymmdd10.;
cards;
001 08-07-01 100 1
001 08-07-02 50  0
001 08-07-06 30  0
001 08-07-10 100  1
001 08-07-15 60  0
002 08-07-01 100 1
002 08-07-05 50 1
002 08-07-06 30 0
;
run;

数据中ID为商品编号，Date为采购或卖出时间，Num为采购或卖出数量，Act识别是采购还是卖出。

希望实现以下目标：按先采购（Act＝1）先卖出（Act＝0）的原则，求每次卖出商品（Act＝0）的加权平均库存时间，这里很麻烦的一个问题是有些卖出对应着几个采购，比如第5行卖出的60中有20来自第一行的采购，另外40则来自第4行的采购。

第2行的库存时间为1天，第3行为5天，第5行为14×20/60＋5×40/60,第8行库存时间为5天

shiyiming

死猪头

shiyiming

敬仰啊，谢谢～～

shiyiming

我喜欢单纯一点
[code:2r5gcp62]
data ahuige;
input ID$ Date yymmdd8. Num Act;
format Date yymmdd10.;
cards;
001 08-07-01 100 1
001 08-07-02 50 0
001 08-07-06 30 0
001 08-07-10 100 1
001 08-07-15 60 0
002 08-07-01 100 1
002 08-07-05 50 1
002 08-07-06 30 0
;
run;

data  buy(rename=(date=indate) drop=i num) sell(drop=i num);
    set ahuige;
    if act=1 then
        do i=1 to num;
        output buy;
        end;
    if act=0 then
        do i=1 to num;
        output sell;
        end;
run;

data temp;
    merge buy sell;
    by id;
    output;
    date=.;
run;

proc sql;
    create table final as
    select distinct id, date, act, avg(date-indate) as avgstore
    from temp( where=(act=0))
    group by id,date ;
[/code:2r5gcp62]

shiyiming

谢谢ahuige
你这个程序是比较简洁，不过如果数据很多的话（比如有100万以上的买记录和卖记录，每个记录的数量num又很大），产生的BUY和SELL会极其大吧，这个问题有没有办法克服掉呢？

shiyiming

很多情况下要代码简洁就有可能要有效率的代价。如果是有100万条等对运算效率的死要求，写法又会不一样了。
不过你的新问题有点含糊，
平均分到没卖完的每组，然后相当于是按各组买的日期分别来算？如果不能平均分怎么办？

现在忙，有时间以后再给你看。

shiyiming

DRAFT东西写了，但没时间QC，明天看看再发出来。

shiyiming

恩，多谢多谢，我等～～

shiyiming

[code:2utn7j7u]data store sellandgift;
    infile datalines truncover;
    input ID$ Date yymmdd8. Num Act Note ;
    format Date yymmdd10.;
    if act=1 and note=. then output store;
    else  output sellandgift;
    datalines;
001 08-07-01 100 1 .
001 08-07-02 10 0 .
001 08-07-06 30 0 .
001 08-07-10 100 1 .
001 08-07-15 60 0 .
001 08-07-16 100 1 .
001 08-07-18 160 1 0
001 08-07-22 100 0 .
001 08-07-25 220 0 .
002 08-07-01 100 1 .
003 08-07-02 100 1 .
003 08-07-03 50 1 .
003 08-07-04 100 0 .
003 08-07-08 30 0 .
003 08-07-09 20 0 .
;
run;
%let offsetid=0;
%macro updateStore(num,note,id,selldate) ;
  %let offsetid=%eval(&offsetid+1);
  %if &note=. %then
    %do;
    data store(keep=id date num act note) avg(keep=id date num diff);
        set store ;
        by id date;
        retain number &num diff 0;
        if id=&id then
          do;
          offnum=min(num,number);
          num=num-offnum;
          number=number-offnum;
          diff=diff+(&selldate-date)*offnum;
          end;
        if num>0 then output store;
        if id=&id and last.id then do;diff=diff/#date=&selldate;output avg; end;
    run;

     proc append base=allavg
          data=avg;
     run;
    %end;
  %if &note=0 %then
    %do;
    proc sql;
      create table store(drop=oldnum ) as
      select *,oldnum+oldnum/sum(oldnum)*&num*(id="&id") as num
      from store(rename=(num=oldnum))
      group by id
      order by id,date;
    %end;
%mend;

data _null_;
  set sellandgift;
  length command $200;
  command=compress('%updateStore('||num||','||note||','||id||','||date||')');
  call execute(command);
run;[/code:2utn7j7u]
我把供应商送的货变成160了，不然都不够卖

shiyiming

谢谢！
你这个是求的每次卖出的平均库存时间，我现在想求的是每次进货的平均库存时间，对着你的程序改了半天，也没成，感觉小改达不到目标，看来学习SAS路还长啊 <!-- s:( --><img src="{SMILIES_PATH}/icon_sad.gif" alt=":(" title="Sad" /><!-- s:( -->  <!-- s:( --><img src="{SMILIES_PATH}/icon_sad.gif" alt=":(" title="Sad" /><!-- s:( -->  <!-- s:( --><img src="{SMILIES_PATH}/icon_sad.gif" alt=":(" title="Sad" /><!-- s:( -->