SAS 統計程序實作
CONTENTS
By DR. Yang , Yi-Chiang
2005.2.1
2011/12/10 1
前言
TITLE ‘ ’;
LIBNAME .... ' C:\ .. \ .. ..' ;
file : import data .xls ->.sas7bdat
tools : options-> preference->
result-> html-> printer
2011/12/10 2
資料管理 1
DATA ..;
SET . ;
/*改造*/
LABEL .... = ' '
.... = ' ' ;
2011/12/10 3
資料管理 2
/* If FORMAT then do before data */
PROC FORMAT;
VALUE A 1=‘ ’ 2=‘ ’ ; VALUE B 3=‘’ 4= ‘’;
DATA C; SET D;
FORMAT E A. F B. ;
2011/12/10 4
資料管理 3
PROC CONTENTS DATA=..VARNUM;
RUN;
2011/12/10 5
PROC MEANS
描述性統計值之計算
以PROC MEANS 來描述資料的特性
• Example AR_11_1
Describe of data characteristics
N MIN MAX RANGE MEAN CSS STD
STDERR SKEWNESS KURTOSIS ;
2011/12/10 6
PROC MEANS
描述性統計值之計算
以PROC MEANS 執行子集資料的描述
Example AR_11_2
Describe Subsets of data
2011/12/10 7
PROC MEANS
描述性統計值之計算
以PROC MEANS 執行一組樣本的t檢定
Example AR_11_3
One sample t-test
2011/12/10 8
PROC MEANS
描述性統計值之計算
以PROC MEANS 執行配對樣本的t檢定
Example AR_11_4
Paired-sample t-test
2011/12/10 9
PROC MEANS
描述性統計值之計算
在PROC MEANS 下界定WEIGHT的指令
Example AR_11_5
How to use WEIGHT statement
2011/12/10 10
PROC MEANS
描述性統計值之計算
在PROC MEANS 下建立輸出資料檔
Example AR_11_6
Creating an output data set
2011/12/10 11
PROC MEANS
描述性統計值之計算
在PROC MEANS 一個分類變項在輸出資
料檔上
Example AR_11_7
One CLASS variable used in an
output data set
_TYPE_ = 0, 1
2011/12/10 12
PROC MEANS
描述性統計值之計算
在PROC MEANS 兩個分類變項在輸出資
料檔上
Example AR_11_8
Two CLASS variables used in an
output data set
_TYPE_ = 0, 1, 2, 3
2011/12/10 13
PROC MEANS
描述性統計值之計算
PROC MEANS 在數個子集是由兩個分類
變項分割資料檔上(如上例 _TYPE_=3)
Example AR_11_8A
Correlating reading and punctuation
averages for six subgroups
2011/12/10 14
PROC MEANS
描述性統計值之計算
PROC MEANS 程序 總結
2011/12/10 15
PROC UNIVARIATE
描述性統計值之計算&繪圖
異於 PROC MEANS
計算 百分位數
圖表 次數分配表
檢定 常態性 符號 符號等級
2011/12/10 16
PROC UNIVARIATE
描述性統計值之計算&繪圖
用PROC UNIVARIATE程序來描述資料的
特性
• Example AR_12_1
Describe data characteristics
動差值 基本統計值 標準數的統計檢定
四分位數&百分位數 極端值
2011/12/10 17
PROC UNIVARIATE
描述性統計值之計算&繪圖
用PROC UNIVARIATE程序來產生次數分
配表
• Example AR_12_2
Frequency tabulation
動差值 基本統計值 標準數的統計檢定
四分位數&百分位數 極端值
次數分配表
2011/12/10 18
PROC UNIVARIATE
描述性統計值之計算&繪圖
用PROC UNIVARIATE程序來產生
莖葉圖 盒狀圖 常態機率圖
• Example AR_12_3
Stem Leaf
Boxplot
Normal probability plot
2011/12/10 19
PROC UNIVARIATE
描述性統計值之計算&繪圖
用PROC UNIVARIATE程序來判斷是否呈
常態分配
• Example AR_12_4
Tests for Normality
2011/12/10 20
PROC UNIVARIATE
描述性統計值之計算&繪圖
用PROC UNIVARIATE程序來從小樣本推
論母體的中位數
• Example AR_12_5
The Sign Test
2011/12/10 21
PROC UNIVARIATE
描述性統計值之計算&繪圖
用PROC UNIVARIATE程序來從小樣本推
論母體的中位數
• Example AR_12_6
Signed Rank Test
2011/12/10 22
PROC UNIVARIATE
描述性統計值之計算&繪圖
用PROC UNIVARIATE程序來執行數據的
加權
• Example AR_12_7
The WEIGHT statement
2011/12/10 23
PROC UNIVARIATE
描述性統計值之計算&繪圖
用PROC UNIVARIATE程序來執行資料子
集之描述性統計分析
• Example AR_12_8
Comprehensive descriptive analysis
of subgroup
2011/12/10 24
PROC UNIVARIATE
描述性統計值之計算&繪圖
在PROC UNIVARIATE程序內建立輸出資
料檔
• Example AR_12_9
Print an output data set
2011/12/10 25
PROC UNIVARIATE
描述性統計值之計算&繪圖
PROC UNIVARIATE程序 總結
2011/12/10 26
PROC TTEST
2母體平均數是否相等之統計檢定
在PROC TTEST程序內執行兩組獨立樣本
之 t 檢定
• Example AR_15_1
Independent-samples t-test of
gender difference on a var test
2011/12/10 27
PROC TTEST
2母體平均數是否相等之統計檢定
在PROC TTEST程序內執行兩組獨立樣本
之 t 檢定 同時比較數個變項
• Example AR_15_1A
Independent-samples t-test of
gender difference on Several var
test
2011/12/10 28
PROC TTEST
2母體平均數是否相等之統計檢定
在PROC TTEST程序內執行兩組獨立樣本
之 t 檢定 同時比較數個觀測體
• Example AR_15_2
Independent-samples t-test of
gender difference on Several Obs
test
2011/12/10 29
PROC TTEST
2母體平均數是否相等之統計檢定
在PROC TTEST程序內執行配對樣本之 t
檢定
• Example AR_15_3
Paired-samples t-test of
2 Var difference on test
Comparison of standardized Var
scores
2011/12/10 30
PROC TTEST
2母體平均數是否相等之統計檢定
在PROC TTEST程序內執行一組樣本之 t
檢定
• Example AR_15_4
One-samples t-test
Comparison of H0=scores
2011/12/10 31
PROC TTEST
2母體平均數是否相等之統計檢定
PROC TTEST程序 總結
2011/12/10 32
PROC FREQ
類別數據之分析
在PROC FREQ程序內執行一個一元次數表
• Example AR_14_1
One-way frequency tables
2011/12/10 33
PROC FREQ
類別數據之分析
在PROC FREQ程序內執行多個一元次數表
• Example AR_14_1A
One-way frequency tables
2011/12/10 34
PROC FREQ
類別數據之分析
在PROC FREQ程序內執行二元次數表
• Example AR_14_2
Two-way frequency tables
Cross tabular freq table
2011/12/10 35
PROC FREQ
類別數據之分析
在PROC FREQ程序內執行三元次數表
• Example AR_14_3
Three-way frequency tables
eg:
TABLES course*grade*race
Controlling for course=X202
Table of grade by race
2011/12/10 36
PROC FREQ
類別數據之分析
用PROC CHART程序來驗證三元的次數表
• Example AR_14_4
Verify Three-way table
2011/12/10 37
PROC FREQ
類別數據之分析
用PROC FREQ程序來鑑定兩個類別變項是
否有關係
• Example AR_14_5
EXPECTED CHISQ MEASURES
the relationship between two
categorical variables
2011/12/10 38
PROC FREQ
類別數據之分析
用PROC FREQ程序如何來處理2X2次數表
中的小樣本
• Example AR_14_6
Fisher's Exact Test
EXACT CHISQUARE EXPECTED
2011/12/10 39
PROC FREQ
類別數據之分析
用PROC FREQ程序如何來處理2X2次數表
中的一致性
• Example AR_14_7
McNemar's Test
Simple Kappa Coefficient
2011/12/10 40
PROC FREQ
類別數據之分析
用PROC FREQ程序如何來處理三元次數表
及其檢定
• Example AR_14_8
( CMH ) Cochran-Mantel-
Haenszel Statistics (Based on Ta
ble Scores) Breslow-Day Test for
Homogeneity of the Odds Ratios
2011/12/10 41
PROC FREQ
類別數據之分析
PROC FREQ程序 總結
2011/12/10 42
PROC CHART
資料之圖示法
在PROC CHART程序內利用縱軸圖及橫軸
圖來計算次數的分布
• Example AR_13_1
Frequency counts by vertical and
horizontal bar charts
2011/12/10 43
PROC CHART
資料之圖示法
在PROC CHART程序內執行縱軸百分比圖
• Example AR_13_2
Vertical percent chart based on two
methods of reading
2011/12/10 44
PROC CHART
資料之圖示法
在PROC CHART程序內執行累積次數的縱
軸圖
• Example AR_13_3
Vertical chart for cumulative
frequency
2011/12/10 45
PROC CHART
資料之圖示法
在PROC CHART程序內執行累積百分比的
縱軸圖
• Example AR_13_4
Vertical chart for cumulative
percent
2011/12/10 46
PROC CHART
資料之圖示法
在PROC CHART程序內執行根據平均值繪
製縱軸圖
• Example AR_13_5
Vertical chart based on averages
2011/12/10 47
PROC CHART
資料之圖示法
在PROC CHART程序內執行用特定的符號
畫次數分布圖
• Example AR_13_6
frequency bar chart with symbols
2011/12/10 48
PROC CHARROUPT
資料之圖示法
在PROC CHART程序內執行用特定的符號
畫次數分布圖
• Example AR_13_7
A frequency bar chart with symbols
and averages
/SUBGROUP=..
2011/12/10 49
PROC CHART
資料之圖示法
在PROC CHART程序內執行用MEANS程
序來驗證CHART產生的圖形
• Example AR_13_8
Verifying the PROC CHART results
by PROC MEANS
2011/12/10 50
PROC CHART
資料之圖示法
在PROC CHART程序內執行並列的縱軸圖
• Example AR_13_9
Vertical charts side by side
/GROUP=..
2011/12/10 51
PROC CHART
資料之圖示法
在PROC CHART程序內執行三維立體直方
圖
• Example AR_13_10
3D block chart
2011/12/10 52
PROC CHART
資料之圖示法
在PROC CHART程序內執行圓形圖
• Example AR_13_11
Pie chart for averages
2011/12/10 53
PROC CHART
資料之圖示法
在PROC CHART程序內執行星形圖
• Example AR_13_12
A Star chart
2011/12/10 54
PROC CHART
資料之圖示法
在PROC CHART程序內執行複雜的星形圖
• Example AR_13_13
Another star chart with more spikes
2011/12/10 55
PROC CHART
資料之圖示法
PROC CHART 程序 總結
2011/12/10 56
PROC PLOT
資料所呈現之趨勢圖
在PROC PLOT程序內描繪兩個連續變項的
關係
• Example AR_18_1
Plotting the relationship between
2 VAR
2011/12/10 57
PROC PLOT
資料所呈現之趨勢圖
在PROC PLOT程序內自訂座標軸的間隔並
添加參考線
• Example AR_18_2
Defining intervals on an axis and
imposing reference lines
2011/12/10 58
PROC PLOT
資料所呈現之趨勢圖
在PROC PLOT程序內改變繪圖的符號
• Example AR_18_3
Changing the plotting symbol
2011/12/10 59
PROC PLOT
資料所呈現之趨勢圖
在PROC PLOT程序內改變座標軸的長度
• Example AR_18_4
Manipulating the length of each axis
2011/12/10 60
PROC PLOT
資料所呈現之趨勢圖
在PROC PLOT程序內改變座標軸的長度
• Example AR_18_4A
Manipulating the length of each axis
2011/12/10 61
PROC PLOT
資料所呈現之趨勢圖
在PROC PLOT程序內重疊繪圖
• Example AR_18_5
Presenting several plots on the
same page
2011/12/10 62
PROC PLOT
資料所呈現之趨勢圖
在PROC PLOT程序內重疊繪圖
• Example AR_18_5A
Presenting several plots on the
same page
2011/12/10 63
PROC PLOT
資料所呈現之趨勢圖
在PROC PLOT程序內將字元標在座標軸上
• Example AR_18_6
Using characters as tick points
on an axis
2011/12/10 64
PROC PLOT
資料所呈現之趨勢圖
在PROC PLOT程序內繪三維的輪廓圖軸上
• Example AR_18_7
A 3-D contour plot
2011/12/10 65
PROC PLOT
資料所呈現之趨勢圖
PROC PLOT程序 總結
2011/12/10 66
PROC GLM
變異數分析
在PROC GLM程序內執行一因子的變異數
分析 先復習PROC MEANS
• Example AR_16_0
PROC MEANS mean std max min;
CLASS indep1; VAR score;
Average drinks and variability
2011/12/10 67
PROC GLM
變異數分析
在PROC ANOVA 程序內執行一因子的
ANOVA 檢定
• Example AR_16_1
Tukey's Studentized Range (HSD)
Test for score
Bonferroni (Dunn) t Tests for
score
2011/12/10 68
PROC GLM
變異數分析
在PROC GLM程序內執行一因子的ANOVA
F 檢定
• Example AR_16_2
The ANOVA Procedure t Tests
(LSD)
Tukey's Studentized Range (HSD)
Test for cdvalue
2011/12/10 69
PROC GLM
變異數分析
在PROC GLM程序內執行二因子的變異數
分析
• Example AR_16_3
Bonferroni (Dunn) t Tests for score
2011/12/10 70
PROC GLM
變異數分析
在PROC GLM程序內執行二因子的變異數
分析
• Example AR_16_3
Bonferroni (Dunn) t Tests for score
2011/12/10 71
PROC GLM
變異數分析
在PROC GLM程序內執行以細格平均數之
點圖來確認無交互作用的存在
• Example AR_16_4
Confirming no interaction with a
plot of cell means
2011/12/10 72
PROC GLM
變異數分析
在PROC GLM程序內執行隨機集區設計
• Example AR_16_5
Sidak t Tests
2011/12/10 73
PROC GLM
變異數分析
在PROC GLM程序內執行拉丁方格設計
• Example AR_16_6
Latin-square design
2011/12/10 74
PROC GLM
變異數分析
在PROC GLM程序內執行拉丁方格設計中
將交互作用併入殘差
• Example AR_16_7
Latin-square design
2011/12/10 75
PROC GLM
變異數分析
在PROC GLM程序內執行重複觀察的實驗
設計
Example AR_16_7 ~ 10
VERY DIFFICULTY !!!
2011/12/10 略 76
PROC GLM
變異數分析
在PROC GLM程序內執行重複觀察的實驗
設計 共變異數分析 ……> 很難
Example AR_16_7 ~ 12
~ VERY DIFFICULTY !!!
2011/12/10 略 77
PROC GLM
變異數分析
PROC GLM程序總結
2011/12/10 78
PROC NPAR1WAY
無母數檢定
在PROC NPAR1WAY程序內執行針對兩組
樣本之中位數進行無母數檢定
• Example AR_17_1
Nonparametric test of two group
means
WILCOXON 2-sample test
2011/12/10 79
PROC NPAR1WAY
無母數檢定
在PROC NPAR1WAY程序內執行VW檢定
來比較兩組樣本之代表值
• Example AR_17_2
Comparison of two sample typical
scores by Van der Waerden test
2011/12/10 80
PROC NPAR1WAY
無母數檢定
在PROC NPAR1WAY程序內執行無母數
單因子變異數分析
• Example AR_17_3
Nonparametric one-way analysis of
variance by ranks
2011/12/10 81
PROC NPAR1WAY
無母數檢定
PROC NPAR1WAY程序 總結
2011/12/10 82
PROC CORR
相關分析
在PROC CORR程序內建立量化變項間關
係強弱的四個指標
• Example AR_19_1
PEARSON SPEARMAN KENDALL
HOEFFDING Four ways to quantify
the strength of a relationship
2011/12/10 83
PROC CORR
相關分析
在PROC CORR程序內用PEARSON簡單相
關係數比較男女間的差異
Example AR_19_2
Compare women & men in term of
PEARSON r
2011/12/10 84
PROC CORR
相關分析
在PROC CORR程序內如何存取相關矩陣
Example AR_19_3
Correlationation Matrix
2011/12/10 85
PROC CORR
相關分析
在PROC CORR程序內如何作偏相關分析
Example AR_19_4
Partial Correlation between..
PEARSON SPEARMAN
2011/12/10 86
PROC CORR
相關分析
在PROC CORR程序內如何作表示測驗信
度的CRONBACH 阿爾法係數
Example AR_19_5
Cronbach Coefficient Alpha
2011/12/10 87
PROC CORR
相關分析
PROC CORR程序 總結
2011/12/10 88
PROC REG
迴歸分析 導出數學方程式
以PROC REG 執行簡單迴歸模型
Example AR_20_1
Simple Regression Analysis
2011/12/10 89
PROC REG
迴歸分析 導出數學方程式
以PROC REG 執行含兩個預測變項的迴歸
分析
Example AR_20_2
2 predictor model for Regression
Analysis
2011/12/10 90
PROC REG
迴歸分析 導出數學方程式
以PROC REG 在迴歸分析中執行類別預測
變項
Example AR_20_3
treat Categorical Predictor in
Regression Analysis
2011/12/10 91
PROC REG
迴歸分析 導出數學方程式
以PROC REG 在迴歸分析中執行最佳變項
的選擇
Example AR_20_4
Selecting the best set of predictors
2011/12/10 92
PROC REG
迴歸分析 導出數學方程式
以PROC REG 在迴歸分析中執行殘差分析
Example AR_20_5A
Checking the independence
assumption
2011/12/10 93
PROC REG
迴歸分析 導出數學方程式
以PROC REG 在迴歸分析中執行殘差分析
Example AR_20_5B
Checking the equal variance
assumption
2011/12/10 94
PROC REG
迴歸分析 導出數學方程式
以PROC REG 在迴歸分析中執行殘差分析
Example AR_20_5C
Checking the normality assumption
2011/12/10 95
PROC REG
迴歸分析 導出數學方程式
以PROC REG 在迴歸分析中使用相關矩陣
作為輸入資料檔
Example AR_20_6
Using a correlation matrix as input
2011/12/10 96
PROC REG
迴歸分析 導出數學方程式
PROC REG 程序 總結
2011/12/10 97
PROC LOGISTIC
邏輯斯複迴歸分析 導出數學方程式
用PROC LOGISTIC作反應變項為二分變
項或序位變項之迴歸分析 勝算比
Example AR_21_1
金屬條分析 2 Covariates & R/N
logist(p)
2011/12/10 98
PROC LOGISTIC
邏輯斯複迴歸分析 導出數學方程式
用PROC LOGISTIC作反應變項為二分變
項或序位變項之迴歸分析 勝算比
Example AR_21_2
癌症復發率 6 Covariates & 01 logist(p)
2011/12/10 99
PROC LOGISTIC
邏輯斯複迴歸分析 導出數學方程式
用PROC LOGISTIC作反應變項為二分變
項或序位變項之迴歸分析 勝算比
Example AR_21_3
癌症復發率的反向迴歸分析
2011/12/10 100
PROC LOGISTIC
邏輯斯複迴歸分析 導出數學方程式
用PROC LOGISTIC作反應變項為二分變
項或序位變項之迴歸分析 勝算比
Example AR_21_4
四肢皮膚收縮的研究分析
2011/12/10 101
PROC LOGISTIC
邏輯斯複迴歸分析 導出數學方程式
PROC LOGISTIC 程序 總結
2011/12/10 102