/* Lesson 10-02b */
/* File Name = les1002b.sas 01/11/18 */
options nocenter linesize=78 pagesize=30; : 画面出力サイズ指定
proc printto print = 'Kougi/les1002_Results.txt' new; : 計算結果をファイル出力
data dragons;
infile 'Kougi/Dragons.csv'
firstobs=4 dlm=',' dsd missover encoding=sjis termstr=crlf; : ファイル特性
input No $ Name : $20. Pos $ DoBY DoBM DoBD : 読み込む変数(項目)のリスト
shintyou taijyuu tou $ da $ bikou $; : 複数行にわたって
proc print data=dragons(obs=20); : 先頭20ケースを表示
run;
proc means data=dragons; : 平均値等を算出
run;
proc freq data=dragons; : 頻度を算出。個々&クロス集計
table DoBY DoBM tou da tou*da;
run;
proc chart data=dragons; : ヒストグラムを描く
hbar shintyou; : 水平棒グラフ。変量を指定。
vbar shintyou; : 垂直棒グラフ。変量を指定。
hbar DoBY;
vbar DoBY;
run;
proc univariate data=dragons plot; : 基礎統計量の算出。plot オプションを指定。
var shintyou DoBY; : 指定した変量について
run;
: Stem and Leaf を表示させるには
ods listing; : 以下の plot オプション用に指定。
ods graphics off; : 以下の plot オプション用に指定。
proc univariate data=dragons plot; : 基礎統計量の算出。plot オプションを指定。
var shintyou DoBY; : 指定した変量について
run;
Thursday, January 11, 2018 11:25:25 AM 55
s
h t
i a
n i b
N D D D t j i
O a P o o o y y t k
b N m o B B B o u o d o
s o e s Y M D u u u a u
1 80 森 繁和 監督 1954 11 18 181 81
2 11 小笠原 慎之莉 投手 1997 10 8 180 83 左 左
3 12 田島 慎二 投手 1989 12 21 181 84 右 右
4 13 岩瀬 仁紀 投手 1974 11 10 181 85 左 左
5 14 佐藤 優 投手 1993 6 29 187 85 右 左
6 16 又吉 克樹 投手 1990 11 4 180 74 右 右
7 17 柳 裕也 投手 1994 4 22 180 83 右 右
8 18 鈴木 翔太 投手 1995 6 16 183 74 右 右
9 19 吉見 一起 投手 1984 9 19 182 91 右 右
10 20 野村 亮介 投手 1993 7 9 187 85 右 右
11 21 岡田 俊哉 投手 1991 12 5 179 65 左 左
12 22 大野 雄大 投手 1988 9 26 183 78 左 左
13 24 福谷 浩司 投手 1991 1 9 183 90 右 右
14 25 武藤 祐太 投手 1989 6 14 178 85 右 右
15 28 丸山 泰資 投手 1995 2 5 176 79 右 右
16 29 山井 大介 投手 1978 5 10 179 82 右 右
17 30 阿知羅 拓馬 投手 1992 11 20 190 95 右 右
18 33 祖父江 大輔 投手 1987 8 11 175 78 右 左
Thursday, January 11, 2018 11:25:25 AM 56
Obs No Name Pos DoBY DoBM DoBD shintyou taijyuu tou da bikou
19 34 福 敬登 投手 1992 6 16 178 87 左 左
20 41 浅尾 拓也 投手 1984 10 22 182 75 右 右
Thursday, January 11, 2018 11:25:25 AM 57
The MEANS Procedure
Variable N Mean Std Dev Minimum Maximum
------------------------------------------------------------------------------
DoBY 81 1989.11 6.4226163 1954.00 1998.00
DoBM 81 6.9629630 3.5122800 1.0000000 12.0000000
DoBD 81 14.5308642 8.1257098 1.0000000 30.0000000
shintyou 81 180.5061728 5.0599493 167.0000000 200.0000000
taijyuu 81 82.7407407 8.7045646 65.0000000 125.0000000
------------------------------------------------------------------------------
Thursday, January 11, 2018 11:25:25 AM 58
The FREQ Procedure
Cumulative Cumulative
DoBY Frequency Percent Frequency Percent
---------------------------------------------------------
1954 1 1.23 1 1.23
1974 1 1.23 2 2.47
1977 2 2.47 4 4.94
1978 2 2.47 6 7.41
1981 2 2.47 8 9.88
1983 1 1.23 9 11.11
1984 5 6.17 14 17.28
1985 3 3.70 17 20.99
1986 2 2.47 19 23.46
1987 3 3.70 22 27.16
1988 7 8.64 29 35.80
1989 7 8.64 36 44.44
1990 4 4.94 40 49.38
1991 9 11.11 49 60.49
1992 7 8.64 56 69.14
1993 8 9.88 64 79.01
1994 6 7.41 70 86.42
1995 4 4.94 74 91.36
1996 4 4.94 78 96.30
1997 1 1.23 79 97.53
1998 2 2.47 81 100.00
Thursday, January 11, 2018 11:25:25 AM 59
The FREQ Procedure
Cumulative Cumulative
DoBM Frequency Percent Frequency Percent
---------------------------------------------------------
1 6 7.41 6 7.41
2 7 8.64 13 16.05
3 5 6.17 18 22.22
4 4 4.94 22 27.16
5 7 8.64 29 35.80
6 5 6.17 34 41.98
7 10 12.35 44 54.32
8 5 6.17 49 60.49
9 7 8.64 56 69.14
10 7 8.64 63 77.78
11 11 13.58 74 91.36
12 7 8.64 81 100.00
Thursday, January 11, 2018 11:25:25 AM 60
The FREQ Procedure
Cumulative Cumulative
tou Frequency Percent Frequency Percent
--------------------------------------------------------
右 63 78.75 63 78.75
左 17 21.25 80 100.00
Frequency Missing = 1
Cumulative Cumulative
da Frequency Percent Frequency Percent
-----------------------------------------------------------
右 44 55.00 44 55.00
左 35 43.75 79 98.75
左右 1 1.25 80 100.00
Frequency Missing = 1
Thursday, January 11, 2018 11:25:25 AM 61
The FREQ Procedure
Table of tou by da
tou da
Frequency|
Percent |
Row Pct |
Col Pct |右 |左 |左右 | Total
---------+--------+--------+--------+
右 | 44 | 18 | 1 | 63
| 55.00 | 22.50 | 1.25 | 78.75
| 69.84 | 28.57 | 1.59 |
| 100.00 | 51.43 | 100.00 |
---------+--------+--------+--------+
左 | 0 | 17 | 0 | 17
| 0.00 | 21.25 | 0.00 | 21.25
| 0.00 | 100.00 | 0.00 |
| 0.00 | 48.57 | 0.00 |
---------+--------+--------+--------+
Total 44 35 1 80
55.00 43.75 1.25 100.00
Frequency Missing = 1
Thursday, January 11, 2018 11:25:25 AM 62
shintyou Cum. Cum.
Midpoint Freq Freq Percent Percent
|
165 |* 1 1 1.23 1.23
|
170 |** 3 4 3.70 4.94
|
175 |******** 15 19 18.52 23.46
|
180 |******************* 37 56 45.68 69.14
|
185 |*********** 21 77 25.93 95.06
|
190 |* 2 79 2.47 97.53
|
195 |* 1 80 1.23 98.77
|
200 |* 1 81 1.23 100.00
|
-----+----+----+----
10 20 30
Frequency
Thursday, January 11, 2018 11:25:25 AM 63
Frequency
| *****
| *****
| *****
| *****
30 + *****
| *****
| *****
| *****
| ***** *****
20 + ***** *****
| ***** *****
| ***** ***** *****
| ***** ***** *****
| ***** ***** *****
10 + ***** ***** *****
| ***** ***** *****
| ***** ***** *****
| ***** ***** ***** *****
| ***** ***** ***** ***** ***** ***** ***** *****
--------------------------------------------------------------------
165 170 175 180 185 190 195 200
shintyou Midpoint
Thursday, January 11, 2018 11:25:25 AM 64
DoBY Cum. Cum.
Midpoint Freq Freq Percent Percent
|
1956 |* 1 1 1.23 1.23
|
1962 | 0 1 0.00 1.23
|
1968 | 0 1 0.00 1.23
|
1974 |* 1 2 1.23 2.47
|
1980 |*** 6 8 7.41 9.88
|
1986 |*********** 21 29 25.93 35.80
|
1992 |********************* 41 70 50.62 86.42
|
1998 |****** 11 81 13.58 100.00
|
-----+----+----+----+-
10 20 30 40
Frequency
Thursday, January 11, 2018 11:25:25 AM 65
Frequency
| *****
40 + *****
| *****
| *****
| *****
| *****
30 + *****
| *****
| *****
| *****
| ***** *****
20 + ***** *****
| ***** *****
| ***** *****
| ***** *****
| ***** ***** *****
10 + ***** ***** *****
| ***** ***** *****
| ***** ***** ***** *****
| ***** ***** ***** *****
| ***** ***** ***** ***** ***** *****
--------------------------------------------------------------------
1956 1962 1968 1974 1980 1986 1992 1998
DoBY Midpoint
<<< 中略:P66-68には以下のP72-74と同じ出力が得られる >>>
<<< 中略:P69-71には以下のP77-79と同じ出力が得られる >>>
Thursday, January 11, 2018 11:25:26 AM 72
The UNIVARIATE Procedure
Variable: shintyou
Moments
N 81 Sum Weights 81
Mean 180.506173 Sum Observations 14621
Std Deviation 5.05994925 Variance 25.6030864
Skewness 0.48774577 Kurtosis 2.40338502
Uncorrected SS 2641229 Corrected SS 2048.24691
Coeff Variation 2.80320012 Std Error Mean 0.56221658
Basic Statistical Measures
Location Variability
Mean 180.5062 Std Deviation 5.05995
Median 180.0000 Variance 25.60309
Mode 180.0000 Range 33.00000
Interquartile Range 5.00000
Thursday, January 11, 2018 11:25:26 AM 73
The UNIVARIATE Procedure
Variable: shintyou
Tests for Location: Mu0=0
Test -Statistic- -----p Value------
Student's t t 321.0616 Pr > |t| <.0001
Sign M 40.5 Pr >= |M| <.0001
Signed Rank S 1660.5 Pr >= |S| <.0001
Quantiles (Definition 5)
Level Quantile
100% Max 200
99% 200
95% 187
90% 186
75% Q3 183
50% Median 180
25% Q1 178
10% 174
5% 173
1% 167
0% Min 167
Thursday, January 11, 2018 11:25:26 AM 74
The UNIVARIATE Procedure
Variable: shintyou
Extreme Observations
----Lowest---- ----Highest---
Value Obs Value Obs
167 30 187 51
170 60 189 25
171 68 190 17
172 81 193 78
173 75 200 24
Thursday, January 11, 2018 11:25:26 AM 75
The UNIVARIATE Procedure
Variable: shintyou
Stem Leaf # Boxplot
200 0 1 *
198
196
194
192 0 1 0
190 0 1 |
188 0 1 |
186 00000 5 |
184 0000000 7 |
182 00000000000000000 17 +-----+
180 0000000000000000000 19 *--+--*
178 0000000000 10 +-----+
176 0000000 7 |
174 000000 6 |
172 000 3 |
170 00 2 0
168
166 0 1 0
----+----+----+----+
Thursday, January 11, 2018 11:25:26 AM 76
The UNIVARIATE Procedure
Variable: shintyou
Normal Probability Plot
201+ *
|
|
|
| * +++
| *++++
| ++*+
| +*****
| ++****
| *******
| ******
| ****++
| ****++
| ****++
| ***++
| *++*+
| ++++
167++*
+----+----+----+----+----+----+----+----+----+----+
-2 -1 0 +1 +2
Thursday, January 11, 2018 11:25:26 AM 77
The UNIVARIATE Procedure
Variable: DoBY
Moments
N 81 Sum Weights 81
Mean 1989.11111 Sum Observations 161118
Std Deviation 6.42261629 Variance 41.25
Skewness -2.4128665 Kurtosis 10.3141417
Uncorrected SS 320484904 Corrected SS 3300
Coeff Variation 0.32288876 Std Error Mean 0.71362403
Basic Statistical Measures
Location Variability
Mean 1989.111 Std Deviation 6.42262
Median 1991.000 Variance 41.25000
Mode 1991.000 Range 44.00000
Interquartile Range 6.00000
Thursday, January 11, 2018 11:25:26 AM 78
The UNIVARIATE Procedure
Variable: DoBY
Tests for Location: Mu0=0
Test -Statistic- -----p Value------
Student's t t 2787.338 Pr > |t| <.0001
Sign M 40.5 Pr >= |M| <.0001
Signed Rank S 1660.5 Pr >= |S| <.0001
Quantiles (Definition 5)
Level Quantile
100% Max 1998
99% 1998
95% 1996
90% 1995
75% Q3 1993
50% Median 1991
25% Q1 1987
10% 1983
5% 1978
1% 1954
0% Min 1954
Thursday, January 11, 2018 11:25:26 AM 79
The UNIVARIATE Procedure
Variable: DoBY
Extreme Observations
----Lowest---- ----Highest---
Value Obs Value Obs
1954 1 1996 78
1974 4 1996 79
1977 44 1997 2
1977 22 1998 26
1978 47 1998 50
Thursday, January 11, 2018 11:25:26 AM 80
The UNIVARIATE Procedure
Variable: DoBY
Stem Leaf # Boxplot
199 55556666788 11 |
199 0000111111111222222233333333444444 34 +-----+
198 5556677788888889999999 22 +--+--+
198 11344444 8 |
197 7788 4 0
197 4 1 0
196
196
195
195 4 1 *
----+----+----+----+----+----+----
Multiply Stem.Leaf by 10**+1
Thursday, January 11, 2018 11:25:26 AM 81
The UNIVARIATE Procedure
Variable: DoBY
Normal Probability Plot
1997.5+ +******** * *
| *************
| ********+++
1982.5+ *****+++
| +++****
|++++*
1967.5+
|
|
1952.5+ *
+----+----+----+----+----+----+----+----+----+----+
-2 -1 0 +1 +2
分布形状と統計量
/* Lesson 11-02 */
/* File Name = les1102.sas 01/11/18 */
options nocenter linesize=78 pagesize=30;
proc printto print = 'Kougi/les1102_Results.txt' new;
data pref;
infile 'Kougi/Prefectures.csv'
firstobs=5 dlm=',' dsd missover encoding=sjis termstr=crlf;
input Pref : $16. CCtiy : $16. Jinkou Menseki Mitsudo SJinkou
Shi Ku Tyou Son;
proc print data=pref(obs=20);
run;
proc means data=pref;
run;
proc chart data=pref;
hbar Jinkou Mitsudo;
vbar Jinkou Mitsudo;
run;
proc univariate data=pref plot;
var Jinkou Mitsudo;
run;
ods listing;
ods graphics off;
proc univariate data=pref plot;
var Jinkou Mitsudo;
run;