2019年10月29日 火曜日 09時02分06秒 27 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 148.9 . . J 60000 . 6 F 149.0 45.0 . G 60000 . 7 F 150.0 46.0 86 40000 . 8 F 151.0 45.0 . J 20000 docomo 5000 9 F 151.0 50.0 . G 60000 J-PHONE . 10 F 151.7 41.5 80 J 35000 . 2019年10月29日 火曜日 09時02分06秒 28 REG プロシジャ モデル : MODEL1 従属変数 : taijyuu 読み込んだオブザベーション数 399 使用されたオブザベーション数 350 欠損値を含むオブザベーション数 49 分散分析 要因 自由度 平方和 平均平方 F 値 Pr > F Model 1 14602 14602 325.81 <.0001 Error 348 15596 44.81605 Corrected Total 349 30198 Root MSE 6.69448 R2 乗 0.4835 従属変数の平均 58.95800 調整済み R2 乗 0.4821 変動係数 11.35466 2019年10月29日 火曜日 09時02分06秒 29 REG プロシジャ モデル : MODEL1 従属変数 : taijyuu パラメータの推定 パラメータ 変数 自由度 推定値 標準誤差 t 値 Pr > |t| Intercept 1 -78.49464 7.62341 -10.30 <.0001 shintyou 1 0.81370 0.04508 18.05 <.0001 2019年10月29日 火曜日 09時02分07秒 30 Obs sex shintyou taijyuu kyoui jitaku kodukai carryer tsuuwa pred1 resid1 1 F 145.0 38.0 . J 10000 . 39.4924 -1.4924 2 F 146.7 41.0 85 J 10000 Vodafone 6000 40.8757 0.1243 3 F 148.0 42.0 . J 50000 . 41.9335 0.0665 4 F 148.0 43.0 80 J 50000 DoCoMo 4000 41.9335 1.0665 5 F 148.9 . . J 60000 . 42.6658 . 6 F 149.0 45.0 . G 60000 . 42.7472 2.2528 7 F 150.0 46.0 86 40000 . 43.5609 2.4391 8 F 151.0 45.0 . J 20000 docomo 5000 44.3746 0.6254 9 F 151.0 50.0 . G 60000 J-PHONE . 44.3746 5.6254 10 F 151.7 41.5 80 J 35000 . 44.9442 -3.4442 11 F 152.0 35.0 77 J 60000 DoCoMo 2000 45.1883 -10.1883 12 F 152.0 43.0 . J 20000 au 3500 45.1883 -2.1883 13 F 152.0 44.0 . 45000 DoCoMo 4000 45.1883 -1.1883 14 F 153.0 41.0 . J 125000 No . 46.0020 -5.0020 15 F 153.0 42.0 . G 0 Vodafone 1000 46.0020 -4.0020 2019年10月29日 火曜日 09時02分07秒 31 Plot of taijyuu*shintyou. Legend: A = 1 obs, B = 2 obs, etc. taijyuu | 100 + A | A | A A | A A 80 + A A AAAA A | B A A B A A | ACB EEB C BBA F A A B | B B BC AHC BBB FDCB CEA B A 60 + A AA BB AACDABJFBHIFADCDC ACA AA | A B EC CA FCIHACDF EBB EB A | A A CDAC DE DEB CFCBA AA A A A | A AAA BABA BAB AB ABA A 40 + A AA A B B | A | | 20 + | --+------------+------------+------------+------------+------------+-- 140 150 160 170 180 190 shintyou NOTE: 49 obs が欠損値です。 2019年10月29日 火曜日 09時02分07秒 32 Plot of pred1*taijyuu. Legend: A = 1 obs, B = 2 obs, etc. 80 + 予 | 測 | 値 | A | A A B A t | A DABCC G A A a | A BB A CAA B A A A A A i | A BABBBFCCAFDCB A AC j | A BBDBCIJBBB ECD A A y 60 + AAABHCDIFDBF CBA A A A u | AAEDCHCDBBA D AA u | BEACABEAAA ABA A | BB ICAB BB | DAGDAAACAAA | BA B CB A A A | B ABABAA A | A AAAB A | AA A 40 + A A ---+-------------+-------------+-------------+-------------+-- 20 40 60 80 100 taijyuu NOTE: 49 obs が欠損値です。 2019年10月29日 火曜日 09時02分07秒 33 Plot of resid1*pred1. Legend: A = 1 obs, B = 2 obs, etc. 40 + | | A | A A | | A A 残 20 + A A 差 | B A A | A A A A B | A A A AA AABDB A A A | A A B ABAA BBAEAACAA AA A | AA A A A CAA A AADAGABABEAB AF 0 +-------------A-A-B--AABABDACAABEA-D-BDEAFGFDCCBCCA--A---A------------ | AAAA AABCF AEBB DFEBGEBFDCCB BE A AA | A BA B CA C ABDAABCBABBBBABBA A | A DA AC BD AA AA A | A A | A -20 + --+------------+------------+------------+------------+------------+-- 30 40 50 60 70 80 予測値 taijyuu NOTE: 49 obs が欠損値です。 2019年10月29日 火曜日 09時02分07秒 34 Plot of resid1*shintyou. Legend: A = 1 obs, B = 2 obs, etc. 40 + | | A | A A | | A A 残 20 + A A 差 | B A A | A A A AAA | A A A A A AAB DB A A A | A A B ABA A BB AEA ACA A A A A | AA A A A CA A AAAD AFB BAB EAB A F 0 +--------A-AB---A-AAB-BDAC-AAB-EA--DBDCBAFFAFDBACBCC-A---A---A-------- | AAA A B B DE AEB BCGE BGDABFD CCB BE A A A | A BA B ABA CABCAAAB CBA BBBB ACA A A | A DA A C B DA A A AA | A A | A -20 + --+------------+------------+------------+------------+------------+-- 140 150 160 170 180 190 shintyou NOTE: 49 obs が欠損値です。 2019年10月29日 火曜日 09時02分07秒 35 Plot of resid1*taijyuu. Legend: A = 1 obs, B = 2 obs, etc. 40 + | | A | A A | | A A 残 20 + A A 差 | AA B | A A BB | A AA AA ADE A A A | A BABAAAI FBA B | B AAAAD ACDEDDCDGA 0 +---------------A-BAACAGCDDAHDEJMDAGBAAA---------------------- | BAADAECJDENCFLDCBD B | CA DCBBEBBDDCFAAAA | A ADB DCB ABA | A A | A -20 + ---+-------------+-------------+-------------+-------------+-- 20 40 60 80 100 taijyuu NOTE: 49 obs が欠損値です。 2019年10月29日 火曜日 09時02分07秒 36 UNIVARIATE プロシジャ 変数 : resid1 (残差) モーメント N 350 重み変数の合計 350 平均 0 合計 0 標準偏差 6.68488112 分散 44.6876356 歪度 1.32528304 尖度 3.64226422 無修正平方和 15595.9848 修正済平方和 15595.9848 変動係数 . 平均の標準誤差 0.35732193 基本統計量 位置 ばらつき 平均 0.00000 標準偏差 6.68488 中央値 -0.99468 分散 44.68764 最頻値 0.16506 範囲 48.81370 四分位範囲 7.24519 2019年10月29日 火曜日 09時02分07秒 37 UNIVARIATE プロシジャ 変数 : resid1 (残差) 位置の検定 H0: Mu0=0 検定 -統計量- ------p 値------- Student の t 検定 t 0 Pr > |t| 1.0000 符号検定 M -22 Pr >= |M| 0.0214 符号付順位検定 S -3526.5 Pr >= |S| 0.0626 正規性の検定 検定 --統計量--- ------p 値------- Shapiro-Wilk W 0.924907 Pr < W <0.0001 Kolmogorov-Smirnov D 0.102661 Pr > D <0.0100 Cramer-von Mises W-Sq 0.870412 Pr > W-Sq <0.0050 Anderson-Darling A-Sq 5.095018 Pr > A-Sq <0.0050 2019年10月29日 火曜日 09時02分07秒 38 UNIVARIATE プロシジャ 変数 : resid1 (残差) 分位点 (定義 5) 水準 分位点 100% 最大値 33.65543 99% 22.23358 95% 11.65543 90% 8.38562 75% Q3 3.09654 50% 中央値 -0.99468 25% Q1 -4.14864 10% -7.49661 5% -8.78568 1% -11.46235 0% 最小値 -15.15827 2019年10月29日 火曜日 09時02分07秒 39 UNIVARIATE プロシジャ 変数 : resid1 (残差) 極値 ----最小値---- ----最大値---- 値 Obs 値 Obs -15.1583 350 21.7239 281 -13.9035 296 22.2336 129 -13.3788 370 28.6554 347 -11.4623 242 29.2347 190 -10.9035 298 33.6554 348 欠損値 ---パーセント--- 欠損値 カウント 全体 欠損値 . 49 12.28 100.00