内容发布更新时间 : 2024/12/22 20:56:14星期一 下面是文章的全部内容请认真阅读。

SAS上机练习题（二）

1、某农村地区1999年14岁女孩的身高资料列于表1。

表1 谋农村地区1999年14岁女孩身高资料（cm）

142.3 148.8 142.7 144.4 144.7 145.1 143.3 154.2 152.3 142.7 137.9 143.9 141.2 139.3 145.8 142.2 137.9 141.2 150.6 142.7 142.4 141.5 141.9 147.9 125.8 139.9 148.9 154.9 145.7 140.8 148.8 147.8 146.7 132.7 149.7 154.0 158.2 138.2 149.8 151.1 140.5 143.4 152.9 147.5 147.7 162.6 141.6 143.6 144.0 150.6 150.8 147.9 136.9 146.5 130.6 142.5 149.0 145.4 139.5 148.9 141.8 148.1 145.4 134.6 130.5 145.2 146.2 146.4 142.4 137.1 144.0 129.4 142.8 132.1 141.8 143.3 143.8 134.7 147.1 140.9 142.5 146.6 135.5 146.8 156.3 150.0 147.3 142.9 141.4 134.7 146.6 134.5 135.1 141.9 142.1 138.1 134.9 146.7 138.5 139.6 148.8 150.3 140.7 143.5 140.2 143.6 138.7 138.9 143.5 139.9 133.1 145.9 139.2 137.4 142.3 160.9 137.7 142.9 126.8 问题：（1）计算均数、中位数；计算均数的95%可信区间；

（2）计算四分位间距、标准差、变异系数；计算标准误； （3）请进行正态性检验； （4）观察频数分布情况；

参考程序：

156.6 151.3 139.6 140.1 138.9 144.5 141.4 137.4 138.5 139.2 134.4

data a;

input x@@; cards;

142.3 148.8 142.7 144.4 144.7 145.1 143.3 154.2 137.9 143.9 141.2 139.3 145.8 142.2 137.9 141.2 142.4 141.5 141.9 147.9 125.8 139.9 148.9 154.9 148.8 147.8 146.7 132.7 149.7 154.0 158.2 138.2 140.5 143.4 152.9 147.5 147.7 162.6 141.6 143.6 150.8 147.9 136.9 146.5 130.6 142.5 149.0 145.4 141.8 148.1 145.4 134.6 130.5 145.2 146.2 146.4 144.0 129.4 142.8 132.1 141.8 143.3 143.8 134.7 142.5 146.6 135.5 146.8 156.3 150.0 147.3 142.9 146.6 134.5 135.1 141.9 142.1 138.1 134.9 146.7 148.8 150.3 140.7 143.5 140.2 143.6 138.7 138.9 133.1 145.9 139.2 137.4 142.3 160.9 137.7 142.9 ;

proc means data=a n mean median clm qrange std cv stderr; var x;

proc univariate data=a normal; histogram x; var x; run;

精选

152.3 150.6 145.7 149.8 144.0 139.5 142.4 147.1 141.4 138.5 143.5 126.8

142.7 142.7 140.8 151.1 150.6 148.9 137.1 140.9 134.7 139.6 139.9 156.6 151.3 139.6 140.1 138.9 144.5 141.4 137.4 138.5 139.2 134.4

2、某研究员欲分析急性染毒对肝脏功能的影响，将40只小鼠分为两组，雄雌各半，试验组进行急性染毒试验，染毒后2小时测定血液中的ALT，整理的结果见下表。 （1） 比较实验组与对照组的ALT有无差异；

ALT(丙氨酸转氨酶) NO 1 2 3 4 5 6 7 8 9 10

参考程序：

male control Treaetment 38.4 30.6 25.5 29.9 29.0 33.6 39.4 42.5 46.7 49.3 57.8 62.4 57.5 49.9 55.6 47.3 50.1 48.7 57.9 60.4 Female control Treaetment 39.7 34.6 40.5 44.3 41.2 25.6 30.6 37.6 39.5 25.1 55.7 63.8 59.9 49.7 48.8 51.3 53.8 42.6 61.7 65.3 data a2;

do sex='M','F'; do grp='C','T'; input alt@@; output; end; end; cards;

38.4 57.8 39.7 30.6 62.4 34.6 25.5 57.5 40.5 29.9 49.9 44.3 29.0 55.6 41.2 33.6 47.3 25.6 39.4 50.1 30.6 42.5 48.7 37.6 46.7 57.9 39.5 49.3 60.4 25.1 ;

