求助各位高手啊 - 经管之家

求助各位高手啊

小鳄鱼a

7674

收藏 2015-05-26

悬赏 1000 个论坛币已解决

想说的几点是：
id有很多个，需要把每个id在每个t时刻的p值都计算出来

多谢了

最佳答案

jl60156 查看完整内容

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

全部回复

jl60156

2015-5-26 23:21:23

proc sql;
create table have1 as
select sum(n) as n, *
from have
group by id, year,half, g;
create table have2 as
select  distinct *
from have1
order by id, year,half, g ;
quit;

proc transpose data=have2 out=trans(drop=_name_);
by id year half;
id g;
var n;
run;
proc sql;
  select distinct id, count(distinct id) into :ids separated by ' ', :idn
  from trans
quit;

%macro want();
%do id=1 %to &idn;
%let idk=%scan(&ids,&id,%str( ));

ods select nlevels;
ods output nlevels=nlev&id;
proc freq data=trans(where=(id=&idk) drop=year half) levels;
run;
ods output close;
proc sql noprint;
  select strip(tablevar) ,strip(tablevar) into :vars separated by ' ',
   :varsd separated by ','
  from nlev&id
  where nnonmisslevels ne 0 and tablevar ne 'id';
quit;
%let nvar=%sysfunc(countw(&vars,' '));
data t(keep=&vars id year half);
set trans;
%do j=1 %to &nvar;
   %let k=%scan(&vars,&j,%str( ));
if missing(&k) then &k=0;
if half=1 then &k=-(&k);
%end;
where id ne &idk;
run;
proc sql;
create table t&id as
select sum(sum(&varsd)) as p,&idk as id,year
from t
group by  year
order by  year;
quit;
%end;
data tall;
set %do id=1 %to &idn; t&id %end;;
run;
proc sort data=have;
by id year;
run;

proc sort data=tall;
by id year;
run;
data want;
merge have(in=a) tall;
by id year;
if a;
run;
%mend;

%want;

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

小鳄鱼a

2015-5-29 10:50:04

在线等啊

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

akalius

2015-5-29 21:15:05

很想帮楼主但是我也不会啊

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

小鳄鱼a

2015-5-29 23:37:36

jl60156 发表于 2015-5-29 23:27
proc sql;
create table have1 as
select sum(n) as n, *

非常感谢，不过显示着一些错误呢

234
235
236
237  proc sql;
238       create table have1 as
239                select sum(n) as n, *
240                from have
241                group by id, year,half, g;
NOTE: 查询要求将汇总统计量与原始的数据重新合并。
WARNING: 变量 n 已经存在于文件“WORK.HAVE1”中。
NOTE: 表 WORK.HAVE1 创建完成，有 75 行，5 列。

242       create table have2 as
243                select  distinct *
244                from have1
245                         order by id, year,half, g ;
NOTE: 表 WORK.HAVE2 创建完成，有 74 行，5 列。

246  quit;
NOTE: “PROCEDURE SQL”所用时间（总处理时间）:
   实际时间       0.01 秒
   CPU 时间       0.01 秒

247
248  proc transpose data=have2 out=trans(drop=_name_);
249       by id year half;
250       id g;
251       var n;
252  run;

NOTE: 从数据集 WORK.HAVE2. 读取了 74 个观测
NOTE: 数据集 WORK.TRANS 有 16 个观测和 13 个变量。
NOTE: “PROCEDURE TRANSPOSE”所用时间（总处理时间）:
   实际时间       0.01 秒
   CPU 时间       0.01 秒

