2018年10月26日 金曜日 18時15分46秒 19 Obs sex shintyou taijyuu kyoui jitaku kodukai carryer tsuuwa 1 F 145.0 38.0 . J 10000 . 2 F 146.7 41.0 85 J 10000 Vodafone 6000 3 F 148.0 42.0 . J 50000 . 4 F 148.0 43.0 80 J 50000 DoCoMo 4000 5 F 149.0 45.0 . G 60000 . 6 F 150.0 46.0 86 40000 . 7 F 151.0 45.0 . J 20000 docomo 5000 8 F 151.0 50.0 . G 60000 J-PHONE . 9 F 151.7 41.5 80 J 35000 . 10 F 152.0 35.0 77 J 60000 DoCoMo 2000 2018年10月26日 金曜日 18時15分46秒 20 CORR プロシジャ 5 変数 : shintyou taijyuu kyoui kodukai tsuuwa 単純統計量 変数 N 平均 標準偏差 合計 shintyou 319 168.52508 7.98682 53760 taijyuu 319 58.19530 8.34648 18564 kyoui 110 86.10909 7.01812 9472 kodukai 301 48271 48841 14529500 tsuuwa 132 6742 4470 890002 単純統計量 変数 最小値 最大値 shintyou 145.00000 185.00000 taijyuu 35.00000 79.00000 kyoui 56.00000 110.00000 kodukai 0 300000 tsuuwa 0 30000 2018年10月26日 金曜日 18時15分46秒 21 CORR プロシジャ Pearson の相関係数 H0: Rho=0 に対する Prob > |r| オブザベーション数 shintyou taijyuu kyoui kodukai tsuuwa shintyou 1.00000 0.72972 0.36345 0.06590 -0.05960 <.0001 <.0001 0.2544 0.4972 319 319 110 301 132 taijyuu 0.72972 1.00000 0.60461 0.01665 -0.04543 <.0001 <.0001 0.7736 0.6050 319 319 110 301 132 kyoui 0.36345 0.60461 1.00000 -0.08466 -0.17722 <.0001 <.0001 0.3882 0.2940 110 110 110 106 37 kodukai 0.06590 0.01665 -0.08466 1.00000 0.26949 0.2544 0.7736 0.3882 0.0021 301 301 106 301 128 tsuuwa -0.05960 -0.04543 -0.17722 0.26949 1.00000 0.4972 0.6050 0.2940 0.0021 132 132 37 128 132 2018年10月26日 金曜日 18時15分46秒 22 REG プロシジャ モデル : MODEL1 従属変数 : taijyuu 読み込んだオブザベーション数 319 使用されたオブザベーション数 319 分散分析 要因 自由度 平方和 平均平方 F 値 Pr > F Model 1 11796 11796 361.06 <.0001 Error 317 10357 32.67126 Corrected Total 318 22153 Root MSE 5.71588 R2 乗 0.5325 従属変数の平均 58.19530 調整済み R2 乗 0.5310 変動係数 9.82189 2018年10月26日 金曜日 18時15分46秒 23 REG プロシジャ モデル : MODEL1 従属変数 : taijyuu パラメータの推定 パラメータ 変数 自由度 推定値 標準誤差 t 値 Pr > |t| Intercept 1 -70.31845 6.77088 -10.39 <.0001 shintyou 1 0.76258 0.04013 19.00 <.0001 2018年10月26日 金曜日 18時15分48秒 24 Obs sex shintyou taijyuu kyoui jitaku kodukai carryer tsuuwa pred1 resid1 1 F 145.0 38.0 . J 10000 . 40.2556 -2.2556 2 F 146.7 41.0 85 J 10000 Vodafone 6000 41.5519 -0.5519 3 F 148.0 42.0 . J 50000 . 42.5433 -0.5433 4 F 148.0 43.0 80 J 50000 DoCoMo 4000 42.5433 0.4567 5 F 149.0 45.0 . G 60000 . 43.3059 1.6941 6 F 150.0 46.0 86 40000 . 44.0685 1.9315 7 F 151.0 45.0 . J 20000 docomo 5000 44.8310 0.1690 8 F 151.0 50.0 . G 60000 J-PHONE . 44.8310 5.1690 9 F 151.7 41.5 80 J 35000 . 45.3648 -3.8648 10 F 152.0 35.0 77 J 60000 DoCoMo 2000 45.5936 -10.5936 11 F 152.0 43.0 . J 20000 au 3500 45.5936 -2.5936 12 F 152.0 44.0 . 45000 DoCoMo 4000 45.5936 -1.5936 13 F 153.0 41.0 . J 125000 No . 46.3562 -5.3562 14 F 153.0 42.0 . G 0 Vodafone 1000 46.3562 -4.3562 15 F 153.0 46.5 87 G 10000 . 46.3562 0.1438 2018年10月26日 金曜日 18時15分48秒 25 : taijyuu*shintyou. A=1, B=2, ... taijyuu | 100 + | | | 80 + A A AAAA | B A A B A A | AAB EDB C AAA F A A B | B B BC AGC BBB DDCB CEA A A 60 + A AA BB AACCABIFBHIFACCDC AAA AA | A B EC AA ECHGACDF EBB DB A | A A CDAC DE DDB CFBAA AA A A | A AAA BABA BAB AB ABA A 40 + A AA A B B | A | | 20 + | --+------------+------------+------------+------------+------------+-- 140 150 160 170 180 190 shintyou 2018年10月26日 金曜日 18時15分48秒 26 : pred1*taijyuu. A=1, B=2, ... 70 + A A B A 予 | B AA A A 測 | A AA BB CA F A A 値 | A A ACA A ABBB C | A B A BA DA ABA BABB A A A A t | A BB BBAC AIDAEBB B D BBA B A a 60 + A AAA BEAD H BCDAACA A AA A A i | AAAB DACA B CA BA A A A j | AD A CBC CBAB B B y | C B AA C A A B A A u | AAB C EA A B B u | B AB CC ACA C A AA 50 + A AA A BA A A | B AA C ABA A | A A A C A A | A A AA A A | AA AA | A 40 + A ---+------------+------------+------------+------------+------------+-- 30 40 50 60 70 80 taijyuu 2018年10月26日 金曜日 18時15分48秒 27 : resid1*pred1. A=1, B=2, ... 40 + | | | | | A 残 20 + A 差 | B A | A A AA AAA | A A A AB AABDA A | A A AA ABA BCABAACBA AA A | AA A A A CA AAABAGDABAEB BA F 0 +----A-B--AABABDACAAB-DA-DBEDBHCGICBB-BBB--A--A--------------- | A BAA AABCF ADBACEECDDDCDCBAC AEAA BA | A BA B ABA CDAB ABCA BBA CAAA A | A AA AA BB A | A A | -20 + ---+-------------+-------------+-------------+-------------+-- 40 50 60 70 80 予測値 taijyuu 2018年10月26日 金曜日 18時15分48秒 28 : resid1*shintyou. A=1, B=2, ... 40 + | | | | | A 残 20 + A 差 | B A | A A A AAAA | A A A A B AAB DA A | A A B ABA BC ABA ACB A AA A | AA A A A CA AAAB AFE BAA EBBA F 0 +----------AB---A-AAB-BDAC-AAB-DA--DBECBAHCAFICAABBB-B---A---A-------- | A AAA A B B DE ADB ACEE CDDADBD CBD AEA A B A | A BA B ABA CDAAA AB CA BC C AAA A | A AA A A B B A | A A | -20 + --+------------+------------+------------+------------+------------+-- 140 150 160 170 180 190 shintyou 2018年10月26日 金曜日 18時15分48秒 29 : resid1*taijyuu. A=1, B=2, ... 40 + | | | | | A 残 20 + A 差 | A A A | A A AAC | A A A B A A ACB B A | A B AB BA FB D BA A AA | AA AA AA C A B AFDCAEACD F 0 +---------------A-AAA-AC-D-CC--DDA-FBDF-BNDBFDA-BCA---A--A------------ | A AAAA BBAB CCE DDE GHBD AIB ACC C B | A BA CABF A ACCA CABB ACB | A AA B BB A | A A | -20 + --+------------+------------+------------+------------+------------+-- 30 40 50 60 70 80 taijyuu 2018年10月26日 金曜日 18時15分48秒 30 : resid1*pred1. A=1, B=2, ... 40 + | | | | | A 残 20 + A 差 | B A | A A AA AAA | A A A AB AABDA A | A A AA ABA BCABAACBA AA A | AA A A A CA AAABAGDABAEB BA F 0 +----A-B--AABABDACAAB-DA-DBEDBHCGICBB-BBB--A--A--------------- | A BAA AABCF ADBACEECDDDCDCBAC AEAA BA | A BA B ABA CDAB ABCA BBA CAAA A | A AA AA BB A | A A | -20 + ---+-------------+-------------+-------------+-------------+-- 40 50 60 70 80 予測値 taijyuu 2018年10月26日 金曜日 18時15分48秒 31 : resid1*shintyou. A=1, B=2, ... 40 + | | | | | A 残 20 + A 差 | B A | A A A AAAA | A A A A B AAB DA A | A A B ABA BC ABA ACB A AA A | AA A A A CA AAAB AFE BAA EBBA F 0 +----------AB---A-AAB-BDAC-AAB-DA--DBECBAHCAFICAABBB-B---A---A-------- | A AAA A B B DE ADB ACEE CDDADBD CBD AEA A B A | A BA B ABA CDAAA AB CA BC C AAA A | A AA A A B B A | A A | -20 + --+------------+------------+------------+------------+------------+-- 140 150 160 170 180 190 shintyou 2018年10月26日 金曜日 18時15分48秒 32 : resid1*taijyuu. A=1, B=2, ... 40 + | | | | | A 残 20 + A 差 | A A A | A A AAC | A A A B A A ACB B A | A B AB BA FB D BA A AA | AA AA AA C A B AFDCAEACD F 0 +---------------A-AAA-AC-D-CC--DDA-FBDF-BNDBFDA-BCA---A--A------------ | A AAAA BBAB CCE DDE GHBD AIB ACC C B | A BA CABF A ACCA CABB ACB | A AA B BB A | A A | -20 + --+------------+------------+------------+------------+------------+-- 30 40 50 60 70 80 taijyuu 2018年10月26日 金曜日 18時15分48秒 33 UNIVARIATE プロシジャ 変数 : resid1 (残差) モーメント N 319 重み変数の合計 319 平均 0 合計 0 標準偏差 5.70688374 分散 32.568522 歪度 0.7505943 尖度 0.92289711 無修正平方和 10356.79 修正済平方和 10356.79 変動係数 . 平均の標準誤差 0.31952415 基本統計量 位置 ばらつき 平均 0.00000 標準偏差 5.70688 中央値 -0.79488 分散 32.56852 最頻値 -1.60778 範囲 35.62579 四分位範囲 6.91353 Note: 3個の最頻値があります(度数: 5)。表では最頻値のなかで最も小さな値を表示します。 2018年10月26日 金曜日 18時15分48秒 34 UNIVARIATE プロシジャ 変数 : resid1 (残差) 位置の検定 H0: Mu0=0 検定 -統計量- ------p 値------- Student の t 検定 t 0 Pr > |t| 1.0000 符号検定 M -18.5 Pr >= |M| 0.0437 符号付順位検定 S -2130 Pr >= |S| 0.1968 正規性の検定 検定 --統計量--- ------p 値------- Shapiro-Wilk W 0.968069 Pr < W <0.0001 Kolmogorov-Smirnov D 0.080749 Pr > D <0.0100 Cramer-von Mises W-Sq 0.472513 Pr > W-Sq <0.0050 Anderson-Darling A-Sq 2.699726 Pr > A-Sq <0.0050 2018年10月26日 金曜日 18時15分48秒 35 UNIVARIATE プロシジャ 変数 : resid1 (残差) 分位点 (定義 5) 水準 分位点 100% 最大値 22.492854 99% 16.205116 95% 10.867060 90% 7.940485 75% Q3 3.118648 50% 中央値 -0.794884 25% Q1 -3.794884 10% -6.694249 5% -8.269725 1% -10.593614 0% 最小値 -13.132940 2018年10月26日 金曜日 18時15分48秒 36 UNIVARIATE プロシジャ 変数 : resid1 (残差) 極値 ----最小値---- ----最大値---- 値 Obs 値 Obs -13.1329 246 15.8671 263 -12.3271 308 16.2051 138 -10.8452 194 17.2051 139 -10.5936 10 18.6800 171 -10.1329 248 22.4929 97