SAS编程求斜率？

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

data a;
input time$16. intensity;
cards;
06OCT14:07:58:00 0.43
06OCT14:07:58:01 0.45
06OCT14:07:58:02 0.49
06OCT14:07:58:03 0.51
06OCT14:07:58:04 0.53
06OCT14:07:58:05 0.53
06OCT14:07:58:06 0.54
06OCT14:07:58:07 0.34
06OCT14:07:58:08 0.54
;
run;
找出intensity几个峰值，并求相邻最大值与最小值之间的斜率，即(intensity2-intensity1)/(time2-time1),SAS如何编程，恳请高手指教！

扫码加我 拉你入群

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

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

分享0 收藏1 回帖

关键词：SAS编程 intensity cards Input Inte 最大值 如何

1. data a;
   
2. input time datetime16. intensity;
   
3. format time datetime16.;
   
4. cards;
   
5. 06OCT14:07:58:00 0.43
   
6. 06OCT14:07:58:01 0.45
   
7. 06OCT14:07:58:02 0.49
   
8. 06OCT14:07:58:03 0.51
   
9. 06OCT14:07:58:04 0.53
   
10. 06OCT14:07:58:05 0.53
    
11. 06OCT14:07:58:06 0.54
    
12. 06OCT14:07:58:07 0.34
    
13. 06OCT14:07:58:08 0.54
    
14. ;
    
15. run;
    
16. data b;
    
17. set a;
    
18. intensity1=lag(intensity);
    
19. if intensity1^=.;
    
20. run;
    
21. data c;
    
22. set b;
    
23. time=lag(time);
    
24. intensity2=intensity;
    
25. if time^=.;
    
26. drop intensity;
    
27. run;
    
28. proc sql;
    
29. create table d as select b.*,c.intensity2 from b,c where b.time=c.time;
    
30. quit;
    
31. data e;
    
32. set d;
    
33. if (intensity-intensity1)>0 and (intensity-intensity2)>0 then flag="极大值点";
    
34. if (intensity-intensity1)<0 and (intensity-intensity2)<0 then flag="极小值点";
    
35. if flag^='';
    
36. run;
    
37. data f;
    
38. set e;
    
39. retain int t;
    
40. if _n_=1 then do;int=intensity;t=time;end;
    
41. else do;t=time-t;int=intensity-int;end;
    
42. k=int/t;
    
43. if _n_=1 then k=.;
    
44. keep time intensity flag k;
    
45. run;
    
46. proc sql;
    
47. create table g as select a.*,f.flag,f.k from a left join f on a.time=f.time;
    
48. quit;

复制代码

mingfeng07 发表于 2014-11-25 00:57

您这么晚还帮助我解决问题，非常感动也非常感谢！我受益匪浅！不过可以用lead和lag函数是不是更简洁一些，我写的程序如下：
data a;
input obs_time datetime16. intensity;
format obs_time datetime16.;
cards;
06OCT14:07:58:00 0.43
06OCT14:07:58:01 0.45
06OCT14:07:58:02 0.49
06OCT14:07:58:03 0.51
06OCT14:07:58:04 0.53
06OCT14:07:58:05 0.53
06OCT14:07:58:06 0.54
06OCT14:07:58:07 0.34
06OCT14:07:58:08 0.54
;
run;
proc expand data=a out=b method=none;

convert intensity = lead1 / transformout = (lead 1);
convert intensity = lag1 / transformout = (lag 1);
run;

data c;
set b;
if (intensity-lag1)>0 and (intensity-lead1)>0 then flag="极大值点";
if (intensity-lag1)<0 and (intensity-lead1)<0 then flag="极小值点";
if flag^='';
run;
proc sort data=c out=d;
by obs_time;
run;
data e;
set d;
slp=dif(intensity)/dif(obs_time);
run;
欢迎提宝贵建议！

mingfeng07 发表于 2014-11-25 00:57

您这么晚还帮助我解决问题，非常感动也非常感谢！我受益匪浅！不过可以用lead和lag函数是不是更简洁一些，我写的程序如下：
data a;
input obs_time datetime16. intensity;
format obs_time datetime16.;
cards;
06OCT14:07:58:00 0.43
06OCT14:07:58:01 0.45
06OCT14:07:58:02 0.49
06OCT14:07:58:03 0.51
06OCT14:07:58:04 0.53
06OCT14:07:58:05 0.53
06OCT14:07:58:06 0.54
06OCT14:07:58:07 0.34
06OCT14:07:58:08 0.54
;
run;
proc expand data=a out=b method=none;

convert intensity = lead1 / transformout = (lead 1);
convert intensity = lag1 / transformout = (lag 1);
run;

data c;
set b;
if (intensity-lag1)>0 and (intensity-lead1)>0 then flag="极大值点";
if (intensity-lag1)<0 and (intensity-lead1)<0 then flag="极小值点";
if flag^='';
run;
proc sort data=c out=d;
by obs_time;
run;
data e;
set d;
slp=dif(intensity)/dif(obs_time);
run;
欢迎提宝贵建议！

data a;
input obs_time datetime16. intensity;
format obs_time datetime16.;
cards;
06OCT14:07:58:00 0.43
06OCT14:07:58:01 0.45
06OCT14:07:58:02 0.49
06OCT14:07:58:03 0.51
06OCT14:07:58:04 0.53
06OCT14:07:58:05 0.53
06OCT14:07:58:06 0.54
06OCT14:07:58:07 0.34
06OCT14:07:58:08 0.54
;
run;
proc expand data=a out=b method=none;

convert intensity = lead1 / transformout = (lead 1);
convert intensity = lag1 / transformout = (lag 1);
run;

data c;
set b;
if (intensity-lag1)>0 and (intensity-lead1)>0 then flag="极大值点";
if (intensity-lag1)<0 and (intensity-lead1)<0 then flag="极小值点";
if flag^='';
run;
proc sort data=c out=d;
by obs_time;
run;
data e;
set d;
slp=dif(intensity)/dif(obs_time);
run;

gotobest 发表于 2014-11-25 11:30
data a;
input obs_time datetime16. intensity;
format obs_time datetime16.;

这个不错

gotobest 发表于 2014-11-25 11:29
您这么晚还帮助我解决问题，非常感动也非常感谢！我受益匪浅！不过可以用lead和lag函数是不是更简洁一些， ...

恩，你这个采用了expand过程步更方便，而且使用了dif函数更有效率。