253  proc sql;
254 select distinct id, count(distinct id) into :ids separated by ' ', :idn
255 from trans
256  quit;
NOTE: 查询要求将汇总统计量与原始的数据重新合并。
257
258  %macro want();
259  %do id=1 %to &idn;
260  %let idk=%scan(&ids,&id,%str( ));
261
262  ods select nlevels;
263  ods output nlevels=nlev&id;
264  proc freq data=trans(where=(id=&idk) drop=year half) levels;
265  run;
266  ods output close;
267  proc sql noprint;
268 select strip(tablevar) ,strip(tablevar) into :vars separated by ' ',
269                                                                                                 :varsd separated
269!  by ','
270 from nlev&id
271 where nnonmisslevels ne 0 and tablevar ne 'id';
272  quit;
273  %let nvar=%sysfunc(countw(&vars,' '));
274  data t(keep=&vars id year half);
275       set trans;
276    %do j=1 %to &nvar;
277       %let k=%scan(&vars,&j,%str( ));
278                         if missing(&k) then &k=0;
279                         if half=1 then &k=-(&k);
280    %end;
281       where id ne &idk;
282  run;
283  proc sql;
284       create table t&id as
285                select sum(sum(&varsd)) as p,&idk as id,year
286                from t
287                group by  year
288                         order by  year;
289  quit;
290  %end;
291  data tall;
292       set %do id=1 %to &idn; t&id %end;;
293  run;
294  proc sort data=have;
295       by id year;
296  run;
297
298  proc sort data=tall;
299       by id year;
300  run;
301  data want;
302       merge have(in=a) tall;
303       by id year;
304       if a;
305  run;
306  %mend;
307
308  %want;
WARNING: 未创建输出“nlevels”。请确保输出对象名称、标签或路径拼写正确。
      还需确保使用了相应的过程选项来生成请求的输出对象。如，确保没有使用 NOPRINT 选项。
WARNING: 未创建输出“nlevels”。请确保输出对象名称、标签或路径拼写正确。
      还需确保使用了相应的过程选项来生成请求的输出对象。如，确保没有使用 NOPRINT 选项。
WARNING: 当前的 ODS SELECT/EXCLUDE/OUTPUT 语句被清除，因为检测到 PROC 步的结尾。
      其原因可能包括交互式过程的无终止（键入 quit; 结束过程）和不带输出的运行组。
NOTE: “PROCEDURE SQL”所用时间（总处理时间）:
   实际时间       0.06 秒
   CPU 时间       0.06 秒

NOTE: 从数据集 WORK.TRANS. 读取了 4 个观测
   WHERE id=1;
NOTE: “PROCEDURE FREQ”所用时间（总处理时间）:
   实际时间       0.11 秒
   CPU 时间       0.03 秒

ERROR: 文件“WORK.NLEV1.DATA”不存在。
NOTE: 由于出错，SAS 系统停止处理该步。
NOTE: “PROCEDURE SQL”所用时间（总处理时间）:
   实际时间       0.00 秒
   CPU 时间       0.01 秒

WARNING: 没有解析符号引用 VARS。

NOTE: 由调用宏“WANT”生成行。
4    data t(keep=&vars id year half);       set trans;
               -
               214
               23
WARNING: 没有解析符号引用 VARS。
WARNING: 没有解析符号引用 VARS。
WARNING: 没有解析符号引用 VARS。
NOTE: 由宏变量 K 生成行。
1 &vars
   -
   22
WARNING: 没有解析符号引用 VARS。
NOTE: 由宏变量 K 生成行。
1    &vars
   -
   180
WARNING: 没有解析符号引用 VARS。
NOTE: 由宏变量 K 生成行。
1    &vars
   -
   180
WARNING: 没有解析符号引用 VARS。
WARNING: 没有解析符号引用 VARS。

ERROR 214-322: 变量名 & 无效。

ERROR 23-7: 对于“KEEP”选项值无效。

ERROR 22-322: 语法错误，期望下列之一: 名称, 带引号的字符串, 数值常数, 日期时间常数, 缺失值, INPUT, PUT.

ERROR 180-322: 语句无效或未按正确顺序使用。

NOTE: 由于出错，SAS 系统停止处理该步。
NOTE: “DATA 语句”所用时间（总处理时间）:
   实际时间       0.05 秒
   CPU 时间       0.04 秒

NOTE: 由调用宏“WANT”生成行。
8                                                    create table t&id as                select sum(sum(&varsd)) as
                                                                                                         -
                                                                                                         22
8  !  p,&idk as id,year                from t                group by  year                      order by  year;
8  ! quit;
WARNING: 没有解析符号引用 VARSD。
ERROR 22-322: 语法错误，期望下列之一: 名称, 带引号的字符串, 数值常数, 日期时间常数, 缺失值, ), BTRIM, INPUT, PUT,
            SUBSTRING, USER.

