/* Lesson 13-1 */
/* File Name = les1301.sas 01/19/06 */
data gakusei;
infile 'all05be.prn'
firstobs=2;
input sex $ shintyou taijyuu kyoui
jitaku $ kodukai carryer $ tsuuwa;
proc print data=gakusei(obs=10);
run;
proc plot data=gakusei; : 散布図
plot shintyou*taijyuu; : 元の変量のプロット
run; :
proc princomp cov data=gakusei out=outprin; : 主成分分析(分散共分散行列)
var shintyou taijyuu; : 2変量
run; :
proc print data=outprin(obs=15); : 結果の出力
run; :
proc plot data=outprin; : 散布図
plot prin2*prin1/vref=0 href=0; : 主成分得点のプロット
run; :
proc sort data=outprin; : 説明のためにソートしてみる
by prin1; : 第一主成分で
run; :
proc print data=outprin; : 体重がややが効いていることの確認
run; :
SAS システム 2 15:12 Friday, January 13, 2006 プロット : SHINTYOU*TAIJYUU. 凡例: A = 1 OBS, B = 2 OBS, ... (NOTE: 46 オブザベーションが欠損値です.) SHINTYOU | 200 + | | A B A 180 + A BECFDDBEA B B A A A | CAGELHTPJHCFECB BB | AFAGIIFECBDEAA AA A A 160 + ADDGDIEDBABB | A ECAEDDA A A | A BAA 140 + ---+-----------+-----------+-----------+-----------+-- 20 40 60 80 100 TAIJYUU SAS システム 3 15:12 Friday, January 13, 2006 Principal Component Analysis 297 Observations 2 Variables Simple Statistics SHINTYOU TAIJYUU Mean 168.5636364 58.52457912 StD 8.1049800 9.31737163 SAS システム 4 15:12 Friday, January 13, 2006 Principal Component Analysis Covariance Matrix SHINTYOU TAIJYUU SHINTYOU 65.69070025 53.24049140 TAIJYUU 53.24049140 86.81341410 Total Variance = 152.50411434 Eigenvalues of the Covariance Matrix Eigenvalue Difference Proportion Cumulative PRIN1 130.530 108.556 0.855911 0.85591 PRIN2 21.974 . 0.144089 1.00000 SAS システム 5 15:12 Friday, January 13, 2006 Principal Component Analysis Eigenvectors PRIN1 PRIN2 SHINTYOU 0.634595 0.772845 TAIJYUU 0.772845 -.634595 SAS システム 6 15:12 Friday, January 13, 2006 S H T K C I A J O A T N I K I D R S P P T J Y T U R U R R O S Y Y O A K Y U I I B E O U U K A E W N N S X U U I U I R A 1 2 1 F 145.0 38.0 . J 10000 . -30.8157 -5.18625 2 F 146.7 41.0 85 J 10000 Vodafone 6000 -27.4183 -5.77620 <中略> SAS システム 8 15:12 Friday, January 13, 2006 プロット : PRIN2*PRIN1. 凡例: A = 1 OBS, B = 2 OBS, ... (NOTE: 46 オブザベーションが欠損値です.) PRIN2 | | 10 + A B A BA | AA C BCACACBAEA AAA | A BBABCAEDDFGDBMFCABE A 0 +-------------B-CABADC-JABEDEEFDGDACDBADA-A-----A--------- | AAABBCCE ACCA A B CFEB EC A A | A AAAAA AA AA A AABBABAB A BA AACA -10 + A A | C A | | | | A A -20 + | A A ---+------------+------------+------------+------------+-- -40 -20 0 20 40 PRIN1 SAS システム 9 15:12 Friday, January 13, 2006 S H T K C I A J O A T N I K I D R S P P T J Y T U R U R R O S Y Y O A K Y U I I B E O U U K A E W N N S X U U I U I R A 1 2 1 F 148.9 . . J 60000 . . . 2 F 153.0 . . G 120000 DoCoMo 200 . . 3 F 155.0 . . J 20000 . . . <中略> SAS システム 46 15:12 Friday, January 13, 2006 S H T K C I A J O A T N I K I D R S P P T J Y T U R U R R O S Y Y O A K Y U I I B E O U U K A E W N N S X U U I U I R A 1 2 310 M 168.0 75.0 . G 150000 . 12.3753 -10.8908 311 M 179.0 66.0 . 30000 . 12.4002 3.3218 312 M 173.0 71.0 100 G 0 . 12.4569 -4.4882 313 M 178.0 67.0 . J 0 . 12.5385 1.9144 314 M 172.0 72.0 89 G 150000 . 12.5951 -5.8957 315 M 177.0 68.0 . G 80000 . 12.6767 0.5070 316 M 182.0 64.0 . G 0 . 12.7583 6.9096 317 M 165.0 78.0 . G 0 2098 12.7900 -15.1131 318 M 170.0 74.0 90 J 0 . 12.8716 -8.7105 319 M 175.0 70.0 95 G 50000 8000 12.9532 -2.3079 320 M 178.0 68.0 . J 100000 DoCoMo 4000 13.3113 1.2798 321 M 184.0 65.0 . G 140000 au 10000 14.8003 7.8206 322 M 175.0 74.0 . J 0 . 16.0446 -4.8463 323 M 180.0 70.0 94 G 70000 au 5000 16.1262 1.5563 324 M 180.0 70.0 . J 40000 au 4000 16.1262 1.5563 325 M 180.0 70.0 . . . 16.1262 1.5563 326 M 180.0 70.0 . J 40000 DoCoMo 6500 16.1262 1.5563 327 M 178.7 71.2 95 0 . 16.2286 -0.2099 328 M 173.5 76.5 . G 100000 . 17.0248 -7.5921 329 M 184.0 68.0 85 30000 . 17.1189 5.9169 330 M 182.0 70.0 90 G 100000 . 17.3954 3.1020 331 M 185.0 68.0 93 J 0 . 17.7535 6.6897 332 M 175.0 77.0 95 G 130000 . 18.3631 -6.7501 333 M 179.1 74.2 . 0 au 4000 18.8010 -1.8046 334 M 175.0 79.0 . J 0 No 0 19.9088 -8.0193 335 M 176.5 78.0 96 J 10000 . 20.0879 -6.2254 336 M 177.0 78.0 . J 40000 . 20.4052 -5.8390 337 M 181.5 74.5 . G 120000 au 3000 20.5559 -0.1401 338 M 178.0 78.0 110 G 50000 . 21.0398 -5.0661 339 M 169.3 88.5 94 J 0 . 23.6337 -18.4531 340 M 186.0 82.0 . J 0 . 29.2079 -1.4218 341 M 182.0 90.0 100 J 40000 . 32.8523 -9.5899 342 M 178.0 95.0 . 1000 No . 34.1781 -15.8543 343 M 178.0 100.0 112 G 60000 . 38.0423 -19.0272
/* Lesson 13-2 */ /* File Name = les1302.sas 01/19/06 */ data gakusei; infile 'all05be.prn' firstobs=2; input sex $ shintyou taijyuu kyoui jitaku $ kodukai carryer $ tsuuwa; proc print data=gakusei(obs=10); run; proc princomp cov data=gakusei out=outprin; : 主成分分析(分散共分散行列) var shintyou taijyuu kyoui; : 3変量 run; : proc print data=outprin(obs=15); : 結果の出力 run; : proc plot data=outprin; : 散布図 plot prin2*prin1/vref=0 href=0; : 主成分得点のプロット plot prin3*prin2/vref=0 href=0; : plot prin3*prin1/vref=0 href=0; : run; :
SAS システム 3
15:12 Friday, January 13, 2006
Principal Component Analysis
108 Observations
3 Variables
Simple Statistics
SHINTYOU TAIJYUU KYOUI
Mean 167.2231481 58.51296296 86.41666667
StD 8.7885398 10.76789181 7.54069954
SAS システム 4
15:12 Friday, January 13, 2006
Principal Component Analysis
Covariance Matrix
SHINTYOU TAIJYUU KYOUI
SHINTYOU 77.2384311 69.6319401 25.1687695
TAIJYUU 69.6319401 115.9474939 53.7151090
KYOUI 25.1687695 53.7151090 56.8621495
SAS システム 5
15:12 Friday, January 13, 2006
Principal Component Analysis
Total Variance = 250.04807459
Eigenvalues of the Covariance Matrix
Eigenvalue Difference Proportion Cumulative
PRIN1 193.624 153.258 0.774349 0.77435
PRIN2 40.366 24.309 0.161434 0.93578
PRIN3 16.058 . 0.064218 1.00000
SAS システム 6
15:12 Friday, January 13, 2006
Principal Component Analysis
Eigenvectors
PRIN1 PRIN2 PRIN3
SHINTYOU 0.533237 -.663184 0.525210
TAIJYUU 0.749414 0.082283 -.656969
KYOUI 0.392475 0.743919 0.540876
SAS システム 7
15:12 Friday, January 13, 2006
S
H T K C
I A J O A T
N I K I D R S P P P
T J Y T U R U R R R
O S Y Y O A K Y U I I I
B E O U U K A E W N N N
S X U U I U I R A 1 2 3
1 F 145.0 38.0 . J 10000 . . . .
2 F 146.7 41.0 85 J 10000 Vodafone 6000 -24.6242 11.1157 -0.03973
3 F 148.0 42.0 . J 50000 . . . .
4 F 148.0 43.0 80 J 50000 DoCoMo 4000 -24.3945 6.6985 -3.37527
5 F 148.9 . . J 60000 . . . .
6 F 149.0 45.0 . G 60000 . . . .
7 F 150.0 46.0 86 40000 . -18.7249 10.0825 -1.05050
8 F 151.0 50.0 . G 60000 J-PHONE . . . .
9 F 151.7 41.5 80 J 35000 . -23.5457 4.1213 -0.44655
10 F 152.0 35.0 77 J 60000 DoCoMo 2000 -29.4343 1.1558 2.35868
<中略>
SAS システム 9
15:12 Friday, January 13, 2006
プロット : PRIN2*PRIN1. 凡例: A = 1 OBS, B = 2 OBS, ...
(NOTE: 235 オブザベーションが欠損値です.)
PRIN2 | |
20 + |
| A A A | A A A A
| AA AAACBB ABA| A AC A A
0 +--------A----A-ABAEABD-CAAACADE-BACAAB-AA-------A----------------
| A AA B B AAB BABD AA
| A | B A A
-20 + |
| A |
| |
-40 + |
---+-----------+-----------+-----------+-----------+-----------+--
-40 -20 0 20 40 60
PRIN1
SAS システム 10
15:12 Friday, January 13, 2006
プロット : PRIN3*PRIN2. 凡例: A = 1 OBS, B = 2 OBS, ...
(NOTE: 235 オブザベーションが欠損値です.)
PRIN3 | |
10 + |
| A A AB A| A A A
| A A A BA D EACC AABD A A
0 +----------------------A---A-A--B--A--FBEAEBCBA-A-AAA-------------
| A A A A CBA CBA AA
| A A AA A | B A A
-10 + | A
| A | A
| |
-20 + |
---+-----------+-----------+-----------+-----------+-----------+--
-30 -20 -10 0 10 20
PRIN2
SAS システム 11
15:12 Friday, January 13, 2006
プロット : PRIN3*PRIN1. 凡例: A = 1 OBS, B = 2 OBS, ...
(NOTE: 235 オブザベーションが欠損値です.)
PRIN3 | |
10 + |
| A |AA AC A A
| A A BACA BAA|ACCCBAA A D
0 +-----------AA--BCBDBBC-AAAABB-FA-A--A----------------------------
| A A AA A BA | B BAC AA
| A A AA A A A A A
-10 + A |
| A | A
| |
-20 + |
---+-----------+-----------+-----------+-----------+-----------+--
-40 -20 0 20 40 60
PRIN1
/* Lesson 13-3 */
/* File Name = les1303.sas 01/19/06 */
data gakusei;
infile 'all05be.prn'
firstobs=2;
input sex $ shintyou taijyuu kyoui
jitaku $ kodukai carryer $ tsuuwa;
proc print data=gakusei(obs=10);
run; :
proc princomp data=gakusei out=outprin; : 相関係数を使って
var shintyou taijyuu kyoui; :
run; :
proc print data=outprin(obs=15);
run;
proc plot data=outprin;
plot prin2*prin1/vref=0 href=0;
plot prin3*prin2/vref=0 href=0;
plot prin3*prin1/vref=0 href=0;
run;
SAS システム 3
15:12 Friday, January 13, 2006
Principal Component Analysis
108 Observations
3 Variables
Simple Statistics
SHINTYOU TAIJYUU KYOUI
Mean 167.2231481 58.51296296 86.41666667
StD 8.7885398 10.76789181 7.54069954
SAS システム 4
15:12 Friday, January 13, 2006
Principal Component Analysis
Correlation Matrix
SHINTYOU TAIJYUU KYOUI
SHINTYOU 1.0000 0.7358 0.3798
TAIJYUU 0.7358 1.0000 0.6615
KYOUI 0.3798 0.6615 1.0000
SAS システム 5
15:12 Friday, January 13, 2006
Principal Component Analysis
Eigenvalues of the Correlation Matrix
Eigenvalue Difference Proportion Cumulative
PRIN1 2.19675 1.57351 0.732250 0.73225
PRIN2 0.62324 0.44322 0.207746 0.94000
PRIN3 0.18001 . 0.060004 1.00000
SAS システム 6
15:12 Friday, January 13, 2006
Principal Component Analysis
Eigenvectors
PRIN1 PRIN2 PRIN3
SHINTYOU 0.559835 -.651126 0.512465
TAIJYUU 0.637080 -.057219 -.768671
KYOUI 0.529824 0.756810 0.382786
[注意] データによっては解釈が困難なことも有り得る。
[参考] 「J:\コンピュータによる統計解析6(林 篤裕)\」に以下のデータを置いておく。
[注意] データによっては解釈が困難なことも有り得る。
また、自分の思い付かない結果を含んでいることもある。
なお、今まで紹介していた私のメールアドレスは実は講義用のものであった。 今後、もし統計に関して何か疑問に出会い、私に連絡・相談してみたいと思った時は、 以下のアドレスを使ってください。
皆さんの期待に応えられたか心許無い部分もありますが、半年間ご苦労様でした。