SAS程序编辑.doc

《SAS程序编辑.doc》由会员分享，可在线阅读，更多相关《SAS程序编辑.doc（7页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、1. 已知正常人乙酰胆碱脂酶的平均值为1.44单位，现测得13例慢性气管炎患者的乙酰胆碱脂酶分别为：1.50，2.19，2.32，2.41，2.11，2.54，2.20，2.36，1.42，2.17，1.84，1.96，2.39 问：慢性气管炎患者与正常人的平均乙酰胆碱脂酶之间的差别有无显著性意义?答案：data lian1;input x;d=x-1.44;cards;1.5 2.19 2.32 2.41 2.11 2.54 2.2 2.36 1.42 2.17 1.84 1.96 2.39;proc univariate normal;var d;proc means mean stder

2、r t prt;var d;run;2.某研究室用甲乙两种血清学方法检查401例确诊的鼻癌患者结果如下表两种血清学方法检查结果甲法乙 法合计+3611037183139合计36941410推断两种血清学方法有无差别。答案：data lian2;do r=1 to 2;do c=1 to 2;input f;output;end;end;cards;361 10 8 31;proc freq;tables r*c/chisq agree;weight f;run ; 3.某医院为了探索导致手术切口感染的原因，收集了305例手术的情况，手术时间小于等于5小时有242例，其中13例感染，手术时间大于

3、5小时有63例，其中7例感染，问两种手术感染率的差别有无统计意义？4.从八窝大白鼠中分别选出同性别、体重相近的两只，喂以水解蛋白和酪蛋白的饲料，四周后测定其体重增加情况，结果如下：窝编号 1 2 3 4 5 6 7 8酪蛋白组 82 86 74 78 82 76 73 90水解蛋白组 12 28 29 28 24 38 21 37比较两种饲料对大白鼠体重增加量的差别。答案：data lian3;input x1 x2;d=x1-x2;cards;82 12 86 28 74 29 78 28 82 24 76 38 73 21 90 37;proc univariate normal;var

4、d;proc means mean stderr t prt;var d;run;5.为研究血清唾液酸含量对慢性胃部疾患的诊断价值，测定了慢性胃炎和溃疡病患者的血清唾液酸含量，结果如下：慢性胃炎患者 45.68 47.20 47.50 44.39 47.80 48.11 46.45 47.01 45.44 溃疡病患者 44.04 45.12 44.07 46.77 46.28 48.56 44.87 48.58 47.32试对两组血清唾液酸含量进行比较。答案：data lian_5;input x g;cards;45.68 1 47.20 1 47.50 1 44.39 1 47.80 1

5、48.11 1 46.45 1 47.01 1 45.44 144.04 2 45.12 2 44.07 2 46.77 2 46.28 2 48.56 2 44.87 2 48.58 2 47.32 2;proc univariate normal;var x; by g;run;proc ttest; class g; var x;run;6.为研究血清唾液酸含量对慢性胃部疾患的诊断价值，测定了正常人、慢性胃炎和溃疡病患者的血清唾液酸含量，结果如下：正常人 43.1 41.46 42.35 41.42 42.01 43.12 42.05 42.20 42.24 41.16慢性胃炎患者 45

6、.68 47.20 47.50 44.39 47.80 48.11 46.45 47.01 45.44 溃疡病患者 44.04 45.12 44.07 46.77 46.28 48.56 44.87 48.58 47.32试对三组血清唾液酸含量进行比较。答案data lian13;do g=1 to 3;input n;do b=1 to n;input x;output;end;end;cards;1043.10 41.46 42.35 41.42 42.01 43.12 42.05 42.20 42.24 41.16945.68 47.20 47.50 44.39 47.80 48.11

7、46.45 47.01 45.44944.04 45.12 44.07 46.77 46.28 48.56 44.87 48.58 47.32;proc univariate normal;var x; by g;proc anova;class g; model x=g;means g/snk hovtest;means g;proc npar1way wilcoxon;var x; class g;run;7. 某地区不同民族包虫感染率调查结果如下，试比较不同民族包虫感染率有无差别。民族阳性数阴性数汉族46317624回族471273哈萨克族53819维吾尔族33829答案：data li

8、an_7;do r= 1 to 4;do c= 1 to 2;input f;output;end;end;cards;463 17624 47 1273 53 819 33 829;proc freq;tables r*c/chisq expected exact;weight f;run;8.某研究者采用随机区组设计进行实验，比较三种抗癌药物对小白鼠肉瘤抑瘤效果，先将15只染有肉瘤小白鼠按体重大小配成5个区组，每个区组内3只小白鼠随机接受三种抗癌药物，以肉瘤的重量为指标，试验结果见表。问三种不同的药物的抑瘤效果有无差别？ 不同药物作用后小白鼠肉瘤重量（g）区组A药B药C药10.820.65

9、0.5120.730.540.2330.430.340.2840.410.210.3150.680.430.24答案：data lian_14;do g= 1 to 3;do b= 1 to 5;input x;output;end;end;cards;A法 B法 + + 142 38 12 680.82 0.65 0.510.73 0.54 0.230.43 0.34 0.280.41 0.21 0.310.68 0.43 0.24;proc univariate normal;var x; by g;proc anova;class g ;model x= g ;means g/snk h

10、ovtest;means g;run;9. 201例就诊者有无右下腹压痛与手术证实是否患了阑尾炎之关系-有无右下 例数 -腹压痛 手术证实情况 阑尾炎 非阑尾炎 合计-有 93 27 120无 15 66 81 -答案：data lian_9;do r= 1 to 2;do c= 1 to 2;input f;output;end; end;cards;93 27 15 66;proc freq;tables r*c/chisq expected exact;weight f;run;10. .有260份血清样品，每份样品一分为二，用两种不同的免疫学检测方法检验类风湿因子，其结果如表所示。试问