NOTE: 由于出错，SAS 系统停止处理该步。
NOTE: “PROCEDURE SQL”所用时间（总处理时间）:
   实际时间       0.01 秒
   CPU 时间       0.01 秒

NOTE: 数据集 WORK.NLEV2 有 11 个观测和 5 个变量。
NOTE: 从数据集 WORK.TRANS. 读取了 6 个观测
   WHERE id=2;
NOTE: “PROCEDURE FREQ”所用时间（总处理时间）:
   实际时间       0.01 秒
   CPU 时间       0.01 秒

NOTE: “PROCEDURE SQL”所用时间（总处理时间）:
   实际时间       0.00 秒
   CPU 时间       0.00 秒

NOTE: 从数据集 WORK.TRANS. 读取了 10 个观测
   WHERE id not = 2;
NOTE: 数据集 WORK.T 有 10 个观测和 10 个变量。
NOTE: “DATA 语句”所用时间（总处理时间）:
   实际时间       0.01 秒
   CPU 时间       0.01 秒

NOTE: 表 WORK.T2 创建完成，有 3 行，3 列。

NOTE: “PROCEDURE SQL”所用时间（总处理时间）:
   实际时间       0.01 秒
   CPU 时间       0.01 秒

NOTE: 数据集 WORK.NLEV3 有 11 个观测和 5 个变量。
NOTE: 从数据集 WORK.TRANS. 读取了 6 个观测
   WHERE id=3;
NOTE: “PROCEDURE FREQ”所用时间（总处理时间）:
   实际时间       0.01 秒
   CPU 时间       0.01 秒

NOTE: “PROCEDURE SQL”所用时间（总处理时间）:
   实际时间       0.00 秒
   CPU 时间       0.00 秒

NOTE: 从数据集 WORK.TRANS. 读取了 10 个观测
   WHERE id not = 3;
NOTE: 数据集 WORK.T 有 10 个观测和 12 个变量。
NOTE: “DATA 语句”所用时间（总处理时间）:
   实际时间       0.00 秒
   CPU 时间       0.00 秒

NOTE: 表 WORK.T3 创建完成，有 3 行，3 列。

NOTE: “PROCEDURE SQL”所用时间（总处理时间）:
   实际时间       0.01 秒
   CPU 时间       0.01 秒

ERROR: 文件“WORK.T1.DATA”不存在。

NOTE: 由于出错，SAS 系统停止处理该步。
WARNING: 数据集 WORK.TALL 可能不完整。该步停止时，共有 0 个观测和 3 个变量。
WARNING: 数据集 WORK.TALL 由于该步已停止，而没有被替换。
NOTE: “DATA 语句”所用时间（总处理时间）:
   实际时间       0.01 秒
   CPU 时间       0.01 秒

NOTE: 从数据集 WORK.HAVE. 读取了 75 个观测
NOTE: 数据集 WORK.HAVE 有 75 个观测和 5 个变量。
NOTE: “PROCEDURE SORT”所用时间（总处理时间）:
   实际时间       0.01 秒
   CPU 时间       0.01 秒

NOTE: 输入数据集为空。
NOTE: 数据集 WORK.TALL 有 0 个观测和 3 个变量。
NOTE: “PROCEDURE SORT”所用时间（总处理时间）:
   实际时间       0.01 秒
   CPU 时间       0.01 秒

NOTE: 从数据集 WORK.HAVE. 读取了 75 个观测
NOTE: 从数据集 WORK.TALL. 读取了 0 个观测
NOTE: 数据集 WORK.WANT 有 75 个观测和 6 个变量。
NOTE: “DATA 语句”所用时间（总处理时间）:
   实际时间       0.01 秒
   CPU 时间       0.00 秒

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

jl60156

2015-5-30 00:30:44

小鳄鱼a 发表于 2015-5-29 23:37
非常感谢，不过显示着一些错误呢

try again

proc sql;
      create table have1 as
            select sum(n) as n, *
            from have
            group by id, year,half, g;
      create table have2 as
            select  distinct *
            from have1
                     order by id, year,half, g ;
