======= regress3-pref-death-heart-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 84.8 45.2 162.78 45.98 61.82 54.98 40.824 140.740 108.0 50.2 164.17 38.25 67.77 58.15 39.704 141.153 98.1 47.5 143.72 47.10 50.34 46.28 38.269 140.872 83.8 42.9 149.12 44.81 54.86 49.45 39.719 140.102 76.6 44.5 152.52 35.93 64.07 52.52 38.240 140.364 86.5 43.6 168.62 50.25 63.59 54.77 37.750 140.468 100.5 49.3 147.78 37.27 68.41 42.10 36.342 140.447 91.0 47.5 159.85 43.14 67.71 49.00 36.565 139.883 96.2 52.3 156.94 43.04 66.16 47.74 36.391 139.060 88.2 49.2 156.92 40.22 57.90 58.80 35.857 139.649 93.5 53.0 168.54 51.27 61.51 55.77 35.605 140.123 89.8 49.3 173.47 57.34 54.00 62.13 35.689 139.692 79.5 43.3 160.05 52.44 49.23 58.38 35.448 139.642 76.0 42.7 176.53 55.26 58.36 62.92 37.903 139.023 79.3 39.6 155.10 44.88 48.34 61.87 36.695 137.211 68.3 36.5 141.20 45.28 43.52 52.40 36.595 136.626 86.1 46.5 152.20 45.25 52.75 54.20 36.065 136.222 77.5 39.7 141.48 44.62 49.12 47.74 35.664 138.568 89.8 50.9 172.58 44.56 68.19 59.84 36.651 138.181 74.4 35.6 186.98 46.16 83.25 57.57 35.391 136.722 90.5 50.4 160.18 44.70 53.50 61.98 34.977 138.383 80.7 41.3 154.26 44.94 56.48 52.83 35.180 136.907 87.4 49.5 157.98 47.27 49.77 60.94 34.730 136.509 88.3 48.5 157.21 51.44 56.02 49.75 35.005 135.869 71.9 41.0 174.58 53.25 55.48 65.85 35.021 135.756 81.8 47.4 187.69 58.19 59.21 70.29 34.686 135.520 86.9 49.7 171.01 53.94 48.44 68.62 34.691 135.183 75.8 44.9 190.92 58.37 52.34 80.22 34.685 135.833 91.3 52.6 191.39 60.67 61.38 69.35 34.227 135.167 101.0 52.2 166.42 50.28 50.15 66.00 35.504 134.238 89.0 37.9 172.04 50.20 61.60 60.25 35.472 133.050 79.0 42.5 173.31 42.63 68.34 62.35 34.662 133.934 73.7 45.8 162.49 51.71 51.07 59.71 34.397 132.460 79.7 44.0 179.72 56.89 59.28 63.55 34.186 131.471 87.7 46.5 169.12 47.34 60.57 61.21 34.066 134.559 84.0 44.1 173.35 52.63 60.34 60.38 34.340 134.043 87.9 45.9 162.06 53.84 47.73 60.49 33.842 132.766 100.0 55.5 176.05 51.13 64.37 60.55 33.560 133.531 89.7 44.5 148.66 37.65 58.38 52.62 33.607 130.418 66.8 35.6 152.39 48.99 50.18 53.22 33.249 130.300 77.3 39.4 151.73 46.80 57.28 47.65 32.745 129.874 79.2 45.7 169.22 44.78 63.26 61.18 32.790 130.742 73.1 40.6 157.08 47.52 60.50 49.06 33.238 131.613 80.2 42.9 166.73 49.92 60.65 56.16 31.911 131.424 82.4 42.1 144.44 37.61 63.49 43.34 31.560 130.558 80.3 43.7 175.01 42.25 80.22 52.53 26.212 127.681 71.2 40.0 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=(b,e) // 年齢調整死亡率【男】の【経度】トランス脂肪酸【計】よる重回帰 ,run,d=(b,e) // 年齢調整死亡率【女】の【経度】トランス脂肪酸【計】よる重回帰 U,@"*,c=(g,h) // 年齢調整死亡率【男】の【液体トランス】【固体トランス】よる回帰係数U V,@"*,d=(g,h) // 年齢調整死亡率【女】の【液体トランス】【固体トランス】よる回帰係数V M,@"*,c=(b,g,h) // 年齢調整死亡率【男】の【経度】【液体トランス】【固体トランス】よる回帰係数M N,@"*,d=(b,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*b-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*b-l) // 【経度】と定数による説明部分を除いた死亡率 ?D=min(D) // 【経度】と定数による説明部分を除いた死亡率の最小値 D=(D-?D) // 【経度】と定数による説明部分を除いた死亡率の最小値からの乖離（定数シフト） n=(n,-?D) // 【液体トランス】【固体トランス】の回帰係数と定数シフトによる新係数ベクトル .................... =pr*(?C,m,?D,n) // 数値プリント =pr*(c,i*b+j,c-i*b-j,C,d,k*b+l,d-k*b-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,b,*,P // 縦軸c,横軸b,回帰線*,個体識別P 年齢調整死亡率【男】　【経度】 c,g,*,P // 縦軸c,横軸g,回帰線*,個体識別P 年齢調整死亡率【男】　【液体トランス】 c,h,*,P // 縦軸c,横軸h,回帰線*,個体識別P 年齢調整死亡率【男】　【固体トランス】 d,b,*,P // 縦軸d,横軸b,回帰線*,個体識別P 年齢調整死亡率【女】　【経度】 d,g,*,P // 縦軸d,横軸g,回帰線*,個体識別P 年齢調整死亡率【女】　【液体トランス】 d,h,*,P // 縦軸d,横軸h,回帰線*,個体識別P 年齢調整死亡率【女】　【固体トランス】 $$ // 終了セクション