常用函数简介

关于常用函数简介SAS函数SAS系统提供的SAS函数比一般高级语言的标准函数多得多，尤其是概率函数、分位数函数、样本统计函数、随机数函数。

如此丰富的SAS函数,对用户统计分析计算以及编写SAS程序带来极大的方便。第2页,共49页，2024年2月25日，星期天SAS函数定义SAS函数是一个子程序，对0个或几个参数（自变量）返回一个结果值。SAS函数的形式：函数名（参数1，参数2，…）例如：

MEAN(x,y,z,a)MIN(SUM(x,y,z),1000)第3页,共49页，2024年2月25日，星期天SAS函数中的参数1）参数类型：变量名；常数；函数；表达式。例1函数及其参数类型举例。第4页,共49页，2024年2月25日，星期天2）参数个数：有些函数不需要自变量，有些只需要一个自变量，有些函数有多个自变量。3）参数表示法：当函数有多个参数时，必须用逗号分隔开。4）参数名缩写方法：函数名(OF变量名1…变量名n)

第5页,共49页，2024年2月25日，星期天参数的表示法举例以逗号(，)分隔参数,如sum(x1,x2,x3,x4,x5,x6,x7,x8)去掉逗号(，)，最前面加上“of”,如sum(ofx1x2x3x4x5x6x7x8)去掉逗号(，)，最前面加上“of”,变量简写：sum(ofx1-x8)第6页,共49页，2024年2月25日，星期天例2正确与错误的表示方法例3数组元素自变量的缩写方法arrayy{10}y1-y10;/*定义数组y{10}*/x=sum(ofy{*});z=sum(ofy1-y10);例中，两个SUM函数的表示法是等价的。第7页,共49页，2024年2月25日，星期天函数结果通常函数的结果由其自变量的属性决定：自变量是字符时结果变量为字符；自变量是数值时结果变量为数值。PUT函数是一个例外，它不管自变量是什么类型，其结果总是字符值。第8页,共49页，2024年2月25日，星期天用PUT语句在LOG窗口显示函数值例4显示概率值和分位数。data;Y=probnorm(1.96);/*标准正态分布Z小于1.96的概率*/putY;q1=tinv(.95,2);/*自由度为2的t分布的0.95分位数*/q2=tinv(.95,2,3);

/*自由为2，非中心参数为3的t分布的0.95分位数*/putq1=q2=;run;0.9750021049q1=2.9199855804q2=13.894376071第9页,共49页，2024年2月25日，星期天利用HELP查找SAS系统提供的各种函数索引中键入：SASFunctionHelp下拉菜单SAS帮助与文件第10页,共49页，2024年2月25日，星期天SAS函数种类在SAS系统9.1版本,SAS函数分为二十五种类型,分别是算术函数、截取函数、数学函数、概率和密度函数、分位数函数、非中心函数、样本统计函数、随机数函数、财政金融函数、日期和时间函数、换算函数和特殊函数等。在SAS系统6.12版本，SAS函数分为十七种类型共有178个,SAS系统8.1版本，SAS函数分为二十二种类型共有319个函数。第11页,共49页，2024年2月25日，星期天第12页,共49页，2024年2月25日，星期天1.日期时间函数简介第13页,共49页，2024年2月25日，星期天第14页,共49页，2024年2月25日，星期天例1.1：计算两个日期之间相隔的天数data

_null_;sdate='01jan2012'd;edate='01jan2013'd;actual=datdif(sdate,edate,'act/act');days360=datdif(sdate,edate,'30/360');putactual=days360=;run;data

_null_;actual=datdif('01jan2002'd,'01jan2003'd,'act/act');days360=datdif('01jan2002'd,'01jan2003'd,'30/360');putactual=days360=;run;第15页,共49页，2024年2月25日，星期天例1.2:以日为单位计算当前日期的天数。data;x=date();y=today();formatxyyymmdd10.;putx=y=;run;第16页,共49页，2024年2月25日，星期天例1.3:将时间值换算为以秒计数的计算方法。data;hms=hms(21,50,51);puthms=;x=21*60*60+50*60+51;putx=;mdy=mdy(08,18,2001);putmdy=;formatmdyyymmdd10.;run;第17页,共49页，2024年2月25日，星期天练习1计算以下数据表中日期值所在的年季月以及处于某月的第几天、某周的第几天。当日期值date=’14JUL2007’d时,试问这一天是星期几?