quit;

proc transpose data=have2 out=trans(drop=_name_);
      by id year half;
      id g;
      var n;
run;
proc sql;
  select distinct id, count(distinct id) into :ids separated by ' ', :idn
  from trans
quit;
ods output nlevels=dummy;
proc freq data=trans( drop=year half) levels;
run;
%macro want();
%do id=1 %to &idn;
%let idk=%scan(&ids,&id,%str( ));

ods select nlevels;
ods output nlevels=nlev&id;
proc freq data=trans(where=(id=&idk) drop=year half) levels;
run;
ods output close;
proc sql noprint;
  select strip(tablevar) ,strip(tablevar) into :vars separated by ' ',
                                                                                             :varsd separated by ','
  from nlev&id
  where nnonmisslevels ne 0 and tablevar ne 'id';
quit;
%let nvar=%sysfunc(countw(&vars,' '));
data t(keep=&vars id year half);
      set trans;
%do j=1 %to &nvar;
   %let k=%scan(&vars,&j,%str( ));
                     if missing(&k) then &k=0;
                     if half=1 then &k=-(&k);
%end;
      where id ne &idk;
run;
proc sql;
      create table t&id as
            select sum(sum(&varsd)) as p,&idk as id,year
            from t
            group by  year
                     order by  year;
quit;
%end;
data tall;
      set %do id=1 %to &idn; t&id %end;;
run;
proc sort data=have;
      by id year;
run;

proc sort data=tall;
      by id year;
run;
data want;
      merge have(in=a) tall;
      by id year;
      if a;
run;
%mend;

%want;

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

点击查看更多内容…

小鳄鱼a

2015-5-30 08:02:28

jl60156 发表于 2015-5-30 00:30
try again

proc sql;

可能我表述的意思不够清晰。

结果显示每个id每年的值都一样，但是应该每个时期的值，也就是每半年的值一样。half是半年的意思，1代表上半年，2代表下半年。
id=1从1997年上半年到1998年下半年。
id=2从1996年上半年到1998年下半年。
id=3从1997年上半年到1999年下半年。

我说的更加准确一些
p值得计算思路是，看看其他的id在t+1时刻对该id在t时期所拥有g的有多少增减，以此来作为id在t+1时刻的p值。t代表半年的。从这个思路来讲，id1,id3对于id2在1996年下半年是有影响的，以此来作为id2在1997年上半年的p值。
因为id3的时间段一直到1999年下半年，而id1和id2到1998年下半年，所以id3对于id2和id1在1998年下半年的p肯定有影响，但是id1和id2对id3在1998年下年的p就没有影响了，因为它们没有1999年上半年的数值。

麻烦您了

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

小鳄鱼a

2015-5-30 08:04:53

jl60156 发表于 2015-5-30 00:30
try again

proc sql;

换句话而讲我不是用year作为t时期的间隔而是用year half作为时期间隔

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

jl60156

2015-6-1 10:04:32

不知道是不是你要的,最好给个例子.比如id=1的所有相对应p值相怎么算

proc sql;
      create table have1 as
            select sum(n) as n, *
            from have
            group by id, year,half, g;
      create table have2 as
            select  distinct *
            from have1
                     order by id, year,half, g ;
quit;
data have2;
set have2;
t=year*10+half*5;
run;

proc transpose data=have2 out=trans(drop=_name_);
      by id t;
      id g;
      var n;
run;
proc sql;
  select distinct id, count(distinct id) into :ids separated by ' ', :idn
  from trans
quit;
ods output nlevels=dummy;
proc freq data=trans levels;
run;
%macro want();
%do id=1 %to &idn;
%let idk=%scan(&ids,&id,%str( ));

proc sql;
  select distinct t, count(distinct t) into :ts separated by ' ', :tn
  from trans
  where id=&idk;
quit;
%do t=1 %to &tn;
%let tk=%scan(&ts,&t,%str( ));
ods select nlevels;
ods output nlevels=nlev&id&t;
proc freq data=trans(where=(id=&idk and t=&tk) ) levels;
run;
ods output close;
proc sql noprint;
  select strip(tablevar) ,strip(tablevar) into :vars separated by ' ',
                                                                                             :varsd separated by ','
  from nlev&id&t
  where nnonmisslevels ne 0 and tablevar not in ('id','t');