proc anova data=a2; class sex grp;

model alt=sex grp; means grp; run; quit;

55.7 63.8 59.9 49.7 48.8 51.3 53.8 42.6 61.7 65.3

精选

3、24名志愿者随机分成两组，每组12人，接受降胆固醇试验，甲组为特殊饮食组，乙组为药物治疗组。受试者试验前后各测量一次血清胆固醇（mmol/L），数据见下表：

甲组

受试者 1 2 3 4 5 6 7 8 9 10 11 12

试验前 6.11 6.81 6.48 7.59 6.42 6.94 9.17 7.33 6.94 7.67 8.15 6.60

试验后 6.00 6.83 4.49 7.28 6.30 6.64 8.42 7.00 6.58 7.22 6.57 6.17

受试者 1 2 3 4 5 6 7 8 9 10 11 12

乙组 试验前 6.90 6.40 6.48 7.00 6.53 6.70 9.10 7.31 6.96 6.81 8.16 6.98

试验后 6.93 6.35 6.41 7.10 6.41 6.68 9.05 6.83 6.91 6.73 7.65 6.52

问题：

（1）试建立能够满足以下统计分析的SAS数据集。

（2）试判断两组受试者试验前血清胆固醇水平是否相等？ （3）分别判断两种降胆固醇措施是否有效？ （4）试判断两种降胆固醇措施的效果是否相同？ 参与程序：

data a3;

do grp='甲组','乙组';

input id before after @@; cha=before-after; output; end; cards; 1 6.11 6.00 1 6.90 2 6.81 6.83 2 6.40 3 6.48 4.49 3 6.48 4 7.59 7.28 4 7.00 5 6.42 6.30 5 6.53 6 6.94 6.64 6 6.70 7 9.17 8.42 7 9.10 8 7.33 7.00 8 7.31 9 6.94 6.58 9 6.96 10 7.67 7.22 10 6.81 11 8.15 6.57 11 8.16 12 6.60 6.17 12 6.98 ;

6.93 6.35 6.41 7.10 6.41 6.68 9.05 6.83 6.91 6.73 7.65 6.52

精选

proc ttest data=a3; class grp; var before;

proc ttest data=a3; paired before*after; where grp='甲组'; run;

proc ttest data=a3; paired before*after; where grp='乙组'; run;

proc ttest data=a3; class grp; var cha; run;

4、某职业病防治所对30名矿工分别测定血清铜蓝蛋白含量（μmol/L），资料如下。问各期血清铜蓝蛋白含量的测定结果有无差别？

疾病分期 0期 I期 II期

测定结果

8.0 9.0 5.8 6.3 5.4 8.5 5.6 5.4 5.5 7.2 5.6 8.5 4.3 11.0 9.0 6.7 9.0 10.5 7.7 7.7

11.3 7.0 9.5 8.5 9.6 10.8 9.0 12.6 13.9 6.5

请用秩和检验完成分析。

参考程序：

data a4; input x@@;

if _n_<=11 then grp='0期'; else if _n_<=20 then grp='1期'; else grp='2期'; cards;

8.0 9.0 5.8 6.3 5.4 8.5 5.6 5.4 5.5 7.2 5.6 8.5 4.3 11.0 9.0 6.7 9.0 10.5 7.7 7.7

11.3 7.0 9.5 8.5 9.6 10.8 9.0 12.6 13.9 6.5 ;

proc npar1way wilcoxon; class grp; var x; run;

精选

教材上的练习题

(卫生统计学，方积乾主编，第六版)

t检验 p121页，1，3，4题 第一题

方法一： data a1; input x@@; cards;

20.99 20.41 20.10 20.00 20.19 22.60 20.99 20.41 20.00 23.00 22.00 ;

proc ttest data=a1 h0=20.70; var x; run; 方法二： data a1; input x@@; d=x-20.70; cards;

20.99 20.41 20.10 20.00 20.19 22.60 20.99 20.41 20.00 23.00 22.00 ;

proc ttest data=a1 ; var d; run;

第三题

data a3; input x@@;

if _n_<=10 then grp=1;else grp=2; cards;

6.2 3.7 5.8 2.7 3.9 6.1 6.7 7.8 3.8 6.9 8.5 6.8 11.3 9.4 9.3 7.3 5.6 7.9 7.2 8.2 ; run;

proc ttest data=a3 ; class grp; var x; run;

精选