sas函数与call子程序.ppt

第3章 SAS函数与CALL子程序;SAS函数定义 ;函数用于组合表达式可以简化编程和统计计算;用函数处理多变量可以简化程序;函数自变量;自变量个数： 有些函数不需要自变量，有些只需要一个自变量，有些函数有多个自变量。SAS9.1中，自变量个数不能超过32767。 如果函数的自变量是表达式，在函数被调用之前，所有表达式自变量首先计算。 自变量表示法： 当函数有多个自变量时，必须用逗号分隔开。 自变量名缩写方法： 函数名(OF变量名1…变量名n ) ;例3.5 正确与错误的表示方法。 例3.6 ???组元素自变量的缩写方法。 array y{10} y1-y10; /*定义数组y{10} */ x=sum (of y{*}); z=sum (of y1-y10); 例中，两个SUM函数的表示法是等价的。 ;函数结果;用PUT语句在LOG窗口显示函数值;日期时间函数;应用举例 ;例3.10 以日为单位计算当前日期的天数。 data; x=date (); y=today(); put x= y=; run; 例中，DATE()和TODAY()结果一样，都是返回当天的天数。当然，也可以用一定的日期格式表示它们。 data; x=date (); y=today(); format x y yymmdd10.; put x= y=; run;;例3.12 计算日期值所在的年季月以及处于某月的第几天。 data a; set ResDat.stk000001 (obs=10); keep date year qtr month day; year=year (date); qtr=qtr (date); month=month (date); day=day (date); proc print noobs; run;;例3.14 将日期时间值换算为以秒计数的计算方法。 data; mdy=mdy (08,18,2001); put mdy=; format mdy yymmdd10.; hms=hms (21,50,51); put hms=; x=21*60*60+50*60+51; put x=; dhms=dhms (date (), 21,50,51); put dhms=; dhms=dhms (15263,21,50,51); put dhms=; y=15263*24*3600+x; put y=; run;;概率分布函数;标准正态分布 ;样本统计函数 ;均值 ;求和 ;方差 ;随机数函数 ;正态分布 ;例3.28 用函数NORMAL产生正态分布随机数。 data RV; retain _seed_ 0; mu=0; sigma=1; do _i_ = 1 to 1000; normal1 = 0 + 1 * normal(_seed_);/*均值为mu,标准差为sigma*/ output; end; drop _seed_ _i_; run;;均匀分布 ;例3.30 用函数UNIFORM产生均匀分布随机数。 data RV1; retain _seed_ 0; a=-1; b=2; do _i_ = 1 to 1000; uniform1 = a + (b - a) * uniform(_seed_); /*区间[a, b]上的均匀分布*/ output; end; drop _seed_ _i_; run; ;随机数函数自变量SEED ;SAS CALL 子程序 ;随机数子程序 ;例3.38 使用随机函数产生两个随机数变量属于同一个随机数流。 data RV; retain seed1 seed2 161321804; do I=1 to 5; x1= ranuni(seed1); x2= ranuni(seed2); output; end; options nocenter; proc print; run; 例中，SEED1和SEED2的初值相同，但X1和X2的第一个观测值却不相同。因为，SEED2的值在这里不起作用，X2的第一个值并不是由SEED2产生的，而是产生第一个X1后一个新种子的结果。所以，这里产生的两个随机数变量属于同一个随机数流。 ;例3.39 使用CALL子程序产生两个独立的种子流和随机数流。 data RV; retain seed3 seed4 135279821; do I=1 to 5; call ranuni(seed3, X3); call ranuni(seed4, X4); output; end; proc print; run; 例中，产生两个独立的种子流和随机数流。因为SEED3=SEE