quit;
%let nvar=%sysfunc(countw(&vars,' '));
data t(keep=&vars id t);
      set trans;
%do j=1 %to &nvar;
   %let k=%scan(&vars,&j,%str( ));
                     if missing(&k) then &k=0;
                     if t=&tk then &k=-(&k);
%end;
where id ne &idk and &tk<=t<=%eval(&tk+5);
run;
proc sql;
      create table t&t as
            select sum(sum(&varsd)) as p,&idk as id,%eval(&tk+5) as t
            from t;
quit;
%end;
data tid&id;
      set %do t=1 %to &tn; t&t %end;;
run;
%end;
data tall;
      set %do id=1 %to &idn;  tid&id %end;;
year=floor(t/10)-(mod(t,10)=0);
half=mod(t/5,2)+(mod(t/5,2)=0)*2;
run;
proc sort data=have;
      by id year half;
run;

proc sort data=tall;
      by id year half;
run;
data want;
      merge have(in=a) tall(drop=t);
      by id year half;
      if a;
run;
%mend;

%want;

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

小鳄鱼a

2015-6-1 15:01:28

jl60156 发表于 2015-6-1 10:04
不知道是不是你要的,最好给个例子.比如id=1的所有相对应p值相怎么算

您好请看看sheet2 ，我回复有个地方讲错了，表中也重新做了说明，这个表更加简单直观一些

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

小鳄鱼a

2015-6-1 15:02:45

jl60156 发表于 2015-6-1 10:04
不知道是不是你要的,最好给个例子.比如id=1的所有相对应p值相怎么算

请看sheet2

附件列表

求助文件.xlsx

大小:17.55 KB

只需: 1000 个论坛币马上下载

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

jl60156

2015-6-1 21:27:02

小鳄鱼a 发表于 2015-6-1 15:02
请看sheet2

如果没理解错,最后发的就是你要的

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

小鳄鱼a

2015-6-1 21:33:51

jl60156 发表于 2015-6-1 10:04
不知道是不是你要的,最好给个例子.比如id=1的所有相对应p值相怎么算

di3在1999 1以后（包括1999 1）是没有取值的.因为id1 或 id2截止到1998 2 ，而考察id3 1999 1的p取值应该是基于其他id在1998 2到 1999 1在这段时间对于id3在1998 2 持有的abcd的增减情况。

而程序中是 1999 1 是有取值的

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

jl60156

2015-6-1 23:06:20

小鳄鱼a 发表于 2015-6-1 21:33
di3在1999 1以后（包括1999 1）是没有取值的.因为id1 或 id2截止到1998 2 ，而考察id3 1999 1的p取值应该 ...

proc sql;
      create table have1 as
            select sum(n) as n, *
            from have
            group by id, year,half, g;
      create table have2 as
            select  distinct *
            from have1
                     order by id, year,half, g ;
quit;
data have2;
      set have2;
      t=year*10+half*5;
run;

proc transpose data=have2 out=trans(drop=_name_);
      by id t;
      id g;
      var n;
run;
proc sql;
  select distinct id, count(distinct id) into :ids separated by ' ', :idn
  from trans
quit;
ods output nlevels=dummy;
proc freq data=trans levels;
run;
%macro want();
%do id=1 %to &idn;
%let idk=%scan(&ids,&id,%str( ));

proc sql;
  select distinct t, count(distinct t) into :ts separated by ' ', :tn
  from trans
  where id=&idk;
quit;
%do t=1 %to &tn;
%let tk=%scan(&ts,&t,%str( ));
ods select nlevels;
ods output nlevels=nlev&id&t;
proc freq data=trans(where=(id=&idk and t=&tk) ) levels;
run;
ods output close;
proc sql noprint;
  select strip(tablevar) ,strip(tablevar) into :vars separated by ' ',
                                                                                             :varsd separated by ','
  from nlev&id&t
  where nnonmisslevels ne 0 and tablevar not in ('id','t');
