统计计算第一次作业

第一次作业
1. 以下数据是调查15个学生是否学过四门课程（记为A、B、C、D）的资料，0表示没有学过，1表示已学过。

0 0 0 1 m 1 1 0 0 f
0 0 1 0 m 0 0 0 1 f
0 1 1 0 f 1 1 0 0 m
1 0 1 1 m 0 0 0 0 m
0 1 1 0 f 1 1 1 1 f
1 0 0 1 f 0 0 1 1 m
1 0 1 1 m 1 0 1 1 f
1 0 1 0 f
(1)对这组资料创建一SAS数据集，然后按男女分别生成两个SAS数据集。

(2)建立一新变量NEW：先把是否学过4门课程的资料看成4 个元素的向量，并作为一个4位二进制数，然后转换为十进制数
作为NEW的值。

比如：（1011）可化为1*23+0*22+1*21+1*20=11。

然后创建一个包含新变量和四门课得分的SAS数据集；并按NEW值从小到大的次序排列后输出。

(3)请统计一下15个学生中学过A、B、C、D四门课的各有多少人，并把结果存为SAS数据集后并打印输出。

程序:
data m f;/*将男女生的数据分别输入到数据集m、f*/
input A B C D S$@@;
if S='m'then output m;
if S='f'then output f;
cards;
0 0 0 1 m 1 1 0 0 f
0 0 1 0 m 0 0 0 1 f
0 1 1 0 f 1 1 0 0 m
1 0 1 1 m 0 0 0 0 m
0 1 1 0 f 1 1 1 1 f
1 0 0 1 f 0 0 1 1 m
1 0 1 1 m 1 0 1 1 f
1 0 1 0 f
;
run;
data a;/*产生数据集放NEW变量*/
set m f;
NEW=A*2**3+B*2**2+C*2**1+D*2**0;
proc sort data=a;by NEW;/*按NEW变量排序*/
proc print;
run;
data b;/*PA,PB,PC,PD分别表示分别通过A,B,C,D四门课的累积人数*/
set m f;
PA+A;PB+B;PC+C;PD+D;
drop A B C D;
proc print;
run;
(2)
(3)
则A B C D 各科通过人数分别为：8，5，9，8个人
2. 试用概率函数计算正态分布N（3，9）的分布函数F（X）。

X F（X）
0 0.158655254
0.1 0.166855348
0.2 0.175323945
0.3 0.184060125
0.4 0.193062337
0.5 0.202328381
0.6 0.211855399
0.7 0.221639864
0.8 0.231677575
0.9 0.241963652
1 0.252492538
1.1 0.263257995
1.2 0.274253118
1.3 0.285470336
1.4 0.296901429
1.5 0.308537539
1.6 0.320369191
1.7 0.332386313
1.8 0.344578258
程序：
data a;
do x=0to3by0.1; y=probnorm((x-3)/3); file print;put y; end;
run;
3.试用分位数函数计算F
程序：
data a;
do ndf=1to5;
do ddf=1to10;
do a=0.05,0.025,0.005;
F=finv(1-a,ndf,ddf);
file print;put F;
end;end;end;
run;
4. 以下四行数据是四次记录下来的数据，每次从4 个厂家（用A，B，C，D代表）生产的同类食品中，每个厂家抽取三包测量每包的数量。

每一行开头三个数值是厂家A生产，以后依次是B、C、D生产的。

72 74 69 61 61 65 62 65 70 85 76 61
67 52 62 60 55 59 64 65 64 67 72 60
57 66 72 72 43 43 63 66 72 56 75 92
57 56 78 60 63 58 61 79 68 73 86 71
（1）请按数据形式输入后生成外部数据文件DB25.TXT.
（2）用DATA步调入外部文件DB25.TXT(用INFILE和INPUT语句)生成包括变量:TIME（次数）,N(包的序
号),FACTORY（厂家）和NUMBER（数量）的SAS数据集
（3）请按以下要求输出(2)产生的SAS数据集:
1.只包含变量TIME，FACTORY和NUMBER；
2.对每个厂家的资料分别输出报告,要求每包食品的数量按小到大的次序排好.
（4）计算每个厂家的12包食品的平均数量，标准差，极差和总和等描述统计量；然后存贮为一个SAS数据集（包含4 个观测）。

程序：
data test;/*输入数据到test数据集*/
do time=1to4;
do factory='A','B','C','D';
do n=1to3;
input number@@;
output;
end;
end;
end;
cards;
72 74 69 61 61 65 62 65 70 85 76 61
67 52 62 60 55 59 64 65 64 67 72 60
57 66 72 72 43 43 63 66 72 56 75 92
57 56 78 60 63 58 61 79 68 73 86 71
;
data_null_;/*将test数据集导出到外部文件DB25.TXT*/
set test;
file'D:\DB25.TXT';
put number;
data test2;/*导入外部文件的数据到test2数据集*/
infile'D:\DB25.TXT';
do time=1to4;
do factory='A','B','C','D';
do n=1to3;
input number @;
output;
end;
end;
end;
proc print;var time factory number; run;/*打印test2数据集中的数据，包含变量 time factory number*/
data a b c d;/*将A,B,C,D厂的数据分别输入到数据集a,b,c,d*/
set test2;
if factory='A'then output a;
if factory='B'then output b;
if factory='C'then output c;
if factory='D'then output d;
proc sort data=a;by number;proc print;
proc sort data=b;by number;proc print;
proc sort data=c;by number;proc print;
proc sort data=d;by number;proc print;
run;
proc summary data=a;var number;/*计算描述统计量*/
output out=aa mean=mean std=std range=range sum=sum; run;
proc summary data=b;var number;
output out=bb mean=mean std=std range=range sum=sum; run;
proc summary data=c;var number;
output out=cc mean=mean std=std range=range sum=sum; run;
proc summary data=d;var number;
output out=dd mean=mean std=std range=range sum=sum; run;
data factory;/*建立工厂名数据集*/
do factory='A','B','C','D';output;end;run;
data sta;/*将各个工厂的描述信息汇总到数据集sta*/
set aa bb cc dd;
data last;
merge factory sta;
proc print;run;
(3)1.
(3)2.A厂：
B厂
C厂
D厂
（4）
5、书上第28页，1.4题
说明：前四题给出SAS程序，并打印出来；第五题要求手写，直接写在打印出来的纸上，不要用本子写。