SAS水平测试，动脑子一起看看吧，高手入

zhengbo8 发表于 2013-6-19 16:25
正则表达式匹配出数值后再操作。

正则表达式？

我是新手，想要学习SAS，刚从图书馆借了几本书回来，却发现书上的SAS版本是9.1.2和v8.2，而且还必须在xp或者vista环境下运行，我的是win7 64位。想问大侠们9.3版本和以上的9.1.2和v8.2有什么区别？若我用9.1.2版本的书来自学可以么（因为图书馆关门了，所以暂时只能用这2个版本的书）

kuhasu 发表于 2013-7-16 21:30
正则表达式？

这个问题，按照你自己提出的思路RUN_MACRO in FCMP,也就是五六行的程序。
京剧

神的境界是不知道1+1等于几

jingju11 发表于 2013-7-20 11:55
这个问题，按照你自己提出的思路RUN_MACRO in FCMP,也就是五六行的程序。
京剧

几年前大胖猫给了段run macro的，五六行搞不定，你试试看
最新的dosubl()搞出来的也奖励1000好了。

boe 发表于 2013-7-20 12:11
神的境界是不知道1+1等于几

a bug of sas

yatming 发表于 2013-2-27 11:35
Although STRPARSER is too obtuse. Please have a try.

你的代码有点小缺陷
当括号是用于计算，而不是用于函数时，你的代码会出现问题
如（a+b）**c这种形式，我这边测试结果好像不对
我在你的基础上，添加了一点，使得它可以处理上述形式的算式了

1. proc fcmp outlib=work.funcs.test;
   
2.     function strparser(expression $) $ 100;
   
3.         length x y $100;
   
4.         x = ''; y = ''; m = 0; n = 0;
   
5.         l = length(expression);
   
6.         do i = 1 to l;
   
7.             b = substr(expression, i, 1);
   
8.             if prxmatch('/[a-z]/', b) then
   
9.             do;*1;
   
10.                 y = cats(y, b);
    
11.             end;*1;
    
12.             else if b='(' then do;*3;
    
13.                 if y ne '' then do;*2;
    
14.                     y = cats('%sysfunc(', cats(y, b));
    
15.                     x = cats(x, y);
    
16.                     y = '';
    
17.                     m = m+1;
    
18.                 end ;*2;
    
19.                                 else do ;
    
20.                  y = cats('%eval(', cats(y, b));
    
21.                     x = cats(x, y);
    
22.                     y = '';
    
23.                     m = m+1;
    
24.                                         end;
    
25.                 n = n + 1;
    
26.             end;*3;
    
27.             else if b = ')' then do;*5;
    
28.                 n = n - 1;
    
29.                 if m = n+1 then do;*4;
    
30.                     x = cats(x, cats(cats(ifc(y='','',cats('&',y)),b),')'));
    
31.                     y = '';
    
32.                     m = m - 1;
    
33.                 end;*4;
    
34.             end;*5;
    
35.             else do;*6;
    
36.                 x = cats(x, cats(ifc(y='','',cats('&',y)),b));
    
37.                 y='';
    
38.             end;*6;
    
39.         end;
    
40.         x = cats(x, ifc(y='','',cats('&',y)));
    
41.         return (x);
    
42.     endsub;
    
43. run;
    
44. 
    
45. options cmplib=work.funcs;
    
46. 
    
47. data _null_;
    
48.     length x $ 1000;
    
49.     input a b c d t:$100.;
    
50.     array y _numeric_;
    
51.     do over y;
    
52.         call symput(lowcase(vname(y)),y);
    
53.     end;
    
54.     x = cats('%sysevalf(',strparser(compress(lowcase(t))),')');
    
55.     call symput('formula', x);
    
56.     result = input(resolve('&formula'),best.);
    
57.     put result;
    
58.     cards;
    
59.     1 2 3 4 a+cos(sum(c,exp(d+log(b))))
    
60.     10 15 20 25 sin(a)+cos(sum(b,exp(c+log(d))))
    
61.     5 6 7 8 b/d
    
62.         1 2 3 4 a*cos(b)
    
63.         1 2 3 4 c*(a+b)**c
    
64. ;
    
65. run;
    
66. data tmp;
    
67.   a=(1+2)**3;
    
68.   put a;
    
69. run;

复制代码

蓝水澈 发表于 2013-7-19 15:51
我是新手，想要学习SAS，刚从图书馆借了几本书回来，却发现书上的SAS版本是9.1.2和v8.2，而且还必须在xp或者 ...

没什么问题。教材一般不涉及新版本更新取消/变化的函数和功能，都是比较基础的。

有论坛币送！

kuhasu 发表于 2013-7-26 20:26
几年前大胖猫给了段run macro的，五六行搞不定，你试试看
最新的dosubl()搞出来的也奖励1000好了 ...

不好意思.自己作了测试,虽然程序短些,但是效率太差,不好意思现丑:
大概是:
data _null_; call symputx('y', %sysfunc(dequote(x))); run;
把这个放到宏里,然后利用
PROC FCMP; RC =RAN_MACRO(...);
我把程序检验了以下,运行一个F =SIN(COS(SUM(1,2,3)**2+1)); 10万次,需要6S,但是直接运行原方程,只需要大概1S.
京剧