[医学]Ch5 SAS基本统计过程一 mean uni ttest npar.ppt

* * * * * * * * * * * * * * * * * * * 语句说明 PROC NPAR1WAY [options]； CLASS variable; VAR variables; 这里，前两句(下划线部分)是必须的语句，后两句是可选的。 PROC NPAR1WAY [options]; 这里的选项options有较多选择，下面介绍两个常用的选项： Data=SAS-data-set：指定分析数据集名。 Wilcoxon：对数据或wilcoxon的秩次得分进行秩分析。对于两组，为wilcoxon秩和检验； 对于多组，为Kruskal-Wallis H检验。 2) Class variable; 这个语句是必须的，它指定一个且只能指定一个分类变量。该变量用来识别数据中的各个类。Class变量可以数值变量，也可以是字符变量。 3) VAR variables; 该语句命名要分析的响应变量或自变量。如果省略，则过程分析数据集中所有数值变量，除了Class变量以外。 例4-6. (两样本比较的秩和检验) 10名健康人与7名铅作业工人血铅测得值如下： 健康人 5 5 5 6 12 13 15 21 21 31 铅作业 17 17 18 20 34 43 44 试作两样本比较的秩和检验。 Data ex; Input x@@; If _n_=10 then g=1; Else g=2; Cards; 5 5 5 6 12 13 15 21 21 31 17 17 18 20 34 43 44 ; PROC NPAR1WAY wilcoxon; CLASS g; VAR x; RUN; 建立分析的SAS数据集，并产生分组变量 进行秩和检验 指明 分析变量x及分组变量g Wilcoxon检验 (两组比较) Z=2.2039, p=0.0275 Kruskal-Wallis检验 (两组及多组比较 ) Chi-Square=5.0754 p=0.0243 例4-7. 某医生用中药治疗胃癌，对梗塞型胃癌和溃疡型胃癌治疗效果见下表： 治愈 显效 有效 无效 梗塞型 6 7 4 1 溃疡型 2 10 8 5 试用两组比较的秩和检验比较其疗效。 首先应产生相应的SAS数据集 下列程序中： f：频数，row：分组(第几行)，col：结果(第几列) Data ex5; input f row col @@; do I=1 to f; output; end; cards; 6 1 1 7 1 2 4 1 3 1 1 4 2 2 1 10 2 2 8 2 3 5 2 4 ; PROC NPAR1WAY data=ex5 wilcoxon; CLASS row; VAR col; RUN; 注意： CLASS语句中一定是“分组变量”，VAR语句中是“结果变量” 我们用UNIVARIATE过程来处理配对资料的秩和检验问题。 例4-8. 某医院测定10名受试者针刺治疗前后痛阈的变化，结果如下表： 针刺前 600 600 685 1050 900 1125 1400 750 1000 1500 针刺后 610 700 575 600 600 1425 1350 825 800 1440 试作配对资料的秩和检验。 2. 配对秩和检验 Data ex3; Input x1 x2; d=x1-x2; Cards; 600 610 600 700 685 575 1050 600 900 600 1125 1425 1400 1350 750 825 1000 800 1500 1440 ; PROC UNIVARIATE data=ex3; VAR d; RUN; Tests?for?Location:?Mu0=0 Test Statistic p Value Students t t 1.015508 Pr |t| 0.3364 Sign M 1 Pr = |M| 0.7539 Signed Rank S 9 Pr = |S| 0.3906 配对秩和检验结果： S=9, p=0.3906 * * * * * * * * * * * * * * * * * * * * * Slide * 彭斌 第四章 SAS基本统计过程 我们曾经