====== regress3-pref-death-cerebrovascular-disease-trans-fatty-acids === ========= 脳血管疾患年齢調整死亡率の ３説明変数 ========= 緯度・【液体食品】【固体食品】トランス脂肪酸摂取量による回帰 ========================================================================== $$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 ====================== $$v // 変量分析セクション $a // 変量記号の割り当て a,aaa // 緯度 b,bbb // 経度 c,ccc // 脳血管疾患年齢調整死亡率・【男】 d,ddd // 脳血管疾患年齢調整死亡率・【女】 e,eee // トランス脂肪酸摂取量【計】 f,fff // 【牛由来食品】トランス脂肪酸摂取量 g,ggg // 【液体食品】トランス脂肪酸摂取量 h,hhh // 【固体食品】トランス脂肪酸摂取量 ---------------- $d // 表示範囲 all // 全範囲 ---------------- $t // 変数変換 P=:ci(a) // 個体識別文字列P作成 =pr*(a,b,c,d,e,f,g,h,P) // 数値プリント ---------------- $r // 回帰分析 ,run,a=(b,c,d,e,f,g,h) // 全変数間の単相関係数行列を計測 ,run,c=(a,e) // 年齢調整死亡率【男】の【経度】トランス脂肪酸【計】よる重回帰 ,run,d=(a,e) // 年齢調整死亡率【女】の【経度】トランス脂肪酸【計】よる重回帰 U,@"*,c=(g,h) // 年齢調整死亡率【男】の【液体トランス】【固体トランス】よる回帰係数U V,@"*,d=(g,h) // 年齢調整死亡率【女】の【液体トランス】【固体トランス】よる回帰係数V M,@"*,c=(a,g,h) // 年齢調整死亡率【男】の【緯度】【液体トランス】【固体トランス】よる回帰係数M N,@"*,d=(a,g,h) // 年齢調整死亡率【女】の【緯度】【液体トランス】【固体トランス】よる回帰係数N ---------------- $t // 変数変換 i=@."(M)1. // 回帰係数ベクトルMの第１説明変数の係数のみのスカラー　定数指示1.の.は必須 j=@."(M)4. // 回帰係数ベクトルMの定数項のみのスカラー　定数指示4.の.は必須 m=@."(M)2.,3. // 回帰係数ベクトルMの第２・第３説明変数の係数ベクトル抽出　定数数字後の.は必須 =pr*(M,i,j,m) // 数値プリント C=(c-i*a-j) // 【緯度】と定数による説明部分を除いた死亡率 ?C=min(C) // 【緯度】と定数による説明部分を除いた死亡率の最小値 C=(C-?C) // 【緯度】と定数による説明部分を除いた死亡率の最小値からの乖離（定数シフト） m=(m,-?C) // 【液体トランス】【固体トランス】の回帰係数と定数シフトによる新係数ベクトル .................... k=@."(N)1. // 回帰係数ベクトルNの第１説明変数の係数のみのスカラー　定数指示1.の.は必須 l=@."(N)4. // 回帰係数ベクトルNの定数項のみのスカラー　定数指示4.の.は必須 n=@."(N)2.,3. // 回帰係数ベクトルNの第２・第３説明変数の係数ベクトル抽出　定数数字後の.は必須 =pr*(N,k,l,n) // 数値プリント D=(d-k*a-l) // 【緯度】と定数による説明部分を除いた死亡率 ?D=min(D) // 【緯度】と定数による説明部分を除いた死亡率の最小値 D=(D-?D) // 【緯度】と定数による説明部分を除いた死亡率の最小値からの乖離（定数シフト） n=(n,-?D) // 【液体トランス】【固体トランス】の回帰係数と定数シフトによる新係数ベクトル .................... =pr*(?C,m,?D,n) // 数値プリント =pr*(c,i*a+j,c-i*a-j,C,d,k*a+l,d-k*a-l,D) // 数値プリント =============== $$g // グラフセクション ------------ ３次元図　------------------------------------- $3 // ３次元図　年齢調整死亡率 c,g,h,P=c,U,* // 縦軸c,横軸g,奥行軸h,個体識別P=バブル変量c,関数U,合成用保存*【男】 d,g,h,P=d,V,* // 縦軸d,横軸g,奥行軸h,個体識別P=バブル変量d,関数V,合成用保存*【女】 // 合成 $3 // ３次元図 【緯度差】考慮済み死亡率格差 C,g,h,P=c,m,* // 縦軸C,横軸g,奥行軸h,個体識別P=バブル変量c,関数m,合成用保存*【男】 D,g,h,P=d,n,* // 縦軸D,横軸g,奥行軸h,個体識別P=バブル変量d,関数n,合成用保存*【女】 // 合成 ------------ 散布（相関）図　------------------------------------- $c // 散布図 c,a,*,P // 縦軸c,横軸a,回帰線*,個体識別P 年齢調整死亡率【男】　【緯度】 c,g,*,P // 縦軸c,横軸g,回帰線*,個体識別P 年齢調整死亡率【男】　【液体トランス】 c,h,*,P // 縦軸c,横軸h,回帰線*,個体識別P 年齢調整死亡率【男】　【固体トランス】 d,a,*,P // 縦軸d,横軸a,回帰線*,個体識別P 年齢調整死亡率【女】　【緯度】 d,g,*,P // 縦軸d,横軸g,回帰線*,個体識別P 年齢調整死亡率【女】　【液体トランス】 d,h,*,P // 縦軸d,横軸h,回帰線*,個体識別P 年齢調整死亡率【女】　【固体トランス】 $$ // 終了セクション