quit;
%let nvar=%sysfunc(countw(&vars,' '));
data t(keep=&vars id t);
      set trans;
%do j=1 %to &nvar;
   %let k=%scan(&vars,&j,%str( ));
                     if missing(&k) then &k=0;
                     if t=&tk then &k=-(&k);
%end;
where id ne &idk and &tk<=t<=%eval(&tk+5);
run;
proc sql;
  select  count(distinct t) into :maxt
  from t;
quit;
%if &maxt eq 2 %then %do;
proc sql;
      create table t&t as
            select sum(sum(&varsd)) as p,&idk as id,%eval(&tk+5) as t
            from t;
quit;
%end;
%else %do;
data t&t;
run;
%end;
%end;
data tid&id;
      set %do t=1 %to &tn; t&t %end;;
run;
%end;
data tall;
      set %do id=1 %to &idn;  tid&id %end;;
            year=floor(t/10)-(mod(t,10)=0);
            half=mod(t/5,2)+(mod(t/5,2)=0)*2;
run;
proc sort data=have;
      by id year half;
run;

proc sort data=tall;
      by id year half;
run;
data want;
      merge have(in=a) tall(drop=t);
      by id year half;
      if a;
run;
%mend;

%want;

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

小鳄鱼a

2015-6-2 22:48:58

jl60156 发表于 2015-6-1 23:06
proc sql;
create table have1 as
select sum(n) as n, *

非常感谢您的帮助有时间的话请再帮我稍微修改一下程序修改完后程序再贴在这即可我再发一个300论坛币的帖子您去领取就可以总体跟例子一样  只是在处理的时候对于每只g要除以所赋值id所持有的g的n    举例来说

Id1 在 1997 1持有  abcdef  ,  分别持有 10 9 10 11 12 13
id2在 1997 1到1997 2拥有的g与id1的交集为abcdef, 那么  a相对id1持有的变化比例（0-11）/10，b ( 0-12)/9 c （17-13）/10 d (18-14)/11  e (19-15)/12 f (0-16)/13  然后求和
Id3在 1997 1到1997 2拥有的g与id1的交集为bcdef, 那么  b  （0-12）/9 c (0-13)/10 d (0-14)/11  e( 0-15)/12 f (18-0)/13  然后求和

加总后，赋给 id1  1997 2

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

jl60156

2015-6-3 09:35:01

小鳄鱼a 发表于 2015-6-2 22:48
非常感谢您的帮助有时间的话请再帮我稍微修改一下程序修改完后程序再贴在这即可我再发一个300论 ...

proc sql;
      create table have1 as
            select sum(n) as n, *
            from have
            group by id, year,half, g;
      create table have2 as
            select  distinct *
            from have1
                     order by year,half,id, g ;
quit;
data have2;
      set have2;
      t=year*10+half*5;
run;

proc transpose data=have2 out=trans(drop=_name_);
      by t id;
      id g;
      var n;
run;
proc sql;
  select distinct id, count(distinct id) into :ids separated by ' ', :idn
  from trans
quit;
ods output nlevels=dummy;
proc freq data=trans levels;
run;
%macro want();
%do id=1 %to &idn;
%let idk=%scan(&ids,&id,%str( ));

proc sql;
  select distinct t, count(distinct t) into :ts separated by ' ', :tn
  from trans
  where id=&idk;
quit;
%do t=1 %to &tn;
%let tk=%scan(&ts,&t,%str( ));
ods select nlevels;
ods output nlevels=nlev&id&t;
proc freq data=trans(where=(id=&idk and t=&tk) ) levels;
run;
ods output close;
proc sql noprint;
  select strip(tablevar) ,strip(tablevar) into :vars separated by ' ', :varsd separated by ','
  from nlev&id&t
  where nnonmisslevels ne 0 and tablevar not in ('id','t');
quit;
%let nvar=%sysfunc(countw(&vars,' '));
data factor(keep=&vars t %do j=1 %to &nvar;
   %let k=%scan(&vars,&j,%str( ));
            rename=(&k=_&k)
%end;);
      set trans;
