==== regress4-pref-death-cerebrovascular-disease-shortage-of-doctors === ========= 脳血管疾患年齢調整死亡率の ４説明変数 ========= 経度・医師不足度・【液体食品】【固体食品】トランスによる回帰 ========================================================================== $$u // ユーザデータ・セクション $c // クロスセクションデータ属性コマンド 0001.00,0047.00,aaa // ケース始点,終点番号 ,★第１系列名 ,bbb // 空白で同一ケース範囲,★第２系列名 ,ccc // 空白で同一ケース範囲,★第３系列名 ,ddd // 空白で同一ケース範囲,★第４系列名 ,eee // 空白で同一ケース範囲,★第５系列名 ,fff // 空白で同一ケース範囲,★第６系列名 ,ggg // 空白で同一ケース範囲,★第７系列名 ,hhh // 空白で同一ケース範囲,★第８系列名 $d // データ入力指示コマンド ctype // ケース毎に読むタイプ ------- ユーザ自身が文字・数値データを ------- テキストファイルまたはExcelシートからコピー＆ペーストする． ------- ユーザデータの各行の末尾にも「//」を挿入してコメント文を記述できる． ------- ★ユーザ文字・数値データをこの行直後にペーストする ★↓↓★ 43.064 141.347 62.7 34.6 162.78 45.98 61.82 54.98 40.824 140.740 84.0 45.3 164.17 38.25 67.77 58.15 39.704 141.153 81.4 44.7 143.72 47.10 50.34 46.28 38.269 140.872 71.6 39.4 149.12 44.81 54.86 49.45 39.719 140.102 76.3 39.5 152.52 35.93 64.07 52.52 38.240 140.364 66.5 39.7 168.62 50.25 63.59 54.77 37.750 140.468 73.7 42.0 147.78 37.27 68.41 42.10 36.342 140.447 72.0 44.6 159.85 43.14 67.71 49.00 36.565 139.883 79.3 46.4 156.94 43.04 66.16 47.74 36.391 139.060 66.3 43.3 156.92 40.22 57.90 58.80 35.857 139.649 64.0 37.3 168.54 51.27 61.51 55.77 35.605 140.123 62.8 37.9 173.47 57.34 54.00 62.13 35.689 139.692 59.5 35.9 160.05 52.44 49.23 58.38 35.448 139.642 58.9 33.9 176.53 55.26 58.36 62.92 37.903 139.023 69.8 37.2 155.10 44.88 48.34 61.87 36.695 137.211 65.1 35.5 141.20 45.28 43.52 52.40 36.595 136.626 53.4 33.6 152.20 45.25 52.75 54.20 36.065 136.222 56.7 35.1 141.48 44.62 49.12 47.74 35.664 138.568 56.1 33.6 172.58 44.56 68.19 59.84 36.651 138.181 68.8 41.5 186.98 46.16 83.25 57.57 35.391 136.722 57.8 37.0 160.18 44.70 53.50 61.98 34.977 138.383 66.0 39.4 154.26 44.94 56.48 52.83 35.180 136.907 59.5 38.0 157.98 47.27 49.77 60.94 34.730 136.509 61.8 36.9 157.21 51.44 56.02 49.75 35.005 135.869 52.0 33.4 174.58 53.25 55.48 65.85 35.021 135.756 58.7 32.4 187.69 58.19 59.21 70.29 34.686 135.520 53.2 31.5 171.01 53.94 48.44 68.62 34.691 135.183 54.3 32.1 190.92 58.37 52.34 80.22 34.685 135.833 49.6 29.2 191.39 60.67 61.38 69.35 34.227 135.167 49.6 33.6 166.42 50.28 50.15 66.00 35.504 134.238 65.9 37.6 172.04 50.20 61.60 60.25 35.472 133.050 54.1 30.2 173.31 42.63 68.34 62.35 34.662 133.934 61.9 34.7 162.49 51.71 51.07 59.71 34.397 132.460 55.4 31.9 179.72 56.89 59.28 63.55 34.186 131.471 64.3 39.1 169.12 47.34 60.57 61.21 34.066 134.559 58.7 36.3 173.35 52.63 60.34 60.38 34.340 134.043 52.6 31.9 162.06 53.84 47.73 60.49 33.842 132.766 61.3 32.4 176.05 51.13 64.37 60.55 33.560 133.531 72.0 37.2 148.66 37.65 58.38 52.62 33.607 130.418 55.0 30.4 152.39 48.99 50.18 53.22 33.249 130.300 60.6 31.0 151.73 46.80 57.28 47.65 32.745 129.874 60.7 34.1 169.22 44.78 63.26 61.18 32.790 130.742 58.4 32.5 157.08 47.52 60.50 49.06 33.238 131.613 59.4 36.4 166.73 49.92 60.65 56.16 31.911 131.424 64.0 37.5 144.44 37.61 63.49 43.34 31.560 130.558 70.7 42.1 175.01 42.25 80.22 52.53 26.212 127.681 51.9 23.1 147.40 39.63 60.80 46.97 --------------------------------------------------------- $c // クロスセクションデータ属性コマンド 0001.00,0047.00,rrr // ケース始点,終点番号 ,★第９系列名 $d // データ入力指示コマンド ctype // ケース毎に読むタイプ ------- ユーザ自身が文字・数値データを ------- テキストファイルまたはExcelシートからコピー＆ペーストする． ------- ユーザデータの各行の末尾にも「//」を挿入してコメント文を記述できる． ------- ★ユーザ文字・数値データをこの行直後にペーストする ★↓↓★ 1.13 1.32 1.40 1.15 1.20 1.24 1.23 1.15 1.17 1.19 1.10 1.12 1.08 1.10 1.22 1.17 1.11 1.18 1.29 1.18 1.24 1.21 1.11 1.20 1.22 1.12 1.09 1.13 1.16 1.15 1.19 1.28 1.12 1.15 1.14 1.22 1.19 1.17 1.24 1.11 1.12 1.12 1.13 1.26 1.18 1.21 1.18 ====================== $$v // 変量分析セクション $a // 変量記号の割り当て a,aaa // 緯度 b,bbb // 経度 c,ccc // 脳血管疾患疾患年齢調整死亡率・【男】 d,ddd // 脳血管疾患疾患年齢調整死亡率・【女】 e,eee // トランス脂肪酸摂取量【計】 f,fff // 【牛由来食品】トランス脂肪酸摂取量 g,ggg // 【液体食品】トランス脂肪酸摂取量 h,hhh // 【固体食品】トランス脂肪酸摂取量 r,rrr // 【医師不足度】 ---------------- $d // 表示範囲 all // 全範囲 ---------------- $t // 変数変換 P=:ci(a) // 個体識別文字列P作成 =pr*(a,b,c,d,e,f,g,h,r,P) // 素データの数値プリント C=z..(c) // 年齢調整死亡率【男】の標準化 D=z..(d) // 年齢調整死亡率【女】の標準化 z=(C,D) // 年齢調整死亡率【男】標準化と【女】の標準化の連結[z] ......... a=(a,a) // 【緯度】の【男】【女】分の連結 b=(b,b) // 【経度】の【男】【女】分の連結 v=(c,d) // 年齢調整死亡率【男】【女】の素データ連結[v]　三次元図のバブル変量として利用 e=(e,e) // トランス脂肪酸【計】の【男】【女】分の連結 f=(f,f) // 【牛由来食品】トランス脂肪酸の【男】【女】分の連結 g=(g,g) // 【液体食品】トランス脂肪酸の【男】【女】分の連結 h=(h,h) // 【固体食品】トランス脂肪酸の【男】【女】分の連結 r=(r,r) // 【医師不足度】の【男】【女】分の連結 P=(P,P) // 【印字】文字列の【男】【女】分の連結 P,nam,:ci,P=(P,P) // 【印字】連結変量Pが文字系列であることを示す変量名（先頭４文字:ci,）に変更 =pr*(a,b,v,z,e,f,g,h,r,P) // 数値プリント ---------------- $r // 回帰分析 ,run,a=(b,z,e,f,g,h,r) // 全変数間の単相関係数行列を計測 F,@"*,z=(b,r) // 死亡率【男】【女】標準化の【経度】【医師不足度】よる回帰係数F M,@"*,z=(b,e,r) // 死亡率【男】【女】標準化の【経度】【トランス計】【医師不足度】よる回帰係数M N,@"*,z=(b,r,g,h) // 死亡率標準化の【経度】【医師不足度】【液体トランス】【固体トランス】よる回帰係数N ---------------- $t // 変数変換 i=@."(M)1. // 回帰係数ベクトルMの第１説明変数の係数のみのスカラー　定数指示1.の.は必須 j=@."(M)4. // 回帰係数ベクトルMの定数項のみのスカラー　定数指示4.の.は必須 m=@."(M)2.,3. // 回帰係数ベクトルMの第２・第３説明変数の係数ベクトル抽出　定数数字後の.は必須 =pr*(M,i,j,m) // 数値プリント Z=(z-i*b-j) // 【経度】と定数による説明部分を除いた死亡率標準化データ ?Z=min(Z) // 【経度】と定数による説明部分を除いた死亡率標準化データの最小値［?Z] Z=(Z-?Z) // 【経度】と定数による説明部分を除いた死亡率標準化の最小値からの乖離（定数シフト） m=(m,-?Z) // 【トランス計】【医師不足度】の回帰係数と定数シフトによる新係数ベクトル =pr*(?Z,m) // 数値プリント =pr*(z,i*b+j,z-i*b-j,Z,P) // 数値プリント .................... k=@."(N)1.,2. // 回帰係数ベクトルNの第１・第２説明変数の係数ベクトル抽出　定数数字後の.は必須 l=@."(N)5. // 回帰係数ベクトルNの定数項のみのスカラー　定数指示5.の.は必須 n=@."(N)3.,4. // 回帰係数ベクトルNの第３・第４説明変数の係数ベクトル抽出　定数数字後の.は必須 =pr*(N,k,l,n) // 数値プリント w=(b,r)"(k,l) // ベクトル内積計算 [w]【経度】と【医師不足度】と定数による説明部分 Y=(z-w) // 【経度】と【医師不足度】と定数による説明部分[w]を除いた死亡率標準化データ ?Y=min(Y) // 【経度】【医師不足度】と定数による説明部分を除いた死亡率標準化の最小値［?Y] Y=(Y-?Y) // 【経度】【医師不足度】と定数による説明部分を除いた死亡率標準化の最小値からの乖離（定数シフト） n=(n,-?Y) // 【液体トランス】【固体トランス】の回帰係数と定数シフトによる新係数ベクトル =pr*(?Y,n) // 数値プリント =pr*(z,w,z-w,Y,P) // 数値プリント =============== $$g // グラフセクション ------------ ３次元図　------------------------------------- $3 // ３次元図　 ........ [z] 死亡率【男】標準化＋【女】標準化の連結データ z,b,r,P=v,F // 縦軸c,横軸b,奥行軸r,個体識別P=バブル変量v,関数F ........ [Z] 【経度差】考慮済み死亡率格差 Z,e,r,P=v,m // 縦軸Z,横軸e,奥行軸r,個体識別P=バブル変量v,関数m ........ [Y] 【経度差】+【医師不足度】考慮済み死亡率格差 Y,g,h,P=v,n // 縦軸Y,横軸g,奥行軸h,個体識別P=バブル変量v,関数n ------------ 散布（相関）図　------------------------------------- $c // 散布図 z 死亡率【男】標準化＋【女】標準化の連結データ z,a,*,P // 縦軸z,横軸a,回帰線*,個体識別P 　【緯度】 z,b,*,P // 縦軸z,横軸b,回帰線*,個体識別P 【経度】 z,e,*,P // 縦軸z,横軸e,回帰線*,個体識別P 【トランス計】 z,f,*,P // 縦軸z,横軸f,回帰線*,個体識別P 【牛由来トランス】 z,g,*,P // 縦軸z,横軸g,回帰線*,個体識別P 【液体トランス】 z,h,*,P // 縦軸z,横軸h,回帰線*,個体識別P 【固体トランス】 z,r,*,P // 縦軸z,横軸r,回帰線*,個体識別P 【医師不足度】 $$ // 終了セクション