实验一SAS系统及数据集创建

1、实验一 SAS 系统及数据集创建一、实验目的:了解 SAS 逻辑库组成,创建。学习 SAS 数据集(数据文件,数 据视窗的命名,创建。学习数据集导入导出,变量定义,变量选择,观测选择 等操作,缺失值插值。二、 实验内容(1用户用 SAS 数据步(DA TA STEP创建一个数据集的方法, DA TA 程序步的主要步骤 : Data 所要创建的数据集名 ;Infile 读取的外部文件名 < FIRSTOBS=开始读入的行 >< OBS=结束行 > ;Input 变量 1 读入模式 变量 2 读入模式 ;Run ;(2 确定变量的读入模式共有四种:column 模式Inpu

2、t 变量 1 < $> 开始列 <-结束列 ><.小数位 > 变量 2 ;例 1:data testd1;infile 'e:dmydirtestdata.dat' firstobs=2 obs=4;input code $1-8 sex $9 bdata $10-16 income 17-23 .2; run;上面的程序的数据步在临时库 WORK 中创建了 SAS 数据集 work.testd1, 而过程步 PRINT 则显示该数据集。formatted 模式格式一:Input <指针控制 > 变量 1 输入格式 <|&g

3、t; ;格式二:Input <指针控制 > (变量表 (输入格式表 <|> ;例 2:Data testd2 ;Infile e:dmydirtestdata.dat ;Input code $ 1-8 sex $ 9 ;Input bdate $ 10-16 income $17-23 .2 ;Proc printdata=testd2 ;Run; list 模式格式一:Input <指针控制 > 变量 1 <$><&> <|> ;格式二:Input <指针控制 > 变量 1 <:|&|

4、> <输入格式 > <|> ;例 3:Data testd5 ;Infile e:dmydirtestdata2.txt dsd;Input name & $15. age address $40. pc ;Proc print data=testd5 ;Run ; 例 4:data A600001(label="邯郸钢铁 "input date yymmdd10. oppr hipr lopr clpr vol sums name $ hstocd $; format date yymmdd10. ; /*设定变量的输出格式 */la

5、bel Date="日期 |Date" /*给变量加标签 */label Oppr="开盘价 |Open Price"cards; /*以下是数据表 */2001-1-16 8.08 8.4 8.05 8.35 45571 37683964 邯郸钢铁 6000012001-1-17 8.39 8.5 8.28 8.29 26345 22076348 邯郸钢铁 6000012001-1-18 8.32 8.4 8.21 8.26 17854 14844434 邯郸钢铁 600001 4、增加变量、选择变量和观测创建新数据集例 5:Data class2 ;

6、Input id test1-test5 ;average=mean(test1,test2,test3,test4,test5;test4=test1+test2+test3+test4+test5;Cards ;980801 100 100 100 100 100980802 90 100 90 100 90980803 81 82 83 84 85Proc print data=class2 ;Run ;Input id test1-test5 ; 输入数据average=mean(test1,test2,test3,test4,test5; 求 test1到 test5的平均值test

7、4=test1+test2+test3+test4+test5; 例 6:Data class3 (drop=test1 test2 test3 test4 test5 ;Set class2 ;Proc print data=class3 ;Run;创建 class3 数据集,数据来源是 class2 中剔除 test1 test2 test3 test4 test5 例 7:Data class4 ;Set class2 ;if average>=80 and (100=test1 or 100=test2 or 100=test3 or 100=test4 or 100=test5

8、 ;Proc printdata=class4 ;Run ;创建数据集 class3, 数据来源是 class2中均值 >=80, 且 test 1test 5中有一值 为 100. 三、习题1、下列为 1975-1980年夏威夷岛莫纳罗亚火山每月释放 CO2的数据,请将它生 成数据集 H1,并剔除释放量小于 330的观测,生成新的数据集 H2。330.45 330.97 331.64 332.87 333.61 333.55331.90 330.05 328.58 328.31 329.41 330.63 331.63 332.46 333.36 334.45 334.82 334.3

9、2 333.05 330.87 329.24 328.87 330.18 331.50 332.81 333.23 334.55 335.82 336.44 335.99 334.65 332.41 331.32 330.73 332.05 333.53 334.66 335.07 336.33 337.39 337.65 337.57 336.25 334.39 332.44 332.25 333.59 334.76 335.89 336.44 337.63 338.54 339.06 338.95 337.41 335.71 333.68 333.69 335.05 336.53 337.

10、81 338.16 339.88 340.57 341.19 340.87 339.25 337.19 335.49 336.63 337.74 338.36 Data H1 ;Input id test1-test6;Cards ;330.45 330.97 331.64 332.87 333.61 333.55 331.90 330.05 328.58 328.31 329.41 330.63 331.63 332.46 333.36 334.45 334.82 334.32 333.05 330.87 329.24 328.87 330.18 331.50 332.81 333.23 3

11、34.55 335.82 336.44 335.99 334.65 332.41 331.32 330.73 332.05 333.53 334.66 335.07 336.33 337.39 337.65 337.57 336.25 334.39 332.44 332.25 333.59 334.76 335.89 336.44 337.63 338.54 339.06 338.95 337.41 335.71 333.68 333.69 335.05 336.53 337.81 338.16 339.88 340.57 341.19 340.87 339.25 337.19 335.49

12、336.63 337.74 338.36 Proc print data =H1 ;Run ;运行结果:co2time330.451975/1/1330.971975/2/1331.641975/3/1332.871975/4/1333.611975/5/1333.551975/6/1331.91975/7/1330.051975/8/1330.631975/12/1331.631976/1/1332.461976/2/1333.361976/3/1334.451976/4/1334.821976/5/1334.321976/6/1333.05 330.87 330.18 331.5 332.

13、81 333.23 334.55 335.82 336.44 335.99 334.65 332.41 331.32 330.73 332.05 333.53 334.66 335.07 336.33 337.39 337.65 337.57 336.25 334.39 332.44 332.25 333.59 334.76 335.89 336.44 337.63 338.54 339.06 338.95 337.41 335.71 333.68 333.69 335.05 336.53 337.81 338.16 1976/7/1 1976/8/1 1976/11/1 1976/12/1

14、1977/1/1 1977/2/1 1977/3/1 1977/4/1 1977/5/1 1977/6/1 1977/7/1 1977/8/1 1977/9/1 1977/10/1 1977/11/1 1977/12/1 1978/1/1 1978/2/1 1978/3/1 1978/4/1 1978/5/1 1978/6/1 1978/7/1 1978/8/1 1978/9/1 1978/10/1 1978/11/1 1978/12/1 1979/1/1 1979/2/1 1979/3/1 1979/4/1 1979/5/1 1979/6/1 1979/7/1 1979/8/1 1979/9

15、/1 1979/10/1 1979/11/1 1979/12/1 1980/1/1 1980/2/1 339.88 340.57 341.19 340.87 339.25 337.19 335.49 336.63 337.74 338.36 1980/3/1 1980/4/1 1980/5/1 1980/6/1 1980/7/1 1980/8/1 1980/9/1 1980/10/1 1980/11/1 1980/12/1 2、1945 年-1950 年费城月度降雨量数据如下（单位：mm） 。 生成数据集 H3。 330.45 330.97 331.64 332.87 333.61 333.5

16、5 331.90 330.05 328.58 328.31 329.41 330.63 331.63 332.46 333.36 334.45 334.82 334.32 333.05 330.87 329.24 328.87 330.18 331.50 332.81 333.23 334.55 335.82 336.44 335.99 334.65 332.41 331.32 330.73 332.05 333.53 334.66 335.07 336.33 337.39 337.65 337.57 336.25 334.39 332.44 332.25 333.59 334.76 335.

17、89 336.44 337.63 338.54 339.06 338.95 337.41 335.71 333.68 333.69 335.05 336.53 337.81 338.16 339.88 340.57 341.19 340.87 339.25 337.19 335.49 336.63 337.74 338.36 data h3; input hight; time=intnx('month','01jan1945'd,_n_-1; format time mmddyy10; cards; 330.45 330.97 331.64 332.87 33

18、3.61 333.55 331.90 330.05 328.58 328.31 329.41 330.63 331.63 332.46 333.36 334.45 334.82 334.32 333.05 330.87 329.24 328.87 330.18 331.50 332.81 333.23 334.55 335.82 336.44 335.99 334.65 332.41 331.32 330.73 332.05 333.53 334.66 335.07 336.33 337.39 337.65 337.57 336.25 334.39 332.44 332.25 333.59 3

19、34.76 335.89 336.44 337.63 338.54 339.06 338.95 337.41 335.71 333.68 333.69 335.05 336.53 337.81 338.16 339.88 340.57 341.19 340.87 339.25 337.19 335.49 336.63 337.74 338.36 ; proc print data=h3; run; hight 330.45 time -5478 330.97 331.64 332.87 333.61 333.55 331.9 330.05 328.58 328.31 329.41 330.63 331.63 332.46 333.36 334.45 334.82 334.32 333.05 330.87 329.24 328.87 330.18 331.5 332.81 333.23 334.55 335.82 336.44 335.99 334.65 332.41 331.32 330.73 332.05 333.53 334.66 335.07 336.33 33