if id eq &idk and t=&tk then do;
         t=&tk;output;
         t=%eval(&tk+5);output;
         end;
run;
data t(keep=&vars id t);
      merge trans(in=in where=(id ne &idk and &tk<=t<=%eval(&tk+5))) factor;
                  by t;
%do j=1 %to &nvar;
   %let k=%scan(&vars,&j,%str( ));
                     if missing(&k) then &k=0;
                     if t=&tk then &k=-(&k);
                                                            &k=&k/_&k;
%end;
         if in;
run;
proc sql;
  select  count(distinct t) into :maxt
  from t;
quit;
%if &maxt eq 2 %then %do;
proc sql;
      create table t&t as
            select sum(sum(&varsd)) as p,&idk as id,%eval(&tk+5) as t
            from t;
quit;
%end;
%else %do;
data t&t;
run;
%end;
%end;
data tid&id;
      set %do t=1 %to &tn; t&t %end;;
run;
%end;
data tall;
      set %do id=1 %to &idn;  tid&id %end;;
            year=floor(t/10)-(mod(t,10)=0);
            half=mod(t/5,2)+(mod(t/5,2)=0)*2;
run;
proc sort data=have;
      by id year half;
run;

proc sort data=tall;
      by id year half;
run;
data want;
      merge have(in=a) tall(drop=t);
      by id year half;
      if a;
run;
%mend;

%want;

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

小鳄鱼a

2015-6-3 11:39:48

jl60156 发表于 2015-6-3 09:35
proc sql;
create table have1 as
select sum(n) as n, *

您好，我后来想了想  ，改进的那个指标可能也不大合理  ，  麻烦您如果有空再帮我修改一下  ，跟之前也差不多  。

Id1 在 1997 1持有  abcdef  ,  分别持有 10 9 10 11 12 13 到 1997  2  abcdef的变动为 a  (14-10) b (8-9) c  (9-10) d  (10-11 ) e  (0-12) f  ( 0-13)  对于这些变动，取绝对值后  ，再带入下面的运算

id2在 1997 1到1997 2拥有的g与id1的交集为abcdef, 那么  a相对id1持有的变化比例（0-11）/4，b ( 0-12)/1 c （17-13）/1  d (18-14)/1  e (19-15)/12 f (0-16)/13  r然后求和
Id3在 1997 1到1997 2拥有的g与id1的交集为bcdef, 那么  b  （0-12）/1 c (0-13)/1 d (0-14)/1  e( 0-15)/12 f (18-0)/13  然后求和

加总后，赋给 id1  1997 2

您

我建个600论坛币的帖子，一会给您贴出网址，一点谢意。您把代码贴到那个帖子上就行。

https://bbs.pinggu.org/forum.php? ... =3742107&extra=

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

小鳄鱼a

2015-6-7 19:54:27

jl60156 发表于 2015-6-3 09:35
proc sql;
create table have1 as
select sum(n) as n, *

您好，可能又要麻烦您了，您若有时间的话，请在14楼那个算法的基础上，也就是发表在（发表于 2015-6-1 23:06:20 ）只是用绝对数值想减的基础上，做如下修改。

现在新引入一个变量m 。上面的做法是采用当期的n减去上一期的n。这次的做法概括为用当期的n减去上一期的m。

举例说明

Id1 在 1997 1持有  abcdef  ,

id2在 1997 1到1997 2拥有的g与id1的交集为abcdef, 那么  持有的变化分别是a:0-5(以前是0-11)，b  0-6（以前是0-12） c 17-7（以前17-13）  d 18-8 （以前18-14） e 19-9（以前19-15） f 0-10（以前0-16）然后求和

Id3在 1997 1到1997 2拥有的g与id1的交集为bcdef, 计算方法同样如id2

这次我开个500论坛币的帖子您写好以后跟我说一下然后我开个帖子       如果还有时间最好修改一下另外两种方法的程序    当然如果修改的部分一样请直接跟我说就可以了

非常感谢

加个图片您随便加一个m变量就可以了

附件列表

1.jpg

原图尺寸 85.87 KB

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

栏目导航

最佳答案

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

分享

扫码加好友，拉您进群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群