11、这两种免疫学检验的结果是否相互独立？两种检测方法的检出率是否相同。 表 两种血清免疫学检测结果答案：data lian_10;do r = 1 to 2;do c = 1 to 2;input f;output;end;end;cards;142 38 12 68;proc freq;tables r*c/chisq agree;weight f;run;11.某人测量了10个家庭兄弟和姐妹的平均身长，试根据本资料的数据对兄弟与姐妹间身长进行相关分析。家庭编号 1 2 3 4 5 6 7 8 9 10 合计兄弟（厘米） 65 66 66 67 68 70 70 71 71 72 686姐妹（厘

12、22.5米） 59 62 65 63 64 65 65 62 69 66 64012为研究国产四类新药阿卡波糖胶囊的降血糖效果，某医院用40名型糖尿病病人进行同期随机对照试验。试验者将这些病人随机等分到试验组（用阿卡波糖胶囊）和对照组（用拜唐苹胶囊），分别测得试验开始前和8周后的空腹血糖，算得空腹血糖下降值见下表，能否认为该国产四类新药阿卡波胶囊与拜唐苹胶囊对空腹血糖的降糖效果不同？试验组-0.70-5.602.002.800.703.504.005.807.10-0.50(n1=20)2.50-1.601.703.000.404.504.602.506.00-1.40对照组3.706.505

13、.005.200.800.200.603.406.60-1.10(n2=20)6.003.802.001.602.002.201.203.101.70-2.00答案：/*data ex3_7; input x ; if _n_21 then c=1; else c=2;其实就是分组cards;-0.70 -5.60 2.00 2.80 0.70 3.50 4.00 5.80 7.10 -0.50 2.50 -1.60 1.70 3.00 0.40 4.50 4.60 2.50 6.00 -1.40 3.70 6.50 5.00 5.20 0.80 0.20 0.60 3.40 6.60 -1.

14、10 6.00 3.80 2.00 1.60 2.00 2.20 1.20 3.10 1.70 -2.00;proc ttest; var x; class c;run; */data bb;do c=1 to 2;input n;do j=1 to n;input x; output;end; end;cards;20-0.70 -5.60 2.00 2.80 0.70 3.50 4.00 5.80 7.10 -0.50 2.50 -1.60 1.70 3.00 0.40 4.50 4.60 2.50 6.00 -1.4020 3.70 6.50 5.00 5.20 0.80 0.20 0.

15、60 3.40 6.60 -1.10 6.00 3.80 2.00 1.60 2.00 2.20 1.20 3.10 1.70 -2.00;proc ttest; var x; class c;run;两组等级资料的秩和检验data p193;do c=1 to 4;do g=1 to 4;input f;do k=1 to f;y=c;output;end;end;end;cards;0 3 5 32 5 7 59 5 3 36 2 2 0;proc npar1way wilcoxon;var y;class g;run;p205随机区组资料的秩和检验data k;do a= 1 to 3;

16、do b =1 to 8;input x;output;end; end;cards;5.3 5.5 5.1 5.2 5.3 5.2 5.6 5.56.2 6.4 4.8 5.0 6.0 4.9 6.7 6.25.8 6.0 5.5 5.5 5.6 5.4 6.0 5.9;proc sort out=ss;by b;proc rank data=ss;var x;by b;proc print;proc freq;tables a*b*x/scores=rank cmh2;run;p206线性相关data pli11_1;input x y;cards;1.1 141.2 131.0 150.

17、9 151.2 131.1 140.9 160.6 171.0 140.9 161.1 150.9 161.1 141.0 150.7 17;proc gplot;plot y*x;proc corr;var x y;run;p211秩相关data p211;input x y;cards;121 3138 2165 1310 0426 2540 2740 01060 01260 01290 01438 32004 0;proc corr spearman;var x y;run;p224直线回归data p224;input x y;cards;1.1 141.2 131.0 150.9 1

18、51.2 131.1 140.9 160.6 171.0 140.9 161.1 150.9 161.1 141.0 150.7 17;proc reg;model y=x;output out=ss p=p r=r; proc plot;plot r*p;run;等级资料秩相关data lian;do a=1 to 3;do b=1 to 3;input f;output;end;end;cards;5 32 2030 38 1015 10 10;proc corr spearman;var a b;7freq f;run;多元回归DATA a1;/*建立SAS数据集*/INPUT x1-x

19、4 y;DATALINES;1300 20.0 80 0.45 0.0661444 23.0 57 0.50 0.076786 26.5 64 1.50 0.0011652 23.0 84 0.40 0.1701756 29.5 72 0.90 0.1561754 30.0 76 0.80 0.1201200 22.5 69 1.80 0.0401500 21.8 77 0.60 0.1201200 27.0 58 1.70 0.1001476 27.0 65 0.65 0.1291820 22.0 83 0.40 0.1351436 28.0 68 2.00 0.099948 22.5 69

20、 2.00 0.0051440 21.5 79 2.40 0.0111084 28.5 59 3.00 0.0031844 26.0 73 1.00 0.1401116 35.0 92 2.80 0.0391656 20.0 83 1.45 0.0591536 23.0 57 1.50 0.087960 24.8 67 1.50 0.0391784 23.3 83 0.90 0.2221496 27.0 65 0.65 0.1451060 26.0 58 1.83 0.0291436 28.0 68 2.00 0.099;PROC REG ;MODEL y=x1-x4 ;PROC corr ;var x1-x4 y;proc corr;var x1 y;partial x2-x4;proc reg;model y=x1-x4/selection=stepwise sle=0.05 sls=0.05;RUN;