我有一组交易数据,以毫秒为单位进行观察但很少见。为了使它成为一个合理的时间序列,我创建了一个宏来构建一个网格(每秒一个观察点)。现在我希望我的数据集适合网格(总结在那一秒发生的卷并使用最后的报价/交易价格。 我非常感谢任何想法。请参阅下面的代码。
最佳
%macro makeGrid;
proc sort data=data.inputdataset; by id date milliseconds descending type order; run;
options nomprint;
data data.outputgrid (keep=id date miliseconds type order);
set data.inputdataset;
by id date;
if first.date;
Type="grid"; Order=0;
%do i=((9*60*60)+(10*60)) %to ((16*60*60)+(50*60));
milliseconds=&i; output;
%end;
run;
options mprint;
%mend makeGrid;
非常绝望的尝试是这样的:
data data.merged;
merge data.outputgrid data.inputdataset;
by id date milliseconds;
if first.id then do;
if milliseconds collate = sum(volume)
run;
网格:
ID date time price volume
ABC 01/01/15 0801000
ABC 01/01/15 0802000
ABC 01/01/15 0803000
ABC 01/01/15 0804000
数据示例:
ID date time price volume
ABC 01/01/15 0800004 1,55 100
ABC 01/01/15 0800110 1,56 200
ABC 01/01/15 0800152 1,52 300
ABC 01/01/15 0800210 1,51 400
ABC 01/01/15 0800352 1,50 200
ABC 01/01/15 0800384 1,51 400
输出:
ID date time price volume
ABC 01/01/15 0801000 1,55 100
ABC 01/01/15 0802000 1,52 500
ABC 01/01/15 0803000 1,51 400
ABC 01/01/15 0804000 1,51 600
如果在那一秒没有交易,则交易量应为0,并定价最后一秒的价格。
编辑:
所以,感谢Chris我没有宏制作网格,但是proc SQL没有按照我的意愿工作。
data grid1;
do seconds = '09:10't to '16:50't ;
output ;
end ;
/* id=.;*/
/* date=.;*/
/* format ric $12. ;*/
/* format Date best12. ;*/
run;
proc sql ;
create table want as
select a.id, a.date, a.miliseconds, sum(a.volume) as Volume
from have a
left join
grid1 b on a.id = b.id
and a.date = b.date
and a.miliseconds = floor(b.seconds)
group by a.id, a.date, a.miliseconds ;
quit;
run;
因此,找不到ID和日期,所以我在数据步骤中创建了它们。但后来id不是同一种格式,我试图修复但失败了。 通过匹配秒网格线之间的毫秒来合并has数据集时,似乎我有一个错误。 我错过了什么?总是感谢投入! 最好的
答案 0 :(得分:2)
首先,为什么要创建一秒钟的宏?
do time = '09:10't to '16:50't ; output ; end ;
假设您的milliseconds
只是小数点后3位的时间,只需相应地加入floor(time)
,然后进行汇总。
proc sql ; create table want as select a.id, a.date, a.time, sum(b.volume) as Volume from grid a left join inputdataset b on a.id = b.id and a.date = b.date and a.time = floor(b.time) group by a.id, a.date, a.time ; quit ;
然后,您可以对此数据集进行后期处理,将缺失的卷设置为零或相应地保留之前的价格。
答案 1 :(得分:0)
首先,谢谢大家的答案。 我想怎么做。可能有一种更简单的方法,但对我来说,这已经足够了:
data data.test;
set data.taq;
next_full_sec = int((miliseconds-1) / 1000) + 1;
run;
data data.grid;
do seconds = 1 to 86400;
output;
end;
run;
proc summary data=data.test;
by ric date time next_full_sec;
var volume;
output out=data.int (drop=_type_) sum=;
run;
proc sort data=data.test (keep=ric date) out=data.id_date nodupkey;
by ric date time;
run;
proc sql;
create table data.testseconds as
select ric, date, time, seconds
from data.id_date, data.grid
order by ric, date, time, seconds
;
quit;
data data.testout;
merge
data.testseconds (in=a rename=(fullseconds=next_full_sec))
data.int (in=b)
;
by ric date time next_full_sec;
if a;
if not b
then do;
_freq_ = 0;
volume = 0;
end;
run;
proc print noobs;
run;