实验报告参考答案

实验报告

实验项目名称SAS语言入门实验室机房四

所属课程名称统计软件

实验类型综合型

实验日期2013年3月

班级

学号

姓名

成绩

【实验方案设计】

实验（1）：计算标准正态分布在x=-3,-2,-1,0,1,2,3时的分布函数F(x)和密度函数f(x)的值；通过该实验请描述SAS编程的流程，此外另举两个常用分布与分布函数有关的5个函数，制表给出。

1

（1）

Data ex1_1;

Input x@@;

fa=CDF('norm',X);

fb=pdf('norm',x);

Datalines;

-3 -2 -1 0 1 2 3

;

Run;

proc print;

run;

（2）

data ex1_2; input x @@;

fx1=CDF('exponential',x); fx2=CDF('logistic',x); fc=PDF('norm',x); fd=PMF('norm',x); fe=LOGPDF('norm',x); ff=LOGPMF('norm',x); datalines ; -3 -2 -1 0 1 2 3 ; run ;

proc print data =ex1_2; run ;

实验（2）：设随机变量T ～t(n),计算t 分布的分位数)(n t p ，其中n=1,5,10,20;p=0.1,0.9,0.95和0.975;通过该实验简述分位数函数与分布函数之间的关系，并在SAS 中举例验证说明。 2. （1）

data exp2_1; input p@@; t1=tinv(p,1); t5=tinv(p,5); t10=tinv(p,10); t20=tinv(p,20); datalines ;

0.1 0.9 0.95 0.975 ; run ;

proc print data =exp2_1; run ;

（2）以正态分布为例

data exp2_2;

input x y;

a=probit(x);

b=Probnorm(y);

datalines;

0.975 1.95996

;

run;

proc print data=exp2_2;

run;

可见，分布函数与分位数函数是反函数关系

实验（3）：当日期值date=’14JUL2007’d时，试问这一天是星期几？此外如何该日期不采用输出格式，输出结果为多少？列举日期的四种格式以及SAS中如何实现？

3.

（1）

data exp3_1;

x='14JUL2007'd;

y=weekday('14JUL2007'd);

put y=;

put x=;

format x yymmdd.6;

run;

proc print;

run;

（2）

data exp3_2;

x='14JUL2007'd;

y=weekday('14JUL2007'd);

put y=;

put x=;

run;

proc print;

run;

(3)

日期格式SAS

10/21/50 MMDDYY8.

10/21/1950 MMDDYY10. 21OCT50 DATE7.

21OCT1950 DATA9. October 21, 1950 WORDDATE.

实验（4）：SAS数据步给定一些变量赋值如下：

code=’VLC’; size=’M’; v1=2; v2=3; v3=7; v4=.;

请验证下列SAS语句的输出结果。

a.（v1+v2-v3）/3

b. v3-v2/v1

c. v1*v2-v3

d. v2*v3/v1

e. v1**2+v2**2

f. code=’VLC’

g. code=’VLC’|size=’M’ h. code=’VLC’ & size=’M’

i. code=’VLC’ & v4^=. j. (v3=.)+(v2=3)

k. v1+v2+v3 ^= 12 i. code=’VLC’ | (size=’M’ & v1=3) m. 3

通过该实验给出SAS中逻辑值与数值之间的变换规则。

4.

data exp4_1;

code='vlc';

size='m';

v1=2;

v2=3;

v3=7;

v4=.;

a=(v1+v2-v3)/3;

b=v3-v2/v1;

c=v1*v2-v3;

d=v2*v3/v1;

e=v1**2+v2**2;

f=(code='vlc');

g=(code='vlc' & size='m');

k=(v1+v2+v3^=12);

j=(v3=.)+(v2=3);

m=(3

put a= b= c= d= e= f= g= k= j= m=;

run;

proc print;

run;

3. 数值与逻辑值的转换规则

逻辑判断语句输出的值只有0,1两种。当去掉比较运算符及逻辑运算符时，仅仅只为数值运算，此时输出的是数值计算结果。当数值运算式中添上比较与逻辑运算符，便能变成逻辑语句，输出逻辑运算值。

实验（5）：给定某航空月度乘客数据集air，请计算每月乘客数较上月增加的百分比；计算逐月乘客数较去年同期增加的百分比。通过该实验掌握SAS如何计算同比和环比，并给出一阶差分和滞后一阶两个概念的区别。此外，请将air数据集放入自己的逻辑库中。

5.

libname data 'd:\dataset\dst';

data shiyan5;

set data.air_2;

pct1=dif(air)/lag(air);

pct2=dif12(air)/lag12(air);

run;

proc print data=shiyan5;

id date;

var air pct1 pct2;

run;

一阶差分是当前期与前一期的值之差，滞后一阶是前一期的值。【小结】

附件：

实验报告说明

1．实验项目名称：要用最简练的语言反映实验的内容。

2．实验类型：一般需说明是验证型实验还是设计型实验、综合型实验或其他实验。

3．实验目的与要求：目的要明确，要抓住重点。

4．实验原理：简要说明本实验项目所涉及的理论知识。

5．实验环境：实验用的软硬件环境（配置）。

6．实验方案设计（思路、步骤和方法等）：这是实验报告极其重要的内容，概括整个实验过程。

对于验证型实验，要写明依据何种原理、何仲操作方法进行实验，并写明需要经过哪几个步骤。

对于设计型和综合型实验，在上述内容基础上还应该画出流程图、设计思路和设计方法，再配以相应的文字说明。

7．实验过程（实验中涉及的记录、数据、分析）：写明具体上述实验方案的具体实施，包括实验过程中的记录、数据和相应的分析。

8．结论（结果）：即根据实验过程中所见到的现象和测得的数据，做出结论。

9．小结：对本次实验的心得体会、思考和建议。

10．指导教师评语及成绩：指导教师依据学生的实际报告内容，用简练语言给出本次实验报告的评价和价值。