第四章SAS系统的基本统计分析功能

数据集PULSE是对15 数据集PULSE是对15个试验者服用某种减缓心跳 PULSE是对15个试验者服用某种减缓心跳 例4.7 药物前后每分钟心跳次数所作的记录（PRE为服 药物前后每分钟心跳次数所作的记录（PRE为服 药前，POST为服药后），考察该种药物的疗效 为服药后），考察该种药物的疗效。 药前，POST为服药后），考察该种药物的疗效。 data tempulse; set sunny.pulse; d=pre-post; d=pre-post; run; run; maxdec= proc means mean stderr prt maxdec=4; d; var d; run; run; 输出结果
编程
分析员 应用
例4.6
proc npar1way data=sunny.bclass data=sunny.bclass wilcoxon; wilcoxon; class sex; var age; run; run;
输出结果
配对样本的均值检验 设( Xi ,Yi ) 为第 i 个样本单位实验前后的对比测量值，且 个样本单位实验前后的对比测量值， D = Xi Yi , ( i =1,2,..., n) ， 为D 总体的均值，则检验假设： i i 总体的均值，则检验假设：
输出结果: 输出结果:
用分析员应用作置信区间 分析员应用作置信区间 （调用方法：Solutions 调用方法： 步骤
Analysis Analyst ）
选定数据集 Statistics Hypothesis tests OneOne-Sample t-test for a mean ttest interval 选定分析变量
用SAS/INSIGHT作置信区间 SAS/INSIGHT作置信区间 （调用方法：Solutions 调用方法：
Analysis Interactive Data Analysis ）
步骤
选定数据集 Analyze Distribution Output 选定分析变量 basic confidence intervals
总体均值的区间估计 （给定显著性水平α） 例4.2 若例4.1 若例4.1中，新生总体IQ成绩的标准差未知，其 4.1中 新生总体IQ成绩的标准差未知， IQ成绩的标准差未知 样本标准差的计算结果为5 样本标准差的计算结果为5分，在显著性水平为 0.05条件下 试建立IQ平均成绩的置信区间。 条件下， IQ平均成绩的置信区间 0.05条件下，试建立IQ平均成绩的置信区间。 data IQSCORE2; IQSCORE2; n=25;mean=120 s=5 alpha=0.05 25;mean=120; 0.05; n=25;mean=120; s=5; alpha=0.05; (t分布的左侧 t=tinv((1 alpha/2),n- (t分布的左侧p分位数) t=tinv((1- alpha/2),n-1); 分布的左侧p分位数) lc=meanlc=mean-t*s/sqrt(n); uc=mean+t*s/sqrt(n); print; proc print; run; run;
第四章 SAS系统的基本统计分析功能 SAS系统的基本统计分析功能
SAS
★ §4.1
参数估计 §4.2 假设检验 §4.3 相关分析 §4.4 回归分析
参数估计
指利用所计算样本的均值 x 、方差 s 、 比率 p 等样本指标来估计总体的均值 、 2 σ 、比率 P等总体指标。 等总体指标。 方差
2
总体方差的区间估计
（给定显著性水平 α ）
2 2
( n 1) s
χ12α
2
2
( n 1) s ≤σ ≤ 2
χα
2
例4.3
某种零件的直径服从正态分布，随机抽出100 某种零件的直径服从正态分布，随机抽出100 个进行测量得到方差 s2 = 9，在显著性水平为 0.01条件下 试建立总体方差的置信区间。 条件下， 0.01条件下，试建立总体方差的置信区间。 data E3; E3; n=100; s2=9; alpha=0.01; n=100; s2=9 alpha=0.01; chi1=cinv(alpha/2,n- (卡方分布的左侧p分位数) chi1=cinv(alpha/2,n-1);卡方分布的左侧p分位数) chi2=cinv((1 alpha/2),nchi2=cinv((1-alpha/2),n-1); lc=(nlc=(n-1)*s2/chi2; uc=(n-1)*s2/chi1; uc=(nproc print;run; print;run;
方差检验 输出结果
若样本呈非正态分布 若样本呈非正态分布，则t检验无效，可使用非参 非正态分布， 检验无效， 无效 数检验中的Wilcoxon秩和检验方法 Wilcoxon秩和检验方法检验两独立样本 数检验中的Wilcoxon秩和检验方法检验两独立样本 的中心位置是否相同。 的中心位置是否相同。 PROC NPAR1WAY DATA=数据集 Wilcoxon； DATA=数据集 Wilcoxon； CLASS分类变量名 分类变量名； CLASS分类变量名； VAR分析变量名列 分析变量名列； VAR分析变量名列； RUN； RUN； 选定数据集 Statistics ANOVA OneNonparametric One-Way ANOVA Wilcoxon. 选定分析变量 Tests
两独立样本的均值检验 设 X1,..., Xn1 和 Y1,...,Yn2 分别是来自两个正态分布总 2 的独立样本，检验假设： 体 N ( 1,σ12 ) 和 N ( 2 ,σ2 ) 的独立样本，检验假设：
H0 : 1 = 2
H1 : 1 ≠ 2
用编程作两独立样本的均值检验 编程作两独立样本的均值检验 TTEST过程的一般形式： TTEST过程的一般形式： 过程的一般形式 PROC TTEST DATA=数据集； DATA=数据集 数据集； CLASS分类变量名 分类变量名； CLASS分类变量名； VAR分析变量名列； VAR分析变量名列； 分析变量名列 RUN； RUN；
PROC UNIVARIATE DATA= SUNNY.BCLASS VAR HEIGHT ; RUN; NORMAL; NORMAL;
输出结果（部分） 输出结果（部分）
用SAS/INSIGHT作正态性检验 步骤 选定数据集 Analyze Output 选定分析变量 用分析员应用作正态性检验 分析员应用作正态性检验 步骤 选定数据集 Distributions Statistics Descriptive Fit Normal. 选定分析变量 Distribution Test For Normality.
信下限、区间和上限，缺省为置信区间）。 信下限、区间和上限，缺省为置信区间）。
例4.5 proc capability data=sunny.bclass data=sunny.bclass noprint; noprint; intervals height weight/method=4,6 weight/method=4 alpha=0.10; alpha=0.10; intervals height/method=4,6 height/method=4 type=lower; type=lower; run; run;
例4.4
proc means data=sunny.bclass clm; mean var clm; var height weight; run; run;
输出结果: 输出结果:
对于正态分布总体，当已知样本的原始数据时， 对于正态分布总体，当已知样本的原始数据时， 可以利用Capability过程的Intervals语句 Capability过程 语句计算总 可以利用Capability过程的Intervals语句计算总 Means过程 体均值和方差的置信区间，或者利用Means过程的 体均值和方差的置信区间，或者利用Means过程的 CLM统计量计算总体均值95 的置信区间。 统计量计算总体均值95﹪ CLM统计量计算总体均值95﹪的置信区间。
CAPABILITY过程的一般形式： CAPABILITY过程的一般形式： 过程的一般形式 PROC CAPABILITY DATA=数据集； DATA=数据集 数据集； INTERVALS 变量名列/选项； 变量名列/选项； RUN； RUN； INTERVALS语句的选项： INTERVALS语句的选项： 语句的选项 Method= Method=1|2|3|4|5|6 设定区间估计统计量（4为均值，6为 设定区间估计统计量（ 为均值，
步骤
选定数据集 Statistics Hypothesis tests TwoTwo-Sample t-test for means tOK. 选定分析变量和分类变量
输出结果
？
（作均值检验前需先检验两样本方差是否相等） 作均值检验前需先检验两样本方差是否相等）
检验两独立样本的方差是否相等： 检验两独立样本的方差是否相等： 的方差是否相等 步骤 选定数据集 Statistics Hypothesis tests TwoTwo-Sample t-test for variances tOK. 选定分析变量和分类变量
H0 : = 0
编程
H1 : ≠ 0
先构造每对观测值的差值变量，然后利用MEAN 先构造每对观测值的差值变量，然后利用MEAN UNIVARIATE过程检验该变量均值是否为零 过程检验该变量均值是否为零. 或UNIVARIATE过程检验该变量均值是否为零.
选定数据集 Statistics Hypothesis tests 分析员 TwoTwo-Sample paired t-test for means t应用 OK. 选定分析变量
输出结果: 输出结果:
第四章 SAS系统的基本统计分析功能 SAS系统的基本统计分析功能SASຫໍສະໝຸດ ★ §4.1 ★ §4.2