Date1991-01-021991-01-031991-01-041991-01-051991-01-071991-01-081991-01-091991-01-101991-01-111991-01-12第18页,共49页，2024年2月25日，星期天解答：dataa;InputdateYYMMDD10.;Cards;1991-01-021991-01-031991-01-041991-01-051991-01-071991-01-081991-01-091991-01-101991-01-111991-01-12;Procprint;FormatdateYYMMDD10.;run;Datab;Seta;year=year(date);qtr=qtr(date);month=month(date);day=day(date);Weekday=weekday(date);procprintnoobs;formatdateYYMMDD10.;Run;第19页,共49页，2024年2月25日，星期天解答：data

_null_;weekday=weekday('12JUL2007'd);putweekday=;run;第20页,共49页，2024年2月25日，星期天2.数学函数ABS(x)：求x的绝对值Mod(x,y)：求x/y的余数MAX(x1,x2,x3,…xn）：求最大值MIN(x1,x2,x3,…xn）：求最小值SQRT(x)：求x的算术平方根LOG(x)：求x的自然对数EXP(x)：求指数函数exSIN(x)、COS(x)、TAN(x)：正弦、余弦、正切函数…….第21页,共49页，2024年2月25日，星期天例2：求-314的绝对值，31除以4的商的余数以及-314、-87、0、1127这4个数的最大值data

_null_;x=abs(-314);y=mod(31,4);z=max(-314,-87,0,1127);putx=y=z=;run;第22页,共49页，2024年2月25日，星期天练习2ex2.1计算0、1.57、3.14对应的正弦函数值，并将结果显示在输出窗口。ex2.2计算e的0次幂,并计算1的自然对数。第23页,共49页，2024年2月25日，星期天解答：Data;a=sin(0);b=sin(1.57);c=sin(3.14);Fileprint;Puta=6.2b=6.2c=6.2;Run;Data;a=sin(0);b=sin(1.57);c=sin(3.14);procprint;formata6.2b6.2c6.2;Run;Data;X=exp(0);Y=log(1);Putxy;Run;第24页,共49页，2024年2月25日，星期天3.数组函数第25页,共49页，2024年2月25日，星期天数组函数应用举例data;arraymult{5,10,2}mult1-mult100;a=DIM1(MULT);b=DIM2(MULT);c=DIM3(MULT);putabc;run;例3：产生一个5×10×2的三维数组，并提取各维度上元素的个数。第26页,共49页，2024年2月25日，星期天练习3产生一个行变量名为2:6，列变量名为4:13的二维数组，输出各维度上的上界和下界。第27页,共49页，2024年2月25日，星期天解答：data;arraymult{2:6,4:13}mult1-mult50;a=HBOUND(MULT);b=HBOUND2(MULT);c=lBOUND(MULT);d=lBOUND(MULT,2);putabcd;run;第28页,共49页，2024年2月25日，星期天4.样本统计量函数第29页,共49页，2024年2月25日，星期天简单统计量函数应用举例例4：请求出2、4、6、3、1这五个数的最大值、最小值、算术平均值、标准差、变异系数。datatest4;inputx1-x5;cards;24631;run;datasimple_statistics;settest1;max_x=max(ofx1-x5);min_x=min(ofx1-x5);mean_x=mean(ofx1-x5);std_x=std(ofx1-x5);cv_x=cv(ofx1-x5);run;odshtml;procprintdata=simple_statistics;varmax_xmin_xmean_xstd_xcv_x;run;odshtmlclose;第30页,共49页，2024年2月25日，星期天用简便的方法实现例1中目的datasimple_statistics;max_x=max(2,4,6,3,1);min_x=min(2,4,6,3,1);mean_x=mean(2,4,6,3,1);std_x=std(2,4,6,3,1);cv_x=cv(2,4,6,3,1);run;odshtml;procprintdata=simple_statistics;varmax_xmin_xmean_xstd_xcv_x;run;odshtmlclose;第31页,共49页，2024年2月25日，星期天练习4某农村地区1998年随机抽查测量10名13岁女孩的身高，资料如下，试求全距、均值、标准差、标准误、变异系数。145.1143.3153.2152.3150.6148.1144.0137.1158.2146.7第32页,共49页，2024年2月25日，星期天解答：dataaa;inputx1-x10;cards;145.1143.3153.2152.3150.6148.1144.0137.1158.2146.7;run;databb;setaa;a=range(ofx1-x10);b=mean(ofx1-x10);c=std(ofx1-x10);d=cv(ofx1-x10);file

print;putabcd;run;第33页,共49页，2024年2月25日，星期天5.截取函数简介第34页,共49页，2024年2月25日，星期天SAS截取函数应用举例例5：请把STD=1.92354和CV=0.601106用3个不同的截取函数(ceil、int、round)进行运算，并分析这些函数的作用。datatest5;STD=1.92354;CV=0.601106;a1=ceil(STD);b1=ceil(CV);a2=int(STD);b2=int(CV);a3=round(STD,0.001);b3=round(CV,0.001);odshtml;procprintdata=test2;Var

a1-a3b1-b3;run;odshtmlclose;第35页,共49页，2024年2月25日，星期天练习5求小于等于31.77的最小整数；将1984.314保留2位小数。第36页,共49页，2024年2月25日，星期天解答：data;x=floor(31.07);y=round(1984.314,0.01);file

print;putxy;run;第37页,共49页，2024年2月25日，星期天6.SAS概率函数第38页,共49页，2024年2月25日，星期天SAS概率函数应用举例例6：试求出服从标准正态分布的随机变量z取值分别小于等于-2.576、-1.960、-1.645、-1.282、0、1.282、1.645、1.960、2.576所对应的概率。datatest3;a1=round(probnorm(-2.576),0.001);a2=round(probnorm(-1.960),0.001);a3=round(probnorm(-1.645),0.001);a4=round(probnorm(-1.282),0.001);a5=round(probnorm(0),0.001);a6=round(probnorm(1.282),0.001);a7=round(probnorm(1.645),0.001);a8=round(probnorm(1.960),0.001);a9=round(probnorm(2.576),0.001);Formata1-a95.3;odshtml;procprintdata=test3;Var

a1-a9;run;odshtmlclose;第39页,共49页，2024年2月25日，星期天练习6设随机变量X服从自由度为3的中心卡方分布，计算X≥7.94的概率。第40页,共49页，2024年2月25日，星期天解答data;p=1-probchi(7.94,3);putp=;run;第41页,共49页，2024年2月25日，星期天7.SAS概率分布下分位数函数第42页,共49页，2024年2月25日，星期天例7：试求t临界值表。optionsLS=95PS=55;dataabc;arrayT(50,9);arraydf(50)df1-df50;g=1;k=1;dod=1to40,50,60,70,80,90,100,200,

500,1000,100000000;df(k)=d;k+1;w=1;doalpha=0.50,0.20,0.10,0.05,0.02,

0.01,0.005,0.002,0.001;p=1-alpha/2;b=tinv(p,d);b=round(b,0.0001);T(g,w)=b;

w=w+1;Output;

end;

g=g+1;

end;Title;

title't临界值表';odshtmlstyle=journal;fileprint;put#1@3‘dfalpha:’

@14‘0.50’@23‘0.20’

@31'0.10'@39‘0.05’

@48‘0.02’@57‘0.01’

@66'0.005'@76'0.002'@86'0.001';put#2'______________________';doL=1to50;c=L+2;

df(50)='9999';put#c@1df(L)4.0@13T(L,1)6.4@21T(L,2)6.4@29T(L,3)6.4@37T(L,4)7.4@46T(L,5)7.4@55T(L,6)7.4@64T(L,7)8.4@74T(L,8)8.4@84T(L,9)8.4;end;run;odshtmlclose;footnote;footnote'note:df=9999standsfordf=100000000';第43页,共49页，2024年2月25日，星期天练习7设随机变量Ｔ~t(df),计算t分布的分位数tp(df),其中自由度和对应概率值分别为：n=1,5,10,20;p=0.10,0.90,0.95,