2018年10月26日 金曜日 18時16分21秒 52 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時16分21秒 53 REG プロシジャ モデル : MODEL1 従属変数 : taijyuu 読み込んだオブザベーション数 325 使用されたオブザベーション数 114 欠損値を含むオブザベーション数 211 分散分析 要因 自由度 平方和 平均平方 F 値 Pr > F Model 2 8070.70705 4035.35353 85.10 <.0001 Error 111 5263.40733 47.41808 Corrected Total 113 13334 Root MSE 6.88608 R2 乗 0.6053 従属変数の平均 58.79298 調整済み R2 乗 0.5982 変動係数 11.71242 2018年10月26日 金曜日 18時16分21秒 54 REG プロシジャ モデル : MODEL1 従属変数 : taijyuu パラメータの推定 パラメータ 変数 自由度 推定値 標準誤差 t 値 Pr > |t| Intercept 1 -106.30023 12.75197 -8.34 <.0001 shintyou 1 0.80655 0.07854 10.27 <.0001 kyoui 1 0.34947 0.08192 4.27 <.0001 2018年10月26日 金曜日 18時16分22秒 55 Obs sex shintyou taijyuu kyoui jitaku kodukai carryer tsuuwa pred1 resid1 1 F 145.0 38.0 . J 10000 . . . 2 F 146.7 41.0 85 J 10000 Vodafone 6000 41.7256 -0.72559 3 F 148.0 42.0 . J 50000 . . . 4 F 148.0 43.0 80 J 50000 DoCoMo 4000 41.0267 1.97328 5 F 149.0 45.0 . G 60000 . . . 6 F 150.0 46.0 86 40000 . 44.7367 1.26333 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 . 44.0109 -2.51095 10 F 152.0 35.0 77 J 60000 DoCoMo 2000 43.2045 -8.20449 11 F 152.0 43.0 . J 20000 au 3500 . . 12 F 152.0 44.0 . 45000 DoCoMo 4000 . . 13 F 153.0 41.0 . J 125000 No . . . 14 F 153.0 42.0 . G 0 Vodafone 1000 . . 15 F 153.0 46.5 87 G 10000 . 47.5058 -1.00578 2018年10月26日 金曜日 18時16分22秒 56 : taijyuu*shintyou. A=1, B=2, ... 100 + A | | A | A | A | A A 75 + A A | A A BAA A A A A | BB A A BAAA A A A taijyuu | A A A CA C AA B | B A AAAA C AA AA A A A A | A A AA A BABA AAA A 50 + AC B CD AAB B | A A A B A A BA A | AA A A | | A | 25 + --+------------+------------+------------+------------+------------+-- 140 150 160 170 180 190 shintyou 2018年10月26日 金曜日 18時16分22秒 57 : taijyuu*kyoui. A=1, B=2, ... 100 + A | | A | A | A | A A 75 + A A | BD AB A | A A CAAAAB B A taijyuu | A A AGA B A | BB ADBDB A | A A A AAA CC AB 50 + A A DD G B | B DAA BA | B B | | A | 25 + ---+-------------+-------------+-------------+-------------+-- 40 60 80 100 120 kyoui 2018年10月26日 金曜日 18時16分22秒 58 : taijyuu*pred1. A=1, B=2, ... 100 + A | | A | A | A | A A 75 + A A | A AB A B AAA | A B AA ABAA A A A A taijyuu | A A A BABB AA B | AAA A A B BAAAA B B A | A A A AA A CABA A A 50 + A B B CBAABAAAB | AA C AA ABA | AA A A | | A | 25 + ---+-------------+-------------+-------------+-------------+-- 40 50 60 70 80 予測値 taijyuu 2018年10月26日 金曜日 18時16分22秒 59 : resid1*pred1. A=1, B=2, ... 40 + | | A | | A | A 残 20 + 差 | | A A | A A A A A | A A A A A BA | A A A A A A A B A BA B A 0 +----A---AA---B-B-AAA---A-C--A-A-ABAAB-AAA----A--------------- | A B AABBAABAAAAA A A AAAB AA A AAB | A A AB B B A A A B B A | A A A | A | -20 + ---+-------------+-------------+-------------+-------------+-- 40 50 60 70 80 予測値 taijyuu 2018年10月26日 金曜日 18時16分22秒 60 : resid1*shintyou. A=1, B=2, ... 40 + | | A | | A | A 残 20 + 差 | | A A | A A A AA | A A B A A A A | A BA A A BA A BA AA A 0 +----------A---A---A-AAA-A--BA--A--AAAA-AEA-A--AA-A-A----------------- | A B B DC AAB BA A B BAA A A B A A | A A AA CB A A A ABA A | A A A | A | -20 + --+------------+------------+------------+------------+------------+-- 140 150 160 170 180 190 shintyou 2018年10月26日 金曜日 18時16分22秒 61 : resid1*kyoui. A=1, B=2, ... 40 + | | A | | A | A 残 20 + 差 | | A A | A AA B | A AA A BB | A BAAAB AAA B B A 0 +-----------------------A---C--A-ABFBBE--AB-----A------------- | EBCAGAEDA AB | AAABBA D AC A | BA | A | -20 + ---+-------------+-------------+-------------+-------------+-- 40 60 80 100 120 kyoui 2018年10月26日 金曜日 18時16分22秒 62 : resid1*taijyuu. A=1, B=2, ... 40 + | | A | | A | A 残 20 + 差 | | A A | A A A AA | A AA A BB | A AA A A B CAABA A 0 +-----------------A-AAACA-B-C-ACBD-C---A---------------------- | A CADCD AB BCBA B AB | A A ABB B A AC B A | A A A | A | -20 + ---+-------------+-------------+-------------+-------------+-- 20 40 60 80 100 taijyuu 2018年10月26日 金曜日 18時16分22秒 63 UNIVARIATE プロシジャ 変数 : resid1 (残差) モーメント N 114 重み変数の合計 114 平均 0 合計 0 標準偏差 6.8248682 分散 46.578826 歪度 2.02681339 尖度 7.21141767 無修正平方和 5263.40733 修正済平方和 5263.40733 変動係数 . 平均の標準誤差 0.63920747 基本統計量 位置 ばらつき 平均 0.00000 標準偏差 6.82487 中央値 -1.33487 分散 46.57883 最頻値 -8.24352 範囲 47.63020 四分位範囲 5.97826 Note: 4個の最頻値があります(度数: 2)。表では最頻値のなかで最も小さな値を表示します。 2018年10月26日 金曜日 18時16分22秒 64 UNIVARIATE プロシジャ 変数 : resid1 (残差) 位置の検定 H0: Mu0=0 検定 -統計量- ------p 値------- Student の t 検定 t 0 Pr > |t| 1.0000 符号検定 M -14 Pr >= |M| 0.0111 符号付順位検定 S -517.5 Pr >= |S| 0.1442 正規性の検定 検定 --統計量--- ------p 値------- Shapiro-Wilk W 0.852435 Pr < W <0.0001 Kolmogorov-Smirnov D 0.132665 Pr > D <0.0100 Cramer-von Mises W-Sq 0.563789 Pr > W-Sq <0.0050 Anderson-Darling A-Sq 3.36241 Pr > A-Sq <0.0050 2018年10月26日 金曜日 18時16分23秒 65 UNIVARIATE プロシジャ 変数 : resid1 (残差) 分位点 (定義 5) 水準 分位点 100% 最大値 34.69175 99% 25.40118 95% 10.63949 90% 7.54104 75% Q3 2.34154 50% 中央値 -1.33487 25% Q1 -3.63672 10% -7.12493 5% -8.24352 1% -10.89564 0% 最小値 -12.93845 2018年10月26日 金曜日 18時16分23秒 66 UNIVARIATE プロシジャ 変数 : resid1 (残差) 極値 -----最小値---- ----最大値---- 値 Obs 値 Obs -12.93845 312 11.7345 171 -10.89564 316 14.5612 319 -8.56575 65 23.5937 293 -8.35741 266 25.4012 150 -8.24352 100 34.6917 233 欠損値 ---パーセント--- 欠損値 カウント 全体 欠損値 . 211 64.92 100.00